倍增求最近公共祖先(LCA)

传送门:P3379 【模板】最近公共祖先(LCA)

#include<bits/stdc++.h>
#define M 514514
using namespace std;

int n,m,s,tot=0,nex[2*M],first[2*M],to[2*M],dep[M],f[M][21];

inline int read(){
	int x=0;char c=getchar();
	while(c<'0'||c>'9'){c=getchar();}
	while(c>='0' and c<='9'){x=(x<<3)+(x<<1)+c-'0';c=getchar();}
	return x;
}

inline void add(int x,int y){
	nex[++tot]=first[x];first[x]=tot;to[tot]=y;
	nex[++tot]=first[y];first[y]=tot;to[tot]=x;
}

inline void dfs(int u,int fa){
	dep[u]=dep[fa]+1;
	for(int i=1;i<=20;i++)
		f[u][i]=f[f[u][i-1]][i-1];//递推 
	for(int e=first[u];e;e=nex[e]){
		int v=to[e];
		if(v==fa)continue;
		f[v][0]=u;
		dfs(v,u);
	}
}

inline int lca(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	for(int i=20;i>=0;i--){
		if(dep[f[x][i]]>=dep[y])x=f[x][i];
		if(x==y)return x;
	}
	for(int i=20;i>=0;i--){
		if(f[x][i]!=f[y][i]){
			x=f[x][i];
			y=f[y][i];
		}
	}
	return f[x][0];
}

int main(){
	n=read(),m=read(),s=read();
	for(int i=1;i<=n-1;i++){
		int x,y;
		x=read(),y=read();
		add(x,y);
	}
	dfs(s,0);
	for(int i=1;i<=m;i++){
		int a,b;
		a=read(),b=read();
		cout<<lca(a,b)<<'\n';
	}
	return 0;
}		
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值