本文原发布时间:\(\texttt{2022-05-21 14:11:52}\)

简介

最经公共祖先 \(\operatorname{LCA}(a,b)=c\),指的是在一棵树上节点 \(a\) 与 \(b\) 之间离两个点最近的一个点 \(c\),使得 \(c\) 是它们的祖先。

比如说,下面这棵树:

那么 \(\operatorname{LCA}(6,7)=2\),注意,不是 \(1\)。

例题(P3379 【模板】最近公共祖先

题目描述

如题,给定一棵有根多叉树,请求出指定两个点直接最近的公共祖先。

输入格式

第一行包含三个正整数 \(N,M,S\),分别表示树的结点个数、询问的个数和树根结点的序号。

接下来 \(N-1\) 行每行包含两个正整数 \(x, y\),表示 \(x\) 结点和 \(y\) 结点之间有一条直接连接的边(数据保证可以构成树)。

接下来 \(M\) 行每行包含两个正整数 \(a, b\),表示询问 \(a\) 结点和 \(b\) 结点的最近公共祖先。

输出格式

输出包含 \(M\) 行,每行包含一个正整数,依次为每一个询问的结果。

提示说明

对于 \(30\%\) 的数据,\(N\leq 10\),\(M\leq 10\)。

对于 \(70\%\) 的数据,\(N\leq 10000\),\(M\leq 10000\)。

对于 \(100\%\) 的数据,\(N\leq 500000\),\(M\leq 500000\)。

2021/10/4 数据更新 @fstqwq:应要求加了两组数据卡掉了暴力跳。

求法

0.暴力法

思路很简单。

首先预处理深度,询问时先统一深度,然后逐级上推即可。

最差时间复杂度 \(O(NM)\) 只能通过 \(70\%\) 的数据。

1.倍增法

倍增法本质上是对暴力法的优化。

首先,我们不必一个一个跳,可以指数级跳跃,效率大幅度提升。

具体如下:

首先,预处理节点 \(i\) 的第 \(2^{j}\) 级祖先 \(\operatorname{FA}[i][j]\),当然,还有深度。

询问时,先调整深度,然后倍增跳跃即可。

时间复杂度 \(O(N+M\log N)\),可以通过所有数据。

代码

我用的是链式前向星存图。

#include <bits/stdc++.h>
using namespace std; struct edge{
int nxt,to,w;
} tree[500005*2];
int head[500005],ec;
void add(int from,int to,int weight){
tree[++ec].nxt=head[from];
tree[ec].to=to;
tree[ec].w=weight;
head[from]=ec;
} int fa[500005][55],deep[500005],lg2[500005]; inline void initlog(int n){
lg2[1]=0;
lg2[2]=1;
for(int i=3;i<=n;i++){
lg2[i]=lg2[i>>1]+1;
}
} void dfs(int now,int parent){
fa[now][0]=parent;
deep[now]=deep[parent]+1;
for(int i=1;(1<<i)<=deep[now];i++){
fa[now][i]=fa[fa[now][i-1]][i-1];
}
for(int i=head[now];i;i=tree[i].nxt){
if(tree[i].to != parent){
dfs(tree[i].to,now);
}
}
} int lca(int x,int y){
if(deep[x]<deep[y]){
swap(x,y);
}
while(deep[x]!=deep[y]){
x=fa[x][lg2[deep[x]-deep[y]-1]];
}
if(x==y){
return x;
}
for(int k=lg2[deep[x]];k>=0;k--){
if(fa[x][k]!=fa[y][k]){
x=fa[x][k];
y=fa[y][k];
}
}
return fa[x][0];
} int n,m,s; int main(){
cin>>n>>m>>s;
for(int i=1,u,v;i<=n-1;i++){
cin>>u>>v;
add(u,v,114514);
add(v,u,1919810);
}
initlog(n);
dfs(s,0);
for(int i=1,arcka,akioi;i<=m;i++){
cin>>arcka>>akioi;
cout<<lca(arcka,akioi)<<endl;
}
return 0;
}

AC in Luogu

Tarjan

埋坑。

例题

(原创)Game2:简单树上问题

题目描述

给你一个 \(N\) 顶点的树,有 \(M\) 个询问,每次询问给出两个节点 \(u,v\) 求这两个点之间的最短距离。

输入格式

第一行为 \(N,M\)

接下来 \(N-1\) 行,每行三个整数 \(u,v,w\),表示一条从 \(u\) 到 \(v\) 的边,请忽略 \(w\)。

然后 \(M\) 行,每行两个整数 \(u,v\)。

输出格式

对于每一个询问,输出一个整数,表示答案。

提示

\(N,M = 100000,1 \le u,v,w \le N\)

时间限制 \(600\operatorname{ms}\),空间限制 \(50\operatorname{MB}\)。

造数据程序

from cyaron import *
n = 100000
m=n
for i in range(1, 25):
test = IO(file_prefix="data/test",data_id=i)
test.input_writeln(n,m)
graph=Graph.tree(n)
test.input_writeln(graph)
for i in range(m):
test.input_writeln(randint(1,n),randint(1,n))
test.output_gen("C:/Users/stu01/Documents/hh.exe")

思路

这道题可以直接用倍增LCA解决,因为:

\[\operatorname{Dis}(u,v)=\operatorname{Deep}(u)+\operatorname{Deep}(v)-2 \times \operatorname{Deep}(\operatorname{LCA}(u,v))
\]

代码

#include <bits/stdc++.h>
using namespace std; struct edge{
int nxt,to,w;
} tree[500005*2];
int head[500005],ec;
void add(int from,int to,int weight){
tree[++ec].nxt=head[from];
tree[ec].to=to;
tree[ec].w=weight;
head[from]=ec;
} int fa[500005][55],deep[500005],lg2[500005]; inline void initlog(int n){
lg2[1]=0;
lg2[2]=1;
for(int i=3;i<=n;i++){
lg2[i]=lg2[i>>1]+1;
}
} void dfs(int now,int parent){
fa[now][0]=parent;
deep[now]=deep[parent]+1;
for(int i=1;(1<<i)<=deep[now];i++){
fa[now][i]=fa[fa[now][i-1]][i-1];
}
for(int i=head[now];i;i=tree[i].nxt){
if(tree[i].to != parent){
dfs(tree[i].to,now);
}
}
} int lca(int x,int y){
if(deep[x]<deep[y]){
swap(x,y);
}
while(deep[x]!=deep[y]){
x=fa[x][lg2[deep[x]-deep[y]-1]];
}
if(x==y){
return x;
}
for(int k=lg2[deep[x]];k>=0;k--){
if(fa[x][k]!=fa[y][k]){
x=fa[x][k];
y=fa[y][k];
}
}
return fa[x][0];
} int n,m,s=1; int main(){
cin>>n>>m;
for(int i=1,u,v,w;i<=n-1;i++){
cin>>u>>v>>w;
add(u,v,114514);
add(v,u,1919810);
}
initlog(n);
dfs(s,0);
for(int i=1,arcka,akioi;i<=m;i++){
cin>>arcka>>akioi;
cout<<deep[arcka]+deep[akioi]-2*deep[lca(arcka,akioi)]<<endl;
}
return 0;
}

P4281 [AHOI2008]紧急集合 / 聚会 | LibreOJ10136. 「一本通 4.4 练习 3」聚会

简要题意

对于一个有根树,有 \(n\) 个节点,有 \(m\) 个询问,每次询问提供三个节点 \(x,y,z\),求一点 \(p\),使得 \(c=\operatorname{Dis}(x,p)+\operatorname{Dis}(y,p)+\operatorname{Dis}(z,p)\) 最小,输出 \(p\) 和 \(c\)。(\(\operatorname{Dis}(u,v)\) 表示 \(u\) 到 \(v\) 的最短路径长度)

对于 \(100\%\) 的数据,\(1\leq x,y,z\leq n\leq 5\times10^5\),\(1\leq m\leq 5\times 10^5\)。

思路

首先先来口胡一个结论,不重复的节点才是最优解。

然后就是算长度了,由于我们是三个点走到一个,那么就是这个:

\[\begin{aligned}
\operatorname{Dis}(x,y,z)=
\operatorname{Deep}(x)+
\operatorname{Deep}(y)+
\operatorname{Deep}(z)-\\
\operatorname{Deep}(\operatorname{LCA}(x,y))-
\operatorname{Deep}(\operatorname{LCA}(x,z))-
\operatorname{Deep}(\operatorname{LCA}(y,z))
\end{aligned}
\]

自己Hand推

时间复杂度 \(O(n+m\log n)\)。注意,本题需要卡常。

代码

#include <bits/stdc++.h>
using namespace std; struct edge{
int nxt,to,w;
} tree[500005*2];
int head[500005],ec;
void add(int from,int to,int weight){
tree[++ec].nxt=head[from];
tree[ec].to=to;
tree[ec].w=weight;
head[from]=ec;
} int fa[500005][55],deep[500005],lg2[500005]; inline void initlog(int n){
lg2[1]=0;
lg2[2]=1;
for(int i=3;i<=n;i++){
lg2[i]=lg2[i>>1]+1;
}
} void dfs(int now,int parent){
fa[now][0]=parent;
deep[now]=deep[parent]+1;
for(int i=1;(1<<i)<=deep[now];i++){
fa[now][i]=fa[fa[now][i-1]][i-1];
}
for(int i=head[now];i;i=tree[i].nxt){
if(tree[i].to != parent){
dfs(tree[i].to,now);
}
}
} int lca(int x,int y){
if(deep[x]<deep[y]){
swap(x,y);
}
while(deep[x]!=deep[y]){
x=fa[x][lg2[deep[x]-deep[y]-1]];
}
if(x==y){
return x;
}
for(int k=lg2[deep[x]];k>=0;k--){
if(fa[x][k]!=fa[y][k]){
x=fa[x][k];
y=fa[y][k];
}
}
return fa[x][0];
} int n,m; inline int cd(int a,int b,int LCA){
return deep[a]+deep[b]-2*deep[LCA];
} int main(){
ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
cin>>n>>m;
for(int i=1,u,v;i<=n-1;i++){
cin>>u>>v;
add(u,v,114514);
add(v,u,1919810);
}
initlog(n);
dfs(1,0);
while(m--){
int x,y,z;
cin>>x>>y>>z;
int l1=lca(x,y),l2=lca(x,z),l3=lca(y,z),ll=0;
if(l1==l2)ll=l3;
else if(l1==l3)ll=l2;
else ll=l1;
cout<<ll<<' ';
cout<<(deep[x]+deep[y]+deep[z]-deep[l1]-deep[l2]-deep[l3])<<'\n';
}
return 0;
}

AC in LibreOJ

AC in Luogu

LCA学习笔记(原洛谷文章)的更多相关文章

  1. 倍增求LCA学习笔记(洛谷 P3379 【模板】最近公共祖先(LCA))

    倍增求\(LCA\) 倍增基础 从字面意思理解,倍增就是"成倍增长". 一般地,此处的增长并非线性地翻倍,而是在预处理时处理长度为\(2^n(n\in \mathbb{N}^+)\ ...

  2. 2-SAT问题学习笔记+例题[洛谷P4792]

    一个不错的2-SAT文章:传送门 问题初入 什么是2-SAT SAT是适定性(Satisfiability)问题的简称 .一般形式为k-适定性问题,简称 k-SAT. 首先,把「2」和「SAT」拆开. ...

  3. dp凸优化/wqs二分学习笔记(洛谷4383 [八省联考2018]林克卡特树lct)

    qwq 安利一个凸优化讲的比较好的博客 https://www.cnblogs.com/Gloid/p/9433783.html 但是他的暴力部分略微有点问题 qwq 我还是详细的讲一下这个题+这个知 ...

  4. 2.《Spring学习笔记-MVC》系列文章,讲解返回json数据的文章共有3篇,分别为:

    转自:https://www.cnblogs.com/ssslinppp/p/4528892.html 个人认为,使用@ResponseBody方式来实现json数据的返回比较方便,推荐使用. 摘要 ...

  5. JavaEE精英进阶课学习笔记《博学谷》

    JavaEE精英进阶课学习笔记<博学谷> 第1章 亿可控系统分析与设计 学习目标 了解物联网应用领域及发展现状 能够说出亿可控的核心功能 能够画出亿可控的系统架构图 能够完成亿可控环境的准 ...

  6. AngularJs学习笔记--Guide教程系列文章索引

    在很久很久以前,一位前辈向我推荐AngularJs.但当时我没有好好学习,仅仅是讲文档浏览了一次.后来觉醒了……于是下定决心好好理解这系列的文档,并意译出来(英文水平不足……不能说是翻译,有些实在是看 ...

  7. 【算法学习】【洛谷】树链剖分 & P3384 【模板】树链剖分 P2146 软件包管理器

    刚学的好玩算法,AC2题,非常开心. 其实很早就有教过,以前以为很难就没有学,现在发现其实很简单也很有用. 更重要的是我很好调试,两题都是几乎一遍过的. 介绍树链剖分前,先确保已经学会以下基本技巧: ...

  8. 3.《Spring学习笔记-MVC》系列文章,讲解返回json数据的文章共有3篇,分别为:

    转自:https://www.cnblogs.com/ssslinppp/p/4528892.html 概述 在文章:<[Spring学习笔记-MVC-3]SpringMVC返回Json数据-方 ...

  9. 1.《Spring学习笔记-MVC》系列文章,讲解返回json数据的文章共有3篇,分别为:

    转自:https://www.cnblogs.com/ssslinppp/p/4528892.html [Spring学习笔记-MVC-3]SpringMVC返回Json数据-方式1:http://w ...

  10. 【算法学习】【洛谷】cdq分治 & P3810 三维偏序

    cdq是何许人也?请参看这篇:https://wenku.baidu.com/view/3b913556fd0a79563d1e7245.html. 在这篇论文中,cdq提出了对修改/询问型问题(Mo ...

随机推荐

  1. Linux实战笔记_CentOS7_无法识别NTFS格式的U盘

    注:因为CentOS 默认不识别NTFS的磁盘格式,所以我们要借助另外一个软件ntfs-3g来挂载.自带的yum源没有这个软件,要用第三方的软件源,比如阿里的epel. #安装ntfs-3g cd / ...

  2. 怎样在GitHub上建立仓库、以及怎样实现分支代码的合并。保姆级别的教程

    GitHub官网地址:https://github.com/ 注意:前提是已经注册了GitHub 文章目录 第一步:创建一个新的仓库 第二步.创建一个分支 第三步.编辑和发布更改的内容 第四步.拉取请 ...

  3. 齐博x1更新了 提供一个部分用户期待已久的功能,修改主题后变为待审

    如下图所示,你可以设置哪些用户组修改主题后,就会把原来已审核通过的主题,变为未审核.适合所有频道.

  4. 时序数据库TDengine 详细安装+集成流程+问题解决

    官方文档:https://docs.taosdata.com/get-started/package/ 点击进入 产品简介 TDengine 是一款高性能.分布式.支持 SQL 的时序数据库 (Dat ...

  5. Day2.1

    HelloWorld 随便新建一个文件夹,存放代码 新建一个java文件 文件后缀名为.java Hello.java 注意:系统可能没有显示文件名后缀,我们需要手动打开 编写代码 public cl ...

  6. 从 QFramework 重新开始

    自从上一篇博文写完后,就进入了繁忙工作状态,直到最近才开始有一点空闲时间. 这次繁忙的状态从 2022 年 10 月 11 日 持续到 2022 年 11 月 5 日. 主要的工作就是 QFramew ...

  7. golang基础语法学习

    1.函数作为一等公民 2.驼峰命名法/大小写决定是否在包外见 3.包名应该是小写的单个单词命名 4. 包名应为其源码的基础名称,如encoding/base64,包名应为base64而不是encodi ...

  8. Winform控件绑定数据

    目录 简介 绑定基类 功能扩展 简单控件绑定 列表控件绑定 绑定BindingList集合 绑定DataTable表格 绑定BindingSource源 表格控件绑定 绑定DataTable 绑定Bi ...

  9. 重新认识下JVM级别的本地缓存框架Guava Cache(2)——深入解读其容量限制与数据淘汰策略

    大家好,又见面了. 本文是笔者作为掘金技术社区签约作者的身份输出的缓存专栏系列内容,将会通过系列专题,讲清楚缓存的方方面面.如果感兴趣,欢迎关注以获取后续更新. 通过<重新认识下JVM级别的本地 ...

  10. Blazor Server完美实现Cookie Authorization and Authentication

    Blazor server-side application用Microsoft.AspNetCore.Identity.EntityFrameworkCore实现Authorization 和 Au ...