1. 程式人生 > 實用技巧 >P2147 [SDOI2008]洞穴勘測(LCT維護樹的連通性)

P2147 [SDOI2008]洞穴勘測(LCT維護樹的連通性)

題目描述

輝輝熱衷於洞穴勘測。

某天,他按照地圖來到了一片被標記為JSZX的洞穴群地區。經過初步勘測,輝輝發現這片區域由n個洞穴(分別編號為1到n)以及若干通道組成,並且每條通道連線了恰好兩個洞穴。假如兩個洞穴可以通過一條或者多條通道按一定順序連線起來,那麼這兩個洞穴就是連通的,按順序連線在一起的這些通道則被稱之為這兩個洞穴之間的一條路徑。 洞穴都十分堅固無法破壞,然而通道不太穩定,時常因為外界影響而發生改變,比如,根據有關儀器的監測結果,123號洞穴和127號洞穴之間有時會出現一條通道,有時這條通道又會因為某種稀奇古怪的原因被毀。

輝輝有一臺監測儀器可以實時將通道的每一次改變狀況在輝輝手邊的終端機上顯示:

如果監測到洞穴u和洞穴v之間出現了一條通道,終端機上會顯示一條指令Connect u v

如果監測到洞穴u和洞穴v之間的通道被毀,終端機上會顯示一條指令Destroy u v

經過長期的艱苦卓絕的手工推算,輝輝發現一個奇怪的現象:無論通道怎麼改變,任意時刻任意兩個洞穴之間至多隻有一條路徑。

因而,輝輝堅信這是由於某種本質規律的支配導致的。因而,輝輝更加夜以繼日地堅守在終端機之前,試圖通過通道的改變情況來研究這條本質規律。 然而,終於有一天,輝輝在堆積成山的演算紙中崩潰了……他把終端機往地面一砸(終端機也足夠堅固無法破壞),轉而求助於你,說道:“你老兄把這程式寫寫吧”。

輝輝希望能隨時通過終端機發出指令Query u v

,向監測儀詢問此時洞穴u和洞穴v是否連通。現在你要為他編寫程式回答每一次詢問。 已知在第一條指令顯示之前,JSZX洞穴群中沒有任何通道存在。

輸入格式

第一行為兩個正整數n和m,分別表示洞穴的個數和終端機上出現過的指令的個數。 以下m行,依次表示終端機上出現的各條指令。每行開頭是一個表示指令種類的字串s("Connect”、”Destroy”或者”Query”,區分大小寫),之後有兩個整數u和v (1≤u, v≤n) 分別表示兩個洞穴的編號。

輸出格式

對每個Query指令,輸出洞穴u和洞穴v是否互相連通:是輸出”Yes”,否則輸出”No”。(不含雙引號)

#include<bits/stdc++.h>
using
namespace std; const int maxn=1e6+100; const int mod=51061; typedef long long ll; ll f[maxn]; ll ch[maxn][2]; ll v[maxn]; ll size[maxn]; ll lm[maxn]; ll la[maxn]; ll s[maxn]; ll st[maxn]; bool r[maxn]; bool nroot (int x) { return ch[f[x]][0]==x||ch[f[x]][1]==x; } void pushup (int x) { s[x]=s[ch[x][0]]+s[ch[x][1]]+v[x]; s[x]%=mod; size[x]=size[ch[x][0]]+size[ch[x][1]]+1; } void rev (int x) { //翻轉 int t=ch[x][0]; ch[x][0]=ch[x][1]; ch[x][1]=t; r[x]^=1; } void push_mul (int x,int c) { // s[x]*=c;s[x]%=mod; v[x]*=c;v[x]%=mod; lm[x]*=c;lm[x]%=mod; la[x]*=c;la[x]%=mod; } void push_add (int x,int c) { // s[x]+=c*size[x];s[x]%=mod; v[x]+=c;v[x]%=mod; la[x]+=c;la[x]%=mod; } void pushdown (int x) { if (lm[x]!=1) push_mul(ch[x][0],lm[x]),push_mul(ch[x][1],lm[x]),lm[x]=1; if (la[x]) push_add(ch[x][0],la[x]),push_add(ch[x][1],la[x]),la[x]=0; if (r[x]) { if (ch[x][0]) rev(ch[x][0]); if (ch[x][1]) rev(ch[x][1]); r[x]=0; } } void rotate (int x) { int y=f[x],z=f[y],k=ch[y][1]==x,w=ch[x][!k]; if (nroot(y)) ch[z][ch[z][1]==y]=x; ch[x][!k]=y; ch[y][k]=w; if (w) f[w]=y; f[y]=x;f[x]=z; pushup(y); } void splay (int x) { int y=x,z=0; st[++z]=y; while (nroot(y)) st[++z]=y=f[y]; while (z) pushdown(st[z--]); while (nroot(x)) { y=f[x];z=f[y]; if (nroot(y)) rotate((ch[y][0]==x)^(ch[z][0]==y)?x:y); rotate(x); } pushup(x); } void access (int x) { for (int y=0;x;x=f[y=x]) splay(x),ch[x][1]=y,pushup(x); } void makeroot(int x) { access(x); splay(x); rev(x); } void split (int x,int y) { makeroot(x); access(y); splay(y); } void link (int x,int y) { makeroot(x); f[x]=y; } void cut (int x,int y) { split(x,y); f[x]=ch[y][0]=0; } int find (int x) { //找樹根 access(x); splay(x); while (ch[x][0]) x=ch[x][0]; return x; } int n,q; int main () { scanf("%d%d",&n,&q); while (q--) { string op; cin>>op; if (op=="Connect") { int x,y; scanf("%d%d",&x,&y); link(x,y); } else if (op=="Destroy") { int x,y; scanf("%d%d",&x,&y); cut(x,y); } else if (op=="Query") { int x,y; scanf("%d%d",&x,&y); if (find(x)==find(y)) printf("Yes\n"); else printf("No\n"); } } }