最近公共祖先 (LCA) 模板

定义

在树中两个节点的最近公共祖先为:两个点的公共祖先里面离根最远的那个。

查询次数较少的情况

可以通过两次 d f s dfs dfs 解决:第一次 d f s dfs dfs 将根节点到其中一个节点的路径上的节点染色,第一次 d f s dfs dfs 遍历到的最后一个染色节点即为最近公共祖先,时间复杂度为 O ( n ) O(n) O(n)。

查询次数较多的情况

可以通过一遍 d f s dfs dfs 预处理求出各节点 i i i 所在的层数 h [ i ] h[i] h[i],和倍增数组 p a r par par : p a r [ i ] [ j ] par[i][j] par[i][j] 为与 i i i 距离为 2 j 2^j 2j 的祖先节点。然后设查询的两个节点为 u u u 和 v v v,设 h [ u ] ≥ [ v ] h[u]\ge[v] h[u]≥[v] ,若 h [ u ] > h [ v ] h[u]>h[v] h[u]>h[v] 则先将 u u u 上移至和 v v v 同一层。之后 u u u 和 v v v 处于同一层, 此时若 u u u 和 v v v 为同一节点,则 v v v 即为最近公共祖先;否则将 u u u 和 v v v 同时尽可能地上移且保证移动后的 u u u 和 v v v 为不同的节点,无法移动之后 u u u 和 v v v 的父节点即为最近公共祖先。预处理的时间复杂度为 O ( n log n ) O(n\text{log}n) O(nlogn),单次查询的时间复杂度为 O ( log n ) O(\text{log}n) O(logn)。求 l c a lca lca 的代码如下:

cpp 复制代码
int lca(int p, int q) {//这里用整数代表节点,若存在相同值的节点则用节点指针
    if (h[p] < h[q])
        swap(p, q);
    for (int dis = h[p] - h[q], i = 0; dis; i++, dis >>= 1)
        if (dis & 1)
            p = par[p][i];
    if (p == q)
        return p;
    for (int i = i0; i >= 0; i--)// i的初值i0可设为ceil(logn)
        if (par[p][i] != par[q][i]) {
            p = par[p][i];
            q = par[q][i];
        }
    return par[p][0];
}

完整的例子包括生成倍增数组等操作可以参考 二叉树的最近公共祖先 这道题的倍增做法:

cpp 复制代码
class Solution {
public:
    TreeNode *lowestCommonAncestor(TreeNode *root, TreeNode *p, TreeNode *q) {
        unordered_map<TreeNode *, vector<TreeNode *>> par;
        unordered_map<TreeNode *, int> h;
        function<void(TreeNode *, TreeNode *, int)> dfs = [&](TreeNode *root, TreeNode *p, int lev) {
            if (!root)
                return;
            h[root] = lev;
            par[root] = vector<TreeNode *>(18);
            par[root][0] = p;
            for (int i = 1; i < 18; i++)
                par[root][i] = par[par[root][i - 1]][i - 1];
            dfs(root->left, root, lev + 1);
            dfs(root->right, root, lev + 1);
        };
        dfs(root, root, 0);
        if (h[p] < h[q])
            swap(p, q);
        for (int dis = h[p] - h[q], i = 0; dis; i++, dis >>= 1)
            if (dis & 1)
                p = par[p][i];
        if (p == q)
            return p;
        for (int i = 17; i >= 0; i--)
            if (par[p][i] != par[q][i]) {
                p = par[p][i];
                q = par[q][i];
            }
        return par[p][0];
    }
};
相关推荐
冷崖2 分钟前
定时器的学习(二)
linux·c++·学习
B站计算机毕业设计之家5 分钟前
深度学习实战:python动物识别分类检测系统 计算机视觉 Django框架 CNN算法 深度学习 卷积神经网络 TensorFlow 毕业设计(建议收藏)✅
python·深度学习·算法·计算机视觉·分类·毕业设计·动物识别
大胆飞猪7 分钟前
高并发内存池日志
c++·项目
And_Ii18 分钟前
LeetCode 3350. 检测相邻递增子数组 II
数据结构·算法·leetcode
想唱rap21 分钟前
C++ string类的使用
开发语言·c++·笔记·算法·新浪微博
JAVA学习通22 分钟前
Replication(下):事务,一致性与共识
大数据·分布式·算法
胖咕噜的稞达鸭22 分钟前
C++中的父继子承(2)多继承菱形继承问题,多继承指针偏移,继承组合分析+高质量习题扫尾继承多态
c语言·开发语言·数据结构·c++·算法·链表·c#
蓝色汪洋25 分钟前
Completed String easy
算法
铭哥的编程日记27 分钟前
贪心算法精选30道编程题 (附有图解和源码)
算法·贪心算法
CoovallyAIHub30 分钟前
顶刊新发!上海交大提出PreCM:即插即用的旋转等变卷积,显著提升分割模型鲁棒性
人工智能·算法·计算机视觉