求公共祖先LCA

最新推荐文章于 2024-07-20 21:06:12 发布

xizi_ghq

最新推荐文章于 2024-07-20 21:06:12 发布

阅读量160

点赞数 1

分类专栏： c++ 算法

本文链接：https://blog.csdn.net/xizi_ghq/article/details/95012418

版权

c++ 同时被 2 个专栏收录

186 篇文章 3 订阅

订阅专栏

算法

166 篇文章 0 订阅

订阅专栏

树剖：

#include<bits/stdc++.h>
using namespace std;
const int N=5e5+10;
struct edge{int to,nex;}e[N*2];
int n,m,R,cnt=1,head[N],d[N],fa[N],size[N],son[N],top[N];
void add(int x,int y){
	e[cnt].to=y;
	e[cnt].nex=head[x];
	head[x]=cnt++;
}
void dfs1(int u,int f,int depth){
	d[u]=depth;
	fa[u]=f;
	size[u]=1;
	for(int x,i=head[u];i;i=e[i].nex){
		x=e[i].to;
		if(x!=f){
			dfs1(x,u,depth+1);
			size[u]+=size[x];
			if(size[son[u]]<size[x])son[u]=x;
		}
	} 
}
void dfs2(int u,int tp){
	top[u]=tp;
	if(son[u])dfs2(son[u],tp);
	for(int x,i=head[u];i;i=e[i].nex){
		x=e[i].to;
		if(x!=fa[u]&&x!=son[u])
			dfs2(x,x);
	}
}
int main(){
	scanf("%d%d%d",&n,&m,&R);
	for(int x,y,i=1;i<n;i++){
		scanf("%d%d",&x,&y);
		add(x,y);
		add(y,x); 
	}
	cnt=0;
	dfs1(R,0,1);
	dfs2(R,R);
	for(int x,y,i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		while(top[x]!=top[y]){
			if(d[top[x]]>=d[top[y]])x=fa[top[x]];
			else y=fa[top[y]];
		}
		printf("%d\n",d[x]<d[y]?x:y);
	}
	return 0;
}

int lca(int x,int y){
    for(;top[x]!=top[y];){
        if(deep[top[x]]<deep[top[y]])
            swap(x,y);
        x=dad[top[x]];
    }
    if(deep[x]>deep[y])
        swap(x,y);
    return x;
}