【BZOJ 2049】[Sdoi2008]Cave 洞穴勘测 LCT模板

觉得LCT是一种特别有意思的东西，%%%%%Tarjan。。。。就去学了下，顺便敲一个版。

学习博客：http://blog.csdn.net/time_major/article/details/51030654

如果splay的基础好的话，看这篇博客应该在半小时内就能搞清楚LCT的大致实现思想和基本代码了，实在讲得太详细，而且版子也很好，除了splay的消除标记部分我觉得那样写的话应该会重复消除把。。反正换了一个。

#include<cstdio>
#include<cstring>
#include<iostream>
#define ls(a) ch[a][0]
#define rs(a) ch[a][1]
#define maxn 200020
using namespace std;
int n,m,ch[maxn][2],fa[maxn],q[maxn],top,flag[maxn];

inline int Q(int u){return u==rs(fa[u]);}
inline bool isrt(int u){return !fa[u]||(u!=ls(fa[u])&&u!=rs(fa[u]));}
void push_down(int u){
    if(!u||!flag[u])return;
    flag[u]=0;swap(ls(u),rs(u));
    if(ls(u))flag[ls(u)]^=1;
    if(rs(u))flag[rs(u)]^=1;
}
void rotate(int u){
    if(u==0)return;                     
    int f=fa[u],ff=fa[f],d=!Q(u);       
    if(!isrt(f))ch[ff][Q(f)]=u;fa[u]=ff;
    if(ch[u][d])fa[ch[u][d]]=f;ch[f][!d]=ch[u][d];
    fa[f]=u;ch[u][d]=f;
}
void splay(int u){
    if(!u)return ;
    top=0;q[++top]=u;
    for(int i=u;!isrt(i);i=fa[i])q[++top]=fa[i];
    for(int i=top;i>=1;i--)push_down(q[i]);
    while(!isrt(u)){
        int f=fa[u];
        if(!isrt(f)){
            if(Q(u)==Q(f))rotate(f);
            else rotate(u);
        }
        rotate(u);
    }
}

void access(int u){
    for(int x=0;u;u=fa[x=u]){
        splay(u);rs(u)=x;
    }
}
void makert(int u){
    access(u);splay(u);flag[u]^=1;
}
void link(int a,int b){
    makert(a),fa[a]=b;
}
void cut(int a,int b){
    makert(a),access(b),splay(b);ls(b)=0,fa[a]=0;
}
void query(int a,int b){
    makert(a),access(b),splay(b);
    while(ls(b))b=ls(b);
    if(a==b)puts("Yes");else puts("No");
}

int main(){
    scanf("%d%d",&n,&m);char s[10];int a,b;
    for(int i=1;i<=m;i++){
        scanf("%s%d%d",s,&a,&b);
        if(s[0]=='Q')query(a,b);
        else if(s[0]=='D')cut(a,b);
        else link(a,b);
    }
    return 0;
}