1. 程式人生 > >Link Cut Tree 動態樹 小結

Link Cut Tree 動態樹 小結

mes idf -- 一個 struct truct 最大 \n 樹鏈剖分

動態樹有些類似 樹鏈剖分+並查集 的思想,是用splay維護的

lct的根是動態的,"輕重鏈"也是動態的,所以並沒有真正的輕重鏈

動態樹的操作核心是把你要把 修改/詢問/... 等等一系列的操作的樹鏈放到一個splay裏,然後用splay根據相對深度大小來維護這個樹鏈

lct利用了splay的神奇性質,通過"認爹不認子"來達到記錄多個子樹的目的

lct的核心,access函數的意義是,在從這個點到它所在聯通塊中 相對深度最小的點 (可以理解為子樹根) 的樹鏈上,打通樹鏈上每個點所在的splay樹,並把需要的那部分樹鏈放到一個splay樹裏,這樣,這顆樹鏈上的所有點就都在同一個splay裏了,而且這顆splay裏只存了這個樹鏈而沒有其他的點

由於根是動態的,所以有了makeroot函數,意義是把x點作為樹根,先用access把需要的一部分樹鏈取出,然後splay到根,發現x點原來是這條樹鏈深度最大的點,現在深度變成了最小,用splay維護一下區間翻轉即可

還有很多細節需要註意,洛谷上神犇的講解很詳細

 1 #include <cstdio>
 2 #include <algorithm>
 3 #include <cstring>
 4 #define root d[0].ch[1]
 5 #define il inline
 6 #define nu 7777
 7 #define inf 500000
 8
#define N 300100 9 using namespace std; 10 11 int n,m,tp; 12 int stk[N]; 13 char str[10]; 14 struct Link_Cut_Tree{ 15 int fa[N],ch[N][2],sum[N],val[N],rev[N]; 16 il int idf(int x){return ch[fa[x]][0]==x?0:1;} 17 il int isroot(int x){return (ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x)?1:0;} 18 il void
con(int x,int ff,int p){fa[x]=ff;ch[ff][p]=x;} 19 il void pushup(int x){sum[x]=sum[ch[x][0]]^sum[ch[x][1]]^val[x];} 20 il void revers(int x){swap(ch[x][0],ch[x][1]),rev[x]^=1;} 21 il void pushdown(int x){ 22 if(rev[x]){ 23 if(ch[x][0]) revers(ch[x][0]); 24 if(ch[x][1]) revers(ch[x][1]); 25 rev[x]=0;}} 26 il void rot(int x){ 27 int y=fa[x];int ff=fa[y];int px=idf(x);int py=idf(y); 28 if(!isroot(y)) ch[ff][py]=x; 29 fa[x]=ff,con(ch[x][px^1],y,px),con(y,x,px^1); 30 pushup(y),pushup(x);} 31 void splay(int x){ 32 int y=x;stk[++tp]=x; 33 while(!isroot(y)){stk[++tp]=fa[y];y=fa[y];} 34 while(tp){pushdown(stk[tp--]);} 35 while(!isroot(x)){ 36 y=fa[x]; 37 if(isroot(y)) rot(x); 38 else if(idf(y)==idf(x)) rot(y),rot(x); 39 else rot(x),rot(x); 40 }} 41 il void access(int x){for(int y=0;x;y=x,x=fa[x])splay(x),ch[x][1]=y,pushup(x);} 42 il void mkroot(int x){access(x),splay(x),revers(x);} 43 il int findrt(int x){ 44 access(x),splay(x); 45 while(ch[x][0])pushdown(x),x=ch[x][0]; 46 return x;} 47 il void split(int x,int y){mkroot(x),access(y),splay(y);} 48 il void link(int x,int y){mkroot(x);if(findrt(y)!=x)fa[x]=y;} 49 il void cut(int x,int y){ 50 mkroot(x); 51 if(findrt(y)==x&&fa[x]==y&&!ch[x][1]) 52 fa[x]=ch[y][0]=0,pushup(y);} 53 }lct; 54 55 int gc() 56 { 57 int rett=0,fh=1;char c=getchar(); 58 while(c<0||c>9){if(c==-)fh=-1;c=getchar();} 59 while(c>=0&&c<=9){rett=(rett<<3)+(rett<<1)+c-0;c=getchar();} 60 return rett*fh; 61 } 62 63 int main() 64 { 65 n=gc(),m=gc(); 66 for(int i=1;i<=n;i++) lct.val[i]=gc(); 67 int fl,x,y; 68 for(int i=1;i<=m;i++) 69 { 70 fl=gc(),x=gc(),y=gc(); 71 if(fl==0) lct.split(x,y),printf("%d\n",lct.sum[y]); 72 if(fl==1) lct.link(x,y); 73 if(fl==2) lct.cut(x,y); 74 if(fl==3) lct.splay(x),lct.val[x]=y,lct.pushup(x); 75 } 76 return 0; 77 }

Link Cut Tree 動態樹 小結