BZOJ 2002 弹飞绵羊 Link-Cut-Tree(LCT)

首先这个序列从右往左捋一下肯定是一棵树。错了不对，是一颗森林。

而第二个操作可以随时改变这棵森林的结构

于是这道题很明显是动态树结构中的Link-Cut-Tree。

关于LCT的概念网上有很多 这里只列出一个关键点

就是如何找到链顶的父节点

其实虽然这个算法的实现方式是对每一条重链(想不出名字来了，就沿用轻重链剖分的名字吧)维护一棵Splay，但实质上每一棵树内所有的Splay树都是连在一起的。

听不懂？没事，我说过的话我自己也听不懂

拿图来说话

其中加粗的边是重链

这里x所在重链和y所在重链存在的Splay树是相连的，x的父节点指向y，而y的子节点并不指向x

即x->fa=y,y->ls=null(或者其它节点)

这样就区分了这两条重链，也保证在Splay的时候由于x->fa->ls!=x，就不会影响到y所在的Splay树，同时也使x所在重链的任意节点Splay后，fa必然指向y

。。。再也不翘语文课了

总之贴代码 实在不懂可以拿这道题的代码调试一下 十分简单明了

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 200010
inline int getint()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m;
int next[maxn],tr[maxn][2],fa[maxn],siz[maxn],st[maxn];
bool rev[maxn];
inline bool isroot(int k)
{
	return tr[fa[k]][0]!=k && tr[fa[k]][1]!=k;
}
void pushup(int x)
{
    siz[x]=siz[tr[x][0]]+siz[tr[x][1]]+1;
}
void pushdown(int x)
{
    int l=tr[x][0],r=tr[x][1];
    if(rev[x])
    {
        rev[x]=0;rev[l]^=1;rev[r]^=1;
        swap(tr[x][0],tr[x][1]);
    }
}
void rotate(int &x)
{
    int l,r,y,z;
    y=fa[x];z=fa[y];
    if(tr[y][0]==x) l=0;
    else l=1; r=l^1;
    if(!isroot(y))
    {
        if(tr[z][0]==y) tr[z][0]=x;
        else tr[z][1]=x;
    }
    fa[x]=z;fa[y]=x;fa[tr[x][r]]=y;
    tr[y][l]=tr[x][r];tr[x][r]=y;
    pushup(y);pushup(x);
}
void splay(int x)
{
    int top=0;st[++top]=x;
    for(int i=x;!isroot(i);i=fa[i])
    {
        st[++top]=fa[i];
    }
    for(int i=top;i;i--) pushdown(st[i]);
    int y,z;
    while(!isroot(x))
    {
        y=fa[x];z=fa[y];
        if(!isroot(y))
        {
            if((tr[y][0]==x)^(tr[z][0]==y)) rotate(x);
            else rotate(y);
        }
        rotate(x);
    }
}
void access(int x)
{
    int t=0;
    while(x)
    {
        splay(x);
        tr[x][1]=t;
        t=x;x=fa[x];
    }
}
void rever(int x)
{
    access(x);splay(x);rev[x]^=1;
}
void link(int x,int y)
{
    rever(x);fa[x]=y;splay(x);
}
void cut(int x,int y)
{
    rever(x);access(y);splay(y);tr[y][0]=fa[x]=0;
}
int main()
{
    n=getint();
	for(int i=1;i<=n;i++)
	{
		int x=getint();
		fa[i]=x+i;
		siz[i]=1;
		if(fa[i]>n+1)fa[i]=n+1;
		next[i]=fa[i];
	}
	siz[n+1]=1;
	m=getint();
	for(int i=1;i<=m;i++)
	{
		int op=getint();
		if(op==1)
		{
			rever(n+1);
			int x=getint();x++;
			access(x);
			splay(x);printf("%d\n",siz[tr[x][0]]);
		}
		else
		{
			int x=getint(),y=getint();x++;
			int t=min(n+1,x+y);
			cut(x,next[x]);link(x,t);next[x]=t;
		}
	}
    return 0;
}