bzoj2049-洞穴勘测(动态树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

样例输入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

题意: 有三种操作，建边，毁边以及查询两个点是否相通，保证任意两个点之间最多只有一条相通的路径，也就是说整个图总是树或森林。

解析: 动态树模板题，另附两个介绍动态树的网址：

https://oi.abcdabcd987.com/summary-of-link-cut-tree/

http://wenku.baidu.com/view/75906f160b4e767f5acfcedb.html

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=10005;
struct lct
{
    lct *fa,*son[2];
    int rev;
    void pushdown() //把延迟更新的标记清除掉
    {
        if(!rev) return;  
        swap(son[0],son[1]); //翻转，交换左右儿子
        son[0]->rev^=1; son[1]->rev^=1; 
        rev=0;
    }
};
struct LCT
{
    lct data[maxn];
    lct *null;
    void init(int Size=maxn-1) //初始化
    {
        null=data; //null指向首元素
        for(int i=0;i<=Size;i++)
            data[i].son[0]=data[i].son[1]=data[i].fa=null;
    }
    bool Same(lct* x,lct* &y) //判断x和x的父亲是否在同一树里
    {
        return (y=x->fa)!=null&&(y->son[0]==x||y->son[1]==x);
    }
    void Rotate(lct* x,int d) //翻转
    {
        lct* y=x->fa; //x的父亲
        y->son[d^1]=x->son[d];
        if(x->son[d]!=null) x->son[d]->fa=y; //x的子节点的父亲指向y
        x->fa=y->fa;  //连接
        if(y->fa->son[0]==y) x->fa->son[0]=x;
        else if(y->fa->son[1]==y) x->fa->son[1]=x;
        x->son[d]=y;
        y->fa=x;
    }
    void Splay(lct* x)
    {
        x->pushdown(); //清除标记
        lct* y;
        while(Same(x,y)) //没有到树的最顶点
        {
            y->pushdown(); 
            x->pushdown();
            Rotate(x,y->son[0]==x); //翻转
        }
    }
    lct* Access(lct* u)  //打通路径
    {
        lct *v=null;
        for(;u!=null;u=u->fa)
        {
            Splay(u);
            u->son[1]=v;
            v=u;
        }
        return v;
    }
    lct* GetRoot(lct* x) //得到根
    {
        for(x=Access(x);x->pushdown(),x->son[0]!=null;x=x->son[0]);
        return x;
    }
    void MakeRoot(lct* x) //使x成为根
    {
        Access(x)->rev^=1;
        Splay(x);
    }
    void Link(lct* x,lct* y) //连接两个点
    {
        MakeRoot(x);
        x->fa=y;
        Access(x);
    }
    void Cut(lct* x,lct* y) //断开两个点
    {
        MakeRoot(x);
        Access(y);
        Splay(y);
        y->son[0]->fa=null;
        y->son[0]=null;
    }
}A;
int N,M;
int main()
{
    scanf("%d%d",&N,&M);
    A.init(N);
    int x,y;
    char S[10];
    while(M--)
    {
        scanf("%s%d%d",S,&x,&y);
        if(S[0]=='Q')
        {
            lct *a=A.GetRoot(A.data+x);
            lct *b=A.GetRoot(A.data+y);
            if(a!=A.null&&a==b) printf("Yes\n");
            else printf("No\n");
        }
        else if(S[0]=='D') A.Cut(A.data+x,A.data+y);
        else A.Link(A.data+x,A.data+y);
    }
    return 0;
}