[学习笔记] Link-Cut-Tree

终于学了 $\text{link-cut-tree}$ ，我是看FlashHu 的博客学的。

概念简述

$\text{lct}$ 的思想主要是把连向某一个儿子的边划分为实边，其它的边划分成虚边。
区别在于虚实是可以动态变化的，因此要使用更高级、更灵活的Splay来维护每一条由若干实边连接而成的实链。

$\text{lct}$的主要性质如下：

• 每一个$\text{splay}$维护的是一条从上到下按在原树中深度严格递增的路径，且中序遍历$\text{splay}$得到的每个点的深度序列严格递增。
• 每个节点包含且仅包含于一个$\text{splay}$中。
• 边分为实边和虚边，实边包含在$\text{splay}$中，而虚边总是由一棵$\text{splay}$指向另一个节点。因为性质2，当某点在原树中有多个儿子时，只能向其中一个儿子拉一条实链，而其它儿子是不能在这个$\text{splay}$中的。那么为了保持树的形状，我们要让到其它儿子的边变为虚边，由对应儿子所属的$\text{splay}$的根节点的父亲指向该点，而从该点并不能直接访问该儿子。

具体操作

结合这一道模板题：【模板】Link Cut Tree，讲一下具体操作

access(x)

$\text{lct}$的核心操作，定义是把$x$到整棵树的根的路径全部改成实边，且$x$要成为这条路径上最深的点。

说白的实现步骤就四个，原谅蒟蒻搞不出图来 ：

• 把当前所指向的父亲转到根
• 修改儿子
• 更新信息
• 切换当前点所指的父亲

void access(int x)
{
	for(int y=0;x;x=par[y=x])
		splay(x),ch[x][1]=y,push_up(x); 
}

makeroot(x)

该操作的定义是把$x$换成大$\text{splay}$的根，我们把$x$到根的路径打通（$access$），然后把$\text{splay(x)}$到根，在$x$上面打上一个区间翻转标记（$flip$），就把深度最大的变成深度最小的，已达到换根的效果。

void flip(int x)
{
	if(!x) return ;
	swap(ch[x][0],ch[x][1]);
	fl[x]^=1;
}
void makeroot(int x)
{
	access(x);splay(x);
	flip(x);
}

findroot(x)

找到$x$在原树上的根，主要用来判断两点之间的连通性，先把$x$到根的路径打通，然后把$\text{splay(x)}$到根，找到$x$最左边的儿子，就是根，注意要把根$\text{splay}$回去，不要改变树的形态，代码如下：

int findroot(int x)
{
	access(x);splay(x);
	while(ch[x][0]) push_down(x),x=ch[x][0];
	splay(x);
	return x;
}

split(x,y)

这个操作就是把$(x,y)$的路径拆出来，我们先把$x$换成根，然后在把$y$到$x$的路径打通，把$y$转到根，那么$y$上面存的信息就是$(x,y)$这一条路径上的，代码如下：

void split(int x,int y)
{
	makeroot(x);
	access(y);splay(y);
}

link(x,y)

连一条$(x,y)$的边，要判断$x,y$的连通性，我们先把$x$给换成根，然后找$y$的$root$即可。

void link(int x,int y)
{
	makeroot(x);
	if(findroot(y)!=x) par[x]=y;
}

cut(x,y)

断一条$(x,y)$的边，要判断$x,y$的连通性，用$\text{link}$类似的方式，这时候要判断$(x,y)$直接相连，先判断$y$的父亲是不是$x$，在判断$y$是否有左儿子（因为$y$是$x$的右儿子，$y$如果有左儿子说明不直接相连）。

void cut(int x,int y)
{
	makeroot(x);
	if(findroot(y)==x && par[y]==x && !ch[y][0])
	{
		par[y]=ch[x][1]=0;
		push_up(x); 
	}
}

所有的操作至此就讲完了，下面附上完整代码。

#include <cstdio>
#include <iostream>
using namespace std;
const int M = 100005;
int read()
{
	int x=0,flag=1;char c;
	while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
	while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return x*flag;
}
int n,m,par[M],ch[M][2],val[M],sum[M],fl[M],st[M];
int nrt(int x)
{
	return ch[par[x]][0]==x || ch[par[x]][1]==x;
}
int chk(int x)
{
	return ch[par[x]][1]==x;
}
void push_up(int x)
{
	if(!x) return ;
	sum[x]=sum[ch[x][0]]^sum[ch[x][1]]^val[x];
}
void flip(int x)
{
	if(!x) return ;
	swap(ch[x][0],ch[x][1]);
	fl[x]^=1;
}
void push_down(int x)
{
	if(!x) return ;
	if(fl[x])
	{
		flip(ch[x][0]);flip(ch[x][1]);
		fl[x]=0;
	}
}
void rotate(int x)
{
	int y=par[x],z=par[y],k=chk(x),w=ch[x][k^1];
	ch[y][k]=w;par[w]=y;
	if(nrt(y)) ch[z][chk(y)]=x;par[x]=z;
	ch[x][k^1]=y;par[y]=x;
	push_up(y);push_up(x);
}
void splay(int x)
{
	int y=x,z=0;
	st[++z]=y;
	while(nrt(y)) st[++z]=y=par[y];
	while(z) push_down(st[z--]);
	while(nrt(x))
	{
		int y=par[x],z=par[y];
		if(nrt(y))
		{
			if(chk(y)==chk(x)) rotate(y);
			else rotate(x);
		}
		rotate(x);
	}
}
void access(int x)
{
	for(int y=0;x;x=par[y=x])
		splay(x),ch[x][1]=y,push_up(x); 
}
void makeroot(int x)
{
	access(x);splay(x);
	flip(x);
}
int findroot(int x)
{
	access(x);splay(x);
	while(ch[x][0]) push_down(x),x=ch[x][0];
	splay(x);
	return x;
}
void split(int x,int y)
{
	makeroot(x);
	access(y);splay(y);
}
void link(int x,int y)
{
	makeroot(x);
	if(findroot(y)!=x) par[x]=y;
}
void cut(int x,int y)
{
	makeroot(x);
	if(findroot(y)==x && par[y]==x && !ch[y][0])
	{
		par[y]=ch[x][1]=0;
		push_up(x); 
	}
}
int main()
{
	n=read();m=read();
	for(int i=1;i<=n;i++) val[i]=read();
	while(m--)
	{
		int op=read(),x=read(),y=read();
		if(op==0)
		{
			split(x,y);
			printf("%d\n",sum[y]);
		}
		if(op==1)
			link(x,y);
		if(op==2)
			cut(x,y);
		if(op==3)
			splay(x),val[x]=y;
	}
}

例题

• 例1：OTOCI，（版题，给luogu传送门）
• 例2：弹飞绵羊（有点意思，给题解）
• 例3：洞穴探测（版题，给luogu传送门）
• 例4：tree（打标记好题）
• 例5：城市旅行（需要推式子的标记题）
• 例6：水管局长（$lct$维护最小生成树）

日后会慢慢补充。
$$Tobecontinued.......$$