求解 LCA の方法
最近公共祖先(LCA)
最近公共祖先简称 LCA(Lowest Common Ancestor)。两个节点的最近公共祖先,就是这两个点的公共祖先里面,离根最远的那个。
-----oi wiki
举个例子
在这张图中,\(5\) 和 \(9\) 的最近公共祖先就是 \(3\),\(9\) 和 \(7\) 的最近公共祖先就是 \(2\)
由于在树上两点间的简单路径是唯一的,并且从上图中我们可以很容易发现两点的 LCA 一定在其的简单路径上,所以我们有时候需要求出某两个点的 LCA
P3379 【模板】最近公共祖先(LCA)
以下所有的测试数据输入输出皆保持一致
暴力跳
最暴力的算法一个一个往上跳,复杂度 \(O(TLE)\)
没写,因为会 T
倍增求 LCA
这个方法是最简单也是最好想的一个方法,而且支持在线的查询,时间复杂度为 \(O(m\log n)\),本人看来是暴力的一个优化
我们维护一个 \(f[i][j]\) 数组表示 \(i\) 点往上跳 \(2^{j}\) 所到达的点的编号,\(dep[i]\) 表示 \(i\) 号点在树中的深度
我们首先要做的就是 DFS 一遍把我们的 \(f\),\(dep\) 数组求出来,然后在查询的时候,按照下面的步骤:
我们首先判断两个点的深度大小,保证传入的 \(x\) 的深度大于 \(y\) 的深度
我们让 \(x\) 点开始往上跳,从大到小枚举 \(f\) 数组的第二维,只要跳完之后深度大于等于 \(y\) 的深度就跳,这样最后我们使 \(x\) 和 \(y\) 来到了同一高度。
我们让 \(x\) 和 \(y\) 一起向上跳,前提是两点跳完后不在同一位置,因为跳的时候可能会跳超,但是两点的编号相同,造成是公共祖先但不是最近的情况,最后我们得到了两点的最近公共祖先下方的两个点 \(x'\) 和 \(y'\),这个时候直接返回其中一个的父节点即可。
注意第二步完成之后需要特判一下,是否两点编号相等,相等直接返回即可,防止两点在同一链上的情况
code:
#include<bits/stdc++.h>
#define N 500010
using namespace std;
struct sb{int t,nex;}e[N<<1];//t存放终点,nex存放下一个边
int head[N],tot;//head存放头节点,tot存放建边数量
int de[N],f[N][22],lg[N];//de存放每一个点的深度,f存放每一个点第2的i次方的祖先,lg存放log值
inline void add(int x,int y)//存边函数
{
e[++tot].t=y;//存放终点
e[tot].nex=head[x];//存放下一个边的头节点
head[x]=tot;//更新
}
inline void dfs(int now,int fa)//深搜函数
{
f[now][0]=fa;//当前节点上一个祖先就是fa
de[now]=de[fa]+1;//当前节点的深度就是祖先的深度加一
for(int i=1;i<=lg[de[now]];i++)//从2的1次方开始循环
f[now][i]=f[f[now][i-1]][i-1];//当前点的第2的i次方个祖先就是从当前的点的第i-1次方祖先的2的i-1次方的祖先
for(int i=head[now];i;i=e[i].nex)//寻找和当前点相连的点
if(e[i].t!=fa)dfs(e[i].t,now);//如果终点不是祖先的话就继续深搜
}
inline int LCA(int x,int y)//寻找最近公共祖先
{
if(de[x]<de[y])swap(x,y);//如果x的深度小于y的深度就把xy的值互换
while(de[x]>de[y])//只要x的深度大于y
x=f[x][lg[de[x]-de[y]]-1];//先让x跳到和y一样的深度
if(x==y)return y;//如果两个点相等说明y是x的祖先,直接返回y
for(int k=lg[de[x]];k>=0;k--)
{
if(f[x][k]==f[y][k])continue;//因为是跳到LCA下面的一层,所以他们一定不相等,
x=f[x][k],y=f[y][k];//替换xy的值继续循环
}
return f[x][0];//返回当前x的祖先即为最近公共祖先
}
int main()
{
int n,m,d;//n个点,m次询问,d为根节点
scanf("%d%d%d",&n,&m,&d);
for(int i=1;i<=n-1;i++)//从第一个到最后一个
{
int x,y;
scanf("%d%d",&x,&y);
add(x,y);//存边
add(y,x);//无向图存两次
}
for(int i=1;i<=n;i++)//预处理log的值
lg[i]=lg[i-1]+(1<<lg[i-1]==i);//计算
dfs(d,0);//深搜处理出f和de的值
for(int i=1;i<=m;i++)//开始询问
{
int x,y;
scanf("%d%d",&x,&y);
printf("%d\n",LCA(x,y));//输出最近公共祖先
}
return 0;//好习惯
}
tarjan 求 LCA
这个算法在 \(n,m\) 差不多的时候是比较快的,因为他复杂度是 \(O(n+m)\) 的,可惜的是只能离线。
我们对于每一个询问的两点之间连一条无向边(并不是在存树的数组里面连边)
在 DFS 的过程中,给走过的节点打上标记,同时维护并查集,如果 \(u\) 节点的这棵子树没搜完,那么 \(fa[u] = u\);,搜完后再更新并查集。
我们假设查询 \(u\) 和 \(v\) 的最近公共祖先,搜到节点 \(u\),如果另一个节点 \(v\) 已经被搜到过了,那么 \(v\) 点的并查集祖先就是 \(u\) 和 \(v\) 的最近公共祖先。
我们来模拟一下用 tarjan 求 LCA 的样例
我们在 DFS 的时候,首先来到 \(4\) 号点,然后我们 \(vis[4]=1\),将两个子节点都合并到 \(4\) 号点上,发现有一组询问 \(4,5\) 但是 \(vis[5]=0\) 所以不能确定 LCA
然后来到 \(2\),将 \(vis[2]=1\) 发现没法往下走,所以不走了,发现有 \((3,2)\) 和 \((1,2)\) 两组询问,但 \(vis[3]=0\),\(vis[1]=0\) 所以不能确定 LCA
然后来到 \(1\),将 \(vis[1]=1\),将子节点合并到 \(1\) 号点,发现有询问 \((1,2)\),这个时候 \(vis[2]=1\) 所以我们能确定 LCA,直接返回 \(fid(2)\),也就是 \(4\)
然后来到 \(3\),将 \(vis[3]=1\),没有子节点,发现有询问 \((3,2)(3,5)\),这时候 \(vis[2]=1,vis[5]=0\),所以能确定第一个询问的 LCA,记录 LCA\((3,2)=fid(2)=4\)
然后来到 \(5\),将 \(vis[5]=1\),发现有询问 \((3,5)(4,5)\),我们发现两个都可以确定,记录 LCA\((3,5)=fid(3)=1\), LCA\((4,5)=fid(4)=4\)
至此算法结束。
有的人可能会问,要是出现下面的情况怎么办?
要是查 \((1,2)\) 的话,不就会返回 \(fid(2)\) 递归下去不就是 \(6\) 了吗?
注意我们是在子树全遍历完才会修改当前点的 \(f\) 的值,所以不会出现这个情况
code:
#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define N 2000100
using namespace std;
int n,m,s,fa[N],vis[N],head[N],cnt,Head[N],cnt1=1;
struct sb{int u,v,next,lca;}e[N],e1[N];
inline void add(int u,int v){e[++cnt].v=v;e[cnt].next=head[u];head[u]=cnt;}
inline void Add(int u,int v){e1[++cnt1].v=v;e1[cnt1].next=Head[u];Head[u]=cnt1;}
inline int fid(int x){if(fa[x]==x)return x;return fa[x]=fid(fa[x]);}
inline int read(){int x=0,f=1;char ch=getchar();while(!isdigit(ch)){f=ch!='-';ch=getchar();}while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}return f?x:-x;}
inline void DFS(int x,int f)
{
vis[x]=1;
for(int i=head[x];i;i=e[i].next)
{
int v=e[i].v;
if(v==f)continue;
if(!vis[v])DFS(v,x),fa[v]=x;
}
for(int i=Head[x];i;i=e1[i].next)
{
int v=e1[i].v;
if(vis[v])e1[i].lca=e1[i^1].lca=fid(v);
}
}
signed main()
{
n=read(),m=read(),s=read();
for(int i=1;i<=n;i++)fa[i]=i;
for(int i=1;i<=n-1;i++)
{
int u=read(),v=read();
add(u,v);
add(v,u);
}
for(int i=1;i<=m;i++)
{
int u=read(),v=read();
Add(u,v);
Add(v,u);
}
DFS(s,0);
for(int i=2;i<=cnt1;i+=2)
cout<<e1[i].lca<<endl;
return 0;
}
/*
6 1 6
6 4
3 1
2 4
5 1
1 4
1 2
*/
树链剖分
以后我会重修的
首先我们就是要进行两遍 dfs 把 \(dep,f,siz,son,top\) 数组给处理出来,他们分别表示当前点的深度,当前点的父节点,以当前点为根的子树大小,当前点的重儿子编号,当前点所在的链的链头元素。
我们处理完之后,我们根据两点的深度大小来往上跳,和倍增一样,我们要防止跳超了,所以跳完以后两个点的 \(top\) 数组相等就说明已经有一个点的链头元素为两个点的 LCA,然后我们直接返回即可。
注意在跳的时候并不用保证谁比谁深度大,因为跳的时候是两个交替着跳,所以没必要
code:
#include <bits/stdc++.h>
#define N 1000100
#define endl '\n'
using namespace std;
struct sb{int u,v,next;}e[N];
int dep[N],fa[N],siz[N],head[N],son[N],rt,n,m,cnt,top[N];
inline void add(int u,int v){e[++cnt].u=u;e[cnt].v=v;e[cnt].next=head[u];head[u]=cnt;}
inline int read(){int x=0,f=1;char ch=getchar();while(!isdigit(ch)){f=ch!='-';ch=getchar();}while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}return f?x:-x;}
inline void dfs1(int x,int deep,int f)
{
dep[x]=deep;fa[x]=f;siz[x]=1;
int maxn=-1;
for(int i=head[x];i;i=e[i].next)
{
if(e[i].v!=f)
{
dfs1(e[i].v,deep+1,x);
siz[x]+=siz[e[i].v];
if(siz[e[i].v]>maxn)maxn=siz[e[i].v],son[x]=e[i].v;
}
}
return ;
}
inline void dfs2(int x,int tp)
{
top[x]=tp;
if(!son[x])return ;
dfs2(son[x],tp);
for(int i=head[x];i;i=e[i].next)
if(e[i].v!=fa[x]&&e[i].v!=son[x])dfs2(e[i].v,e[i].v);
return ;
}
inline int LCA(int x,int y)
{
while(top[x]!=top[y])
{
if(dep[top[x]]>=dep[top[y]])x=fa[top[x]];
else y=fa[top[y]];
}
return dep[x]<dep[y]?x:y;
}
signed main()
{
n=read();m=read(),rt=read();
for(int i=1;i<=n-1;i++)
{
int u=read(),v=read();
add(u,v);add(v,u);
}
dfs1(rt,0,rt);
dfs2(rt,rt);
for(int i=1;i<=m;i++)cout<<LCA(read(),read())<<endl;
return 0;
}
其他算法
由于本人实在太弱所以只整理了前面的三种
其实也可以用 DFS 序来求 LCA,可以去看看魏老师的博客:https://www.luogu.com.cn/blog/AlexWei/leng-men-ke-ji-dfs-xu-qiu-lca
也可以用欧拉序来求 LCA,比如:https://www.luogu.com.cn/blog/shenhy1205/solution-p3379
或者转化为 RMQ 问题求解:https://www.cnblogs.com/cj-xxz/p/11142232.html
或者用 LCT :https://www.luogu.com.cn/blog/surf/solution-p3379
求解 LCA の方法的更多相关文章
- 求解LCA问题的几种方式
求解LCA问题的几种方式 这篇随笔讲解图论中LCA问题(最近公共祖先)的几种求解方式及实现方法.LCA问题属于高级图论,所以希望读者学习过初级图论,知道图的一些基本知识,并懂得深搜算法的实现方式.这样 ...
- 浅谈倍增法求解LCA
Luogu P3379 最近公共祖先 原题展现 题目描述 如题,给定一棵有根多叉树,请求出指定两个点直接最近的公共祖先. 输入格式 第一行包含三个正整数 \(N,M,S\),分别表示树的结点个数.询问 ...
- 快速幂的求解-java方法(int范围之内)
思想就是,将十进制数化成二进制数.其它就是很简单了. 如:2的11次幂,11的二进制位1011,所以2(11) = 2(2(0) + 2(1) + 2(3)); 具体实现步骤,看代码比较简单 impo ...
- 最近公共祖先(LCA)的三种求解方法
转载来自:https://blog.andrewei.info/2015/10/08/e6-9c-80-e8-bf-91-e5-85-ac-e5-85-b1-e7-a5-96-e5-85-88lca- ...
- LCA—倍增法求解
LCA(Least Common Ancestors) 即最近公共祖先,是指在有根树中,找出某两个结点u和v最近的公共祖先. 常见解法一般有三种 这里讲解一种在线算法-倍增 首先我们定义fa[u][j ...
- (转载)LCA问题的Tarjan算法
转载自:Click Here LCA问题(Lowest Common Ancestors,最近公共祖先问题),是指给定一棵有根树T,给出若干个查询LCA(u, v)(通常查询数量较大),每次求树T中两 ...
- [CF 191C]Fools and Roads[LCA Tarjan算法][LCA 与 RMQ问题的转化][LCA ST算法]
参考: 1. 郭华阳 - 算法合集之<RMQ与LCA问题>. 讲得很清楚! 2. http://www.cnblogs.com/lazycal/archive/2012/08/11/263 ...
- 【算法】RMQ LCA 讲课杂记
4月4日,应学弟要求去了次学校给小同学们讲了一堂课,其实讲的挺内容挺杂的,但是目的是引出LCA算法. 现在整理一下当天讲课的主要内容: 开始并没有直接引出LCA问题,而是讲了RMQ(Range Min ...
- 关于LCA
LCA:最近公共祖先 指在有根树中,找出某两个结点u和v最近的公共祖先 如图,5,7的最近公共祖先就是3 接下来,我们来了解如何求解LCA No.1 暴力 首先想到的肯定是暴力,我们搜索,从两个节点一 ...
- Luogu 2680 NOIP 2015 运输计划(树链剖分,LCA,树状数组,树的重心,二分,差分)
Luogu 2680 NOIP 2015 运输计划(树链剖分,LCA,树状数组,树的重心,二分,差分) Description L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之 ...
随机推荐
- centos7 启动Tomcat7时报错:The APR based Apache Tomcat Native library which allows optimal performance in production environments was not found
INFO: The APR based Apache Tomcat Native library which allows optimal performance in production envi ...
- WSL安装Ubuntu 22.04 (2)
1.安装系统环境 1.1.设置软件包源为国内镜像 参考:Ubuntu更换国内镜像源 - 知乎 1.2.更新系统软件包 sudo apt-get update && sudo apt-g ...
- 使用Jmeter进行https接口测试时,如何导入证书?
转载:https://www.cnblogs.com/tester-zhangxiaona/p/12295473.html
- STM32的USART的DMA不定长度收发代码
/* * 函数名:USART1_DMA_Config * 描述 :串口1的DMA 的初始化配置 * 输入 :无 * 输出 : 无 * 调用 :在USART1_INIT()中被调用 */ voi ...
- windows下 mstsc 远程Ubuntu 图形界面2
采用 xrdp-0.9.2 + xorgxrdp-0.2.0 + xfce ,可完美实现远程,且: 1.支持中断后重新连接上一次连接. 2.支持clipbord跨平台复制. 3.xfce桌面,双机打开 ...
- 项目构建node-sass源码报错 SyntaxError:Unexpectedtoken"?"
背景 vue2项目,之前一直构建正常.今天改了代码,构建时报错,报错原因显示编译node-sass源码时出错. 报错信息: Modulebuild failed:/node_modules/node- ...
- 重学c#系列—— explicit、implicit与operator[三十四]
前言 我们都知道operator 可以对我们的操作符进行重写,那么explicit 和 implicit 就是对转换的重写. 正文 explicit 就是强制转换,然后implicit 就是隐式转换. ...
- Spring 注解整理
更多内容,前往IT-BLOG 一.核心容器 [1]@Configuration:告诉 Spring 这是一个配置类(配置类=配置文件)[2]@Bean:给容器中注册一个 Bean :类型为返回值类型, ...
- c++与linux详细计划,精确到每一天(仅80天)
好的,以下是三个月中每一天的学习计划: 第一个月: 第1天:阅读C++教程,熟悉环境,了解基本语法和数据类型 第2天:编写Hello World程序,加深对C++环境的了解 第3天:学习函数的定义和调 ...
- Oracle_用户-授权-角色
Oracle创建用户及表空间 1. 用户 创建用户: sql> create user <用户名> IDENTIFIED BY <用户密码> default tables ...