Link Cut Tree 动态树 小结
Posted guapisolo
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Link Cut Tree 动态树 小结相关的知识,希望对你有一定的参考价值。
动态树有些类似 树链剖分+并查集 的思想,是用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 ",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 动态树 小结的主要内容,如果未能解决你的问题,请参考以下文章