动态树Link-Cut Tree(LCT)

简介

一般的树上带修改问题,树剖(轻重链剖分)就可以搞定了,但是万一有连边、断边之类的东西?Good Game.
我们想一想,什么数据结构能处理移动点的问题?平衡树。
那我们就尝试把平衡树挪到树上,这时我们的树剖方式就叫做实链剖分。

定义与性质

先给定义。
实边:包含在一个splay中的边。一条实边恰好包含在一个splay中。一个节点的儿子只有一个连的是实边,其它都是虚边。
虚边:由一个splay的某个节点拉向另一个splay的根节点的边。
实链:一个splay里面从上到下连续的实边(感性理解一下,我说的不太清楚也不太准确)。
然后给出一些重要性质。
1、LCT由很多splay组成,这些splay之间由虚边连接。每个splay中序遍历得到的点,在原树中都按深度从小到大严格升序。
2、(此处的父亲,代表splay中的父子关系,并不是原树的)对于实边,父节点的儿子就是子节点,子节点的父亲就是父节点。但是对于虚边,子节点肯定是一个splay的根,父节点肯定是另一个splay的某个节点,所以通过父节点不能访问到子节点,但通过子节点能访问到父节点,这就是“认父不认子”。利用该性质可以判断一个点是不是splay的根。

操作

核心操作

access

这是核心中的战斗机,就是把某个点到原树根节点的路径全部拉成一条实链,且这条实链要以该点结尾。我看网上的讲解都写得很简洁,但不知道为什么scarlyw大佬的access函数写得忒复杂,还需要另外两个函数来铺垫,我这里就不提了,因为我没学会她的神奇写法(其实是我没听讲+我太菜了)。
代码

inline void access(int x) {for(register int y=0;x;y=x,x=ft[x]) splay(x),ch[x][1]=y,pushup(x);}

makeroot

把某个点钦定为原树的根节点,需要先access一下,再把这个splay从头到尾翻转一下,这样你会发现该splay中所有的点在原树中的深度都颠倒过来了,但是这条链旁侧的分支毫无影响,于是原先在原树最下方深度最大的钦定点就变成深度最小的原树根节点。是不是很巧妙啊?
代码

inline void makeroot(int x) {access(x),splay(x),pushrev(x);}

findroot

查找该点所在原树的根节点,只需access一下,把该点转为该splay的根,接着一直找左儿子,最后找到的就是原树根节点。
代码

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

其他操作

有了这个,现在很多树上操作都是小事一桩。

连、断边

首先要用findroot操作判断连、断边是否合法。当然,如果没有断边操作,可以并查集,常数更小。
连边只需把一个点makeroot一下,再接到另一个点下面。
断边只需把其中一个点makeroot一下,再把原树中的父亲的儿子和儿子的父亲设为0即可,要pushup。

inline void link(int x,int y) {makeroot(x);if(findroot(y)!=x) ft[x]=y;}
inline void cut(int x,int y) {makeroot(x);if(findroot(y)==x&&ft[y]==x&&!ch[y][0]) ft[y]=ch[x][1]=0,pushup(x);}

提取一条链

此处以查询一条链的点权异或值为例。我们一个点makeroot一下,另一个点access一下,这条链就提出来了。

inline int qxor(int x,int y) {makeroot(x),access(y),splay(y);return val[y];}

FBI Warning

这里我们用的splay与一般写法有点小小的差别。
首先,我们要自行写函数判断一个点是否为splay的根,在旋转函数中你会遇到。
接着,在splay时,要先从splay的根开始沿着这条链一直往下pushdown。为什么?也许是我们旋转时并未pushdown的缘故,但到底怎么回事我不知道,毕竟我没听scarlyw大佬讲课+我太菜了。

代码

究极压行,但fsy大佬压行比我的厉害多了

#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e5+50;
int n,m,a[MAXN],pstk[MAXN];
int ch[MAXN][2],val[MAXN],ft[MAXN],rev[MAXN];
inline int Read()
{
    int x=0,f=1;char ch;
    while(!isdigit(ch=getchar())) if(ch=='-') f=-1;
    while(isdigit(ch)) x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
    return x*f;
}
inline void pushup(int x) {val[x]=a[x]^val[ch[x][0]]^val[ch[x][1]];}
inline void pushrev(int x) {if(x) rev[x]^=1,swap(ch[x][0],ch[x][1]);}
inline void pushdown(int x) {if(rev[x]) pushrev(ch[x][0]),pushrev(ch[x][1]),rev[x]=0;}
inline int isroot(int x) {return ch[ft[x]][0]!=x&&ch[ft[x]][1]!=x;}
inline void rota(int x)
{
    int y=ft[x],z=ft[y],k=(x==ch[y][1]);
    ft[x]=z;if(!isroot(y)) ch[z][y==ch[z][1]]=x;
    ft[ch[x][k^1]]=y,ch[y][k]=ch[x][k^1];
    ft[y]=x,ch[x][k^1]=y;
    pushup(x),pushup(y);
}
inline void splay(int x)
{
    int top=0,y=x,z;
    pstk[++top]=y;
    while(!isroot(y)) pstk[++top]=y=ft[y];
    while(top) pushdown(pstk[top--]);
    while(!isroot(x))
    {
        y=ft[x],z=ft[y];
        if(!isroot(y)) (x==ch[y][1])^(y==ch[z][1])?rota(x):rota(y);
        rota(x);
    }
    pushup(x);
}
inline void access(int x) {for(register int y=0;x;y=x,x=ft[x]) splay(x),ch[x][1]=y,pushup(x);}
inline int findroot(int x) {access(x),splay(x);while(ch[x][0]) x=ch[x][0];splay(x);return x;}
inline void makeroot(int x) {access(x),splay(x),pushrev(x);}
inline void link(int x,int y) {makeroot(x);if(findroot(y)!=x) ft[x]=y;}
inline void cut(int x,int y) {makeroot(x);if(findroot(y)==x&&ft[y]==x&&!ch[y][0]) ft[y]=ch[x][1]=0,pushup(x);}
inline int qxor(int x,int y) {makeroot(x),access(y),splay(y);return val[y];}
int main()
{
    n=Read(),m=Read();
    for(register int i=1;i<=n;++i) a[i]=Read();
    int opt,x,y;
    for(register int i=1;i<=m;++i)
    {
        opt=Read(),x=Read(),y=Read();
        if(opt==0) printf("%d\n",qxor(x,y));
        else if(opt==1) link(x,y);
        else if(opt==2) cut(x,y);
        else splay(x),a[x]=y,pushup(x);
    }
    return 0;
}

原文地址:https://www.cnblogs.com/SKTT1Faker/p/12116315.html

时间: 2024-10-26 03:00:33

动态树Link-Cut Tree(LCT)的相关文章

LuoguP3690 【模板】Link Cut Tree (动态树) 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. 输入输出

P3690 【模板】Link Cut Tree (动态树)

P3690 [模板]Link Cut Tree (动态树) https://www.luogu.org/problemnew/show/P3690 分析: LCT模板 代码: 注意一下cut! 1 #include<cstdio> 2 #include<algorithm> 3 4 using namespace std; 5 6 const int N = 300100; 7 8 int val[N],fa[N],ch[N][2],rev[N],sum[N],st[N],top;

脑洞大开加偏执人格——可持久化treap版的Link Cut Tree

一直没有点动态树这个科技树,因为听说只能用Splay,用Treap的话多一个log.有一天脑洞大开,想到也许Treap也能从底向上Split.仔细思考了一下,发现翻转标记不好写,再仔细思考了一下,发现还是可以写的,只需要实时交换答案二元组里的两棵树,最后在吧提出来的访问节点放回去就行了.本着只学一种平衡树的想法,脑洞大开加偏执人格的开始写可持久化Treap版的Link Cut Tree... 写了才发现,常数硕大啊!!!代码超长啊!!!因为merge是从上到下,split从下到上,pushdow

动态树之link-cut tree

说好的专题... lct的一些概念看论文 杨哲<QTREE解法的一些研究> 简单易懂. 首先不要把lct想象得很难,其实很水的.lct就是很多splay树维护的树... lct的access操作就是在原树中拓展一条点到根的类二叉树出来(用splay来维护) 这里,splay树是按深度作为关键字的,当然,在无向图中(无环)可以任意指定一个点为根,这点要切记(因为在这里操作时,有些操作需要换根,所以一定要理解) link操作就是将点作为这颗类二叉树的根,然后合并 cut操作就是将点作为这颗二叉树的

Link Cut Tree学习笔记

从这里开始 动态树问题和Link Cut Tree 一些定义 access操作 换根操作 link和cut操作 时间复杂度证明 Link Cut Tree维护链上信息 Link Cut Tree维护子树信息 小结 动态树问题和Link Cut Tree 动态树问题是一类要求维护一个有根树森林,支持对树的分割, 合并等操作的问题. Link Cut Tree(林可砍树?简称LCT)是解决这一类问题的一种数据结构. 一些无聊的定义 Link Cut Tree维护的是动态森林中每棵树的任意链剖分. P

HDOJ 题目3966 Aragorn&#39;s Story(Link Cut Tree成段加减点权,查询点权)

Aragorn's Story Time Limit: 10000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 5505    Accepted Submission(s): 1441 Problem Description Our protagonist is the handsome human prince Aragorn comes from The Lor

AC日记——【模板】Link Cut Tree 洛谷 P3690

[模板]Link Cut Tree 思路: LCT模板: 代码: #include <bits/stdc++.h> using namespace std; #define maxn 300005 int n,m,val[maxn]; int top,ch[maxn][2],f[maxn],xr[maxn],q[maxn],rev[maxn]; inline void in(int &now) { int if_z=1;now=0; char Cget=getchar(); while

bzoj2049 [Sdoi2008]Cave 洞穴勘测 link cut tree入门

link cut tree入门题 首先说明本人只会写自底向上的数组版(都说了不写指针.不写自顶向下QAQ……) 突然发现link cut tree不难写... 说一下各个函数作用: bool isroot(int x):判断x是否为所在重链(splay)的根 void down(int x):下放各种标记 void rotate(int x):在x所在重链(splay)中将x旋转到fa[x]的位置上 void splay(int x):在x坐在重链(splay)中将x旋转到根 void acce

link cut tree 入门

鉴于最近写bzoj还有51nod都出现写不动的现象,决定学习一波厉害的算法/数据结构. link cut tree:研究popoqqq那个神ppt. bzoj1036:维护access操作就可以了. #include<cstdio> #include<cstring> #include<cctype> #include<algorithm> #include<queue> using namespace std; #define rep(i,s,

Link Cut Tree(无图慎入)

类似树链剖分(其实直接记住就可以了),提前放代码 1 #include<cstdio> 2 #include<cstdlib> 3 #include<iostream> 4 #include<algorithm> 5 #include<cstring> 6 #include<climits> 7 #include<cmath> 8 #define N (int)(3e5+5) 9 using namespace std;