对于树上区间信息的维护
静态处理当然首选树链剖分
但如果一棵树是动态变化的,树剖便不再适用
于是我们就有了特别高级的——Link Cut Tree
LCT平均复杂度
O
(
l
o
g
N
)
O(logN)
O(logN),但常数贼大
一般跑的比树剖的
O
(
l
o
g
2
N
)
O(log^2N)
O(log2N)还慢
但树剖能做的题LCT 都能做,反过来就不是了
LCT概念与基础
LCT和树剖一样采用了轻重链剖分的思想
但既然LCT支持动态变化
线段树维护肯定不再适用,所以我们就使用更灵活的Splay来维护
(然而其实LCT比Splay早发明来着)
重儿子:这里的可字面理解为比较偏重的儿子吧,即需要先行处理的那一个,每个节点有且只有一个重儿子
重边:每个节点与其重儿子的连边称为重边,对每个节点这样的边有且只有一条
重链:由若干条连续的重边组成的路径
辅助树:
对于森林中每一条重链
我们都用一棵splay来维护,称为辅助树
每一个辅助树维护的是在原树中一条从上到下按深度严格递增的路径
即中序遍历辅助树得到的每个点的深度序列严格递增
LCT Access(x)
access操作是LCT算法的核心
表示打通一条从x到根节点的重路径
比如这样
代码实现其实十分简单
画图模拟一下大概就懂了~~(其实是我不懂怎么描述)~~
void access(int x)
{
int t=0;
while(x)
{
splay(x);
ch[x][1]=t;
update(x);
t=x; x=fa[x];
}
}
LCT Make-root(x)
表示将x结点作为原树的根进行重构
我们先打通从原根到x的重路径
这是x一定是其所在辅助树中深度最大的结点
而原根一定是与x在同一棵辅助树并且深度最小
换言之,在这棵splay中原根在序列最前,x在序列最后
那么我们要使x成为整棵树的根
就直接反转其所在辅助树表示的序列
即将x伸展到所在辅助树根
然后为x打上旋转标记
void mkrt(int x)
{
access(x); splay(x);
lzy[x]^=1;
}
LCT find(x)
表示找到当前x结点所在树的树根
打通从当前根到x的重路径
这时根一定在这个splay的最前面
所以直接把x伸展到其辅助树的根,然后找到x最左的儿子即可
int find(int x)
{
access(x); splay(x);
while(ch[x][0]) x=ch[x][0];
return x;
}
LCT Link(x,y)
表示连接森林中两节点x,y
先找到x,y两节点所在树的树根
如果树根相同直接返回
如果不相同就使x成为其所在树的根
然后将x 父指针指向y
void match()//听说link做函数名在linux下有时玄学CE???
{
int x=read(),y=read();
if(find(x)!=find(y))
mkrt(x),fa[x]=y;
}
LCT Cut(x,y)
表示切断当前x与y的连边
使x成为其所在树的根
这时y一定是x的某个儿子
要怎么确定是哪个呢
我们打通从x(当前根)到y的重路径
此时在辅助树中x一定是y的前驱
所以将y伸展到辅助树的根
这时x一定是y的左儿子
直接断开就好啦
断开后要记得更新信息
void cut()
{
int x=read(),y=read();
mkrt(x);
access(y); splay(y);
if(ch[y][0]!=x) return;
fa[x]=0; ch[y][0]=0;
update(y);
}
LCT isroot(x)
表示判断x是否是其所在辅助树的根
int isrt(int x)
{
return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;
}
LCT应用
洛谷P3690【模板】Link Cut Tree(动态树)
给定n个点以及每个点的权值,要你处理接下来的m个操作。操作有4种。操作从0到3编号。点从1到n编号。
0:后接两个整数(x,y),代表询问从x到y的路径上的点的权值的xor和。保证x到y是联通的。
1:后接两个整数(x,y),代表连接x到y,若x到y已经联通则无需连接。
2:后接两个整数(x,y),代表删除边(x,y),不保证边(x,y)存在。
3:后接两个整数(x,y),代表将点x上的权值变成y。
基本包含了上面所提到的所有操作
用
x
r
[
x
]
xr[x]
xr[x]维护以x为根的子树的异或和
修改操作直接
a
c
c
e
s
s
(
x
)
,
s
p
l
a
y
(
x
)
access(x),splay(x)
access(x),splay(x)
然后修改v[x],再更新xr[x]
查询x到y的异或和
先mkrt(x),再access(y),splay(y)
这时y是其所在辅助树的根
y所在辅助树表示的是x到y的路径经过的结点
因为y是所在辅助树的最后一个,所以直接用其左儿子的xr[]值异或v[y]
#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;
int read()
{
int f=1,x=0;
char ss=getchar();
while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
return f*x;
}
const int maxn=500010;
int n,m;
int v[maxn],xr[maxn];
int ch[maxn][2],fa[maxn];
int lzy[maxn];
int st[maxn];
void update(int p)
{
xr[p]=v[p]^xr[ch[p][0]]^xr[ch[p][1]];
}
void push(int p)
{
if(!lzy[p]) return;
swap(ch[p][0],ch[p][1]);
lzy[ch[p][0]]^=1;
lzy[ch[p][1]]^=1;
lzy[p]=0;
}
int isrt(int x)
{
return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;
}
void rotate(int x)
{
int y=fa[x],z=fa[y];
int t=(ch[y][0]==x);
if(!isrt(y))
{
if(ch[z][0]==y) ch[z][0]=x;
else ch[z][1]=x;
}
fa[y]=x; fa[ch[x][t]]=y; fa[x]=z;
ch[y][t^1]=ch[x][t]; ch[x][t]=y;
update(y); update(x);
}
void splay(int x)
{
int top=0; st[++top]=x;
for(int i=x;!isrt(i);i=fa[i])//要把splay中根到x路径都先下推
st[++top]=fa[i];
while(top) push(st[top--]);
while(!isrt(x))
{
int y=fa[x],z=fa[y];
if(!isrt(y))
{
if((ch[y][0]==x)^(ch[z][0]==y)) rotate(x);
else rotate(y);
}
rotate(x);
}
}
void access(int x)
{
int t=0;
while(x)
{
splay(x);
ch[x][1]=t;
update(x);
t=x; x=fa[x];
}
}
int find(int x)
{
access(x); splay(x);
while(ch[x][0]) x=ch[x][0];
return x;
}
void mkrt(int x)
{
access(x); splay(x);
lzy[x]^=1;
}
void match()
{
int x=read(),y=read();
if(find(x)!=find(y))
mkrt(x),fa[x]=y;
}
void cut()
{
int x=read(),y=read();
mkrt(x);
access(y); splay(y);
if(ch[y][0]!=x) return;
fa[x]=0; ch[y][0]=0;
update(y);
}
void change()
{
int x=read(),y=read();
access(x); splay(x);
v[x]=y;
xr[x]=xr[ch[x][0]]^v[x];
}
void get()
{
int x=read(),y=read();
mkrt(x);
access(y); splay(y);
printf("%d\n",xr[ch[y][0]]^v[y]);
}
int main()
{
n=read();m=read();
for(int i=1;i<=n;++i)
v[i]=read(),xr[i]=v[i];
while(m--)
{
int k=read();
if(k==0) get();
else if(k==1) match();
else if(k==2) cut();
else if(k==3) change();
}
return 0;
}
另一题洛谷P2147 [SDOI2008]洞穴勘测,也是模板题,比上面这个还简单点
其他应用
BZOJ2594 [Wc2006]水管局长数据加强版【LCT】题解
洛谷P4219 [BJOI2014]大融合【LCT】题解 LCT维护子树的应用
洛谷P4234 最小差值生成树【LCT】题解 LCT维护边权