Cod sursa(job #2600340)

Utilizator retrogradLucian Bicsi retrograd Data 12 aprilie 2020 13:53:29
Problema Lowest Common Ancestor Scor 60
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 3.17 kb
#include <bits/stdc++.h>

using namespace std;


struct LinkCut {
  struct Node {
    int p = 0, c[2] = {0, 0}, pp = 0;
    bool flip = 0;
    int val = 0, dp = 0;
  };
  vector<Node> T;
 
  LinkCut(int n) : T(n + 1) {}
 
  // SPLAY TREE OPERATIONS START

  int dir(int x, int y) { return T[x].c[1] == y; }
 
  void set(int x, int d, int y) {
    if (x) T[x].c[d] = y, pull(x);
    if (y) T[y].p = x;
  }
 
  void pull(int x) {
    if (!x) return;
    int &l = T[x].c[0], &r = T[x].c[1];
    T[x].dp = max({T[x].val, T[l].dp, T[r].dp});
  }
 
  void push(int x) {
    if (!x || !T[x].flip) return;
    int &l = T[x].c[0], &r = T[x].c[1];
    swap(l, r); T[l].flip ^= 1; T[r].flip ^= 1;
    T[x].flip = 0;
  }
 
  void rotate(int x, int d) { 
    int y = T[x].p, z = T[y].p, w = T[x].c[d];
    swap(T[x].pp, T[y].pp);
    set(y, !d, w);
    set(x, d, y);
    set(z, dir(z, y), x);
  }
 
  void splay(int x) { 
    for (push(x); T[x].p;) {
      int y = T[x].p, z = T[y].p;
      push(z); push(y); push(x);
      int dx = dir(y, x), dy = dir(z, y);
      if (!z) 
        rotate(x, !dx); 
      else if (dx == dy) 
        rotate(y, !dx), rotate(x, !dx); 
      else
        rotate(x, dy), rotate(x, dx);
    }
  }
 
  // SPLAY TREE OPERATIONS END
 
  void MakeRoot(int u) {
    Access(u);
    int l = T[u].c[0];
    T[l].flip ^= 1;
    swap(T[l].p, T[l].pp);
    set(u, 0, 0);
  }
 
  // int Access(int u) {
  //   int v = u, last = 0;
  //   while (v) {
  //     splay(u); splay(v);
  //     int r = T[v].c[1];
  //     swap(T[r].p, T[r].pp);
  //     set(v, 1, u == v ? 0 : u);
  //     last = v; v = T[v].pp; T[last].pp = 0;
  //   }
  //   return splay(u), last;
  // }

  // int Access(int _u) {
  //   int last = 0;
  //   for (int v = 0, u = _u; u; u = T[v = u].pp) {
  //     last = v;
  //     splay(u); splay(v);
  //     int r = T[u].c[1];
  //     T[v].pp = 0;
  //     swap(T[r].p, T[r].pp);
  //     set(u, 1, v);
  //   }
  //   return splay(_u), last;
  // }

  int Access(int u) {
    int _u = u, v = 0;
    for (; u; u = T[v = u].pp) {
      splay(u); 
      int &r = T[u].c[1];
      swap(T[r].pp, T[r].p);
      swap(T[v].pp, T[v].p);
      r = v;
    }
    return splay(_u), v;
  }

  int LCA(int u, int v) {
    Access(u);
    return Access(v);
  }

  // void Access(int u) {
  //   int _u = u;
  //   for (int v = 0; u; u = T[v = u].pp) {
  //     splay(u); 
  //     int r = T[u].c[1]; T[u].c[1] = v;
  //     swap(T[r].p, T[r].pp); 
  //     swap(T[v].p, T[v].pp);
  //   }
  //   splay(_u);
  // }

  void Link(int u, int v) { 
    assert(!Connected(u, v));
    MakeRoot(v);
    T[v].pp = u;
  }

  void Cut(int u, int v) {
    MakeRoot(u); Access(u); splay(v);
    assert(T[v].pp == u);
    T[v].pp = 0;
  }

  bool Connected(int u, int v) {
    if (u == v) return true;
    MakeRoot(u); Access(v); splay(u);
    return T[v].p == u || T[T[v].p].p == u;
  }

  int GetPath(int u, int v) {
    MakeRoot(u); Access(v); return v;
  }
};

int main() {
  ifstream cin("lca.in");
  ofstream cout("lca.out");

  int n, q; cin >> n >> q;
  LinkCut lc(n);
  for (int i = 2; i <= n; ++i) {
    cin >> lc.T[i].pp;
  }

  for (int i = 0; i < q; ++i) {
    int a, b; cin >> a >> b;
    cout << lc.LCA(a, b) << '\n';
  }

  return 0;
}