BZOJ 2002: [Hnoi2010]Bounce 弹飞绵羊 lct 动态树 splay

http://www.lydsy.com/JudgeOnline/problem.php?id=2002

http://blog.csdn.net/frods/article/details/52244250

本来以为lct是很难的算法没想到这么简单。。。也可能只是这道题简单。

看上去挺暴力的算法却有着很优秀的复杂度，很美丽了orz。

做法就是每次查询把子节点连一条到根的通路，所有的儿子都在父亲右儿子上，再splay查询的点。此时查询的点的左儿子的size就是所求的距离。

这个算法太神奇了orz。明明感觉有些粗笨简单但是透出朴实而容易ac的美，就像约翰内斯·维米尔的画，QAQ太令人愉悦了，QAQ谁会不喜欢splay呢。

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 using namespace std;
 const int maxn=;
 int n,m;
 int kid[maxn][]={},fa[maxn]={},siz[maxn]={};
 inline bool CheckRoot(int x){return kid[fa[x]][]!=x&&kid[fa[x]][]!=x;}
 inline void updata(int x){
     siz[x]=+siz[kid[x][]]+siz[kid[x][]];
 }
 void rotate(int x){
     int y=fa[x];int fy=fa[y];
     int l,r;
     l=kid[y][]==x?:; r=l^;
     if(!CheckRoot(y)){
         if(kid[fy][]==y)kid[fy][]=x;
         else kid[fy][]=x;
     }fa[x]=fy;
     fa[y]=x;
     fa[kid[x][r]]=y;
     kid[y][l]=kid[x][r];kid[x][r]=y;
     updata(y);//updata(x);这里不需要updata(x)因为splay只updata一下就可以惹
 }
 void Splay(int x){
     while(!CheckRoot(x)){
         int y=fa[x];int fy=fa[y];
         if(!CheckRoot(y)){
             if((x==kid[y][])^(y==kid[fy][]))rotate(x);
             else rotate(y);
         }rotate(x);
     }updata(x);
 }
 void Access(int x){
     int y=;
     while(x){
         Splay(x);kid[x][]=y;
         updata(x);
         y=x;x=fa[x];
     }
 }
 void Cut(int x){
     Access(x);Splay(x);
     kid[x][]=fa[kid[x][]]=;
 }
 void Link(int x,int y){
     Cut(x);
     fa[x]=y;
     updata(x);
 }
 int main(){
     scanf("%d",&n);
     int k,x,v;
     for(int i=;i<=n;i++){
         scanf("%d",&v);
         fa[i]=i+v>n?:i+v;
         siz[i]=;
     }
     scanf("%d",&m);
     for(int i=;i<=m;i++){
         scanf("%d%d",&k,&x);x++;
         if(k==){
             Access(x);
             Splay(x);
             printf("%d\n",siz[kid[x][]]+);
         }
         else{
             scanf("%d",&v);
             Link(x,x+v>n?:x+v);
         }
     }
     return ;
 }