-
Notifications
You must be signed in to change notification settings - Fork 21
/
Copy path101875I - I Will Go.cpp
135 lines (102 loc) · 2.42 KB
/
101875I - I Will Go.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
#include<bits/stdc++.h>
using namespace std;
//{
#define si(a) scanf("%d",&a)
#define sii(a,b) scanf("%d %d",&a,&b);
#define siii(a,b,c) scanf("%d %d %d",&a,&b,&c);
#define sl(a) scanf("%lld",&a)
#define sll(a,b) scanf("%lld %lld",&a,&b);
#define slll(a,b,c) scanf("%lld %lld %lld",&a,&b,&c);
#define outi(a) printf("%d\n",a)
#define outii(a,b) printf("%d %d\n",a,b)
#define outis(a) printf(" %d",a)
#define outl(a) printf("%lld\n",a)
#define outll(a,b) printf("%lld %lld\n",a,b)
#define outls(a) printf(" %lld",a)
#define cel(n,k) ((n-1)/k+1)
#define sets(a) memset(a, -1, sizeof(a))
#define clr(a) memset(a, 0, sizeof(a))
#define fr(n) for(int i=0;i<n;i++)
#define fr1(n) for(int i=1;i<=n;i++)
#define frj(n) for(int j=0;j<n;j++)
#define frj1(n) for(int j=1;j<=n;j++)
#define pb push_back
#define all(v) v.begin(),v.end()
#define mp make_pair
#define ff first
#define ss second
#define INF 10000007
#define fastIO() ios_base::sync_with_stdio(false); cin.tie(NULL);
typedef long long i64;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<long long,long long> pll;
//}
const int N= 1e5+5;
int n, m ,a,b, T[N], P[N][22], L[N], par[N];
vector<int> v[N], parents;
void dfs(int u, int p){
T[u]= p;
L[u]= (p>-1) ? L[p]+1 : 0;
for(int s: v[u]){
if(s==T[u]) continue;
dfs(s,u);
}
}
void init_LCA(int n){
sets(P);
fr(n){
P[i][0]= T[i];
}
for(int j=1; (1<<j)<n;j++){
fr(n){
if(P[i][j-1]!=-1)
P[i][j]= P[P[i][j-1]][j-1];
}
}
}
int lca_query(int p, int q)
{
int tmp, log, i;
if (L[p] < L[q])
tmp = p, p = q, q = tmp;
log = 1;
while (1) {
int next = log + 1;
if ((1 << next) > L[p])break;
log++;
}
for (i = log; i >= 0; i--)
if (L[p] - (1 << i) >= L[q])
p = P[p][i];
if (p == q)
return p;
for (i = log; i >= 0; i--)
if (P[p][i] != -1 && P[p][i] != P[q][i])
p = P[p][i], q = P[q][i];
return T[p];
}
main() {
int q;
sii(n, q);
fr(n) par[i]= i;
fr(n) {
si(a);
if(a==-1){
parents.pb(i);
continue;
}
v[i].pb(a);
v[a].pb(i);
}
L[1]= 0;
for(int p: parents)
dfs(p,-1);
init_LCA(n);
while(q--){
sii(a,b);
int lca= lca_query(a,b);
if(b==lca) puts("Yes");
else puts("No");
}
}