2049: [Sdoi2008]Cave 洞穴勘测
Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 9096 Solved: 4351
[Submit][Status][Discuss]
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
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
No
Yes
No
样例输出2 cave.out
Yes
No
思路{
LCT的入门题吧。。。。。思想+代码看了一个晚自习才看懂。。.。。
我终于发现我的Splay之前写得多么丑。。。。。
自己也终于LCT入门辣。。。~(≧▽≦)/~啦啦啦～～～～兴奋。。。。。
说说自己的理解吧：LCT就是以深度为键值用平衡树维护树上的每一条Preferred path！
access(x)就是把从x到根的路径变为关键路径！注意：这时x没有Preferred child！
rever(x)则是先access(x),然后把xSplay到根，然后区间翻转一下，x就成了根了。。。。
link(x,y)将x和y连边，先rever(x),把x变为其树的根，然后把x作为一条Preferred path直接连到y下面！
cut(x,y)删除x和y的连边，首先保证正确性，先rever(x)到根，
再与y连成Preferred pat直接Splay(y),切断即可。
若不rever(x)，cut之后可能还是有问题。。。。
findrt(x)找x所树中的根节点，先access(x),再跳ch[x][0]都可以啦。。。。
反正自己非常弱啊，LCT还有好多啦。。。。。。
}
#include<bits/stdc++.h>
#define RG register
#define il inline
#define N 10010
using namespace std;
int ch[N][2],fa[N],n,Q,st[N];bool lazy[N];
bool isroot(int x){return ch[fa[x]][1]!=x&&ch[fa[x]][0]!=x;}
void down(int x){
if(lazy[x]){
lazy[ch[x][0]]^=1,lazy[ch[x][1]]^=1;
swap(ch[x][0],ch[x][1]);
lazy[x]^=1;
}return;
}
void Rotate(int x){
int y=fa[x],z=fa[y],l,r;
if(ch[y][0]==x)l=0;else l=1;r=l^1;
if(!isroot(y))ch[z][ch[z][1]==y]=x;
fa[x]=z,fa[y]=x;fa[ch[x][r]]=y;
ch[y][l]=ch[x][r],ch[x][r]=y;
}
void Splay(int x){int y=x;
int top(0);st[++top]=x;
while(!isroot(y))st[++top]=fa[y],y=fa[y];
for(int i=top;i;i--)down(st[i]);
while(!isroot(x)){
y=fa[x];int z=fa[y];
if(!isroot(y)){
if(ch[z][0]==y^ch[y][0]==x)Rotate(x);
else Rotate(y);
}Rotate(x);
}
}
void access(int x){int t=0;while(x)Splay(x),ch[x][1]=t,t=x,x=fa[x];}
void rever(int x){access(x),Splay(x),lazy[x]^=1;}
int findrt(int x){access(x),Splay(x);while(ch[x][0])x=ch[x][0];return x;}
void link(int x,int y){rever(x);fa[x]=y;Splay(x);}
void cut(int x,int y){rever(x);access(y);Splay(y);ch[y][0]=fa[x]=0;}
void Query(int x,int y){if(findrt(x)==findrt(y))cout<<"Yes
";else cout<<"No
";}
int main(){
freopen("cave.in","r",stdin);
freopen("cave.out","w",stdout);
char cha[15];scanf("%d%d",&n,&Q);
for(int i=1;i<=Q;++i){
scanf("%s",cha);int x,y;scanf("%d%d",&x,&y);
if(cha[0]=='C')link(x,y);
if(cha[0]=='D')cut(x,y);
if(cha[0]=='Q')Query(x,y);
}
return 0;
}