[BZOJ]2049 [洞穴勘测] LCT[裸]

2049: [Sdoi2008]Cave 洞穴勘测

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

[ Submit][ Status][ Discuss]

几句闲话

           终于学LCT了，早就听说LCT的强大功能...不过话说因为splay自带巨大常数，LCT也很慢啊.不过功能也很强就是咯...感觉splay是个好东西.

题解

           这道题确实是裸得不能再裸的LCT了，你要是懂了LCT的话，你就会知道这道题是大水题，这里暂时没有写LCT的博客，因为还不是很熟练，不敢乱说瞎bb啥的...但是日后一定会填这个坑，比如说以【史上最详细讲解LCT】或【浅谈LCT】或【论LCT】或什么感觉好像很厉害好像国家集训队似的水平的标题出世...回归正题，这道题连边砍边判连通性，你总不能用并查集吧...看到link cut就上link cut tree LCT.初学模板参照hzwer(%%%%%%)，给出一份更精简的代码.(不算读优60行左右).

   代码:

     

#include
#include
using namespace std;
const int maxn=10005;
char ss[10];
bool rev[maxn];
int fa[maxn],c[maxn][2],s[maxn],n,m;
inline const int read(){
    register int f=1,x=0;
    register char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return f*x;
} 
inline bool isroot(int x){
    return c[fa[x]][0]!=x&&c[fa[x]][1]!=x;
}
inline void pushdown(int x){
    int l=c[x][0],r=c[x][1];
    if(rev[x]){
       rev[x]^=1,rev[l]^=1,rev[r]^=1;
       swap(c[x][0],c[x][1]);
	}
}
inline void rotate(int x){
    int y=fa[x],z=fa[y];
    int l=(c[y][0]!=x),r=l^1;
    if(!isroot(y)){
       if(c[z][0]==y) c[z][0]=x;
       else c[z][1]=x;
	}
	fa[x]=z,fa[y]=x,fa[c[x][r]]=y;
	c[y][l]=c[x][r],c[x][r]=y;
}
inline void splay(int x){
	int top=0;s[++top]=x;
    for(int i=x;!isroot(i);i=fa[i]) s[++top]=fa[i];//
    for(int i=top;i;i--) pushdown(s[i]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if(!isroot(y)){
           if(c[y][0]==x^c[z][0]==y) rotate(x);
           else rotate(y);
		}
		rotate(x);
	}
}
inline void access(int x){
    for(int t=0;x;x=fa[x]) splay(x),c[x][1]=t,t=x;
}
inline void make_root(int x){
    access(x),splay(x),rev[x]^=1;
}
inline void cut(int x,int y){
    make_root(x),access(y),splay(y),c[y][0]=fa[x]=0;
}
inline void link(int x,int y){
    make_root(x),fa[x]=y,splay(x);
}
inline int find(int x){
    access(x),splay(x);
    int y=x;
    while(c[y][0]) y=c[y][0];
    return y;
}
int main(){
    n=read(),m=read();
    while(m--){
      scanf("%s",ss);
      int x=read(),y=read();
	  if(ss[0]=='C') link(x,y);
	  if(ss[0]=='D') cut(x,y);
	  if(ss[0]=='Q') puts(find(x)==find(y)?"Yes":"No");
	}
	return 0;
}