[bzoj2049][Sdoi2008]Cave 洞穴勘测

来自FallDream的博客，未经允许，请勿转载，谢谢。

---

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

n<=10000  m<=200000

貌似是一道lct模板题啊qaq   我一开始写了什么哈希一下边然后线段树分治加启发式合并啥的 狂wa不止... 丢人啊。

然后交了以前的lct模板居然T了  我开始怀疑是不是以前的模板是错的.......不说了 抄一个靠谱的模板，以后写lct就用它了。

#include<iostream>
#include<cstdio>
#define MN 10000
using namespace std;
inline int read()
{
    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,fa[MN+5],c[MN+5][2],rev[MN+5],q[MN+5],top;
char op[20];

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

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

void rotate(int x)
{
    int y=fa[x],z=fa[y],l=c[y][1]==x,r=l^1;
    if(!isroot(y)) c[z][c[z][1]==y]=x;
    fa[x]=z;fa[y]=x;fa[c[x][r]]=y;
    c[y][l]=c[x][r];c[x][r]=y;
}

void splay(int x)
{
    q[top=1]=x;
    for (int i=x;!isroot(i);) q[++top]=i=fa[i];
    for(int i=top;i;i--) pushdown(q[i]);
    while(!isroot(x))
    {
        int y=fa[x],z=fa[y];
        if(!isroot(y))
        {
            if(c[y][1]==x^c[z][1]==y) rotate(x);
            else rotate(y);
        }
        rotate(x);
    }    
}

void access(int x)
{
    for(int y=0;x;x=fa[y=x])
       splay(x),c[x][1]=y;
}

void link(int x,int y)
{
    access(x);splay(x);rev[x]^=1;
    access(y);c[y][1]=x;fa[x]=y;
}

void cut(int x,int y)
{
    access(x);splay(x);rev[x]^=1;
    access(y);splay(y);c[y][0]=fa[x]=0;
}

int find(int x)
{
    access(x);splay(x);
    for(;c[x][0];x=c[x][0]);
    return x;
}

int main()
{
    n=read();m=read();
    for(int i=1;i<=m;i++)
    {
        scanf("%s",op);int x=read(),y=read();
        if(op[0]=='C') link(x,y);
        if(op[0]=='D') cut(x,y);
        if(op[0]=='Q') puts(find(x)==find(y)?"Yes":"No");
    }
    return 0;
}