#include <bits/stdc++.h>
#define int long long
using namespace std;
const int sz=400;
int lim,n,R,q,r1,r2,cnt[25001],c[25001][200000/sz+1],c2[25001][200000/sz+1],sum[400001],r[200001],p[200001],tin[200001],tout[200001],t,C[25001],pos[25001];
vector <int> ke[200001],ve[25001],V;
void dfs(int u){
tin[u]=++t;
ve[r[u]].push_back(t);
for (int i=1;i<=lim;i++)
c[r[u]][i]+=cnt[i];
if (C[r[u]]>sz)
cnt[pos[r[u]]]++;
for (int v:ke[u])
dfs(v);
tout[u]=++t;
ve[r[u]].push_back(-t);
if (C[r[u]]>sz)
cnt[pos[r[u]]]--;
}
int32_t main(){
ios_base::sync_with_stdio(NULL);cin.tie(nullptr);
cin >> n >> R >> q;
for (int i=1;i<=n;i++){
if (i>1){
cin >> p[i];
ke[p[i]].push_back(i);
}
cin >> r[i];
C[r[i]]++;
}
for (int i=1;i<=R;i++)
V.push_back(i);
sort(V.begin(),V.end(),[](int i, int j){return C[i]>C[j];});
for (int i=0;i<V.size();i++){
pos[V[i]]=i+1;
if (C[V[i]]>sz)
lim=i+1;
}
dfs(1);
for (int i=1;i<=R;i++)
if (C[i]>sz){
memset(sum,0,sizeof(sum));
for (int j:ve[i])
if (j>0)
sum[j]++;
for (int j=1;j<=n*2;j++)
sum[j]+=sum[j-1];
for (int j=1;j<=n;j++)
c2[r[j]][pos[i]]+=sum[tout[j]]-sum[tin[j]-1];
}
while (q--){
cin >> r1 >> r2;
if (ve[r1].empty()||ve[r2].empty()){
cout << 0 << endl;
continue;
}
if (C[r1]>sz){
cout << c[r2][pos[r1]] << endl;
continue;
}
if (C[r2]>sz){
cout << c2[r1][pos[r2]] << endl;
continue;
}
int j=0,cnt=0,res=0;
for (int i:ve[r2]){
while (j<ve[r1].size()){
if (abs(ve[r1][j])>abs(i))
break;
cnt+=(ve[r1][j]>0)-(ve[r1][j]<0);
j++;
}
if (i>0)
res+=cnt;
}
cout << res << endl;
}
}
Compilation message
regions.cpp: In function 'int32_t main()':
regions.cpp:35:19: warning: comparison of integer expressions of different signedness: 'long long int' and 'std::vector<long long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
35 | for (int i=0;i<V.size();i++){
| ~^~~~~~~~~
regions.cpp:68:21: warning: comparison of integer expressions of different signedness: 'long long int' and 'std::vector<long long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
68 | while (j<ve[r1].size()){
| ~^~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
10840 KB |
Output is correct |
2 |
Correct |
2 ms |
10840 KB |
Output is correct |
3 |
Correct |
2 ms |
10840 KB |
Output is correct |
4 |
Correct |
5 ms |
10840 KB |
Output is correct |
5 |
Correct |
4 ms |
10840 KB |
Output is correct |
6 |
Correct |
10 ms |
10840 KB |
Output is correct |
7 |
Correct |
15 ms |
10840 KB |
Output is correct |
8 |
Correct |
19 ms |
10840 KB |
Output is correct |
9 |
Correct |
28 ms |
11352 KB |
Output is correct |
10 |
Correct |
43 ms |
13412 KB |
Output is correct |
11 |
Correct |
58 ms |
13812 KB |
Output is correct |
12 |
Correct |
88 ms |
14432 KB |
Output is correct |
13 |
Correct |
90 ms |
14124 KB |
Output is correct |
14 |
Correct |
112 ms |
14716 KB |
Output is correct |
15 |
Correct |
170 ms |
18416 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
601 ms |
24088 KB |
Output is correct |
2 |
Correct |
671 ms |
22956 KB |
Output is correct |
3 |
Correct |
1121 ms |
26444 KB |
Output is correct |
4 |
Correct |
147 ms |
14700 KB |
Output is correct |
5 |
Correct |
201 ms |
16976 KB |
Output is correct |
6 |
Correct |
385 ms |
75352 KB |
Output is correct |
7 |
Correct |
631 ms |
101100 KB |
Output is correct |
8 |
Correct |
592 ms |
107736 KB |
Output is correct |
9 |
Correct |
1016 ms |
23052 KB |
Output is correct |
10 |
Runtime error |
81 ms |
131072 KB |
Execution killed with signal 9 |
11 |
Correct |
1705 ms |
22172 KB |
Output is correct |
12 |
Runtime error |
74 ms |
131072 KB |
Execution killed with signal 9 |
13 |
Runtime error |
65 ms |
131072 KB |
Execution killed with signal 9 |
14 |
Runtime error |
114 ms |
131072 KB |
Execution killed with signal 9 |
15 |
Runtime error |
71 ms |
131072 KB |
Execution killed with signal 9 |
16 |
Runtime error |
58 ms |
131072 KB |
Execution killed with signal 9 |
17 |
Runtime error |
83 ms |
131072 KB |
Execution killed with signal 9 |