本文主要是介绍LCT入门笔记,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
LCT是动态树的一种,通过维护实链和虚链来维护所有路径之间的关系(类似于树链剖分)。这样做的目的是为了减少某些链上的修改、查询等操作的复杂度。虽然LCT常数巨大。
学LCT的大部分都会树剖吧?我们都知道树剖维护子树最大的儿子并形成一条重链,由于树剖是静态的,所以可以用线段树来维护。而由于LCT需要维护动态的边,要加边删边。所以需要用更灵活的数据结构来维护,也就是splay(也可以用非旋Treap,但是打的人貌似比较少)。
LCT将所有的边分为实边和虚边,对于所有的实边,都有互相连接的边。而对于虚边就只有儿子连到父亲的边。我们的splay维护所有的实边,而剩下的虚边则可以连接所有的splay。由于我们用splay维护,所以我们可以发现对于所有点x,它的儿子中只有一条实边。每一个Splay维护的是一条从上到下按在原树中深度严格递增的路径,且中序遍历Splay得到的每个点的深度序列严格递增。如图是一棵splay。接下去LCT的讲解就以它为例。(红色为实边)
LCT支持维护并修改链上信息,换根,删边连边,合并分离树等操作。
LCT中有关splay的操作:
void up(int k){ //向上更新数据(本题是求xor值)xyz[k]=xyz[son[k][0]]^xyz[son[k][1]]^val[k];
}
void revers(int k){ //区间翻转 swap(son[k][0],son[k][1]);rev[k]^=1;`}
void down(int k){ //下推标记if(rev[k]){revers(son[k][0]);revers(son[k][1]);rev[k]=0;`}
}
void rotate(int k){if(isroot(k)) return;int f=fa[k],gran=fa[f],d=(son[f][1]==k);if(!isroot(f)) son[gran][son[gran][1]==f]=k; //额外注意if(!isroot(f))语句,此处不判断会引起致命错误son[f][d]=son[k][d^1];if(son[k][d^1]) fa[son[k][d^1]]=f;son[k][d^1]=f;fa[f]=k;fa[k]=gran;up(f);up(k);
}
void splay(int k){int z=0,y=k,gran;sta[++z]=y;while(!isroot(y)) sta[++z]=y=fa[y];while(z) down(sta[z--]);while(!isroot(k)){y=fa[k];gran=fa[y];if(!isroot(y)) rotate((son[y][1]==k)^(son[gran][1]==y)?k:y);rotate(k);}up(k);
}
splay的操作类似于普通的splay,但是又有所不同,比如在进行rotate时我们需要判断点是否为当前splay的根。而在进行splay操作时,必须时一定要从上往下放标记。
Access操作:
access操作是LCT的精髓。access(x)的操作是将当前LCT的根到x的路径上的所有边变为实边,比如下图。
我们将从根到黑点的所有边变为实边,那么就必须将一些本来的实边变为虚边。
那我们对于每个splay都将x旋转到根,则x的连向右儿子的边变为虚边。并将x的右儿子变成上次access的点,然后上提x(将x变为其父亲),继续进行access操作即可。由于splay的性质,右儿子的深度大于点x,而由于我们已经确定了实边,剩下的边都是虚边,所以将右儿子变为虚边即可。(删去x的右儿子,但是fa[右儿子]=x),我们这样一直往上提,就可以完成access操作。
void access(int k){ //打通LCT根到k的路径 int y=0;do{splay(k);son[k][1]=y;up(k);y=k;k=fa[k];}while(k);
}
Makeroot操作:
makeroot(x)操作指的是将点x变为当前LCT的根,那我们先access(x),打通x到根的路径,然后再splay一下,则x就变成了当前splay的根,但是目前的splay是不满足性质的,因为(当前)深度最小的x在树根处,而所有的点都是x的左子树,而实际上所有的点都应该在x的右子树里(因为深度大于x),所以我们直接将整个区间reverse一下,翻转一遍就可以满足splay的性质了。
void makeroot(int k){ //将k变为LCT根 access(k);splay(k);revers(k);
}
Findroot操作:
findroot(x)表示x在的splay的树根。我们先打通x到根的链,然后把x旋转到根,由于我们要找的根深度最小(之前),所以直接一直寻找当前点的左儿子即可。注意在进行findroot时一定要下推标记,不然可能会错。
int findroot(int k){ //找k所在LCT的根 access(k);splay(k);while(son[k][0]) down(k),k=son[k][0];return k;
}
Split操作:
split(x,y)可以处理出x到y的路径。我们makeroot(x),然后再access(y),这样就形成了一个从x到y路径的splay,然后我们将y旋转到根,则y点上的值就是我们要求的路径的值。
void split(int x,int y){ //查询x到y的路径 makeroot(x);access(y);splay(y);
}
Link操作:
link(x,y)表示连边(x,y),我们连x和y之前要先判断一下,如果find(x)==find(y),即x和y在同一个splay里,在同一个联通块里,那么就直接return,不需要连边。否则我们先makeroot(x),然后将x的父亲变成y即可。
void link(int x,int y){ //连边x到y makeroot(x);if(findroot(y)==x) return;fa[x]=y;
}
Cut操作:
cut(x,y)表示删边(x,y)。我们考虑哪些情况不需要cut,首先如果find(x)!=find(y),肯定不需要cut,因为它们根本不在同一个联通块里。我们makeroot(x),那么如果fa[x]!=y那也可以不删,因为不存在x到y的边。那如果x与y有边就一定可以cut吗?不,如果x有右儿子就不能删,因为x有右儿子就会意味着x和y之间有别的边。
void cut(int x,int y){ //删边x到y makeroot(x);if(findroot(y)!=x||fa[x]!=y||son[x][1]) return; //顺序不能打反fa[x]=son[y][0]=0;up(y);
}
这些操作做完之后就完成了最基础的LCT操作,下面是LUOGU上LCT模板题的CODE
#include<bits/stdc++.h>
#define MAXN 300005
using namespace std;
int read(){char c;int x=0,y=1;while(c=getchar(),(c<'0'||c>'9')&&c!='-');if(c=='-') y=-1;else x=c-'0';while(c=getchar(),c>='0'&&c<='9')x=x*10+c-'0';return x*y;
}
int son[MAXN][2],fa[MAXN],val[MAXN],xyz[MAXN],rev[MAXN],sta[MAXN];
void up(int k){xyz[k]=xyz[son[k][0]]^xyz[son[k][1]]^val[k];
}
void revers(int k){ //区间翻转 swap(son[k][0],son[k][1]);rev[k]^=1;
}
void down(int k){if(rev[k]){revers(son[k][0]);revers(son[k][1]);rev[k]=0;}
}
int isroot(int k){ //判断k是否为其splay的根 return (son[fa[k]][0]!=k)&&(son[fa[k]][1]!=k);
}
void rotate(int k){if(isroot(k)) return;int f=fa[k],gran=fa[f],d=(son[f][1]==k);if(!isroot(f)) son[gran][son[gran][1]==f]=k;son[f][d]=son[k][d^1];if(son[k][d^1]) fa[son[k][d^1]]=f;son[k][d^1]=f;fa[f]=k;fa[k]=gran;up(f);up(k);
}
void splay(int k){int z=0,y=k,gran;sta[++z]=y;while(!isroot(y)) sta[++z]=y=fa[y];while(z) down(sta[z--]);while(!isroot(k)){y=fa[k];gran=fa[y];if(!isroot(y)) rotate((son[y][1]==k)^(son[gran][1]==y)?k:y);rotate(k);}up(k);
}
void access(int k){ //打通LCT根到k的路径 int y=0;do{splay(k);son[k][1]=y;up(k);y=k;k=fa[k];}while(k);
}
void makeroot(int k){ //将k变为LCT根 access(k);splay(k);revers(k);
}
int findroot(int k){ //找k所在LCT的根 access(k);splay(k);while(son[k][0]) down(k),k=son[k][0];return k;
}
void split(int x,int y){ //查询x到y的路径 makeroot(x);access(y);splay(y);
}
void link(int x,int y){ //连边x到y makeroot(x);if(findroot(y)==x) return;fa[x]=y;
}
void cut(int x,int y){ //删边x到y makeroot(x);if(findroot(y)!=x||fa[x]!=y||son[x][1]) return;fa[x]=son[y][0]=0;up(y);
}
int n,m;
int main()
{n=read();m=read();for(int i=1;i<=n;i++) val[i]=read();for(int i=1;i<=m;i++){int type=read(),x=read(),y=read();if(type==0){split(x,y);printf("%d\n",xyz[y]);}if(type==1) link(x,y);if(type==2) cut(x,y);if(type==3){splay(x);val[x]=y;}}return 0;
}
这篇关于LCT入门笔记的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!