【BZOJ】2049: [Sdoi2008]Cave 洞穴勘测（lct/并查集）

http://www.lydsy.com/JudgeOnline/problem.php?id=2049

bzoj挂了。。在wikioi提交，，1A～写lct的速度越来越快了～都不用debug～～

新学习了换根操作及link cut的换根操作～很简单的。到时候我开专题再说吧。

（军训终于完了T_T

#include <cstdio>

#include <cstring>

#include <string>

#include <iostream>

#include <cmath>

#include <algorithm>

using namespace std;

#define rep(i, n) for(int i=0; i<(n); ++i)

#define for1(i,a,n) for(int i=(a);i<=(n);++i)

#define for2(i,a,n) for(int i=(a);i<(n);++i)

#define for3(i,a,n) for(int i=(a);i>=(n);--i)

#define for4(i,a,n) for(int i=(a);i>(n);--i)

#define CC(i,a) memset(i,a,sizeof(i))

#define max(a,b) ((a)>(b)?(a):(b))

#define min(a,b) ((a)<(b)?(a):(b))

#define read(a) a=getnum()

#define print(a) printf("%d", a)

#define dbg(x) cout << #x << "=" << x << endl;

#define dbgarr(a, n) for1(i, 0, n) cout << a[i] << " "; cout << endl

inline int getnum() { int ret=0, k=1; char c; for(c=getchar(); c<'0' || c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0' && c<='9'; c=getchar()) ret=ret*10+c-'0'; return ret*k; }

const int N=10005;

int n, m;

struct node* null;

struct node {

	node* ch[2], *fa;

	bool rev;

	bool d() const { return fa->ch[1]==this; }

	void setc(node* c, bool d) { ch[d]=c; c->fa=this; }

	bool check() const { return fa->ch[0]!=this && fa->ch[1]!=this; }

	void pushdown() {

		if(rev) {

			ch[0]->rev^=true;

			ch[1]->rev^=true;

			swap(ch[0], ch[1]);

			rev=false;

		}

	}

}*nd[N];

inline void rot(node* r) {

	node* fa=r->fa; bool d=r->d();

	fa->pushdown(); r->pushdown();

	if(fa->check()) r->fa=fa->fa;

	else fa->fa->setc(r, fa->d());

	fa->setc(r->ch[!d], d);

	r->setc(fa, !d);

}

void fix(node* x) {

	if(!x->check()) fix(x->fa);

	x->pushdown();

}

inline void splay(node* r) {

	fix(r);

	while(!r->check())

		if(r->fa->check()) rot(r);

		else r->d()==r->fa->d()?(rot(r->fa), rot(r)):(rot(r), rot(r));

}

inline node* access(node* x) {

	node* y=null;

	for(; x!=null; y=x, x=x->fa) {

		splay(x);

		x->ch[1]=y;

	}

	return y;

}

inline void mkroot(node* x) {

	access(x)->rev^=true; splay(x);

}

inline void link(node* x, node* y) {

	mkroot(x); x->fa=y;

}

inline void cut(node* x, node* y) {

	mkroot(x);

	access(y); splay(y);

	y->ch[0]->fa=null; y->ch[0]=null;

}

inline node* findrt(node* x) {

	access(x); splay(x);

	while(x->ch[0]!=null) x=x->ch[0];

	return x;

}

int main() {

	read(n); read(m);

	null=new node; null->fa=null->ch[0]=null->ch[1]=null; null->rev=false;

	char s[10];

	int u, v;

	for1(i, 1, n) {

		nd[i]=new node;

		nd[i]->fa=nd[i]->ch[0]=nd[i]->ch[1]=null; nd[i]->rev=false;

	}

	rep(i, m) {

		scanf("%s", s);

		read(u); read(v);

		if(s[0]=='Q')

			findrt(nd[u])==findrt(nd[v])?(printf("Yes\n")):(printf("No\n"));

		else if(s[0]=='C') link(nd[u], nd[v]);

		else cut(nd[u], nd[v]);

	}

	return 0;

}

还有一种奇葩做法。。因为图只要判联通，我们干脆用并查集做，orz（在seter博客看到的http://seter.is-programmer.com/posts/29428.html）orzzzzzzz

link u v 将u提根，然后设置u的父亲为v

cut u v 将u提根，然后将v的父亲设为0（反向了所以是设v的父亲而不是u的父亲）

query u v u提根后直接暴力并查集查找v的父亲是不是u

提根将u的父亲及其边全部反向。

代码十分简洁，速度十分快，挤进了rank10.。orz

#include <cstdio>

using namespace std;

#define read(a) a=getnum()

inline int getnum() { int ret=0; char c; for(c=getchar(); c<'0' || c>'9'; c=getchar()); for(; c>='0' && c<='9'; c=getchar()) ret=ret*10+c-'0'; return ret; }

const int N=10005;

int n, m, u, v, f[N];

inline void sroot(int x) { for(int c=0, fa=f[x]; x; fa=f[x]) { f[x]=c; c=x; x=fa; } }

int main() {

	read(n); read(m); char c;

	for(int i=0; i<m; ++i) {

		for(c=getchar(); c<'A'||c>'Z'; c=getchar());

		read(u); read(v);

		sroot(u);

		if(c=='C') f[u]=v;

		else if(c=='D') f[v]=0;

		else {

			for(; v!=u && v; v=f[v]);

			puts(v==u?"Yes":"No");

		}

	}

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

Source

动态树

巴特西