辉辉热衷于洞穴勘测。某天,他按照地图来到了一片被标记为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
题意:一个森林,有三种操作。
L,u,v: 在u和v之间加一条边。
C,u,v:删去u和v之间的边。
Q,u,v:查询u和v是否在同一棵树里。
思路:简单的LCT。
具体的:
- 关于查找两个节点X、Y是否在同一棵树上:因为splay的根节点是不固定的,也没有代表性,所以需要找“原树”的根节点。 具体的:把Xaceess,然后splay到根,然后就可以向左边找深度第的节点,直到找到原树的根节点。
- 由于有多棵splay树,所以不能仅仅记录root节点,而需要用isroot函数判定X是否是一棵树的根。所以在rotate的时候要先把X和fa[fa[X]]的关系搞定,不然会错。
- splay操作要先pushdown。
#include<cstdio> #include<cstdlib> #include<iostream> #include<algorithm> using namespace std; const int maxn=10010; struct LCT { int rev[maxn],ch[maxn][2],fa[maxn],stc[maxn],cnt; int isroot(int x){ return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x; } void rotate(int x) { int old=fa[x],fold=fa[old]; int opt=(ch[old][1]==x); if(!isroot(old)) ch[fold][ch[fold][1]==old]=x; fa[x]=fold; //x和fold的关系最先判定,不然会和下面冲突。 ch[old][opt]=ch[x][opt^1]; fa[ch[old][opt]]=old; ch[x][opt^1]=old; fa[old]=x; } void splay(int x) { int top=0; stc[++top]=x; for(int i=x;!isroot(i);i=fa[i]) stc[++top]=fa[i]; for(int i=top;i;i--) pushdown(stc[i]); for(int f;!isroot(x);rotate(x)){ f=fa[x]; if(!isroot(f)) rotate((ch[f][1]==x)==(ch[fa[f]][1]==f)?f:x); } } void pushdown(int x) { if(!rev[x]) return ; int l=ch[x][0],r=ch[x][1]; rev[l]^=1; rev[r]^=1; rev[x]=0; swap(ch[x][0],ch[x][1]); } void access(int x) { int rson=0; while(x){ splay(x); ch[x][1]=rson; rson=x; x=fa[x]; } } int find(int x){ access(x); splay(x); while(ch[x][0]) x=ch[x][0]; return x; } int make_root(int x){ access(x); splay(x); rev[x]^=1; } void link(int x,int y){ make_root(x); fa[x]=y; splay(x); } void cut(int x,int y){ make_root(x); access(y); splay(y); fa[x]=ch[y][0]=0; } }S; char ch[10]; int main() { int n,m,x,y; scanf("%d%d",&n,&m); while(m--){ scanf("%*c%s%d%d",&ch,&x,&y); if(ch[0]=='C') S.link(x,y); if(ch[0]=='D') S.cut(x,y); if(ch[0]=='Q'){ if(S.find(x)==S.find(y)) printf("Yes "); else printf("No "); } } }