BZOJ 2049: [Sdoi2008]Cave 洞穴勘测 （动态树入门）

2049: [Sdoi2008]Cave 洞穴勘测

Time Limit: 10 Sec   Memory Limit: 259 MB
Submit: 1528   Solved: 644
[ Submit][ Status]

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

动态树

 

 

很简单的动态树操作，只有cut和link,判断在不在同一颗子树

 

  1 /* ***********************************************
  2 Author        :kuangbin
  3 Created Time  :2013-9-3 22:29:16
  4 File Name     :F:\2013ACM练习\专题学习\动态树-LCT\BZOJ2049洞穴勘测.cpp
  5 ************************************************ */
  6 
  7 #include <stdio.h>
  8 #include <string.h>
  9 #include <iostream>
 10 #include <algorithm>
 11 #include <vector>
 12 #include <queue>
 13 #include <set>
 14 #include <map>
 15 #include <string>
 16 #include <math.h>
 17 #include <stdlib.h>
 18 #include <time.h>
 19 using namespace std;
 20 
 21 //支持三种操作：
 22 //1.合并子树
 23 //2.分裂子树
 24 //3.询问x是否可达y
 25 //
 26 //Access(x) : 其实是把路径x到根打通，并且整条链用splay维护，这颗splay的
 27 //所有节点就是这条路径上的节点
 28 //link(u,v) : 把两颗子树合并为一颗树。先把u置为根，然后把u的父亲置为v.
 29 //cut(u,v)  : 把u和v分离成两颗子树。先把u置为根，然后把y splay到根，把
 30 //y 把左儿子与y分离
 31 const int MAXN = 10010;
 32 int ch[MAXN][2],pre[MAXN];
 33 int rev[MAXN];
 34 bool rt[MAXN];//rt[]表示该点是否为splay的根
 35 void Update_Rev(int r)
 36 {
 37     if(!r)return;
 38     swap(ch[r][0],ch[r][1]);
 39     rev[r] ^= 1;
 40 }
 41 void push_down(int r)
 42 {
 43     if(rev[r])
 44     {
 45         Update_Rev(ch[r][0]);
 46         Update_Rev(ch[r][1]);
 47         rev[r] = 0;
 48     }
 49 }
 50 void push_up(int r)
 51 {
 52     
 53 }
 54 void Rotate(int x)
 55 {
 56     int y = pre[x], kind = ch[y][1]==x;
 57     ch[y][kind] = ch[x][!kind];
 58     pre[ch[y][kind]] = y;
 59     pre[x] = pre[y];
 60     pre[y] = x;
 61     ch[x][!kind] = y;
 62     if(rt[y])
 63         rt[y] = false, rt[x] = true;
 64     else
 65         ch[pre[x]][ch[pre[x]][1]==y] = x;
 66     push_up(y);
 67 }
 68 void P(int r)
 69 {
 70     if(!rt[r])P(pre[r]);
 71     push_down(r);
 72 }
 73 void Splay(int r)
 74 {
 75     P(r);
 76     while( !rt[r] )
 77     {
 78         int f = pre[r], ff = pre[f];
 79         if(rt[f])
 80             Rotate(r);
 81         else if( (ch[ff][1]==f)==(ch[f][1]==f) )
 82             Rotate(f), Rotate(r);
 83         else
 84             Rotate(r), Rotate(r);
 85     }
 86     push_up(r);
 87 }
 88 //Access(x) :打通路径x到根节点
 89 int Access(int x)
 90 {
 91     int y = 0;
 92     do
 93     {
 94         Splay(x);
 95         rt[ch[x][1]] = true, rt[ch[x][1]=y] = false;
 96         push_up(x);
 97         x = pre[y=x];
 98     }
 99     while(x);
100     return y;
101 }
102 //判断是否同根(真实的树，非splay)
103 bool judge(int u,int v)
104 {
105     while(pre[u]) u = pre[u];
106     while(pre[v]) v = pre[v];
107     return u == v;
108 }
109 //使r成为所在树的根
110 void mroot(int r)
111 {
112     Access(r);
113     Splay(r);
114     Update_Rev(r);
115 }
116 void link(int u,int v)
117 {
118     if(judge(u,v))
119     {
120         puts("-1");
121         return;
122     }
123     mroot(u);
124     pre[u] = v;
125 }
126 void cut(int u,int v)
127 {
128     if(u == v || !judge(u,v))
129     {
130         puts("-1");
131         return;
132     }
133     mroot(u);
134     Splay(v);
135     pre[ch[v][0]] = pre[v];
136     pre[v] = 0;
137     rt[ch[v][0]] = true;
138     ch[v][0] = 0;
139     push_up(v);
140 }
141 int main()
142 {
143     //freopen("in.txt","r",stdin);
144     //freopen("out.txt","w",stdout);
145     int n,m;
146     while(scanf("%d%d",&n,&m) == 2)
147     {
148         for(int i = 0;i <= n;i++)
149         {
150             pre[i] = 0;
151             ch[i][0] = ch[i][1] = 0;
152             rev[i] = 0;
153             rt[i] = true;
154         }
155         char op[20];
156         int u,v;
157         while(m--)
158         {
159             scanf("%s%d%d",op,&u,&v);
160             if(op[0] == 'Q')
161             {
162                 if(judge(u,v))printf("Yes\n");
163                 else printf("No\n");
164             }
165             else if(op[0] == 'C')
166                 link(u,v);
167             else cut(u,v);
168         }
169     }
170     return 0;
171 }