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

题目传送门

好吧，有些时候不能盲目地相信板子——LCT的splay和普通的splay有区别的……

我就是因为直接copy板子而一直TLE，还不知道哪里错……

直接给出ZZK的blog，还是ZZK最强啦！！！

总之LCT的最核心的操作就是access，把一棵树分成很多实链，每条实链对应一棵splay，然后处理两个点之间的关系（连接与断开）就变成了splay的操作。

至于时间复杂度，好像也是可以证明是$O((n+m)log_2n)$的。但是我不会……

回到这题，发现就是一道LCT的裸题，直接拿来入门就行了。

附上AC代码：

#include <cstdio>
#include <cctype>
#include <algorithm>
using namespace std;

const int N=1e4+10;
int ti,f[N],ch[N][2],mk[N];

inline char nc(void){
    static char ch[100010],*p1=ch,*p2=ch;
    return p1==p2&&(p2=(p1=ch)+fread(ch,1,100010,stdin),p1==p2)?EOF:*p1++;
}

inline void read(int &a){
    static char c=nc();int f=1;
    for (;!isdigit(c);c=nc()) if (c=='-') f=-1;
    for (a=0;isdigit(c);a=(a<<3)+(a<<1)+c-'0',c=nc());
    return (void)(a*=f);
}

#define isrt(x) ((x)!=ch[f[x]][0]&&(x)!=ch[f[x]][1])
inline void rotate(int x){
    int y=f[x],op=(ch[y][1]==x);
    ch[y][op]=ch[x][op^1];
    if (ch[x][op^1]) f[ch[x][op^1]]=y;
    f[x]=f[y];
    if (!isrt(y)) ch[f[y]][ch[f[y]][1]==y]=x;
    f[y]=x,ch[x][op^1]=y;
}
inline void push(int x){if (mk[x]) mk[x]=0,mk[ch[x][0]]^=1,mk[ch[x][1]]^=1,swap(ch[x][0],ch[x][1]);}
inline void splay(int x){
    static int top,sk[N];sk[top=1]=x;
    for (int i=x; !isrt(i); i=f[i]) sk[++top]=f[i];
    while (top) push(sk[top--]);
    for (int fa=f[x]; !isrt(x); rotate(x),fa=f[x])
        if (!isrt(fa)) rotate((x==ch[fa][1])==(fa==ch[f[fa]][1])?fa:x);
}
inline void access(int x){for (int lst=0; x; lst=x,x=f[x]) splay(x),ch[x][1]=lst;}
inline void root(int x){access(x),splay(x),mk[x]^=1;}
inline void connect(int x,int y){root(x),f[x]=y;}
inline void cut(int x,int y){root(x),access(y),splay(y),f[x]=ch[y][0]=0;}
inline int get_rt(int x){access(x),splay(x);while (ch[x][0]) x=ch[x][0]; return x;}

int main(void){
    for (read(ti),read(ti); ti; --ti){
        int x,y;char c=nc();
        while (c!='Q'&&c!='C'&&c!='D') c=nc();
        read(x),read(y);
        switch (c){
            case 'Q': puts(get_rt(x)==get_rt(y)?"Yes":"No"); break;
            case 'C': connect(x,y); break;
            case 'D': cut(x,y); break;
        }
    }
    return 0;
}