洛谷 P3379 【模板】最近公共祖先(LCA) 做题反思(2024.3.21)

2024-03-21 22:44

本文主要是介绍洛谷 P3379 【模板】最近公共祖先(LCA) 做题反思(2024.3.21),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

【模板】最近公共祖先(LCA)

题目描述

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

输入格式

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

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

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

输出格式

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

样例 #1

样例输入 #1

5 5 4
3 1
2 4
5 1
1 4
2 4
3 2
3 5
1 2
4 5

样例输出 #1

4
4
1
4
4

提示

对于 30 % 30\% 30% 的数据, N ≤ 10 N\leq 10 N10 M ≤ 10 M\leq 10 M10

对于 70 % 70\% 70% 的数据, N ≤ 10000 N\leq 10000 N10000 M ≤ 10000 M\leq 10000 M10000

对于 100 % 100\% 100% 的数据, 1 ≤ N , M ≤ 500000 1 \leq N,M\leq 500000 1N,M500000 1 ≤ x , y , a , b ≤ N 1 \leq x, y,a ,b \leq N 1x,y,a,bN不保证 a ≠ b a \neq b a=b

样例说明:

该树结构如下:

第一次询问: 2 , 4 2, 4 2,4 的最近公共祖先,故为 4 4 4

第二次询问: 3 , 2 3, 2 3,2 的最近公共祖先,故为 4 4 4

第三次询问: 3 , 5 3, 5 3,5 的最近公共祖先,故为 1 1 1

第四次询问: 1 , 2 1, 2 1,2 的最近公共祖先,故为 4 4 4

第五次询问: 4 , 5 4, 5 4,5 的最近公共祖先,故为 4 4 4

故输出依次为 4 , 4 , 1 , 4 , 4 4, 4, 1, 4, 4 4,4,1,4,4

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

最近公共祖先(LCA)的概念
两个节点的最近公共祖先(LCA)就是这两个点的公共祖先里面,离他们最近的那个。
如图所示:
在这里插入图片描述
例如:
9和8的祖先节点有1、5、6,但是离9和8最近的公共祖先是节点6

倍增算法:是最经典的求LCA算法
∙ \bullet dep[u]是存u点的深度;
∙ \bullet fa[u][i]存从u点向上跳 2 i 2^i 2i层的祖先节点。(i=0,1,2,3,…);

例如:
从节点9向上跳 2 0 2^0 20层的祖先是3,跳 2 1 2^1 21层的祖先是6,跳 2 2 2^2 22层的祖先是1,跳 2 3 2^3 23层的祖先是0(即越界)。

1、dfs一遍,创建ST表
倍增递推,fa[u][i]=fa[fa[u][i-1]][i-1]
在这里插入图片描述
假设u是节点9,则层数增加两倍到倒数第四层,再除以2到倒数第二层,找到节点6,依次递归到节点1。

如表所示:

012319
100000
510000
251000
651000
365000
936100

2、利用ST表求LCA
(1)第一阶段,将u,v跳到同一层
设u,v两点的深度之差为y,将y进行二进制拆分,可以将y次游标跳跃优化为“y的二进制表示所有含1的个数”次游标跳跃,一定能跳到同一层。
例如,y=1019(0…01111111011)=512+256+…+8+2+1,不越界则跳,共跳9次到达。
(2)第二阶段,将u,v一起跳到LCA的下一层
从最大的i开始循环尝试,一直尝试到0,最后游标u,v一定能停在LCA的下一层。
两游标会跳512+256+…+8+2=1018层,共跳8次到达LCA的下一层。

代码示例如下:

#include<iostream>
#include<algorithm>
using namespace std;
const int N = 5e5 + 10;
int n, m, s;
int f[N][22], dep[N];
int e[2 * N], ne[2 * N], h[N], idx;
void add(int a, int b) {e[++idx] = b, ne[idx] = h[a], h[a] = idx;
}
void dfs(int u, int fa) {dep[u] = dep[fa] + 1;//计算节点u的深度f[u][0] = fa;//找节点u的父节点for (int i = 1; i <= 20; i++) {//这里从1循环到20的原因是因为n的最大范围是5e5+10,所以2的20次方是一百万,远远大于5e5+10,这样保证数据不会溢出f[u][i] = f[f[u][i - 1]][i - 1];//找到所有的公共祖先}for (int i = h[u]; i; i = ne[i]) {int j = e[i];if (j != fa)dfs(j, u);//如果没有找到父节点就继续搜}
}
int lca(int u, int v) {if (dep[u] < dep[v])swap(u, v);//先跳到同一层for (int i = 20; i >= 0; i--) {if (dep[f[u][i]] >= dep[v]) {//将节点u继续往上跳u = f[u][i];}}if (u == v)return v;//找到最近公共祖先for (int i = 20; i >= 0; i--) {if (f[u][i] != f[v][i]) {//如果没有找到最近公共祖先u = f[u][i];v = f[v][i];}}return f[u][0];
}
int main() {scanf("%d%d%d", &n, &m, &s);for (int i = 1; i < n; i++) {int x, y;scanf("%d%d", &x, &y);add(x, y);add(y, x);}dfs(s, 0);while (m--) {int a, b;scanf("%d%d", &a, &b);printf("%d\n", lca(a, b));}return 0;
}

这篇关于洛谷 P3379 【模板】最近公共祖先(LCA) 做题反思(2024.3.21)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/834026

相关文章

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

C++11第三弹:lambda表达式 | 新的类功能 | 模板的可变参数

🌈个人主页: 南桥几晴秋 🌈C++专栏: 南桥谈C++ 🌈C语言专栏: C语言学习系列 🌈Linux学习专栏: 南桥谈Linux 🌈数据结构学习专栏: 数据结构杂谈 🌈数据库学习专栏: 南桥谈MySQL 🌈Qt学习专栏: 南桥谈Qt 🌈菜鸡代码练习: 练习随想记录 🌈git学习: 南桥谈Git 🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈�

poj1330(LCA最近公共祖先)

题意:求最近公共祖先 思路:之前学习了树链剖分,然后我就用树链剖分的一小部分知识就可以解这个题目了,记录每个结点的fa和depth。然后查找时,每次将depth大的结点往上走直到x = y。 代码如下: #include<iostream>#include<algorithm>#include<stdio.h>#include<math.h>#include<cstring>

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

poj 2104 and hdu 2665 划分树模板入门题

题意: 给一个数组n(1e5)个数,给一个范围(fr, to, k),求这个范围中第k大的数。 解析: 划分树入门。 bing神的模板。 坑爹的地方是把-l 看成了-1........ 一直re。 代码: poj 2104: #include <iostream>#include <cstdio>#include <cstdlib>#include <al

最大流、 最小费用最大流终极版模板

最大流  const int inf = 1000000000 ;const int maxn = 20000 , maxm = 500000 ;struct Edge{int v , f ,next ;Edge(){}Edge(int _v , int _f , int _next):v(_v) ,f(_f),next(_next){}};int sourse , mee

【LabVIEW学习篇 - 21】:DLL与API的调用

文章目录 DLL与API调用DLLAPIDLL的调用 DLL与API调用 LabVIEW虽然已经足够强大,但不同的语言在不同领域都有着自己的优势,为了强强联合,LabVIEW提供了强大的外部程序接口能力,包括DLL、CIN(C语言接口)、ActiveX、.NET、MATLAB等等。通过DLL可以使用户很方便地调用C、C++、C#、VB等编程语言写的程序以及windows自带的大

C++语法知识点合集:11.模板

文章目录 一、非类型模板参数1.非类型模板参数的基本形式2.指针作为非类型模板参数3.引用作为非类型模板参数4.非类型模板参数的限制和陷阱:5.几个问题 二、模板的特化1.概念2.函数模板特化3.类模板特化(1)全特化(2)偏特化(3)类模板特化应用示例 三、模板分离编译1.概念2.模板的分离编译 模版总结 一、非类型模板参数 模板参数分类类型形参与非类型形参 非类型模板