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
4000
, 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

题目传送门

人生中第一发LCT，模版题，也不是我负责的专题，那我就不讲了（以后再慢慢填坑）

代码如下：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
struct node
{
int f,son[2];
bool fz;
node(){fz=false;}
}tr[11000];
void pushdown(int x)
{
tr[x].fz=false;
swap(tr[x].son[0],tr[x].son[1]);
int lc=tr[x].son[0],rc=tr[x].son[1];
tr[lc].fz=1-tr[lc].fz;
tr[rc].fz=1-tr[rc].fz;
}
void rotate(int x,int w)
{
int f=tr[x].f,ff=tr[f].f;
int R,r;

R=f;r=tr[x].son[w];
tr[R].son[1-w]=r;
if(r!=0)tr[r].f=R;

R=ff;r=x;
if(tr[R].son[0]==f)tr[R].son[0]=r;
else if(tr[R].son[1]==f)tr[R].son[1]=r;
tr[r].f=R;

R=x;r=f;
tr[R].son[w]=r;
tr[r].f=R;

}
int tmp[11000];
void splay(int x,int rt)
{
int s=0,i=x;
while(tr[i].f!=0&&(tr[tr[i].f].son[0]==i||tr[tr[i].f].son[1]==i))
{
tmp[++s]=i;
i=tr[i].f;
}
tmp[++s]=i;
while(s!=0)
{
i=tmp[s];s--;
if(tr[i].fz==true) pushdown(i);
}

while(tr[x].f!=rt&&(tr[tr[x].f].son[0]==x||tr[tr[x].f].son[1]==x))//还有虚边啊！
{
int f=tr[x].f,ff=tr[f].f;
if(ff==rt||(tr[ff].son[0]!=f&&tr[ff].son[1]!=f))
{
if(x==tr[f].son[0])rotate(x,1);
else rotate(x,0);
}
else
{
if(tr[f].son[0]==x&&tr[ff].son[0]==f){rotate(f,1);rotate(x,1);}
else if(tr[f].son[1]==x&&tr[ff].son[0]==f){rotate(x,0);rotate(x,1);}
else if(tr[f].son[0]==x&&tr[ff].son[1]==f){rotate(x,1);rotate(x,0);}
else if(tr[f].son[1]==x&&tr[ff].son[1]==f){rotate(f,0);rotate(x,0);}
}
}
}
void access(int x)//访问x
//还记得树剖的重儿子吗？这是令点x到整棵动态树的根这条路径变成偏爱路径(相当于树剖的重链)，这一条路径就是一棵伸展树。
{
int y=0;
while(x!=0)
{
splay(x,0);
tr[x].son[1]=y;
if(y!=0)tr[y].f=x;
y=x;x=tr[x].f;
}
}
void makeroot(int x)//让x成为动态树的根
{
access(x);splay(x,0);//因为是链，splay之后只有左孩子(上面y=0)
tr[x].fz=1-tr[x].fz;//因为要让x成为整棵树的根，所以x的深度要最小(通过翻转实现),为find_root做准备
}
void link(int x,int y)
{//为什么可以直接用makeroot？？不用担心构成环，x和y是不联通的。
makeroot(x);tr[x].f=y;access(x);//删去access是没有影响的，但从定义上说应该加上
}
void cut(int x,int y)
{
makeroot(x);
access(y);splay(y,0);
tr[tr[y].son[0]].f=0;tr[y].son[0]=0;
}
int find_root(int x)//访问完x后，x所属的伸展树的最左端的点就是树的根，因为伸展树实际意义上就是一条链啊！！
{
access(x);splay(x,0);
while(tr[x].son[0]!=0)x=tr[x].son[0];
return x;
}
char ss[10];
int main()
{
int n,m,x,y;
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
scanf("%s%d%d",ss+1,&x,&y);
if(ss[1]=='C')link(x,y);
else if(ss[1]=='D')cut(x,y);
else if(ss[1]=='Q')
{
if(find_root(x)==find_root(y)) printf("Yes\n");
else printf("No\n");
}
}
return 0;
}

by_lmy