关于 LCA (最近公共祖先)

news/2025/3/22 18:08:56/文章来源:https://www.cnblogs.com/bgf0212/p/18786918

对于一棵树上的两个点,他们的所有公共祖先中深度最大的那一个被称为它们的最近公共祖先(LCA)。求 LCA 有很多不同的方法。


倍增

倍增求 LCA,首先需要对树进行 dfs(废话),标记每个节点的直接父亲(\(2^0\) 级祖先)。然后我们就可以利用倍增的思想预处理它的第 \(2^1, 2^2, ... ,2^{\log n}\) 级祖先。然后就可以在线地求两个点(设为 \(x, y\),且有深度关系 \(d_x>d_y\))的 LCA。步骤如下:

  1. 判断 \(x\)\(y\) 是否相等,是则返回 \(x\)
  2. \(x\) 往上跳,直到与 \(y\) 的深度相同;
  3. 重复第 \(1\) 步;
  4. \(x\)\(y\) 往上跳,但是两个点不能相遇;
  5. 返回 \(fa_{x,0}\)

复杂度为 \(O(q \log n)\)。一般来说是够用的。

代码:

int lmt,d[N],f[N][50];
void init()
{for(int i=1;i<=n;i++)l[i]=l[i-1]+(1<<l[i-1]==i);for(int j=1;j<=lmt;j++)for(int i=1;i<=n;i++)f[i][j]=f[f[i][j-1]][j-1];
}
inline int lca(int x,int y)
{if(x==y)return x;if(d[x]<d[y])x^=y^=x^=y;for(int i=lmt;i>=0;i--)if(d[f[x][i]]>=d[y])x=f[x][i];if(x==y)return x;for(int i=lmt;i>=0;i--)if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];return f[x][0];
}

Tarjan

你学会了倍增求 LCA,然后有一天遇到了这道模板题:【模板】最近公共祖先(LCA)。你交了你最爱的倍增算法,但是发现 T 了若干个点。所以你去学习了(也许)更快的 Tarjan 求 LCA。

Tarjan 求 LCA 是一种离线算法。我们先存下来所有的询问(记得特判 \(x=y\))。然后我们对树进行一次 dfs。当一个点到达过但是未回溯时,给这个点打上标记 \(1\),回溯之后打上标记 \(2\),并且在并查集上将其与它的直接父亲连起来。如果我们当前在 \(y\) 点,需要查询的 \(x\) 点已经被标记为 \(2\),我们在并查集上查询 \(x\) 的祖先,就是 \(x,y\) 的 LCA。
image
如图,如果我们当前在 \(5\) 号点,需要查询的是 \(4\) 号点和它的 LCA,这个时候查询到 \(4\) 号点的祖先是 \(2\)\(2\) 号点就是它们的 LCA。反正就是这么求的,证明略。感性理解一下qwq
时间复杂度大概为 \(O(n+q)\),非常的够用。

代码:

struct query{int y,id,nxt;
}q[N*2];
int lq[N],leq=0;//存询问
inline void insertq(int x,int y,int i)
{q[++leq]={y,i,lq[x]};lq[x]=leq;
}
int vis[N],f[N];
int getf(int x)
{if(f[x]==x)return x;return f[x]=getf(f[x]);
}
int ans[N];
void tarjan(int x,int fa)
{vis[x]=1;for(int i=lk[x];i;i=e[i].nxt)if(e[i].y!=fa)tarjan(e[i].y,x);for(int i=lq[x];i;i=q[i].nxt)//遍历询问if(vis[q[i].y]==2||x==q[i].y)ans[q[i].id]=getf(q[i].y);vis[x]=2;f[x]=fa;//标记 2 并连上父亲
}

然后你通过了此题觉得万事大吉,以后再也不用倍增 LCA 了。


转化为 RMQ 问题

由于前一天问了人很好的 dalao 一道题,他对你说 kruskal 重构树,然而你并不知道那是什么,所以你刻苦学习,然后遇到了这道题:最小瓶颈路(加强版)。
注意到此题的 \(q\) 达到了惊人的 \(10^7\),感觉内存不是很妙。然后你找到了一种通过欧拉序把求 LCA 转化为 RMQ 问题的神奇做法。
欧拉序即在 dfs 时到达一个点和从儿子返回这个点时都将此点加入生成的序列。我们令 \(pos_x\) 表示 \(x\) 第一次在欧拉序中出现的位置, \(eun_i\) 表示欧拉序的第 \(i\) 位的点,\(d_x\) 依然表示 \(x\) 的深度。显然 \(x,y\) 的 LCA 是欧拉序上从 \(pos_x\) 走到 \(pos_y\) 经过的所有点中深度最小的,也是欧拉序最小的。然后就很自然地转化为了 RMQ 问题,使用 st 表即可 \(O(1)\) 单次查询。

代码:

int d[N*2],eun[N*4],pos[N*2],now=0;
bool vis[N*2];
void dfs(int x)
{vis[x]=true;eun[++now]=x;pos[x]=now;for(int i=lk[x],y;i;i=nxt[i]){d[y=to[i]]=d[x]+1;dfs(y);eun[++now]=x;}
}
int lg[N*4];
int st[30][N*4],rev[30][N*4];
void init()
{for(int i=2;i<=cnt*2;i++)lg[i]=lg[i>>1]+1;//欧拉序的长度为 2*cnt-1for(int i=1;i<=cnt;i++)//点数为 cnt{if(!vis[getf(i)])d[fa[i]]=1,dfs(fa[i]);}for(int i=1;i<=now;i++)st[0][i]=d[eun[i]],rev[0][i]=eun[i];int lmt=lg[now];for(int j=1;j<=lmt;j++)for(int i=1;i+(1<<j)<=now+1;i++)if(st[j-1][i]<st[j-1][i+(1<<(j-1))])st[j][i]=st[j-1][i],rev[j][i]=rev[j-1][i];elsest[j][i]=st[j-1][i+(1<<(j-1))],rev[j][i]=rev[j-1][i+(1<<(j-1))];
}
inline int lca(int x,int y)
{if(x==y)return 0;x=pos[x];y=pos[y];if(x>y)swap(x,y);int k=lg[y-x+1];return st[k][x]<st[k][y-(1<<k)+1]?rev[k][x]:rev[k][y-(1<<k)+1];
}

树链剖分也可以,但是懒得用,这里已经有非常好的 \(O(1)\) LCA 了(逃走)。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/903554.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

初入博客园

Start 第一篇随笔记录一下对于博客园的便利使用技巧。 自定义皮肤 安装和配置博客园皮肤全称按照了guangzan大佬的说明教程。[1] 里面附带有一有获取音乐链接以及歌词的网站,感觉很有用😍;[2] 贴一个网站头部背景图(苹果初代机-麦金塔Macintosh):待更新 🍀🍀🍀htt…

PCIe总线-RK3588 PCIe平台驱动分析

1.简介 RK3588 PCIe RC和EP使用同一个平台驱动,其主要的作用是解析设备树中的资源、初始化中断、使能电源、初始化PHY、使能时钟和释放复位,然后根据compatible属性初始化RC或者EP驱动。 2.入口 平台驱动的定义如下,当compatible属性为"rockchip,rk3588-pcie",则…

202413350081刁嘉怡博客园2

TASK11 #include <stdio.h>2 #include <stdlib.h>3 #include <time.h>4 5 #define N 56 7 int main() {8 int number;9 int i; 10 11 srand(time(0)); // 以当前系统时间作为随机种子 12 for(i = 0; i < N; ++i) { 13 numbe…

四则运算结对项目

这个作业属于哪个课程 https://edu.cnblogs.com/campus/gdgy/SoftwareEngineeringClassof2023这个作业的目标 结对合作生成一个四则运算题目生成器结对成员 3123004365彭颂华 &&3123004354黄子恒github地址 https://github.com/PShua/zuoye &&https://github.c…

Vue3 vite 集成 sass

一、安装依赖npm install sass-embedded二、配置全局变量 1、新建文件 src/styles/variables.scss 2、配置Vite 修改 vite.config.jscss: {preprocessorOptions: {scss: {additionalData: `@use "@/styles/variables" as *;`,},},},3、测试使用 a、在variables.scss…

在Android studio 里面想要写入文本内容怎么做

在 AndroidManifest.xml 文件中添加写入外部存储的权限。如果是 Android 6.0(API 级别 23)及以上的系统,还需要在运行时请求权限。在 MainActivity 的 writeToFile 方法里调用 FileHelper 类: private void writeToFile(String content) { FileHelper.writeToInternalStora…

iptables 介绍与实战

iptables是Linux内核中用于配置防火墙规则的工具。它基于Netfilter框架,可以对通过网络接口的数据包进行过滤、修改等操作。通过设置一系列规则,iptables能够控制哪些数据包可以进入或离开系统,从而实现网络安全防护等功能。它主要工作在网络层,能够根据数据包的源地址、目…

CF771E题解

CF771E题解很容易设出 \(dp_{i,j}\) 表示第一行选到 \(i\),第二行选到 \(j\) 的方案数 首先考虑部分分。 \(|a_i|\le 1\),那么产生贡献的一个矩阵不会超过 \(2\),那么就没必要考虑 \(|i-j|\ge 4\) 的状态了。证明如下: 不妨设 \(i<j\),那么我与其从 \(dp_{i,j}\to dp_{…

Word目录链接中只选中开头的文字

Word目录链接中只选中开头的文字 Word目录想选第一个字却选中整行 word自动生成目录里选中第一个字不是选整段就是选一行,怎么解决解决方法: 【不行】按上档键Shift+Home键。 【不行】Windows有个很强大的Alt健,按住Alt之后,在目录里面就可以想选哪个选哪个,不会出现再你这…

基于MPPT控制的锂离子电池充电器

基于MPPT的集成电路,可跟踪电源点并对电池充电,最大额定功率为1A,外形尺寸为TP4056。当谈到独立的单电池充电器时,只有一个流行的名字进入我的脑海,那就是我们的多功能锂离子/锂po TP4056电池充电器。这是广泛使用的,有很多功能与电池保护有关。锂电池在市场上很容易获得…

【Linux文件】把/etc/passwd删除了,该怎么办?

场景: 在做渗透测试发现网站存在任意文件删除漏洞,测试删除了/etc/passwd,那么删除后该如何恢复?一、 /etc/passwd文件的作用 /etc/passwd 是 Linux 系统中存储用户账户信息的关键文件,包含用户名、UID、GID、主目录路径和默认 Shell 等。删除 /etc/passwd 仍会导致以下问题…

20244210 实验一《Python程序设计》实验报告

20244210 2024-2025-2 《Python程序设计》实验一报告 课程:《Python程序设计》 班级:2442 姓名: 陈可 学号:20244210 实验教师:王志强 实验日期:2024年3月18日 必修/选修: 公选课 1.实验内容 1.熟悉Python开发环境; 2.练习Python运行、调试技能; 3.编写程序,练习变…