bzoj 2049 SDOI 2008 Cave 洞穴勘测 [Link-Cut Tree]

2049: [Sdoi2008]Cave 洞穴勘测
Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 6796 Solved: 3166

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

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

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

Sample Input
样例输入1 cave.in
200 5
Query 123 127
Connect 123 127
Query 123 127
Destroy 127 123
Query 123 127

样例输入2 cave.in
3 5
Connect 1 2
Connect 3 1
Query 2 3
Destroy 1 3
Query 2 3

Sample Output
样例输出1 cave.out
No
Yes
No

样例输出2 cave.out
Yes
No

HINT
数据说明 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指令将摧毁的是一条存在的通道本题输入、输出规模比较大，建议c\c++选手使用scanf和printf进行I\O操作以免超时

---

只有link和cut两种操作的LCT。。
初写LCT。。。
splay和rotate要稍微改一下，把对to的判断改成isroot(y)
reverse操作的时候不要保持当前结点最新！！
不然TREE_ROOT满天飞。。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
const int INF=0x3f3f3f3f;
const int maxn = 10005;
struct Node
{
    int fa;
    int ch[2];
    int rev;
    bool isroot;
    Node () { fa=ch[0]=ch[1]=rev=0; isroot=true; }
}node[maxn];
int maxnode;
#define fa(x) node[x].fa
#define ch(x,d) node[x].ch[d]
#define rev(x) node[x].rev
#define isroot(x) node[x].isroot
inline void update(int x) {}
inline void pushdown(int x)
{
    if(!x) return;
    if(!rev(x)) return;
    rev(x)^=1; swap(ch(x,0),ch(x,1));
    if(ch(x,0)) rev(ch(x,0))^=1;
    if(ch(x,1)) rev(ch(x,1))^=1;
}
int sta[maxn];
int top;
inline void upgrade(int x)
{
    while(!isroot(x))
    {
        sta[++top]=x;
        x=fa(x);
    }
    sta[++top]=x;
    while(top) pushdown(sta[top--]);
}
inline void rotate(int x)
{
    int y=fa(x),z=fa(y);
    int l=(ch(y,1)==x),r=l^1;
    if(isroot(y)) isroot(y)=false,isroot(x)=true;
    else ch(z,ch(z,1)==y)=x;
    fa(ch(x,r))=y; fa(y)=x; fa(x)=z;
    ch(y,l)=ch(x,r); ch(x,r)=y;
    update(y); update(x);
}
inline void splay(int x)
{
    upgrade(x);
    while(!isroot(x))
    {
        int y=fa(x),z=fa(y);
        if(!isroot(y))
            if(ch(y,0)==x ^ ch(z,0)==y) rotate(x);
            else rotate(y);
        rotate(x);
    }
}
inline int access(int x)
{
    int y=0;
    do
    {
        splay(x);
        isroot(ch(x,1))=true;
        isroot(ch(x,1)=y)=false;
        update(x);
        x=fa(y=x);
    }while(x);
    return y;
}
void makeroot(int x)
{
    int rt=access(x);
    rev(rt)^=1;
}
void cut(int u,int v)
{
    makeroot(u);
    access(v);
    splay(v);
    fa(ch(v,0))=0; isroot(ch(v,0))=true; ch(v,0)=0; 
}
void join(int u,int v)
{
    access(u); splay(u);
    rev(u)^=1;
    fa(u)=v;
}
int find_min(int x)
{
    while(ch(x,0)) x=ch(x,0);
    return x;
}
bool query(int u,int v)
{
    access(u); splay(u);
    int t1=find_min(u);
    access(v); splay(v);
    int t2=find_min(v);
    return t1==t2;
}
int main()
{
    #ifndef ONLINE_JUDGE
    freopen("cave.in","r",stdin);
    freopen("cave.out","w",stdout);
    #endif
    int n,q;
    char cmd[10];
    int x,y;
    scanf("%d%d",&n,&q);
    for(int i=1;i<=q;i++)
    {
        scanf("%s%d%d",cmd,&x,&y);
        if(cmd[0]=='Q')
            if(query(x,y)) printf("Yes\n");
            else printf("No\n");
        else if(cmd[0]=='C') join(x,y);
             else cut(x,y);
    }
    return 0;
}