Submission #229692

# Submission time Handle Problem Language Result Execution time Memory
229692 2020-05-05T22:57:38 Z xiaowuc1 None (KOI16_tree) C++14
100 / 100
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