[Scoi2015]情报传递 主席树

Description
给定一个n个节点的树，有两种操作：
1、求x~y这条链上的节点数和大于c的节点数。
2、i这个节点从下一时刻开始每次操作++。

---

Sample Input
7
0 1 1 2 2 3 3
6
1 4 7 0
2 1
2 4
2 7
1 4 7 1
1 4 7 3

---

Sample Output
5 0
5 2
5 1

---

这道题挺裸的。。。
我们考虑节点的值大于c的情况，设开始累加的时间是i，当前时态为now，可得：
now - i > c，得i < now - c
那么我们就现将每个节点按照它开始累加的时间建树，求一下1~now-c-1区间有多少个数即可，其他的树上搞来搞去你就像COT那样做就好了。

---

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
int _max(int x, int y) {return x > y ? x : y;}

struct node {
    int lc, rc, c;
} t[20 * 210000]; int cnt, rt[210000];
struct edge {
    int x, y, next;
} e[410000]; int len, last[210000];
struct e {
    int opt, l, r, c;
} gg[210000];
int dep[210000], fa[210000][20];

void ins(int x, int y) {
    e[++len].x = x; e[len].y = y;
    e[len].next = last[x]; last[x] = len;
}

void Link(int &u, int l, int r, int p) {
    if(!u) u = ++cnt;
    t[u].c++;
    if(l == r) return ;
    int mid = (l + r) / 2;
    if(p <= mid) Link(t[u].lc, l, mid, p);
    else Link(t[u].rc, mid + 1, r, p);
}

void Merge(int &u1, int u2) {
    if(!u1 || !u2) {u1 = u1 + u2; return ;}
    t[u1].c += t[u2].c;
    Merge(t[u1].lc, t[u2].lc);
    Merge(t[u1].rc, t[u2].rc);
}

int query(int u1, int u2, int u3, int u4, int l, int r, int p) {
    if(l == r) return t[u1].c + t[u2].c - t[u3].c - t[u4].c;
    int mid = (l + r) / 2;
    int c = t[t[u1].lc].c + t[t[u2].lc].c - t[t[u3].lc].c - t[t[u4].lc].c;
    if(p <= mid) return query(t[u1].lc, t[u2].lc, t[u3].lc, t[u4].lc, l, mid, p);
    else return c + query(t[u1].rc, t[u2].rc, t[u3].rc, t[u4].rc, mid + 1, r, p);
}

void dfs(int x) {
    for(int i = 1; i <= 19; i++) fa[x][i] = fa[fa[x][i - 1]][i - 1];
    for(int k = last[x]; k; k = e[k].next) {
        int y = e[k].y;
        if(y != fa[x][0]) {
            fa[y][0] = x;
            dep[y] = dep[x] + 1;
            Merge(rt[y], rt[x]);
            dfs(y);
        }
    }
}

int LCA(int x, int y) {
    if(dep[x] < dep[y]) swap(x, y);
    for(int i = 19; i >= 0; i--)
        if(dep[x] - dep[y] >= (1 << i))
            x = fa[x][i];
    if(x == y) return x;
    for(int i = 19; i >= 0; i--)
        if(fa[x][i] != fa[y][i])
            x = fa[x][i], y = fa[y][i];
    return fa[x][0];
}

int main() {
    int n; scanf("%d", &n);
    for(int i = 1; i <= n; i++) {
        int x; scanf("%d", &x);
        ins(x, i); ins(i, x);
    }
    int m; scanf("%d", &m);
    for(int i = 1; i <= m; i++) {
        scanf("%d", &gg[i].opt);
        if(gg[i].opt == 2) {
            scanf("%d", &gg[i].l);
            Link(rt[gg[i].l], 1, m, i);
        }
        else scanf("%d%d%d", &gg[i].l, &gg[i].r, &gg[i].c);
    }
    dfs(1);
    for(int i = 1; i <= m; i++) {
        if(gg[i].opt == 1) {
            int x = gg[i].l, y = gg[i].r;
            int lca = LCA(x, y);
            printf("%d ", dep[x] + dep[y] - dep[lca] * 2 + 1);
            if(i - gg[i].c <= 2) printf("0\n");
            else printf("%d\n", query(rt[x], rt[y], rt[lca], rt[fa[lca][0]], 1, m, i - gg[i].c - 1));
        }
    }
    return 0;
}