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

相关文章

从入门到精通详解Python虚拟环境完全指南

《从入门到精通详解Python虚拟环境完全指南》Python虚拟环境是一个独立的Python运行环境,它允许你为不同的项目创建隔离的Python环境,下面小编就来和大家详细介绍一下吧... 目录什么是python虚拟环境一、使用venv创建和管理虚拟环境1.1 创建虚拟环境1.2 激活虚拟环境1.3 验证虚

Java List 使用举例(从入门到精通)

《JavaList使用举例(从入门到精通)》本文系统讲解JavaList,涵盖基础概念、核心特性、常用实现(如ArrayList、LinkedList)及性能对比,介绍创建、操作、遍历方法,结合实... 目录一、List 基础概念1.1 什么是 List?1.2 List 的核心特性1.3 List 家族成

Python学习笔记之getattr和hasattr用法示例详解

《Python学习笔记之getattr和hasattr用法示例详解》在Python中,hasattr()、getattr()和setattr()是一组内置函数,用于对对象的属性进行操作和查询,这篇文章... 目录1.getattr用法详解1.1 基本作用1.2 示例1.3 原理2.hasattr用法详解2.

c++日志库log4cplus快速入门小结

《c++日志库log4cplus快速入门小结》文章浏览阅读1.1w次,点赞9次,收藏44次。本文介绍Log4cplus,一种适用于C++的线程安全日志记录API,提供灵活的日志管理和配置控制。文章涵盖... 目录简介日志等级配置文件使用关于初始化使用示例总结参考资料简介log4j 用于Java,log4c

史上最全MybatisPlus从入门到精通

《史上最全MybatisPlus从入门到精通》MyBatis-Plus是MyBatis增强工具,简化开发并提升效率,支持自动映射表名/字段与实体类,提供条件构造器、多种查询方式(等值/范围/模糊/分页... 目录1.简介2.基础篇2.1.通用mapper接口操作2.2.通用service接口操作3.进阶篇3

Python自定义异常的全面指南(入门到实践)

《Python自定义异常的全面指南(入门到实践)》想象你正在开发一个银行系统,用户转账时余额不足,如果直接抛出ValueError,调用方很难区分是金额格式错误还是余额不足,这正是Python自定义异... 目录引言:为什么需要自定义异常一、异常基础:先搞懂python的异常体系1.1 异常是什么?1.2

Python实现Word转PDF全攻略(从入门到实战)

《Python实现Word转PDF全攻略(从入门到实战)》在数字化办公场景中,Word文档的跨平台兼容性始终是个难题,而PDF格式凭借所见即所得的特性,已成为文档分发和归档的标准格式,下面小编就来和大... 目录一、为什么需要python处理Word转PDF?二、主流转换方案对比三、五套实战方案详解方案1:

Spring WebClient从入门到精通

《SpringWebClient从入门到精通》本文详解SpringWebClient非阻塞响应式特性及优势,涵盖核心API、实战应用与性能优化,对比RestTemplate,为微服务通信提供高效解决... 目录一、WebClient 概述1.1 为什么选择 WebClient?1.2 WebClient 与

Spring Boot 与微服务入门实战详细总结

《SpringBoot与微服务入门实战详细总结》本文讲解SpringBoot框架的核心特性如快速构建、自动配置、零XML与微服务架构的定义、演进及优缺点,涵盖开发环境准备和HelloWorld实战... 目录一、Spring Boot 核心概述二、微服务架构详解1. 微服务的定义与演进2. 微服务的优缺点三

从入门到精通详解LangChain加载HTML内容的全攻略

《从入门到精通详解LangChain加载HTML内容的全攻略》这篇文章主要为大家详细介绍了如何用LangChain优雅地处理HTML内容,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录引言:当大语言模型遇见html一、HTML加载器为什么需要专门的HTML加载器核心加载器对比表二