动态树(LCT、Top Tree、ETT)
LCT
Upd:
一个细节:假如我们要修改某个节点的数据,那么要先把它makeroot再修改,改完之后pushup。
LCT是一种维护森林的数据结构,本质是用Splay维护实链剖分。
实链剖分大概是这样的:每个节点往一个儿子连实边,其它的儿子连虚边。
而我们用Splay维护实链剖分后的每一条实链。
因此LCT有一些基本的性质:
\(1.\)每一棵Splay维护树上一条直上直下的实链,且其中序遍历的点的序列的深度递增。
\(2.\)每个节点包含且仅包含于一棵Splay。
\(3.\)实边包含于Splay中,而虚边则连接两棵Splay。对于每个节点而言,它会记录它在Splay中的父亲,而play的根节点的父亲则是原树中这棵Splay代表的链的链顶的父亲。但是每个节点只会记录它在Splay中的左右儿子,并不会记录由虚边连接的儿子。(认父不认子)
下面如果没有特殊说明,我们默认我们所说的树为Splay而非原树。
然后我们先来说几个预备的基本操作:
nrooot
nroot实现判断一个节点是否为该节点所在Splay的根节点。
根据认父不认子的特性,我们只需要判断该节点的父亲是否有它这个儿子即可。
int nroot(int x){return ch[fa[x]][0]==x||ch[fa[x]][1]==x;}
pushrev
pushrev实现把一个子树翻转。
在后面的操作需要用到。
具体为交换左右儿子,并给左右儿子打上翻转标记。
void pushrev(int x){swap(lc,rc),rev[x]^=1;}
pushdown
pushdown实现下放标记。
在后面的操作中我们有翻转整棵子树的标记。
有的题目可能还会有其它的标记。
void pushdown(int x){if(!rev[x])return ;rev[x]=0,pushrev(lc),pushrev(rc);}
pushall
pushall实现把某个点到该点所在Splay根节点路径上的标记全部下放。
在后面的操作需要用到。
可以用栈实现不过函数堆栈更加方便。
void pushall(int x){if(nroot(x))pushall(fa[x]);pushdown(x);}
然后我们来分析一下Splay的基本操作:
这里我们需要支持的基本上就只有rotate和splay两个操作了。
rotate
和一般的Splay没有什么较大的差别。
判断一下当前节点的父亲是否为当前Splay的根,免得破坏认父不认子的特性。
void rotate(int x)
{
int y=fa[x],z=fa[y],k=ch[y][1]==x;
if(nroot(y)) ch[z][ch[z][1]==y]=x;
fa[x]=z,fa[y]=x,fa[ch[x][!k]]=y,ch[y][k]=ch[x][!k],ch[x][!k]=y,pushup(y);
}
splay
和一般的Splay没有什么较大的差别。不过我们现在只需要旋转到根节点了。
不过在splay之前需要pushall一下,典型的查询前pushdown。
void splay(int x)
{
pushall(x);
for(int y;nroot(x);rotate(x)) if(nroot(y=fa[x])) rotate((ch[fa[y]][0]==y)^(ch[y][0]==x)? y:x);
pushup(x);
}
然后就是LCT的基本操作了。
LCT所有的操作都依赖于一个核心操作:access。
access
access实现将原树中某个点到根的路径变为一条实链,单独拿出来做一棵Splay。
假设我们要拉的是\(x\)点。
首先我们把\(x\)Splay一下。
这时\(x\)的左子树中的点都会在这条实链上,而右子树的点都不在。所以我们把\(rc\)置为\(0\)。(注意认父不认子的特性)
然后我们跳到\(y=fa_x\)(根据上面的性质,我们跳到的实际上就是原树中\(x\)当前实链链顶的父亲。),把y\(Splay\)一下。
那么此时\(y\)的左子树中的点还是都在这条实链上,而右子树的点都不在。所以我们把刚才的\(x\)接在\(y\)的右儿子处。
注意因为修改后pushup的原则,我们需要在更新右儿子的时候pushup一下。
这样一直做到原树的根为止即可。
void access(int x){for(int y=0;x;x=fa[y=x])splay(x),rc=y,pushup(x);}
makeroot
相比access,makeroot更进一步地实现了把一个节点\(x\)转成原树中的根节点。
原理还是很简单的。我们先把\(x\)access一下,然后把\(x\)splay到它所在Splay的根节点。
此时因为它是原树中这条链上深度最大的点,所以它没有右儿子。
为了让它变成原树中深度最小的点,我们把它的左右子树交换一下,这样他就没有左儿子了,就变成了原树中深度最小的点。
交换子树可以通过打标记的方法来完成。
void makeroot(int x){access(x),splay(x),pushrev(x);}
findroot
findroot实现找到一个节点\(x\)所在原树中的根节点。
和makeroot类似,我们先access、splay\(x\)。
那么此时\(x\)所在原树中的根就是它左儿子中一直跳左儿子最后到的点。
根据查询前pushdown的原则我们要一边跳左儿子一边pushdown。
Upd:根据xzz的说法这里由于我们pushdown的写法可以不用pushdown。
最后可以顺便把查询到的根节点给splay一下。
int findroot(int x){access(x),splay(x);while(lc)x=lc;return splay(x),x;}
split
split实现把原树中一条链\((x,y)\)单独拿出来做一条实链。
很轻松地,我们先让\(x\)做原树的根,然后拉一条\((x,y)\)的实链出来,再让\(y\)做这条实链的Splay的根。这样我们在过程中通过pushup和pushdown就会让这条实链(这棵Splay)上的信息反映到这棵Splay的根节点\(y\)上。
然后如果我们要查询一条路径的信息,就可以快乐地split然后查\(y\)的信息了。
void split(int x,int y){makeroot(x),access(y),splay(y);}
下面则是一些真正意义上LCT能做而树剖做不了的东西了。
link
link实现在两点\(x,y\)之间新连一条边\((x,y)\)。(如果\(x,y\)连通就不连)
保证\(x,y\)不连通:
直接让\(x\)做\(x\)所在原树的根,然后把\(fa_x\)置为\(y\)即可。
void link(int x,int y){makeroot(x),fa[x]=y;}
不保证\(x,y\)不连通:
还要判一下连通性。即\(x\)成为原树所在的根之后,判断\(y\)所在原树的根是否为\(x\)。
void link(int x,int y){makeroot(x);if(findroot(y)^x)fa[x]=y;}
cut
cut实现断掉边\((x,y)\)。(如果不存在边\(x,y\)就不断)。
保证\((x,y)\)存在:
先把\((x,y)\)给split出来。
然后根据split的写法,此时\(y\)一定是该Splay的根,\(x\)一定是该原树的根,所以\(x\)一定是\(y\)的左儿子,那么我们直接断就完事了。
记得pushup。
void cut(int x,int y){split(x,y),fa[x]=ch[y][0]=0,pushup(y);}
不保证\((x,y)\)存在:
先判断连通性。
由于我们makeroot了,所以此时\(x\)一定为原树的根。
又由于我们findroot了,所以此时\(x\)一定是这条链的Splay的根。
那么如果此时\(y\)是\(x\)的右儿子,\(x\)是\(y\)的父亲,且\(y\)没有左儿子(\(x,y\)之间没有其它深度的点),那么\((x,y)\)就是存在的。
事实上如果满足了\(x\)是\(y\)的父亲,\(y\)没有左儿子,那么\((x,y)\)就是存在的。
因为\(x\)不可能有左儿子。
void cut(int x,int y){makeroot(x);if(findroot(y)==x&&fa[y]==x&&!ch[y][0])fa[y]=rc=0,pushup(x);}
这样我们就成功完成了LCT的基本操作。
然后就是一些技巧性的东西了:
LCT维护e-dcc:
删边似乎不太好做?只考虑加边(虽然树剖也能做)。
假如加入的边的两端不连通就直接加入,否则我们就会形成一个e-dcc,此时缩点即可。
LCT维护v-dcc:
圆方树?溜了溜了。
LCT维护边权:
一般而言我们最基本的思路是把边权放到其深度更大的那一端。
但是由于LCT会改变父子关系所以这东西没办法做了。
因此我么考虑拆点,把一条边拆成一个点,向其两端连边。
然后就可以做了。
LCT维护子树:
这是个大头。
如果会Top Tree的话似乎会轻松很多。
对于某些较容易维护的信息,我们可以考虑开一个新的数组来记录每个点的所有轻儿子所在子树的信息和,那么这个点的子树信息和就是两个重儿子的子树信息和加上所有轻儿子的子树信息和了。
(下面的\(sum\)表示该节点所有轻儿子的子树信息和,\(Sum\)表示该节点子树的信息和)
这样对于pushup操作,我们要多加一点东西。
void pushup(int x){Sum[x]=Sum[lc]+Sum[rc]+sum[x]+1;}
对于改变了轻重关系的操作,我们需要实时维护这个东西。比如access。
void access(int x){for(int y=0;x;x=fa[y=x])splay(x),sum[x]+=Sum[rc],sum[x]-=Sum[rc=y];}
link:
\(y\)多了一个轻儿子,所以要进行更改。
所以\(y\)在Splay中的祖先也要跟着更改,这太麻烦了。
我们把\(y\)旋转到它所在Splay的根,就不用了改\(y\)的祖先了。
保证\(x,y\)不连通:
void link(int x,int y){split(x,y),sum[fa[x]=y]+=Sum[x],pushup(y);}
这里的split是一个偷懒的写法。
不保证\(x,y\)不连通:
void link(int x,int y){makeroot(x);if(findroot(y)^x)sum[fa[x]=y]+=Sum[x],pushup(y);}
cut:
断掉的一定是重边,因为我们pushup了所以并没有特殊的影响。
LCT维护染色连通块:
可以直接开一个LCT,把同色的点连起来,不过很辣鸡。
更加优秀一点的是每个颜色开一个LCT。
还有一些诸如把点的颜色赋给其父边的技巧。
动态树(LCT、Top Tree、ETT)的更多相关文章
- hdu 5398 动态树LCT
GCD Tree Time Limit: 5000/2500 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Su ...
- hdu 5002 (动态树lct)
Tree Time Limit: 16000/8000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submi ...
- HDU 4718 The LCIS on the Tree (动态树LCT)
The LCIS on the Tree Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Oth ...
- [模板] 动态树/LCT
简介 LCT是一种数据结构, 可以维护树的动态加边, 删边, 维护链上信息(满足结合律), 单次操作时间复杂度 \(O(\log n)\).(不会证) 思想类似树链剖分, 因为splay可以换根, 用 ...
- 动态树LCT小结
最开始看动态树不知道找了多少资料,总感觉不能完全理解.但其实理解了就是那么一回事...动态树在某种意思上来说跟树链剖分很相似,都是为了解决序列问题,树链剖分由于树的形态是不变的,所以可以通过预处理节点 ...
- 动态树LCT(Link-cut-tree)总结+模板题+各种题目
一.理解LCT的工作原理 先看一道例题: 让你维护一棵给定的树,需要支持下面两种操作: Change x val: 令x点的点权变为val Query x y: 计算x,y之间的唯一的最短路径的点 ...
- bzoj2049-洞穴勘测(动态树lct模板题)
Description 辉辉热衷于洞穴勘测.某天,他按照地图来到了一片被标记为JSZX的洞穴群地区.经过初步勘测,辉辉发现这片区域由n个洞穴(分别编号为1到n)以及若干通道组成,并且每条通道连接了恰好 ...
- SPOJ OTOCI 动态树 LCT
SPOJ OTOCI 裸的动态树问题. 回顾一下我们对树的认识. 最初,它是一个连通的无向的无环的图,然后我们发现由一个根出发进行BFS 会出现层次分明的树状图形. 然后根据树的递归和层次性质,我们得 ...
- BZOJ 2002: [Hnoi2010]Bounce 弹飞绵羊 (动态树LCT)
2002: [Hnoi2010]Bounce 弹飞绵羊 Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 2843 Solved: 1519[Submi ...
- HDU 5002 Tree(动态树LCT)(2014 ACM/ICPC Asia Regional Anshan Online)
Problem Description You are given a tree with N nodes which are numbered by integers 1..N. Each node ...
随机推荐
- docker命令集锦
sudo docker image ls 查看有哪些image镜像sudo docker run hello 运行image 删除全部containerdocker rm $(docker conta ...
- ZAP-Queries【luogu3455】
题目大意 有不超过\(50000\)个询问,每次询问有多少正整数对\(x\),\(y\),满足\(x\leqslant a\),\(y \leqslant b\),并且\(gcd(x,y)=c\).其 ...
- Android_(自动化)获取手机存储卡的容量
手机上的存储卡是可以随时插拔的,每次插拔时会像操作系统总发送Action广播事件. 使用StatFs文件系统来获取MicroSD存储卡的剩余容量,在使用前先判断是否插入了存储卡,如果不存在则不于计算 ...
- From 7.15 To 7.21
From 7.15 To 7.21 大纲 竞赛 一周七天, 总共做三十五到四十道题吧, 题解要一起写了, 周六之前写不完的话就只能回家补咯 加上考试题总共做了25道题... 还要学计算几何, 生成函数 ...
- Zjoi2010排列计数Perm
这东西还是挺有思想的,道听途说一些东西,问问DuanYue同志,然后自己打表画树推了推,就搞出来了. 首先根据p i>p i/2(向下取整)这种形式,如果线段树学的好的人,一定能看出来,这是在唯 ...
- Vue双向绑定的实现原理系列(一):Object.defineproperty
了解Object.defineProperty() github源码 Object.defineProperty()方法直接在一个对象上定义一个新属性,或者修改一个已经存在的属性, 并返回这个对象. ...
- tomcat 配置域名部署war 项目
第一步把打包好的war包 放到 tomcat目录下的webapps 下 截图: 第二步:找到tomcat目录下的server.xml配置文件 server.xml在conf文件夹下面 编辑server ...
- LeetCode 21. 合并两个有序链表(Merge Two Sorted Lists)
题目描述 将两个有序链表合并为一个新的有序链表并返回.新链表是通过拼接给定的两个链表的所有节点组成的. 示例: 输入:1->2->4, 1->3->4 输出:1->1-& ...
- layui template list
//第一步:编写模版.你可以使用一个script标签存放模板,如: <script id="demo" type="text/html"> < ...
- Oracle 11g的日志路径
Oracle数据库的最常用问题定位日志是alert日志,Oracle数据库的日志文件alert_$ORACLE_SID.log记录了重作日志的转换,数据库启动和关闭,数据库结构的改变,回退段的修改,死 ...