2049: [Sdoi2008]Cave 洞穴勘测
Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 10405 Solved: 5070
[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
HINT
100%的数据满足n≤10000, m≤200000 保证所有Destroy指令将摧毁的是一条存在的通道本题输入、输出规模比较大，建议cc++选手使用scanf和printf进行IO操作以免超时.
分析：LCT的模板题.
维护五种操作：
1.Reverse(x) :先Access(x),然后Splay(x),最后splay区间翻转.
2.Access(x)：从x往上跳到根，每次都找到t的父亲x,splay(x),x的右儿子设为t.
3.Link(x,y)：Reverse(x),x的父亲设置为y,splay(x).
4.Cut(x,y)：Reverse(x),Access(y),splay(y),将y的左子树切断.
5.Find(x)：Access(x),splay(x),从x开始不断往左子树上跳，直到不能跳为止.
发现我以前的splay模板并不能过，于是换了一个代码比较短的模板，比较容易写错. splay那里要倒着更新.
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; const int maxn = 200010; int n,m,son[maxn][20],fa[maxn],rev[maxn],sta[maxn]; void pushdown(int x) { if (rev[x]) { rev[son[x][0]] ^= 1; rev[son[x][1]] ^= 1; rev[x] = 0; swap(son[x][0],son[x][1]); } } bool get(int x) { return son[fa[x]][1] == x; } bool is_root(int x) { return son[fa[x]][0] != x && son[fa[x]][1] != x; } void zhuan(int x) { int y = fa[x]; int z = fa[y]; int temp = get(x); if (!is_root(y)) son[z][son[z][1] == y] = x; fa[x] = z; son[y][temp] = son[x][temp ^ 1]; fa[son[y][temp]] = y; son[x][temp ^ 1] = y; fa[y] = x; } void splay(int x) { int top = 0; sta[++top] = x; int y = x; for (; !is_root(y); y = fa[y]) sta[++top] = fa[y]; for (int i = top; i >= 1; i--) pushdown(sta[i]); int temp; for (;!is_root(x);zhuan(x)) { if (!is_root(temp = fa[x])) { if (get(temp) == get(x)) zhuan(temp); else zhuan(x); } } } void Access(int x) { int t = 0; for (; x; t = x,x = fa[x]) { splay(x); son[x][1] = t; } } 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() { scanf("%d%d",&n,&m); while (m--) { char ch[10]; scanf("%s",ch); int x,y; scanf("%d%d",&x,&y); if (ch[0] == 'Q') { if (find(x) == find(y)) puts("Yes"); else puts("No"); } if (ch[0] == 'C') Link(x,y); if (ch[0] == 'D') Cut(x,y); } return 0; }