当前位置: 代码迷 >> 综合 >> [题解]bzoj2049(SDOI2008)Cave 洞穴勘探
  详细解决方案

[题解]bzoj2049(SDOI2008)Cave 洞穴勘探

热度:53   发布时间:2023-12-22 02:57:05.0

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操作以免超时

Solution

        这一题是LCT的模板题,初学者可以练练手,只需要学会基础的link和cut操作就可以完成。(听说用不压缩路径的并查集在某些数据较水的OJ可以AC)

代码:

#include<cstdio>
#include<algorithm>
using namespace std;const int maxn=10010;
struct node{int fa,ch[2];bool reverse,is_root;
}T[maxn];
int n,m;void pushreverse(int x){if(!x)return;swap(T[x].ch[0],T[x].ch[1]);T[x].reverse^=1;
}
void pushdown(int x){if(T[x].reverse){pushreverse(T[x].ch[0]);pushreverse(T[x].ch[1]);T[x].reverse=false;}
}
int getson(int x){return x==T[T[x].fa].ch[1];
}
void rotate(int x){if(T[x].is_root)return;int k=getson(x),fa=T[x].fa;int fafa=T[fa].fa;T[fa].ch[k]=T[x].ch[k^1];if(T[x].ch[k^1])T[T[x].ch[k^1]].fa=fa;T[x].ch[k^1]=fa;T[fa].fa=x;T[x].fa=fafa;if(!T[fa].is_root)T[fafa].ch[fa==T[fafa].ch[1]]=x;else T[x].is_root=true,T[fa].is_root=false;
}
void push(int x){if(!T[x].is_root)push(T[x].fa);pushdown(x);
}
void Splay(int x){push(x);for(int fa;!T[x].is_root;rotate(x)){if(!T[fa=T[x].fa].is_root){rotate((getson(x)==getson(fa))?fa:x);}}
}
void access(int x){int y=0;do{Splay(x);T[T[x].ch[1]].is_root=true;T[T[x].ch[1]=y].is_root=false;x=T[y=x].fa;}while(x);
}
bool judge(int u,int v){while(T[u].fa)u=T[u].fa;while(T[v].fa)v=T[v].fa;return u==v;
}
void mroot(int x){access(x);Splay(x);pushreverse(x);
}
void link(int u,int v){mroot(u);T[u].fa=v;
}
void cut(int u,int v){mroot(u);Splay(v);T[T[v].ch[0]].fa=T[v].fa;T[T[v].ch[0]].is_root=true;T[v].fa=0;T[v].ch[0]=0;
}
void Work(){char opt[50];int u,v;while(m--){scanf("%s%d%d",opt,&u,&v);if(opt[0]=='Q'){if(judge(u,v))printf("Yes\n");else printf("No\n");}else if(opt[0]=='C'){link(u,v);}else if(opt[0]=='D'){cut(u,v);}}
}int main(){scanf("%d%d\n",&n,&m);for(int i=1;i<=n;i++){T[i].reverse=T[i].fa=T[i].ch[0]=T[i].ch[1]=0;T[i].is_root=true;}Work();return 0;
}