https://www.luogu.org/problemnew/show/P3379
题目如上
关于LCA,利用二进制的特性解决问题,与快速幂有些相似;
须知定理:
任一实数n,有2^k+1 > n >= 2^k 那么n一定等于2^k , 2^k-1, 2^k-2, … 2^1, 2^0,中的任意几个的的和;
例如 9=2^3 + 2^0,
10=2^3 + 2^1;
11= 2^3 + 2^1 + 2^0;
那么求LCA可以不用一个一个往上跳了,利用这一规律便可以快速的找到LCA;
具体步骤如下:
1.将两个节点调整到同一高度
2.判断是否为同一节点
3.若2不成立将x,y一起往上跳,直到他们的父节点为同一个为止
f[i][j] 表示i节点向上跳2^i所到达的节点,先dfs求出f[i][0]即i的父节点,在递推求出其他;
递推公式如下
f[i][j]=f[f[i][j-1]][j-1],意思是i节点向上跳了2^j-1 后再向上跳2^j-1;
代码如下
// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;
vector<int> E[500010];
int deep[500010],f[500010][22];
int n,m,s;
void dfs(int now,int last){
deep[now]=deep[last]+1;
f[now][0]=last;
for(int i=1;i*i<=deep[now];i++) f[now][i]=f[f[now][i-1]][i-1];
for(int i=0;i<E[now].size();i++) if(!deep[E[now][i]]) dfs(E[now][i],now);
}
int Find_Lca(int x,int y){
if(deep[x]<deep[y]) swap(x,y);//保证x较深
for(int i=sqrt(deep[x]-deep[y])+1;i>=0;i--) if(deep[f[x][i]]>=deep[y]) x=f[x][i];//调整到同一高度
if(x==y) return y;//情况2
for(int i=sqrt(deep[x])+1;i>=0;i--) if(f[x][i]!=f[y][i]){x=f[x][i]; y=f[y][i];}
return f[x][0];
}
int read(){
int x=0,f=1;
char ch=getchar();
if(ch=='-') f=-1;
while(ch<'0'||ch>'9') ch=getchar();
while(ch>='0'&&ch<='9') x*=10,x+=ch-'0',ch=getchar();
return x*f;
}
int main(){
n=read(); m=read(); s=read();
for(int i=1;i<=n-1;i++){
int u=read(),v=read();
E[u].push_back(v); E[v].push_back(u);
}
dfs(s,0);//预处理出f[i][0],和深度
for(int i=1;i<=m;i++){
int u=read(),v=read();
printf("%d\n",Find_Lca(u,v));
}
return 0;
}