【BZOJ】洞穴勘测-Link-cut-tree

传送门：BZOJ洞穴勘测
LCT教程：LCT（Link-Cut Tree）详解（蒟蒻自留地）

---

题意

    辉辉热衷于洞穴勘测。某天，他按照地图来到了一片被标记为JSZX的洞穴群地区。经过初步勘测，辉辉发现这片区域由n个洞穴（分别编号为1到n）以及若干通道组成，并且每条通道连接了恰好两个洞穴。假如两个洞穴可以通过一条或者多条通道按一定顺序连接起来，那么这两个洞穴就是连通的，按顺序连接在一起的这些通道则被称之为这两个洞穴之间的一条路径。洞穴都十分坚固无法破坏，然而通道不太稳定，时常因为外界影响而发生改变，比如，根据有关仪器的监测结果，123号洞穴和127号洞穴之间有时会出现一条通道，有时这条通道又会因为某种稀奇古怪的原因被毁。辉辉有一台监测仪器可以实时将通道的每一次改变状况在辉辉手边的终端机上显示：如果监测到洞穴u和洞穴v之间出现了一条通道，终端机上会显示一条指令 Connect u v 如果监测到洞穴u和洞穴v之间的通道被毁，终端机上会显示一条指令 Destroy u v 经过长期的艰苦卓绝的手工推算，辉辉发现一个奇怪的现象：无论通道怎么改变，任意时刻任意两个洞穴之间至多只有一条路径。因而，辉辉坚信这是由于某种本质规律的支配导致的。因而，辉辉更加夜以继日地坚守在终端机之前，试图通过通道的改变情况来研究这条本质规律。然而，终于有一天，辉辉在堆积成山的演算纸中崩溃了……他把终端机往地面一砸（终端机也足够坚固无法破坏），转而求助于你，说道：“你老兄把这程序写写吧”。辉辉希望能随时通过终端机发出指令 Query u v，向监测仪询问此时洞穴u和洞穴v是否连通。现在你要为他编写程序回答每一次询问。已知在第一条指令显示之前，JSZX洞穴群中没有任何通道存在。

---

输入

第一行为两个正整数n和m，分别表示洞穴的个数和终端机上出现过的指令的个数。以下m行，依次表示终端机上出现的各条指令。每行开头是一个表示指令种类的字符串s（”Connect”、”Destroy”或者”Query”，区分大小写），之后有两个整数u和v (1≤u, v≤n且u≠v) 分别表示两个洞穴的编号。

---

输出

对每个Query指令，输出洞穴u和洞穴v是否互相连通：是输出”Yes”，否则输出”No”。（不含双引号）

---

数据范围

10%的数据满足n≤1000, m≤20000 20%的数据满足n≤2000, m≤40000 30%的数据满足n≤3000, m≤60000 40%的数据满足n≤4000, m≤80000 50%的数据满足n≤5000, m≤100000 60%的数据满足n≤6000, m≤120000 70%的数据满足n≤7000, m≤140000 80%的数据满足n≤8000, m≤160000 90%的数据满足n≤9000, m≤180000 100%的数据满足n≤10000, m≤200000 保证所有Destroy指令将摧毁的是一条存在的通道

---

题意

LCT裸题……
只需要link-cut的操作。具体看代码。

---

代码

#include
#include
#include
#include
#include
using namespace std;
const int N=1e6+10;
int n,m,c[N][2],f[N],rev[N];
int now[N];
char s[10];

inline int read(){
    char ch=getchar();int x=0;
    while(ch<'0' || ch>'9') ch=getchar();
    while(ch>='0' && ch<='9') {x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x;
}

inline bool isroot(int k)
{
    return (c[f[k]][0]!=k && c[f[k]][1]!=k);
}

inline void pushdown(int k)
{
    if(rev[k]){
        rev[c[k][0]]^=1;
        rev[c[k][1]]^=1;
        swap(c[k][0],c[k][1]);
        rev[k]=0;
    }
}

inline void rotate(int x)
{
    int y=f[x];int z=f[y];
    int ptr= c[y][1]==x;
    if(!isroot(y)){
        c[z][c[z][1]==y]=x;
    }
    f[x]=z;f[y]=x;f[c[x][ptr^1]]=y;
    c[y][ptr]=c[x][ptr^1];c[x][ptr^1]=y;
}

inline void splay(int x)
{
    int tot=0,p=x;now[++tot]=x;
    for(;!isroot(p);p=f[p]){
        now[++tot]=f[p];
    }
    for(int i=tot;i>0;i--) pushdown(now[i]);
    while(!isroot(x)){
        int y=f[x];int z=f[y];
        if(!isroot(y)){
            if( (c[y][0]==z)^(c[z][0]==y)) rotate(x);
            else rotate(y);
        }
        rotate(x);
    }
}

inline void access(int k)
{
    int t=0;
    while(k){
        splay(k);
        c[k][1]=t;
        t=k;
        k=f[k];
    }
}

inline void rever(int k)
{
    access(k);splay(k);rev[k]^=1;
    //mroot:access->splay->rev[]^=1; 
}

inline void cut(int u,int v)
{
    rever(u);access(v);splay(v);c[v][0]=f[u]=0;
}

inline void link(int u,int v)
{
    rever(u);f[u]=v;splay(u);
}

inline int find(int k)
{
    access(k);splay(k);
    int p=k;
    while(c[p][0]) p=c[p][0];
    return p;
}

int main(){
    n=read();m=read();
    while(m--){
        scanf("%s",s);int u=read(),v=read();
        if(s[0]=='Q'){
            if(find(u)==find(v)) printf("Yes\n");
            else printf("No\n");
        }else if(s[0]=='C') link(u,v);
        else cut(u,v);
    }
    return 0;
}