bzoj2049洞穴勘测(lct模板题,lct详解)

2049: [Sdoi2008]Cave 洞穴勘测

Time Limit: 10 Sec   Memory Limit: 259 MB
Submit: 8231   Solved: 3881
[ 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

样例输入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操作以免超时

Source




对lct有了更进一步的理解,代码中对除splay操作外每一段代码都有详尽的解释,建议和论文配套食用:

                  https://wenku.baidu.com/view/75906f160b4e767f5acfcedb


#include
#define N 200020
using namespace std;
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int n,m;
bool rev[N];
int fa[N],tr[N][2];
bool isroot(int u){return tr[fa[u]][0]!=u&&tr[fa[u]][1]!=u;}
void rotate(int x)
{
	int y=fa[x],z=fa[y];
	int l=(tr[y][1]==x),r;r=l^1;
	if(!isroot(y))
	{
		if(tr[z][0]==y)tr[z][0]=x;
		else tr[z][1]=x;
	}
	fa[x]=z,fa[y]=x,fa[tr[x][r]]=y;
	tr[y][l]=tr[x][r],tr[x][r]=y;
}
// splay的同时进行pushdown,将路径上所有点入栈,自栈顶向下pushdown 
int st[N],top;
void pushdown(int u)
{
	int l=tr[u][0],r=tr[u][1];
	if(rev[u])
	{
		rev[u]^=1,rev[l]^=1,rev[r]^=1;
		swap(tr[u][0],tr[u][1]);
	}
}
void splay(int x)
{
	top=0;st[++top]=x;
	for(int i=x;!isroot(i);i=fa[i])st[++top]=fa[i];
	while(top)pushdown(st[top--]);
	while(!isroot(x))
	{
		int y=fa[x],z=fa[y];
		if(!isroot(y))
		{
			if(tr[y][0]==x^tr[z][0]==y)
				rotate(x);
			else rotate(y);
		}rotate(x);
	}
}
// access 是将x到最浅的根路径变成prefer path,将x'不断splay到当前链的根,再让x'与fa[x']连重边
void access(int x)
{
	int t=0;
	while(x)
	{
		splay(x);
		tr[x][1]=t;
		t=x;x=fa[x];
	}
}
// rever是换根操作,是将x换成原树中的根,前两个操作显然,先access,将路径变为prefer path,
// 再splay,就可以一路到根了。翻转式因为换到根后,父节点变成儿子节点 ,儿子节点变为父节点 
void rever(int x)
{
	access(x);splay(x);rev[x]^=1;
}
// link操作是将原图中不相联通的两棵树上的点链接起来,使两树合并成一棵树
// 若要链接(x,y),则先将x换成它所在树的根,让y成为它的父亲
// (不用在y中有改动,因为默认它不是prefer child) 
void link(int x,int y)
{
	rever(x);fa[x]=y;
}
// cut 操作与link操作相对,是将有连边的两个点(x,y)之间的边断开,形成两颗树
// 论文上的方式是直接将y换为其所在链的splay根,由于link-cut-tree的性质:
//     某一颗splay维护的一条链上y节点的左儿子x一定是它在原树中的父节点
// 由此将y splay到所在链的根后,其左儿子就是x,断开即可
// 但是注意到,上一行仅当x在此时原树中是y父亲成立,所以需要先满足这个条件
// 那么方便起见,就将x换成原树的根,这样x在原树中就一定是y的父亲了
// 但还要注意就是有可能y不是x的prefer son,所以需先access(y),保证y为x的prefer son 
void cut(int x,int y)
{
	rever(x);
	access(y);splay(y);
//  需要注意的是,我们发现这里对y进行的操作与rever操作只缺一个rev[y]^=1,即翻转操作
//  原因是我们只是要将y转到最顶端方便操作,不想改变原图中整棵树的父子关系,勿混淆 
	fa[x]=tr[y][0]=0; 
}
// find 操作是寻根,这里寻的根是原树中的根,根据link-cut-tree的性质(上面已经提到)
// 我们将x转到最顶端(不是换根!!!上面已经提醒过,不改变父子关系),
// 一直向左儿子走(等价于在原图中一直向父节点走,最后就会走到原图的根) 
int find(int x)
{
	access(x);splay(x);
	while(tr[x][0])x=tr[x][0];
	return x; 
}
int main()
{
	char ch[10];int x,y;
	n=read(),m=read();
	for(int i=1;i<=m;i++)
	{
		scanf("%s",ch);
		int x=read(),y=read();
		if(ch[0]=='C')link(x,y);
		else if(ch[0]=='D')cut(x,y);
		else 
		{
			int fx=find(x),fy=find(y);
			if(fx==fy)puts("Yes");
			else puts("No"); 
		}
	}
}


你可能感兴趣的:(高级数据结构)