HNOI 2010 弹飞绵羊LCT

%SAK:blog.csdn.net/frods

动态树要求我们维护一个由若干棵子结点无序的有根树组成的森林。 要求这个数据结构支持对树的分割，合并，对某个点到它的根的路径的某些操作, 以及对某个点的子树进行的某些操作。其中解决问题使用最多的是LCT，这里就主要介绍一下LCT吧。

其实动态树主题思想跟树链剖分差不多，非常频繁的使用了Splay。

类似树链剖分的定义：
PreferredChild：如果结点v的子树中, 最后被访问的结点在子树w中, 这里w是v的儿子, 那么就称w是v的PreferredChild
PreferredEdge：每个点到它的PreferredChild的边称作PreferredEdge
PreferredPath：由PreferredEdge连接成的不可再延伸的路径称为 PreferredPath

一棵LCT
容易得出整棵树就被划分成了若干条 PreferredPath。对每条 PreferredPath， 用这条路上的点的深度作为关键字, 用一棵平衡树来维护它（一般使用Splay，理论上Treap也可以，可是我从没有看见有人这样做过）。然后这棵平衡树就被叫做AuxiliaryTree。我们把AuciliaryTree中深度最小的节点的父亲节点称为PathParent。

Link-Cut Trees 就是将要维护的森林中的每棵树 T 表示为若干个 AuxiliaryTree， 并通过PathParent将这些AuxiliaryTree连接起来的数据结构

图解
^

基本操作
首先是它最基本的操作 access(x)。
一旦我们调用access(x)，那么从点x到根结点的路径就成为一条新的 PreferredPath。 如果路径上的某个节点u不是它的父亲v的PreferredChild， 那么我们要将v的PreferredChild变为u ， 原本包含v的PreferredPath/AuxiliaryTree将不再包含节点v及其之上的部分。

其他操作
-find_root(x)操作-即找到节点x所在树的根节点。
先access(x)，然后将x Splay到所在AuxiliaryTree的根节点。找到这棵AuxiliaryTree最左端的点即可。

-cut(x)操作-即断开x与其父亲节点的边。
先access(x)，然后将x Splay到所在AuxiliaryTree的根节点。断开x和父亲节点的边。

-link(v,w)操作-让v成为w的新的儿子。其中v是一棵树的根结点，并且v和w是不同的两棵树中的结点。
先访问v，然后修改v所属的AuxiliaryTree的PathParent为w，然后再次访问v。

题目大意 ：n个装置编号为0-n-1，每个装置一个权值表示一旦到这个装置上就会往后弹的步数，需要支持：

１）最少跳多少步能从输入的第x号节点弹出去；

２）更改某个装置的权值。

思路：

若可以从i弹到j则连一条从j到i的边，若弹飞了，则连一条从0到i的边。

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

#define For(i,j,k) for(int i = j;i <= k;i++)
#define Forr(i,j,k) for(int i = j;i >= k;i--)
#define Set(i,j) memset(i, j, sizeof(i))

const int N = 200010;

int n, m;

struct Link_Cut_Tree{
    int ch[N][2], sz[N], fa[N];

    bool isroot(int x){
        return x != ch[fa[x]][0] && x != ch[fa[x]][1];
    }

    void pushup(int h){
        sz[h] = sz[ch[h][0]] + sz[ch[h][1]] + 1;
    }

    void rotate(int x){
        int y = fa[x], z = fa[y], c = ch[y][1] == x;
        if(!isroot(y)) ch[z][ch[z][1] == y] = x; fa[x] = z;
        ch[y][c] = ch[x][c ^ 1], fa[ch[y][c]] = y;
        ch[x][c ^ 1] = y, fa[y] = x;
        pushup(y);
    }

    void splay(int x){
        while(!isroot(x)){
            int y = fa[x], z = fa[y];
            if(!isroot(y)){
                if((ch[y][0] == x) ^ (ch[z][0] == y)) rotate(x);
                else rotate(y);
            }
            rotate(x);
        }
        pushup(x);
    }

    void access(int x){
        int t = 0;
        while(x){
            splay(x);
            ch[x][1] = t;
            pushup(x);
            t = x, x = fa[x];
        }
    }

    void cut(int x){
        access(x);
        splay(x);
        ch[x][0] = fa[ch[x][0]] = 0;
    }

    void link(int x, int y){
        cut(x);
        fa[x] = y;
        pushup(x);
    }

}T;

int main(){
    int ty, x, y;
    scanf("%d", &n);
    For(i,1,n) scanf("%d", &x), T.fa[i] = x + i > n ? 0 : x + i, T.sz[i] = 1;
    scanf("%d", &m);
    while(m--){
        scanf("%d%d", &ty, &x);
        ++x;
        if(ty == 1){
            T.access(x);
            T.splay(x);
            printf("%d\n", T.sz[T.ch[x][0]] + 1);
        }
        else scanf("%d", &y), T.link(x, x + y > n ? 0 : x + y);
    }
    return 0;
}