BZOJ2049 SDOI2008 Cave 洞穴勘测 【LCT】

BZOJ2049 SDOI2008 Cave 洞穴勘测

---

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操作，就可以当成Link-Cut Tree模板写，第一次写LCT还是遇到一些困难，这里的Splay作为辅助树和写伸展树时的splay有一些区别，导致调了很久没有看出来。

但是LCT本身的Access，Link，Cut函数的实现都非常容易，是一种极其方便的树上处理方式，而且LCT本身需要注意的细节比较少，而且辅助Splay的使用又比较机械，总体来说是一个较为简单的算法吧

---

#include<bits/stdc++.h>
using namespace std;
const int N=1e4+10,INF=1e9;
inline int read(){
    int ans=0,w=1;char c=getchar();
    while(!isdigit(c)&&c!='-')c=getchar();
    if(c=='-')w=-1,c=getchar();
    while(isdigit(c))ans=ans*10+c-'0',c=getchar();
    return ans*w;
}
int n,m;
int fa[N],son[N][2],st[N];
bool rev[N]={0};
inline bool isroot(int x){
    return son[fa[x]][0]!=x&&son[fa[x]][1]!=x;
}
void pushdown(int x){
    if(rev[x]){
        rev[x]^=1;
        rev[son[x][0]]^=1;
        rev[son[x][1]]^=1;
        swap(son[x][0],son[x][1]);
    }
}
inline bool Son(int t){
    return son[fa[t]][1]==t;
}
inline void rotate(int x){
    int f=fa[x],g=fa[f],a=Son(x),b=a^1;
    if(!isroot(f)){
        if(Son(f))son[g][1]=x;
        else son[g][0]=x;
    }
    fa[x]=g;
    fa[f]=x;son[f][a]=son[x][b];
    fa[son[x][b]]=f;son[x][b]=f;
}
void splay(int x){
    int top=0;st[++top]=x;
    for(int i=x;!isroot(i);i=fa[i])st[++top]=fa[i];
    for(int i=top;i;i--)pushdown(st[i]);
    while(!isroot(x)){
        int f=fa[x];
        if(!isroot(f)){
            if(Son(x)^Son(f))rotate(x);
            else rotate(f);
        }
        rotate(x);
    }
}
void access(int x){
    int t=0;
    while(x){
        splay(x);
        son[x][1]=t;
        t=x;x=fa[x];
    }
}
void reverse(int x){
    access(x);
    splay(x);
    rev[x]^=1;
}
void link(int x,int y){
    reverse(x);
    fa[x]=y;
    splay(x);
}
void cut(int x,int y){
    reverse(x);
    access(y);
    splay(y);
    son[y][0]=fa[x]=0;
}
int find(int x){
    access(x);
    splay(x);
    while(son[x][0])x=son[x][0];
    return x;
}
int main(){
    char ch[10];
    n=read();m=read();
    for(int i=1;i<=m;i++){
        scanf("%s",ch);
        int x=read(),y=read();
        if(ch[0]=='C')link(x,y);
        else if(ch[0]=='D')cut(x,y);
        else{
            if(find(x)==find(y))printf("Yes\n");
            else printf("No\n");
        }
    }
    return 0;
} 

---

第一次学，代码主体思路借鉴hzwer，感谢