最近公共祖先(lca)倍增算法【模板】

P3379 【模板】最近公共祖先(LCA) - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

复制代码
#include<bits/stdc++.h>
#include<cstdio>
using namespace std;
const int N=5e5+100;
const int inf=0x3f3f3f;
int n,m,s;
vector<int>g[N];
int dep[N];//存u点的深度
int fa[N][20];//存从u点向上跳2^i层的祖先节点
void dfs(int u,int father)
{
    dep[u]=dep[father]+1;//该点深度为他的父节点深度加1
    fa[u][0]=father;//当前节点向上跳一步就是他的父节点
    for(int i=1;i<=19;i++)
        fa[u][i]=fa[fa[u][i-1]][i-1];//递推点的所有能跳到的祖先节点,越界为0
    for(auto i:g[u])
        if(i!=father)
            dfs(i,u);
}
int lca(int u,int v)
{
    if(dep[u]<dep[v])//让u点跳
        swap(u,v);
    /*for(int i=19;i>=0;i--)//先让两点到达同一层
    {
        if(dep[fa[u][i]]>=dep[v])//u节点跳到祖先节点后的深度大于v点的深度
            u=fa[u][i];//u点变到其祖先节点处
    }*/
    int dx=dep[u]-dep[v];
    for(int i=0;dx>0;i++,dx=dx/2)
    {
        if(dx&1)
        {
            u=fa[u][i];
        }
    }
    //两种写法,让两点到达同一层
    
    if(u==v)//特判,v点为u的祖先节点
        return v;
    for(int i=19;i>=0;i--)
    {
        if(fa[u][i]!=fa[v][i])//两点跳后的祖先节点的父节点是否相等,不相等i-1
            u=fa[u][i],v=fa[v][i];
        //停在最近公共祖先的下层的儿子节点(此时两点的父节点就是公共最近祖先)
    }
    return fa[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);
        g[x].push_back(y);
        g[y].push_back(x);
    }
    dfs(s,0);//从根节点遍历树
    while(m--)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        printf("%d\n",lca(a,b));
    }
    return 0;
}
相关推荐
倒头就睡的小比特3 天前
算法竞赛C++常用的STL
c++·算法
小羊没烦恼!3 天前
初探性能优化——2个月到4小时的性能提升
java·开发语言·windows·算法·c#
猎头南楼3 天前
知识社区推荐系统实践:新用户冷启动与长短期兴趣建模的挑战 资深推荐算法工程师
人工智能·深度学习·算法·机器学习
旖旎夜光3 天前
力控面试题 01.01: 判定字符是否唯一(位运算) —— 题解
c++·学习·算法·leetcode·力控
wzdark3 天前
大规模并行计算中的负载均衡算法研究4
算法
Because_of_Her13 天前
并查集-听课笔记
笔记·算法·并查集
码流子3 天前
高速公路安全监测实践:碰撞监测预警+物联网底座,从感知到处置的闭环
大数据·人工智能·物联网·算法·架构
another heaven3 天前
【算法/C++ MD5算法能否逆解码?原理、C++实现与同类哈希算法对比】
c++·算法·哈希算法
wzdark3 天前
从算法设计模式看编程思维的抽象能力4
算法
2601_962218613 天前
万象生鲜系统称重自动多退少补算法解决生鲜非标品痛点
大数据·数据库·人工智能·python·算法