Lct浅谈
Lct浅谈
1.对lct的认识
首先要知道$lct$是什么。$lct$的全称为$link-cut-tree$。通过全称可以看出,这个数据结构是维护树上的问题,并且是可以支持连边断边操作。$lct$是一种动态树,可以处理动态森林问题,并支持查询与修改。
2.前置知识
在学习$lct$之前,建议先学习一下树链剖分,以及树链剖分套线段树,因为$lct$的思想与其相似,且树链剖分套线段树维护的是静态树,并且好实现。并且在$lct$的代码实现中我们要用到$splay$这个小小的知识点,所以建议学一下。
3.$lct$的讲解
我们先看树链剖分套线段树,这是先将整棵树构建好,我们直接在构建的树上进行树剖,并在树剖序上建出线段树,直接在线段树上乱搞就好了。我们想怎么能将静态树变成动态树呢?由于线段树是静态的,不能支持插入,所以我们需要改变我们所运用的数据结构,我们将线段树改变成为$splay$。
在$lct$中我们需要讲解八个操作。在讲解这些操作之前,我们需要知道一些名词以及其含义。
1) 辅助树与虚子树:我们在$lct$这个算法中要运用到$splay$,这个$splay$的运用之处与树剖中线段树的运用之处基本相似。$splay$是用来维护这棵树上的链,一个链就是一个$splay$。这些$splay$所构成的东西就叫辅助树。维护辅助树的关键字是当前节点的深度,也就是一个$splay$上的一个节点的右儿子的子树中的所有点的深度都要比当前点的深度小,反之都大。但是我们用$splay$维护了链,这些$splay$怎么连到一起?我们在维护$splay$的时候要记录当前节点的父亲,到整棵$splay$的根节点是没有的父亲的,这时我们将这个根的父亲指到整个链顶在原树中的父亲就好了,这是我们称原树中的链顶的父亲有一个虚子树,这棵虚子树就是整个链的$splay$。
知道上面两个名词之后,我们来讲解一下操作:
1) $Splay$:这个操作就是将当前节点旋转成为当前节点所在$splay$的根。这个操作要配合$rotate$来实现,这两个函数和正常的$splay$没有太大的区别,只是判根不一样,这里的判根是看当前节点的父亲的两个儿子是否都不等于自己,如果是则为根,反之不是。
bool check(int p) {return son[fa[p]][1]==p;}
bool isroot(int p) {return son[fa[p]][0]!=p&&son[fa[p]][1]!=p;}
void rotate(int p)
{
int tmp1=fa[p],tmp2=fa[tmp1],tmp3=check(p);
fa[son[tmp1][tmp3]=son[p][tmp3^1]]=tmp1,fa[p]=tmp2;
if(!isroot(tmp1)) son[tmp2][check(tmp1)]=p;
fa[son[p][tmp3^1]=tmp1]=p,pushup(tmp1),pushup(p);
}
void splay(int p)
{
update(p);
for(int i;i=fa[p],!isroot(p);rotate(p))
if(!isroot(i)) rotate(check(p)==check(i)?i:p);
}
2) $Access$:这个操作是将当前节点到当前节点所在辅助树的根的链变为一颗$splay$,代码不难理解。
void access(int p)
{
for(int t=0;p;t=p,p=fa[p]) splay(p),pushup(p);
}
3) $Makeroot$:这个操作是将当前节点变为当前节点所在辅助树的根。显然这个操作就是将当前节点先$Access$一下,再将当前节点$Splay$一下。由于希望保证时间复杂度,所以我们将当前
void makeroot(int p)
{
access(p),splay(p),swap(son[p][0],son[p][1]),rev[p]^=1;
}
4) $Update$:这个函数的功能在上面已经提到。
void pushdown(int p)
{
if(!rev[p]) return;
swap(son[son[p][0]][0],son[son[p][0]][1]);
swap(son[son[p][1]][0],son[son[p][1]][1]);
rev[son[p][0]]^=1,rev[son[p][1]]^=1,rev[p]=0;
}
void update(int p)
{
if(!isroot(p)) update(fa[p]);pushdown(p);
}
5) $Pushup$:这个操作就是用来统计答案的,类似于线段树上的$pushup$。
void pushup(int p)
{
if(p) size[p]=1+size[son[p][1]]+size[son[p][0]];
}
6) $Link$:这个操作就是连边,将$x,y$连接在一起。我们需要先$Makeroot(x)$,并$Makeroot(y)$。这样我们就可以直接将$fa[x]=y$,这样我们将以$x$为根的辅助树变为以$y$为根的辅助树的虚子树,这样就解决了。
void link(int x,int y)
{
makeroot(x),makeroot(y),fa[x]=y,ord[y]+=size[x];
}
7) $Cut$:这个操作就是断边,将$x,y$所连接的边断开。我们先$Makeroot(x)$,并$Access(y),Splay(y)$,这时$y$的左儿子就是$x$,所以只需要$son[y][0]=fa[x]=0$就好了。
void cut(int x,int y)
{
makeroot(x),access(y),splay(y),son[y][0]=fa[x]=0;
}
8) $Find$:这个操作是用来寻找当前节点所在辅助树中的根是谁。
int find(int x)
{
access(x),splay(x);
while((son[x][0]&&rev[x]==false)||(son[x][1]&&rev[x])) pushdown(x),x=son[x][0];
return x;
}
整合一下,就是下面的模板:
bool check(int p) {return son[fa[p]][1]==p;}
bool isroot(int p) {return son[fa[p]][0]!=p&&son[fa[p]][1]!=p;}
void pushdown(int p)
{
if(!rev[p]) return;
swap(son[son[p][0]][0],son[son[p][0]][1]);
swap(son[son[p][1]][0],son[son[p][1]][1]);
rev[son[p][0]]^=1,rev[son[p][1]]^=1,rev[p]=0;
}
void pushup(int p) {if(p) size[p]=1+size[son[p][1]]+size[son[p][0]];}
void update(int p) {if(!isroot(p)) update(fa[p]);pushdown(p);}
void rotate(int p)
{
int tmp1=fa[p],tmp2=fa[tmp1],tmp3=check(p);
fa[son[tmp1][tmp3]=son[p][tmp3^1]]=tmp1,fa[p]=tmp2;
if(!isroot(tmp1)) son[tmp2][check(tmp1)]=p;
fa[son[p][tmp3^1]=tmp1]=p,pushup(tmp1),pushup(p);
}
void splay(int p)
{
update(p);
for(int i;i=fa[p],!isroot(p);rotate(p))
if(!isroot(i)) rotate(check(p)==check(i)?i:p);
}
void access(int p) {for(int t=0;p;t=p,p=fa[p]) splay(p),son[p][1]=t,pushup(p);}
void makeroot(int p) {access(p),splay(p),swap(son[p][0],son[p][1]),rev[p]^=1;}
void link(int x,int y) {makeroot(x),makeroot(y),fa[x]=y,ord[y]+=size[x];}
void cut(int x,int y) {makeroot(x),access(y),splay(y),son[y][0]=fa[x]=0;}
Lct浅谈的更多相关文章
- 【转】Android Canvas的save(),saveLayer()和restore()浅谈
Android Canvas的save(),saveLayer()和restore()浅谈 时间:2014-12-04 19:35:22 阅读:1445 评论:0 收藏: ...
- [随机化算法] 听天由命?浅谈Simulate Anneal模拟退火算法
Simulate Anneal模拟退火算法,是一种用于得到最优解的随机化算法. 如果可以打一手漂亮的随机化搜索,也许当你面对一筹莫展的神仙题时就有一把趁手的兵器了. 这篇题解将教你什么?SA的基本思路 ...
- 浅谈 Fragment 生命周期
版权声明:本文为博主原创文章,未经博主允许不得转载. 微博:厉圣杰 源码:AndroidDemo/Fragment 文中如有纰漏,欢迎大家留言指出. Fragment 是在 Android 3.0 中 ...
- 浅谈 LayoutInflater
浅谈 LayoutInflater 版权声明:本文为博主原创文章,未经博主允许不得转载. 微博:厉圣杰 源码:AndroidDemo/View 文中如有纰漏,欢迎大家留言指出. 在 Android 的 ...
- 浅谈Java的throw与throws
转载:http://blog.csdn.net/luoweifu/article/details/10721543 我进行了一些加工,不是本人原创但比原博主要更完善~ 浅谈Java异常 以前虽然知道一 ...
- 浅谈SQL注入风险 - 一个Login拿下Server
前两天,带着学生们学习了简单的ASP.NET MVC,通过ADO.NET方式连接数据库,实现增删改查. 可能有一部分学生提前预习过,在我写登录SQL的时候,他们鄙视我说:“老师你这SQL有注入,随便都 ...
- 浅谈WebService的版本兼容性设计
在现在大型的项目或者软件开发中,一般都会有很多种终端, PC端比如Winform.WebForm,移动端,比如各种Native客户端(iOS, Android, WP),Html5等,我们要满足以上所 ...
- 浅谈angular2+ionic2
浅谈angular2+ionic2 前言: 不要用angular的语法去写angular2,有人说二者就像Java和JavaScript的区别. 1. 项目所用:angular2+ionic2 ...
- iOS开发之浅谈MVVM的架构设计与团队协作
今天写这篇博客是想达到抛砖引玉的作用,想与大家交流一下思想,相互学习,博文中有不足之处还望大家批评指正.本篇博客的内容沿袭以往博客的风格,也是以干货为主,偶尔扯扯咸蛋(哈哈~不好好工作又开始发表博客啦 ...
随机推荐
- USACO Section1.4 Arithmetic Progressions 解题报告
ariprog解题报告 —— icedream61 博客园(转载请注明出处)-------------------------------------------------------------- ...
- 每天一个Linux命令(3):ls命令
ls命令用来显示目标列表,在Linux中是使用率较高的命令.ls命令的输出信息可以进行彩色加亮显示,以分区不同类型的文件. 语法 ls(选项)(参数) 选项 -a:显示所有档案及目录(ls内定将档案名 ...
- Linux忘记root密码的解决办法
这里以centos6为例: 第一步:先将系统重新启动,在读秒的时候按下任意键就会出现如下图的菜单界面: 第二步:按下『e』就能够进入grub的编辑模式,如图: 第三步:将光标移动到kernel那一行, ...
- .net 匿名方法
匿名方法 核心就是lambda语法,下面是使用举例: var conn= MySqlHelper.GetConn(); var list=conn.Query<User>("SE ...
- OZ customize windows iamge
1.之前通过修改auto文件和windows.py文件可实现自动安装OS,但是disk bus只能以IDE模式启动,进入系统自动下载运行bat文件,安装cloudbase-init和qga,及清理日志 ...
- 蓝桥杯_风险度量_dfs_无向图两节点间的所有路径
标题:风险度量 X星系的的防卫体系包含 n 个空间站.这 n 个空间站间有 m 条通信链路,构成通信网.两个空间站间可能直接通信,也可能通过其它空间站中转. 对于两个站点x和y (x != y), 如 ...
- linux下的静态库和共享库
转载&&增加: 我们在编写一个C语言程序的时候,经常会遇到好多重复或常用的部分,如果每次都重新编写固然是可以的,不过那样会大大降低工作效率,并且影响代码的可读性,更不利于后期 ...
- shit vue-cli & path bug & baseUrl bug
vue-cli path bug https://cli.vuejs.org/zh/guide/#cli baseUrl bug baseUrl: "././" , https:/ ...
- [UOJ#131][BZOJ4199][NOI2015]品酒大会
[UOJ#131][BZOJ4199][NOI2015]品酒大会 试题描述 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品酒家”和“首席猎手”两个 ...
- Codeforces 498D Traffic Jams in the Land | 线段树
题目大意: 给坐标轴1~n的点,每个点有一个权值,从一个点走到下一个点需要1s,如果当前时间是权值的倍数就要多花1s 给出q组操作,C表示单点修改权值,A表示询问0时刻x出发到y的时间 题解:因为权值 ...