倍增算法 最近公共祖先(LCA)C++

倍增算法(Binary Lifting)是一种用于解决最近公共祖先(LCA)问题的算法。

最近公共祖先问题是指在一棵有根树中,给定两个节点,找到它们的最近公共祖先节点。最近公共祖先节点是指在树中同时是这两个节点的祖先节点,并且离这两个节点的距离最近。

倍增算法的核心思想是在树的深度方向上使用二进制的方式进行跳跃。具体步骤如下:

  1. 预处理:对于每个节点,计算它的2^i级祖先节点,其中i从0开始递增。这可以通过自底向上的动态规划来实现。第一级祖先节点是每个节点的父节点,第二级祖先节点是每个节点的父节点的父节点,以此类推。

  2. 查询:给定两个节点x和y,首先判断它们的深度d,将深度较大的节点向上移动到与另一个节点处于同一深度。然后,从最大的2^i开始,每次向上跳跃2^i个级别,直到两个节点达到同一级别。在跳跃过程中,每次都检查是否达到了公共祖先节点。每个节点跳跃的级别数取决于两个节点的深度差。

倍增算法的时间复杂度为O(log n),其中n为树的节点数。预处理需要O(n log n)的时间来计算每个节点的祖先节点,查询阶段需要O(log n)的时间来计算最近公共祖先节点。空间复杂度为O(n log n),用于存储每个节点的祖先节点。

复制代码
#include <bits/stdc++.h>
#define ll long long
#define endl "\n"
#define KUI ios::sync_with_stdio(0), cin.tie(0), cout.tie(0)
using namespace std;
const int con = 5e5 + 5;
const int mod = 998244353;
int n, m, s, k, dep[con], fa[con][20];
vector<int> v[con];
void dfs(int u, int f)
{
    dep[u] = dep[f] + 1;
    fa[u][0] = f;
    for (int i = 1; i <= 19; i++)
    {
        fa[u][i] = fa[fa[u][i - 1]][i - 1];
    }
    for (int v : v[u])
    {
        if (v != f)
        {
            dfs(v, u);
        }
    }
}
int lca(int u, int v)
{
    if (dep[u] < dep[v])
    {
        swap(u, v);
    }
    for (int i = 19; i >= 0; i--)
    {
        if (dep[fa[u][i]] >= dep[v])
        {
            u = fa[u][i];
        }
    }
    if (u == v)
    {
        return v;
    }
    for (int i = 19; i >= 0; i--)
    {
        if (fa[u][i] != fa[v][i])
        {
            u = fa[u][i], v = fa[v][i];
        }
    }
    return fa[u][0];
}
int main()
{
    KUI;
    cin >> n >> m >> s;
    for (int i = 1; i <= n - 1; i++)
    {
        int a1, a2;
        cin >> a1 >> a2;
        v[a1].push_back(a2);
        v[a2].push_back(a1);
    }
    dfs(s, 0);
    for (int i = 1; i <= m; i++)
    {
        int a1, a2;
        cin >> a1 >> a2;
        cout << lca(a1, a2) << endl;
    }
    return 0;
}
相关推荐
HAPPY酷2 分钟前
std::pair` 与 `std::map` 基础
开发语言·c++·算法
懒神降世2 分钟前
基于iVentoy的PXE服务器的部署实战指南
运维·服务器·开发语言·云原生·vmware·openeuler·iventoy
柒儿吖3 分钟前
基于 lycium 在 OpenHarmony 上交叉编译 cppDES 完整实践
c++·harmonyos
山东布谷网络科技7 分钟前
对标Yalla和Chamet:海外直播语聊APP中多人派对房的关键技术细节
java·开发语言·人工智能·php·语音识别·软件需求·海外电商系统开发
爱搞事的程小猿14 分钟前
qml自定义扩展模块
c++·qt·qml
喜欢吃燃面15 分钟前
基础算法:高精度
开发语言·c++·学习·算法
3GPP仿真实验室15 分钟前
【Matlab源码】6G候选波形:OFDM-IM 增强仿真平台 GIM、MM、IQ
开发语言·网络·matlab
xuxie9916 分钟前
【无标题】
java·开发语言
堕27417 分钟前
java数据结构当中的《Lambda表达式》
java·数据结构·python
叫我一声阿雷吧19 分钟前
【JS实战案例】实现图片懒加载(基础版)原生JS+性能优化,新手可直接复现
开发语言·javascript·性能优化·js图片懒加载