LCT小结
前言
去年其实已经学过\(LCT\)了 ，但因为准备\(noip\)就没做什么题，忘得差不多啦，来补份总结
其实\(LCT\)是可以用\(FHQ\_treap\)实现的，但似乎更慢？？某\(dalao\)测试过的
反正\(LCT\)码量小(稍微压点就\(60\)行以下了)，而且网上大部分\(LCT\)都是用\(splay\)操作的，暂时不建议大家用\(FHQ\_treap\)
注：本文过于简短，只适合做总结观看
性质
\(LCT\)，全称\(Link Cut Tree\)，是一种维护动态树的利器
\(1.\)每个原树节点存在且仅存在于一棵\(splay\)
\(2.\)整个\(LCT\)，是由多棵\(splay\)构成的森林，其中每棵\(splay\)维护的是一条从上至下在原树中深度递增的路径，中序遍历后的点序列的在原树中的深度递增
\(3.LCT\)里的点靠实边和虚边连接，实边在\(splay\)中
虚边是由一棵\(splay\)指向另一棵\(splay\)的某个节点：该\(splay\)中的根指向原树深度最小的点的父亲(根与深度最小的点不等价)
当原树某点有多个儿子时，其中一个儿子存在于同棵\(splay\)且拉实边，其他儿子存在于其他的\(splay\)向该点拉虚边
\(Update(x)\)更新\(x\)某个要维护的值
inline void Update(LL x){
sum[x]=val[x]^sum[son[x][0]]^sum[son[x][1]];
}
\(Notroot(x)(\)不该棵\(splay\)为根\()Splay(x)(\)上旋到该棵\(splay\)的根\(),Rotate(x)\)
这三个函数不建议大家改动，老老实实打就好了，容易出错 且与标准\(splay\)有部分差别
inline bool Notroot(LL x){
return son[fa[x]][0]==x||son[fa[x]][1]==x;
}
inline void Rotate(LL x){
LL y=fa[x],z=fa[y],lz=(son[y][1]==x);;
if(Notroot(y)){
son[z][son[z][1]==y]=x;//**
}
son[y][lz]=son[x][lz^1],fa[son[y][lz]]=y;
son[x][lz^1]=y,fa[y]=x;
fa[x]=z;
Update(y),Update(x);
}
inline void Splay(LL x){
LL y(x),top(0);
sta[++top]=y;//**
while(Notroot(y)) sta[++top]=y=fa[y];//**
while(top) Pushdown(sta[top--]);//**
while(Notroot(x)){
y=fa[x];
if(Notroot(y)){
LL z(fa[y]);
if(((son[y][0]==x)^(son[z][0]==y))==0) Rotate(y);
else Rotate(x);
}
Rotate(x);
}
}
\(Access(x)\)拉链1(维护\(x\)到原树根)
每次把节点\(x\)上旋到根，因为我们是要维护这一条链，当然比\(x\)深度更大的点(右子树)丢到另外的\(splay\)里去
细节：并不等效于拉一条以实根为根的链，故查询链时需注意直接\(Split(x,y)\)再查询\(y\)
inline void Access(LL x){
for(LL y=0;x;y=x,x=fa[x]){
Splay(x),son[x][1]=y,Update(x);
}
}
比如我们\(Access(N)\)，看一下是怎么具体移动的
\(Makeroot(x)\)(换根)
通常我们维护的\((x,y)\)是位于原根的两个不同的子树，可是\(Access(x)\)只能维护\(x\)到原根这条链
一棵树随便提一点上来都还是树，那我们维护引入换根操作通过重构\(LCT\)，使得满足原有性质上换原根
比如我们维护\((x,y)\)这条链，就\(Makeroot(x),Access(y)\)
来看看\(Makeroot(x)\)的具体代码：当我们把\(x\)上旋到根时是没有右子树的(\(x\)深度最大)，而我们既然要换根得让\(x\)变成深度最小的，就交换左右子树
inline void Pushr(LL x){
swap(son[x][0],son[x][1]);r[x]^=1;
}
inline void Makeroot(LL x){
Access(x),Splay(x),Pushr(x);
}
\(Pushdown(x)\)下传标记
inline void Pushdown(LL x){
if(r[x]){
if(son[x][0]){Pushr(son[x][0]);}
if(son[x][1]){Pushr(son[x][1]);}
r[x]=0;
}
}
\(Findroot(x)\)找根
inline LL Findroot(LL x){
Access(x),Splay(x);
while(son[x][0]){
Pushdown(x),x=son[x][0];
}Splay(x);
return x;
}
\(Split(x,y)\)拉链2(维护\(x\)到\(y\))
inline void Split(LL x,LL y){
Makeroot(x),Access(y),Splay(y);
}
\(Link(x,y)\)动态连边
inline void Link(LL x,LL y){
Makeroot(x);
if(Findroot(y)!=x) fa[x]=y;
}
\(Cut(x,y)\)动态删边
inline void Delet(LL x,LL y){
Split(x,y);
if(son[y][0]!=x||son[x][1]) return;//注意这个地方
fa[x]=son[y][0]=0;
Update(y);
}
模板题
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<string>
#include<cmath>
#include<stack>
using namespace std;
typedef int LL;
const LL maxn=1e6;
inline LL Read(){
LL x(0),f(1);char c=getchar();
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9')x=(x<<3)+(x<<1)+c-'0',c=getchar();
return x*f;
}
LL n,m;
LL fa[maxn],son[maxn][2],val[maxn],sum[maxn],r[maxn];
inline bool Notroot(LL x){
return son[fa[x]][0]==x||son[fa[x]][1]==x;
}
inline void Pushr(LL x){
swap(son[x][0],son[x][1]);r[x]^=1;
}
inline void Pushdown(LL x){
if(r[x]){
if(son[x][0]){Pushr(son[x][0]);}
if(son[x][1]){Pushr(son[x][1]);}
r[x]=0;
}
}
inline void Update(LL x){
sum[x]=val[x]^sum[son[x][0]]^sum[son[x][1]];
}
inline void Rotate(LL x){
LL y=fa[x],z=fa[y],lz=(son[y][1]==x);;
if(Notroot(y)){
son[z][son[z][1]==y]=x;//**
}
son[y][lz]=son[x][lz^1],fa[son[y][lz]]=y;
son[x][lz^1]=y,fa[y]=x;
fa[x]=z;
Update(y),Update(x);
}
LL sta[maxn];
inline void Splay(LL x){
LL y(x),top(0);
sta[++top]=y;//**
while(Notroot(y)) sta[++top]=y=fa[y];//**
while(top) Pushdown(sta[top--]);//**
while(Notroot(x)){
y=fa[x];
if(Notroot(y)){
LL z(fa[y]);
if(((son[y][0]==x)^(son[z][0]==y))==0) Rotate(y);
else Rotate(x);
}
Rotate(x);
}
}
inline void Access(LL x){
for(LL y=0;x;y=x,x=fa[x]){
Splay(x),son[x][1]=y,Update(x);
}
}
inline void Makeroot(LL x){
Access(x),Splay(x),Pushr(x);
}
inline LL Findroot(LL x){
Access(x),Splay(x);
while(son[x][0]){
Pushdown(x),x=son[x][0];
}Splay(x);
return x;
}
inline void Split(LL x,LL y){
Makeroot(x),Access(y),Splay(y);
}
inline void Link(LL x,LL y){
Makeroot(x);
if(Findroot(y)!=x) fa[x]=y;
}
inline void Delet(LL x,LL y){
Split(x,y);
if(son[y][0]!=x||son[x][1]) return;
fa[x]=son[y][0]=0;
Update(y);
}
int main(){
n=Read(),m=Read();
for(LL i=1;i<=n;++i) val[i]=Read();
while(m--){
LL op(Read()),x(Read()),y(Read());
if(op==0){Split(x,y),printf("%d\n",sum[y]);}
else if(op==1){Link(x,y);}
else if(op==2){Delet(x,y);}
else{Splay(x),val[x]=y;Update(x);}
}
return 0;
}
