BZOJ 2049 [Sdoi2008]Cave 洞穴勘测

title: ‘BZOJ 2049 [Sdoi2008]Cave 洞穴勘测’
categories: BZOJ
date: 2016-1-14 22:16:00
tags: [LCT]

---

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.txt
200 5
Query 123 127
Connect 123 127
Query 123 127
Destroy 127 123
Query 123 127

output.txt
No
Yes
No

Solution

裸的LCT删边和添边操作。

Code

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>

#define maxn 200000+5

using namespace std;

struct Splay_Tree{
    int ch[2];
    int flip,sz,pf,fa;
}tr[maxn];

char s[8];
int st[maxn];
int n,m,top;

inline int in(){
    int x=0;
    char ch=getchar();
    while(ch<'0' || ch>'9') ch=getchar();
    while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
    return x;
}

void Pushdown(int k){
    if(tr[k].flip){
        tr[k].flip^=1;
        tr[tr[k].ch[0]].flip^=1;
        tr[tr[k].ch[1]].flip^=1;
        swap(tr[k].ch[0],tr[k].ch[1]);
    }
}

void Update(int k){
    tr[k].sz=tr[tr[k].ch[0]].sz+tr[tr[k].ch[1]].sz+1;
}

void Rotate(int &k,int d){
    int t=tr[k].ch[d^1];
    tr[t].pf=tr[k].pf; tr[t].fa=tr[k].fa;
    if(tr[t].ch[d]) tr[tr[t].ch[d]].fa=k;
    tr[k].pf=0; tr[k].fa=t;
    tr[k].ch[d^1]=tr[t].ch[d]; tr[t].ch[d]=k;
    Update(k); Update(t); k=t;
}

void Splay(int &k,int x){
    int d1=(tr[tr[k].ch[0]].sz<x?1:0),t=tr[k].ch[d1];
    if(d1==1) x-=tr[tr[k].ch[0]].sz+1;
    if(x){
        int d2=(tr[tr[t].ch[0]].sz<x?1:0);
        if(d2==1) x-=tr[tr[t].ch[0]].sz+1;
        if(x){
            Splay(tr[t].ch[d2],x);
            if(d1==d2) Rotate(k,d1^1);
            else Rotate(tr[k].ch[d1],d1);
        }
        Rotate(k,d1^1);
    }
}

int rank(int k,int &rt){
    int res=tr[tr[k].ch[0]].sz+1;
    st[++top]=k;
    while(tr[k].fa){
        int p=tr[k].fa;
        st[++top]=p;
        if(tr[p].ch[0]!=k)
            res+=tr[tr[p].ch[0]].sz+1;
        k=p;
    }
    rt=k;
    return res;
}

void Haha(int k){
    int rt,x=rank(k,rt);
    while(top) Pushdown(st[top--]);
    x=rank(k,rt);
    Splay(rt,x);
}

void Cutright(int k){
    if(tr[k].ch[1]){
        tr[tr[k].ch[1]].pf=k;
        tr[tr[k].ch[1]].fa=0;
        tr[k].ch[1]=0;
    }
}

void Access(int k){
    Haha(k); Cutright(k); Update(k);
    while(tr[k].pf){
        int u=tr[k].pf;
        Haha(u); Cutright(u);
        tr[u].ch[1]=k;
        tr[k].pf=0; tr[k].fa=u;
        Update(u); k=u;
    }
}

int Findroot(int k){
    Access(k); Haha(k);
    while(tr[k].ch[0]) k=tr[k].ch[0];
    return k;
}

void Beroot(int k){
    Access(k); Haha(k); tr[k].flip^=1;
}

void Cut(int k,int p){
    Beroot(k); Access(p); Haha(p);
    tr[p].ch[0]=tr[k].pf=tr[k].fa=0;
    Update(p);
}

void Link(int k,int p){
    int rt1=Findroot(k),rt2=Findroot(p);
    if(rt1==rt2) return;  
    Access(k); Beroot(p);
    tr[p].pf=k;
}

int main(){
    n=in(); m=in();
    for(int i=1;i<=n;i++)
        tr[i].sz=1;
    for(int i=1;i<=m;i++){
        int x,y;
        scanf("%s",s);
        x=in(); y=in();
        if(s[0]=='Q'){
            Beroot(x);
            if(Findroot(y)==x) puts("Yes");
            else puts("No");
        }
        else if(s[0]=='C') Link(x,y);
        else Cut(x,y);
    }
    return 0;
}