本文原发布时间:\(\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. win10设置或更改硬盘图标

    1.首先要找到".ico"的图标素材,存放在要更改图标的硬盘根目录,可命名为a.ico.图片素材大家可以去一些素材网找找,一般可以找到很多. 2.在该硬盘分区空白处鼠标右键新建一个 ...

  2. C++11绑定器bind及function机制

    前言 之前在学muduo网络库时,看到陈硕以基于对象编程的方式,大量使用boost库中的bind和function机制,如今,这些概念都已引入至C++11,包含在头文件<functional&g ...

  3. JS逆向实战7-- 某省在线审批网站params 随机生成

    参数分析 我们首先通过抓包 发现这个就是我们所需要的数据 然后我们通过fidder 发起请求 结果: 通过我们反复测试 发现这个params的参数是每次请求中都会变化的 断点查找 我们通过 这个t参数 ...

  4. 小菜鸡的学习笔记---<正则表达式(1)>

    正则表达式学习笔记(1) (纯新手学习笔记,大佬绕路 QAQ) 一.简介 正则表达式就是一种文本模式用来匹配一系列满足特定条件的字符串,可以对比一下数学里面的表达式,比如我们要用一个表达式表示一串数字 ...

  5. 《ASP.NET Core技术内幕与项目实战》精简集-目录

    本系列是杨中科2022年最新作品<ASP.NET Core技术内幕与项目实战>及B站配套视频(强插点赞)的精简集,是一个读书笔记.总结和提炼了主要知识点,遵守代码优先原则,以利于快速复习和 ...

  6. 使用@Transactional注解的方法所在的类获取不到注解的解决方案

    前段时间遇到一个问题,一个service叫做A吧,有多个实现类分别是B,C,D,需要根据前端传的不同参数去匹配不同的实现类,我就自定义了一个注解@OrderDeal放在B,C,D上面,然后匹配前端传的 ...

  7. 质数之和【计算第x个到第y个质数之和】

    题目:质数之和 已知,第一个质数是2,第二个质数是3,第三个质数是5,第四个质数是7,第五个质数是11,第六个质数是13,第七个质数是17,输入两个不相等的正整数a和b,求出第a个质数到第b个质数当中 ...

  8. golang实现一个简单的http代理

    代理是网络中的一项重要的功能,其功能就是代理网络用户去取得网络信息.形象的说:它是网络信息的中转站,对于客户端来说,代理扮演的是服务器的角色,接收请求报文,返回响应报文:对于web服务器来说,代理扮演 ...

  9. 单例模式实现的多种方式、pickle序列化模块、选课系统需求分析等

    目录 单例模式实现的多种方式 方式一: 方式二: 方式三 方式四 pickle序列化模块 选课系统需求分析 功能提炼 选课系统架构设计 三层架构 选课系统目录搭建 选课系统功能搭建 单例模式实现的多种 ...

  10. Guess Next Session

    打开又是一个输入框的界面,点一下下面的看源码 很简短的一个源码 大概意思是如果password等于session[password]就输出flag 直接搜了下session函数的漏洞,发现sessio ...