bzoj 2049 [Sdoi2008]Cave 洞穴勘测 LCT

题目：

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

题意：

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

思路：

LCT 简单题。 LCT 跟树链剖分有一定的相似性，然后就是感觉更复杂了一点，要用 splay 维护

#include 

using namespace std;

const int N = 10000 + 10, INF = 0x3f3f3f3f;

int son[N][2], fat[N], siz[N], rev[N];
int top, stk[N];
bool is_root(int x)
{
    return son[fat[x]][0] != x && son[fat[x]][1] != x;
}
void push_down(int x)
{
    if(rev[x])
    {
        swap(son[x][0], son[x][1]);
        rev[son[x][0]] ^= 1, rev[son[x][1]] ^= 1;
        rev[x] ^= 1;
    }
}
void Rotate(int x)//进入Rotate后再决定何种旋转
{
    int y = fat[x], p = son[y][0] == x;
    son[y][!p] = son[x][p], fat[son[x][p]] = y;
    if(! is_root(y)) son[fat[y]][son[fat[y]][1]==y] = x;
    fat[x] = fat[y];
    son[x][p] = y, fat[y] = x;
}
void splay(int x)
{
    top = 0;
    stk[++top] = x;
    for(int i = x; !is_root(i); i = fat[i]) stk[++top] = fat[i];
    for(int i = top; i >= 1; i--) push_down(stk[i]);
    while(! is_root(x))
    {
        int y = fat[x], z = fat[y];
        if(is_root(y)) Rotate(x);
        else
        {
            if((x == son[y][0]) ^ (y == son[z][0])) Rotate(x), Rotate(x);
            else Rotate(y), Rotate(x);
        }
    }
}
void access(int x)
{
    int y = 0;
    while(x)
    {
        splay(x);
        son[x][1] = y;
        y = x, x = fat[x];
    }
}
void make_root(int x)
{
    access(x); splay(x);
    rev[x] ^= 1;
}
void cut(int x, int y)
{
    make_root(x);
    access(y); splay(y);
    son[y][0] = 0, fat[x] = 0;
}
void link(int x, int y)
{
    make_root(x); fat[x] = y; splay(x);
}
int find_root(int x)
{
    while(son[x][0]) x = son[x][0];
    return x;
}
//先把x作为原树的根，然后对y进行access操作，判断x有没有在y的左子树上的最左节点
bool query(int x, int y)
{
    make_root(x);
    access(y); splay(y);
    if(find_root(y) == x) return true;
    else return false;
}
//直接判断两个点是不是在同一个splay中
//int find_root(int x)
//{
//    access(x); splay(x);
//    while(son[x][0]) x = son[x][0];
//    return x;
//}
//bool query(int x, int y)
//{
//    return find_root(x) == find_root(y);
//}
int main()
{
    int n, m;
    scanf("%d%d", &n, &m);
    memset(fat, 0, sizeof fat);
    char opt[10];
    int x, y;
    for(int i = 1; i <= m; i++)
    {
        scanf(" %s%d%d", opt, &x, &y);
        if(opt[0] == 'C') link(x, y);
        else if(opt[0] == 'D') cut(x, y);
        else puts(query(x, y) ? "Yes" : "No");
    }
    return 0;
}