// #pragma GCC optimize ("Ofast")
// #pragma GCC target ("avx,avx2")
#include <bits/stdc++.h>
#define ll long long
#define vi vector<int>
#define f first
#define s second
#define pb push_back
#define inf ~((long long)1<<63)
using namespace std;
const int mod = 1000000007;
void setio(string name = "") {
ios::sync_with_stdio(false), cin.tie(nullptr);
if(name.length()) {
freopen((name + ".in").c_str(), "r", stdin);
freopen((name + ".out").c_str(), "w", stdout);
}
}
class BIT {
private:
int size;
vector<int> bit;
public:
BIT(int size) : size(size), bit(size + 1) {}
void add(int ind, int val) {
ind++;
for (; ind <= size; ind += ind & -ind) { bit[ind] += val; }
}
int prev(int ind) {
ind++;
int total = 0;
for (; ind > 0; ind -= ind & -ind) { total += bit[ind]; }
return total;
}
int sum(int a, int b){
return prev(b)-(a == 0 ? 0:prev(a-1));
}
};
int N, R, timer = 0;
vector<vector<int>> cnt;
vector<int> big;
vector<int> st, en;
vector<vector<int>> g;
vector<vector<int>> ans;
vector<int> nums;
vector<BIT> sums;
void dfs(int s, vector<int>& curr){
if(big[nums[s]]+1) curr[big[nums[s]]]++;
else for(int i = 0; i<curr.size(); i++) ans[s][i] = curr[i];
st[s] = timer++;
for(int i: g[s]) dfs(i, curr);
if(big[nums[s]]+1) curr[big[nums[s]]]--;
en[s] = timer-1;
}
void getComp(int size){
vector<int> b(size);
for(int i = R-1; i>=0; i--) if(big[i]+1) b[--size] = i;
for(int i = 0; i<b.size(); i++) for(int j: cnt[b[i]])
sums[i].add(st[j], 1);
}
int main() {
setio();
int Q; cin >> N >> R >> Q;
g.resize(N, vector<int>(0));
nums.resize(N);
cin >> nums[0];
for(int i = 1; i<N; i++){
int par; cin >> par;
g[par-1].pb(i);
cin >> nums[i];
}
cnt.resize(R, vector<int>(0));
for(int i = 0; i<N; i++) cnt[--nums[i]].pb(i);
int sqrt = 1;
while(sqrt*sqrt <= N) sqrt++;
big.resize(R, -1);
int p = 0;
for(int i = 0; i<R; i++) if(cnt[i].size() >= sqrt) big[i] = p++;
st.resize(N);
en.resize(N);
ans.resize(N, vector<int>(p));
vector<int> curr(p);
dfs(0, curr);
sums.resize(p, BIT(N));
getComp(p);
BIT b(N);
unordered_map<int, int> mem;
while(Q--> 0){
int up, down; cin >> up >> down; --up, --down;
if(mem.find(up*R+down) != mem.end()) cout << mem[up*R+down] << endl;
else if(big[down]+1){
int ret = 0;
for(int i: cnt[up]) ret+=sums[big[down]].sum(st[i], en[i]);
mem[up*R+down] = ret;
cout << ret << endl;
}else if(big[up]+1){
int ret = 0;
for(int i: cnt[down])
ret+=ans[i][big[up]];
mem[up*R+down] = ret;
cout << ret << endl;
}else{
int ret = 0;
for(int i: cnt[down]) b.add(st[i], 1);
for(int i: cnt[up]) ret+=b.sum(st[i], en[i])-(up == down ? 1:0);
for(int i: cnt[down]) b.add(st[i], -1);
mem[up*R+down] = ret;
cout << ret << endl;
}
}
}
Compilation message
regions.cpp: In function 'void dfs(int, std::vector<int>&)':
regions.cpp:49:26: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
49 | else for(int i = 0; i<curr.size(); i++) ans[s][i] = curr[i];
| ~^~~~~~~~~~~~
regions.cpp: In function 'void getComp(int)':
regions.cpp:58:21: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
58 | for(int i = 0; i<b.size(); i++) for(int j: cnt[b[i]])
| ~^~~~~~~~~
regions.cpp: In function 'int main()':
regions.cpp:78:47: warning: comparison of integer expressions of different signedness: 'std::vector<int>::size_type' {aka 'long unsigned int'} and 'int' [-Wsign-compare]
78 | for(int i = 0; i<R; i++) if(cnt[i].size() >= sqrt) big[i] = p++;
| ~~~~~~~~~~~~~~^~~~~~~
regions.cpp: In function 'void setio(std::string)':
regions.cpp:15:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
15 | freopen((name + ".in").c_str(), "r", stdin);
| ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
regions.cpp:16:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
16 | freopen((name + ".out").c_str(), "w", stdout);
| ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
| # |
결과 |
실행 시간 |
메모리 |
Grader output |
| 1 |
Correct |
1 ms |
208 KB |
Output is correct |
| 2 |
Correct |
1 ms |
208 KB |
Output is correct |
| 3 |
Correct |
2 ms |
208 KB |
Output is correct |
| 4 |
Correct |
4 ms |
336 KB |
Output is correct |
| 5 |
Correct |
7 ms |
336 KB |
Output is correct |
| 6 |
Correct |
21 ms |
480 KB |
Output is correct |
| 7 |
Correct |
32 ms |
532 KB |
Output is correct |
| 8 |
Correct |
31 ms |
568 KB |
Output is correct |
| 9 |
Correct |
25 ms |
1228 KB |
Output is correct |
| 10 |
Correct |
72 ms |
1600 KB |
Output is correct |
| 11 |
Correct |
136 ms |
2168 KB |
Output is correct |
| 12 |
Correct |
151 ms |
2992 KB |
Output is correct |
| 13 |
Correct |
158 ms |
2948 KB |
Output is correct |
| 14 |
Correct |
298 ms |
5076 KB |
Output is correct |
| 15 |
Correct |
374 ms |
8280 KB |
Output is correct |
| # |
결과 |
실행 시간 |
메모리 |
Grader output |
| 1 |
Correct |
1102 ms |
24704 KB |
Output is correct |
| 2 |
Correct |
1341 ms |
17816 KB |
Output is correct |
| 3 |
Correct |
2213 ms |
21600 KB |
Output is correct |
| 4 |
Correct |
376 ms |
4428 KB |
Output is correct |
| 5 |
Correct |
433 ms |
6720 KB |
Output is correct |
| 6 |
Correct |
563 ms |
32336 KB |
Output is correct |
| 7 |
Correct |
776 ms |
36512 KB |
Output is correct |
| 8 |
Correct |
1054 ms |
96824 KB |
Output is correct |
| 9 |
Correct |
3073 ms |
22232 KB |
Output is correct |
| 10 |
Runtime error |
108 ms |
131072 KB |
Execution killed with signal 9 |
| 11 |
Correct |
5527 ms |
28308 KB |
Output is correct |
| 12 |
Correct |
1647 ms |
30540 KB |
Output is correct |
| 13 |
Correct |
2420 ms |
32484 KB |
Output is correct |
| 14 |
Correct |
2623 ms |
62132 KB |
Output is correct |
| 15 |
Correct |
3875 ms |
42116 KB |
Output is correct |
| 16 |
Correct |
4022 ms |
47656 KB |
Output is correct |
| 17 |
Correct |
3940 ms |
78344 KB |
Output is correct |
