# |
Submission time |
Handle |
Problem |
Language |
Result |
Execution time |
Memory |
229692 |
2020-05-05T22:57:38 Z |
xiaowuc1 |
None (KOI16_tree) |
C++14 |
|
481 ms |
51452 KB |
#include <algorithm>
#include <bitset>
#include <cassert>
#include <chrono>
#include <complex>
#include <cstring>
#include <iomanip>
#include <iostream>
#include <map>
#include <queue>
#include <random>
#include <set>
#include <stack>
#include <vector>
using namespace std;
// BEGIN NO SAD
#define rep(i, a, b) for(int i = a; i < (b); ++i)
#define trav(a, x) for(auto& a : x)
#define all(x) x.begin(), x.end()
#define sz(x) (int)(x).size()
#define derr if(1) cerr
typedef vector<int> vi;
// END NO SAD
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<vector<ll>> matrix;
typedef pair<int, pii> state;
const int KOOSAGA_DEPTH = 18;
const int KOOSAGA_SZ = 1 << KOOSAGA_DEPTH;
bool koosagatree[2 * KOOSAGA_SZ];
void upd(int idx) {
idx += KOOSAGA_SZ;
while(idx) {
koosagatree[idx] = true;
idx /= 2;
}
}
bool qry(int lhs, int rhs) {
lhs += KOOSAGA_SZ;
rhs += KOOSAGA_SZ;
while(lhs <= rhs) {
if(lhs == rhs) return koosagatree[lhs];
if(lhs%2) if(koosagatree[lhs++]) return true;
if(rhs%2==0) if(koosagatree[rhs--]) return true;
lhs /= 2;
rhs /= 2;
}
return false;
}
vector<int> children[KOOSAGA_SZ];
int vertextoedge[KOOSAGA_SZ]; // vertextoedge[i] gives koosaga tree ID for (parent of i, i) edge
int treedepth[KOOSAGA_SZ];
int treetop[KOOSAGA_SZ];
int treesz[KOOSAGA_SZ];
int lca[KOOSAGA_SZ][KOOSAGA_DEPTH];
bool _sort_by_sz(int a, int b) {
return treesz[a] > treesz[b];
}
void dfsforhld(int curr, int top, int& id) {
treetop[curr] = top;
if(sz(children[curr]) == 0) return;
sort(all(children[curr]), _sort_by_sz);
vertextoedge[children[curr][0]] = id++;
dfsforhld(children[curr][0], top, id);
for(int i = 1; i < sz(children[curr]); i++) {
vertextoedge[children[curr][i]] = id++;
dfsforhld(children[curr][i], children[curr][i], id);
}
}
void dfsforsz(int curr) {
treesz[curr] = 1;
for(int out: children[curr]) {
treedepth[out] = treedepth[curr] + 1;
dfsforsz(out);
treesz[curr] += treesz[out];
}
}
int getlca(int a, int b) {
if(treedepth[a] < treedepth[b]) swap(a, b);
for(int d = KOOSAGA_DEPTH-1; d >= 0; d--) {
if(treedepth[a] - (1<<d) >= treedepth[b]) a = lca[a][d];
}
assert(treedepth[a] == treedepth[b]);
for(int d = KOOSAGA_DEPTH-1; d >= 0; d--) {
if(lca[a][d] != lca[b][d]) {
a = lca[a][d];
b = lca[b][d];
}
}
while(a != b) {
a = lca[a][0];
b = lca[b][0];
}
return a;
}
bool pathqry(int root, int curr) {
assert(treedepth[root] <= treedepth[curr]);
while(treedepth[curr] > treedepth[root]) {
assert(curr);
if(curr == treetop[curr]) {
if(qry(vertextoedge[curr], vertextoedge[curr])) return false;
curr = lca[curr][0];
continue;
}
if(treedepth[treetop[curr]] >= treedepth[root]) {
// go all the way to top[curr]
int nume = treedepth[curr] - treedepth[treetop[curr]] - 1;
assert(nume >= 0);
if(qry(vertextoedge[curr]-nume, vertextoedge[curr])) return false;
curr = treetop[curr];
continue;
}
// go all the way to root
int nume = treedepth[curr] - treedepth[root] - 1;
assert(nume >= 0);
if(qry(vertextoedge[curr]-nume, vertextoedge[curr])) return false;
return true;
}
return true;
}
bool connected(int a, int b) {
int r = getlca(a, b);
assert(r);
return pathqry(r, a) && pathqry(r, b);
}
void solve() {
int n, q;
cin >> n >> q;
for(int i = 2; i <= n; i++) {
cin >> lca[i][0];
children[lca[i][0]].push_back(i);
}
dfsforsz(1);
for(int d = 1; d < KOOSAGA_DEPTH; d++) {
for(int i = 1; i <= n; i++) {
lca[i][d] = lca[lca[i][d-1]][d-1];
}
}
{
int id = 0;
dfsforhld(1, 1, id);
}
while(q--) {
int a, b, t;
cin >> a >> b >> t;
assert(t == 0 || t == 1);
if(connected(a, b)) {
cout << "YES\n";
if(t==1 && a > 1) upd(vertextoedge[a]);
}
else {
cout << "NO\n";
if(t==1 && b > 1) upd(vertextoedge[b]);
}
}
}
// are there edge cases (N=1?)
// are array sizes proper (scaled by proper constant, for example 2* for koosaga tree)
// integer overflow?
// DS reset properly between test cases
int main() {
ios_base::sync_with_stdio(false);
cin.tie(NULL); cout.tie(NULL);
solve();
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
9 ms |
6784 KB |
Output is correct |
2 |
Correct |
9 ms |
6784 KB |
Output is correct |
3 |
Correct |
9 ms |
6784 KB |
Output is correct |
4 |
Correct |
9 ms |
6784 KB |
Output is correct |
5 |
Correct |
9 ms |
6784 KB |
Output is correct |
6 |
Correct |
8 ms |
6784 KB |
Output is correct |
7 |
Correct |
9 ms |
6784 KB |
Output is correct |
8 |
Correct |
9 ms |
6784 KB |
Output is correct |
9 |
Correct |
9 ms |
6784 KB |
Output is correct |
10 |
Correct |
11 ms |
6784 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
9 ms |
6784 KB |
Output is correct |
2 |
Correct |
9 ms |
6784 KB |
Output is correct |
3 |
Correct |
9 ms |
6784 KB |
Output is correct |
4 |
Correct |
9 ms |
6784 KB |
Output is correct |
5 |
Correct |
9 ms |
6784 KB |
Output is correct |
6 |
Correct |
8 ms |
6784 KB |
Output is correct |
7 |
Correct |
9 ms |
6784 KB |
Output is correct |
8 |
Correct |
9 ms |
6784 KB |
Output is correct |
9 |
Correct |
9 ms |
6784 KB |
Output is correct |
10 |
Correct |
11 ms |
6784 KB |
Output is correct |
11 |
Correct |
9 ms |
6656 KB |
Output is correct |
12 |
Correct |
8 ms |
6656 KB |
Output is correct |
13 |
Correct |
8 ms |
6656 KB |
Output is correct |
14 |
Correct |
9 ms |
6656 KB |
Output is correct |
15 |
Correct |
9 ms |
6656 KB |
Output is correct |
16 |
Correct |
9 ms |
6656 KB |
Output is correct |
17 |
Correct |
9 ms |
6656 KB |
Output is correct |
18 |
Correct |
9 ms |
6656 KB |
Output is correct |
19 |
Correct |
9 ms |
6656 KB |
Output is correct |
20 |
Correct |
9 ms |
6656 KB |
Output is correct |
21 |
Correct |
9 ms |
6656 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
9 ms |
6784 KB |
Output is correct |
2 |
Correct |
9 ms |
6784 KB |
Output is correct |
3 |
Correct |
9 ms |
6784 KB |
Output is correct |
4 |
Correct |
9 ms |
6784 KB |
Output is correct |
5 |
Correct |
9 ms |
6784 KB |
Output is correct |
6 |
Correct |
8 ms |
6784 KB |
Output is correct |
7 |
Correct |
9 ms |
6784 KB |
Output is correct |
8 |
Correct |
9 ms |
6784 KB |
Output is correct |
9 |
Correct |
9 ms |
6784 KB |
Output is correct |
10 |
Correct |
11 ms |
6784 KB |
Output is correct |
11 |
Correct |
9 ms |
6656 KB |
Output is correct |
12 |
Correct |
8 ms |
6656 KB |
Output is correct |
13 |
Correct |
8 ms |
6656 KB |
Output is correct |
14 |
Correct |
9 ms |
6656 KB |
Output is correct |
15 |
Correct |
9 ms |
6656 KB |
Output is correct |
16 |
Correct |
9 ms |
6656 KB |
Output is correct |
17 |
Correct |
9 ms |
6656 KB |
Output is correct |
18 |
Correct |
9 ms |
6656 KB |
Output is correct |
19 |
Correct |
9 ms |
6656 KB |
Output is correct |
20 |
Correct |
9 ms |
6656 KB |
Output is correct |
21 |
Correct |
9 ms |
6656 KB |
Output is correct |
22 |
Correct |
95 ms |
9720 KB |
Output is correct |
23 |
Correct |
98 ms |
9720 KB |
Output is correct |
24 |
Correct |
95 ms |
9720 KB |
Output is correct |
25 |
Correct |
101 ms |
9720 KB |
Output is correct |
26 |
Correct |
102 ms |
9720 KB |
Output is correct |
27 |
Correct |
103 ms |
9720 KB |
Output is correct |
28 |
Correct |
107 ms |
9848 KB |
Output is correct |
29 |
Correct |
106 ms |
9720 KB |
Output is correct |
30 |
Correct |
121 ms |
9848 KB |
Output is correct |
31 |
Correct |
112 ms |
9848 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
9 ms |
6784 KB |
Output is correct |
2 |
Correct |
9 ms |
6784 KB |
Output is correct |
3 |
Correct |
9 ms |
6784 KB |
Output is correct |
4 |
Correct |
9 ms |
6784 KB |
Output is correct |
5 |
Correct |
9 ms |
6784 KB |
Output is correct |
6 |
Correct |
8 ms |
6784 KB |
Output is correct |
7 |
Correct |
9 ms |
6784 KB |
Output is correct |
8 |
Correct |
9 ms |
6784 KB |
Output is correct |
9 |
Correct |
9 ms |
6784 KB |
Output is correct |
10 |
Correct |
11 ms |
6784 KB |
Output is correct |
11 |
Correct |
481 ms |
46900 KB |
Output is correct |
12 |
Correct |
465 ms |
46712 KB |
Output is correct |
13 |
Correct |
460 ms |
46712 KB |
Output is correct |
14 |
Correct |
473 ms |
46876 KB |
Output is correct |
15 |
Correct |
461 ms |
46712 KB |
Output is correct |
16 |
Correct |
477 ms |
46792 KB |
Output is correct |
17 |
Correct |
328 ms |
46844 KB |
Output is correct |
18 |
Correct |
335 ms |
51320 KB |
Output is correct |
19 |
Correct |
335 ms |
51452 KB |
Output is correct |
20 |
Correct |
324 ms |
51324 KB |
Output is correct |
21 |
Correct |
481 ms |
51192 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
9 ms |
6784 KB |
Output is correct |
2 |
Correct |
9 ms |
6784 KB |
Output is correct |
3 |
Correct |
9 ms |
6784 KB |
Output is correct |
4 |
Correct |
9 ms |
6784 KB |
Output is correct |
5 |
Correct |
9 ms |
6784 KB |
Output is correct |
6 |
Correct |
8 ms |
6784 KB |
Output is correct |
7 |
Correct |
9 ms |
6784 KB |
Output is correct |
8 |
Correct |
9 ms |
6784 KB |
Output is correct |
9 |
Correct |
9 ms |
6784 KB |
Output is correct |
10 |
Correct |
11 ms |
6784 KB |
Output is correct |
11 |
Correct |
9 ms |
6656 KB |
Output is correct |
12 |
Correct |
8 ms |
6656 KB |
Output is correct |
13 |
Correct |
8 ms |
6656 KB |
Output is correct |
14 |
Correct |
9 ms |
6656 KB |
Output is correct |
15 |
Correct |
9 ms |
6656 KB |
Output is correct |
16 |
Correct |
9 ms |
6656 KB |
Output is correct |
17 |
Correct |
9 ms |
6656 KB |
Output is correct |
18 |
Correct |
9 ms |
6656 KB |
Output is correct |
19 |
Correct |
9 ms |
6656 KB |
Output is correct |
20 |
Correct |
9 ms |
6656 KB |
Output is correct |
21 |
Correct |
9 ms |
6656 KB |
Output is correct |
22 |
Correct |
95 ms |
9720 KB |
Output is correct |
23 |
Correct |
98 ms |
9720 KB |
Output is correct |
24 |
Correct |
95 ms |
9720 KB |
Output is correct |
25 |
Correct |
101 ms |
9720 KB |
Output is correct |
26 |
Correct |
102 ms |
9720 KB |
Output is correct |
27 |
Correct |
103 ms |
9720 KB |
Output is correct |
28 |
Correct |
107 ms |
9848 KB |
Output is correct |
29 |
Correct |
106 ms |
9720 KB |
Output is correct |
30 |
Correct |
121 ms |
9848 KB |
Output is correct |
31 |
Correct |
112 ms |
9848 KB |
Output is correct |
32 |
Correct |
481 ms |
46900 KB |
Output is correct |
33 |
Correct |
465 ms |
46712 KB |
Output is correct |
34 |
Correct |
460 ms |
46712 KB |
Output is correct |
35 |
Correct |
473 ms |
46876 KB |
Output is correct |
36 |
Correct |
461 ms |
46712 KB |
Output is correct |
37 |
Correct |
477 ms |
46792 KB |
Output is correct |
38 |
Correct |
328 ms |
46844 KB |
Output is correct |
39 |
Correct |
335 ms |
51320 KB |
Output is correct |
40 |
Correct |
335 ms |
51452 KB |
Output is correct |
41 |
Correct |
324 ms |
51324 KB |
Output is correct |
42 |
Correct |
481 ms |
51192 KB |
Output is correct |
43 |
Correct |
263 ms |
30456 KB |
Output is correct |
44 |
Correct |
253 ms |
30584 KB |
Output is correct |
45 |
Correct |
262 ms |
30456 KB |
Output is correct |
46 |
Correct |
266 ms |
30456 KB |
Output is correct |
47 |
Correct |
275 ms |
30548 KB |
Output is correct |
48 |
Correct |
172 ms |
29424 KB |
Output is correct |
49 |
Correct |
174 ms |
29944 KB |
Output is correct |
50 |
Correct |
180 ms |
30072 KB |
Output is correct |
51 |
Correct |
318 ms |
39288 KB |
Output is correct |