树上最近公共祖先(lca 模板)

注意至少fa[N][log2(N)+1]
fa[i][j]:从i开始向上走 2 j 2^j 2j的大小即 i + 2 j i+2^j i+2j<==> i + 2 j − 1 + 2 j − 1 i+2^{j-1}+2^{j-1} i+2j1+2j1可得
fa[i][j]=fa[fa[i][j-1]][j-1]
循环应从log2[N]到0此时向上走由 2 p > ( 2 p − 1 ) 2^p>(2^p-1) 2p>(2p1)可知必定有一点刚好相等

传送门

#include <bits/stdc++.h>
using namespace std;
using ll = long long ;
const int N = 4e4+10;
vector<int>g[N<<1];
int dep[N],fa[N][16];
void bfs(int x)
{
    queue<int>q;
    memset(dep,0x3f,sizeof dep);
    dep[x]=1;
    dep[0]=0;
    q.push(x);
    while(!q.empty()){
        int temp=q.front();q.pop();
        for(auto &t:g[temp]){
            if(dep[t]>dep[temp]+1){
                dep[t]=dep[temp]+1;
                fa[t][0]=temp;
                q.push(t);
                for(int i=1;i<=15;i++)
                    fa[t][i]=fa[fa[t][i-1]][i-1];
            }
        }
    }
}
int lca(int a,int b)
{
    if(dep[a]<dep[b])swap(a,b);
    for(int i=15;i>=0;i--)
        if(dep[fa[a][i]]>=dep[b])
            a=fa[a][i];
    if(a==b)return a;
    for(int i=15;i>=0;i--)
        if(fa[a][i]!=fa[b][i])
        a=fa[a][i],b=fa[b][i];
    return fa[a][0];
}
int main()
{
    ios::sync_with_stdio(false);
    int n;
    cin>>n;
    int root;
    for(int i=0,u,v;i<n;i++){
        cin>>u>>v;
        if(~v) g[u].push_back(v),g[v].push_back(u);
        else root=u;
    }
    bfs(root);
    int q;cin>>q;
    while(q--){
        int l,r;
        cin>>l>>r;
        int k=lca(l,r);
        if(k==l)cout<<"1\n";
        else if(k==r)cout<<"2\n";
        else cout<<"0\n";
    }
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

倾海、

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值