LCT入门笔记

2024-02-04 11:48
文章标签 入门 笔记 lct

本文主要是介绍LCT入门笔记,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

LCT是动态树的一种,通过维护实链和虚链来维护所有路径之间的关系(类似于树链剖分)。这样做的目的是为了减少某些链上的修改、查询等操作的复杂度。虽然LCT常数巨大。
学LCT的大部分都会树剖吧?我们都知道树剖维护子树最大的儿子并形成一条重链,由于树剖是静态的,所以可以用线段树来维护。而由于LCT需要维护动态的边,要加边删边。所以需要用更灵活的数据结构来维护,也就是splay(也可以用非旋Treap,但是打的人貌似比较少)。
LCT将所有的边分为实边和虚边,对于所有的实边,都有互相连接的边。而对于虚边就只有儿子连到父亲的边。我们的splay维护所有的实边,而剩下的虚边则可以连接所有的splay。由于我们用splay维护,所以我们可以发现对于所有点x,它的儿子中只有一条实边。每一个Splay维护的是一条从上到下按在原树中深度严格递增的路径,且中序遍历Splay得到的每个点的深度序列严格递增。如图是一棵splay。接下去LCT的讲解就以它为例。(红色为实边)

这里写图片描述LCT支持维护并修改链上信息,换根,删边连边,合并分离树等操作。


LCT中有关splay的操作:
void up(int k){  //向上更新数据(本题是求xor值)xyz[k]=xyz[son[k][0]]^xyz[son[k][1]]^val[k];
}
void revers(int k){  //区间翻转 swap(son[k][0],son[k][1]);rev[k]^=1;`}
void down(int k){  //下推标记if(rev[k]){revers(son[k][0]);revers(son[k][1]);rev[k]=0;`}
}
void rotate(int k){if(isroot(k)) return;int f=fa[k],gran=fa[f],d=(son[f][1]==k);if(!isroot(f)) son[gran][son[gran][1]==f]=k;  //额外注意if(!isroot(f))语句,此处不判断会引起致命错误son[f][d]=son[k][d^1];if(son[k][d^1]) fa[son[k][d^1]]=f;son[k][d^1]=f;fa[f]=k;fa[k]=gran;up(f);up(k);
}
void splay(int k){int z=0,y=k,gran;sta[++z]=y;while(!isroot(y)) sta[++z]=y=fa[y];while(z) down(sta[z--]);while(!isroot(k)){y=fa[k];gran=fa[y];if(!isroot(y)) rotate((son[y][1]==k)^(son[gran][1]==y)?k:y);rotate(k);}up(k);
}
splay的操作类似于普通的splay,但是又有所不同,比如在进行rotate时我们需要判断点是否为当前splay的根。而在进行splay操作时,必须时一定要从上往下放标记。

Access操作:
access操作是LCT的精髓。access(x)的操作是将当前LCT的根到x的路径上的所有边变为实边,比如下图。

这里写图片描述

我们将从根到黑点的所有边变为实边,那么就必须将一些本来的实边变为虚边。
那我们对于每个splay都将x旋转到根,则x的连向右儿子的边变为虚边。并将x的右儿子变成上次access的点,然后上提x(将x变为其父亲),继续进行access操作即可。由于splay的性质,右儿子的深度大于点x,而由于我们已经确定了实边,剩下的边都是虚边,所以将右儿子变为虚边即可。(删去x的右儿子,但是fa[右儿子]=x),我们这样一直往上提,就可以完成access操作。
void access(int k){  //打通LCT根到k的路径 int y=0;do{splay(k);son[k][1]=y;up(k);y=k;k=fa[k];}while(k);
}

Makeroot操作:
makeroot(x)操作指的是将点x变为当前LCT的根,那我们先access(x),打通x到根的路径,然后再splay一下,则x就变成了当前splay的根,但是目前的splay是不满足性质的,因为(当前)深度最小的x在树根处,而所有的点都是x的左子树,而实际上所有的点都应该在x的右子树里(因为深度大于x),所以我们直接将整个区间reverse一下,翻转一遍就可以满足splay的性质了。
void makeroot(int k){  //将k变为LCT根 access(k);splay(k);revers(k);
}

Findroot操作:
findroot(x)表示x在的splay的树根。我们先打通x到根的链,然后把x旋转到根,由于我们要找的根深度最小(之前),所以直接一直寻找当前点的左儿子即可。注意在进行findroot时一定要下推标记,不然可能会错。
int findroot(int k){ //找k所在LCT的根 access(k);splay(k);while(son[k][0]) down(k),k=son[k][0];return k;
}

Split操作:
split(x,y)可以处理出x到y的路径。我们makeroot(x),然后再access(y),这样就形成了一个从x到y路径的splay,然后我们将y旋转到根,则y点上的值就是我们要求的路径的值。
void split(int x,int y){  //查询x到y的路径 makeroot(x);access(y);splay(y);
}

Link操作:
link(x,y)表示连边(x,y),我们连x和y之前要先判断一下,如果find(x)==find(y),即x和y在同一个splay里,在同一个联通块里,那么就直接return,不需要连边。否则我们先makeroot(x),然后将x的父亲变成y即可。
void link(int x,int y){  //连边x到y makeroot(x);if(findroot(y)==x) return;fa[x]=y;
}

Cut操作:
cut(x,y)表示删边(x,y)。我们考虑哪些情况不需要cut,首先如果find(x)!=find(y),肯定不需要cut,因为它们根本不在同一个联通块里。我们makeroot(x),那么如果fa[x]!=y那也可以不删,因为不存在x到y的边。那如果x与y有边就一定可以cut吗?不,如果x有右儿子就不能删,因为x有右儿子就会意味着x和y之间有别的边。
void cut(int x,int y){   //删边x到y makeroot(x);if(findroot(y)!=x||fa[x]!=y||son[x][1]) return;  //顺序不能打反fa[x]=son[y][0]=0;up(y);
}

这些操作做完之后就完成了最基础的LCT操作,下面是LUOGU上LCT模板题的CODE
#include<bits/stdc++.h>
#define MAXN 300005
using namespace std;
int read(){char c;int x=0,y=1;while(c=getchar(),(c<'0'||c>'9')&&c!='-');if(c=='-') y=-1;else x=c-'0';while(c=getchar(),c>='0'&&c<='9')x=x*10+c-'0';return x*y;
}
int son[MAXN][2],fa[MAXN],val[MAXN],xyz[MAXN],rev[MAXN],sta[MAXN];
void up(int k){xyz[k]=xyz[son[k][0]]^xyz[son[k][1]]^val[k];
}
void revers(int k){  //区间翻转 swap(son[k][0],son[k][1]);rev[k]^=1;
}
void down(int k){if(rev[k]){revers(son[k][0]);revers(son[k][1]);rev[k]=0;}
}
int isroot(int k){  //判断k是否为其splay的根 return (son[fa[k]][0]!=k)&&(son[fa[k]][1]!=k);
}
void rotate(int k){if(isroot(k)) return;int f=fa[k],gran=fa[f],d=(son[f][1]==k);if(!isroot(f)) son[gran][son[gran][1]==f]=k;son[f][d]=son[k][d^1];if(son[k][d^1]) fa[son[k][d^1]]=f;son[k][d^1]=f;fa[f]=k;fa[k]=gran;up(f);up(k);
}
void splay(int k){int z=0,y=k,gran;sta[++z]=y;while(!isroot(y)) sta[++z]=y=fa[y];while(z) down(sta[z--]);while(!isroot(k)){y=fa[k];gran=fa[y];if(!isroot(y)) rotate((son[y][1]==k)^(son[gran][1]==y)?k:y);rotate(k);}up(k);
}
void access(int k){  //打通LCT根到k的路径 int y=0;do{splay(k);son[k][1]=y;up(k);y=k;k=fa[k];}while(k);
}
void makeroot(int k){  //将k变为LCT根 access(k);splay(k);revers(k);
}
int findroot(int k){ //找k所在LCT的根 access(k);splay(k);while(son[k][0]) down(k),k=son[k][0];return k;
}
void split(int x,int y){  //查询x到y的路径 makeroot(x);access(y);splay(y);
}
void link(int x,int y){  //连边x到y makeroot(x);if(findroot(y)==x) return;fa[x]=y;
}
void cut(int x,int y){   //删边x到y makeroot(x);if(findroot(y)!=x||fa[x]!=y||son[x][1]) return;fa[x]=son[y][0]=0;up(y);
}
int n,m;
int main()
{n=read();m=read();for(int i=1;i<=n;i++) val[i]=read();for(int i=1;i<=m;i++){int type=read(),x=read(),y=read();if(type==0){split(x,y);printf("%d\n",xyz[y]);}if(type==1) link(x,y);if(type==2) cut(x,y);if(type==3){splay(x);val[x]=y;}}return 0;
}

这篇关于LCT入门笔记的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/677387

相关文章

Spring Security 从入门到进阶系列教程

Spring Security 入门系列 《保护 Web 应用的安全》 《Spring-Security-入门(一):登录与退出》 《Spring-Security-入门(二):基于数据库验证》 《Spring-Security-入门(三):密码加密》 《Spring-Security-入门(四):自定义-Filter》 《Spring-Security-入门(五):在 Sprin

数论入门整理(updating)

一、gcd lcm 基础中的基础,一般用来处理计算第一步什么的,分数化简之类。 LL gcd(LL a, LL b) { return b ? gcd(b, a % b) : a; } <pre name="code" class="cpp">LL lcm(LL a, LL b){LL c = gcd(a, b);return a / c * b;} 例题:

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

Java 创建图形用户界面(GUI)入门指南(Swing库 JFrame 类)概述

概述 基本概念 Java Swing 的架构 Java Swing 是一个为 Java 设计的 GUI 工具包,是 JAVA 基础类的一部分,基于 Java AWT 构建,提供了一系列轻量级、可定制的图形用户界面(GUI)组件。 与 AWT 相比,Swing 提供了许多比 AWT 更好的屏幕显示元素,更加灵活和可定制,具有更好的跨平台性能。 组件和容器 Java Swing 提供了许多

【IPV6从入门到起飞】5-1 IPV6+Home Assistant(搭建基本环境)

【IPV6从入门到起飞】5-1 IPV6+Home Assistant #搭建基本环境 1 背景2 docker下载 hass3 创建容器4 浏览器访问 hass5 手机APP远程访问hass6 更多玩法 1 背景 既然电脑可以IPV6入站,手机流量可以访问IPV6网络的服务,为什么不在电脑搭建Home Assistant(hass),来控制你的设备呢?@智能家居 @万物互联

系统架构师考试学习笔记第三篇——架构设计高级知识(20)通信系统架构设计理论与实践

本章知识考点:         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识

poj 2104 and hdu 2665 划分树模板入门题

题意: 给一个数组n(1e5)个数,给一个范围(fr, to, k),求这个范围中第k大的数。 解析: 划分树入门。 bing神的模板。 坑爹的地方是把-l 看成了-1........ 一直re。 代码: poj 2104: #include <iostream>#include <cstdio>#include <cstdlib>#include <al

MySQL-CRUD入门1

文章目录 认识配置文件client节点mysql节点mysqld节点 数据的添加(Create)添加一行数据添加多行数据两种添加数据的效率对比 数据的查询(Retrieve)全列查询指定列查询查询中带有表达式关于字面量关于as重命名 临时表引入distinct去重order by 排序关于NULL 认识配置文件 在我们的MySQL服务安装好了之后, 会有一个配置文件, 也就

论文阅读笔记: Segment Anything

文章目录 Segment Anything摘要引言任务模型数据引擎数据集负责任的人工智能 Segment Anything Model图像编码器提示编码器mask解码器解决歧义损失和训练 Segment Anything 论文地址: https://arxiv.org/abs/2304.02643 代码地址:https://github.com/facebookresear

数学建模笔记—— 非线性规划

数学建模笔记—— 非线性规划 非线性规划1. 模型原理1.1 非线性规划的标准型1.2 非线性规划求解的Matlab函数 2. 典型例题3. matlab代码求解3.1 例1 一个简单示例3.2 例2 选址问题1. 第一问 线性规划2. 第二问 非线性规划 非线性规划 非线性规划是一种求解目标函数或约束条件中有一个或几个非线性函数的最优化问题的方法。运筹学的一个重要分支。2