倍增法
fa[i,j]表示从第i个点,向上走2^j步所能走到的节点 0<=j<=logn
depth[i]表示节点i的深度
depth[i]=第i个点到根节点长度+1
步骤
(1)先将两个点跳到同一层
(2)让两个点同时往上一直跳,一直跳到最近公共祖先下一层
预处理: O(logn)
查询: O(logn)
来个模板
给定一棵包含 n 个节点的有根无向树,节点编号互不相同,但不一定是 1∼n。
有 m 个询问,每个询问给出了一对节点的编号 x 和 y,询问 x 与 y 的祖孙关系。
输入格式
输入第一行包括一个整数 表示节点个数;
接下来 n 行每行一对整数 a 和 b,表示 a 和 b 之间有一条无向边。如果 b 是 −1,那么 a 就是树的根;
第 n+2 行是一个整数 m 表示询问个数;
接下来 m 行,每行两个不同的正整数 x 和 y,表示一个询问。
输出格式
对于每一个询问,若 x 是 y 的祖先则输出 1,若 y 是 x 的祖先则输出 2,否则输出 0。
数据范围
1≤n,m≤4×104,
1≤每个节点的编号≤4×104
输入样例:
10
234 -1
12 234
13 234
14 234
15 234
16 234
17 234
18 234
19 234
233 19
5
234 233
233 12
233 13
233 15
233 19
输出样例:
1
0
0
0
2
#include<bits/stdc++.h>
using namespace std;
const int N=4e4+1000;
int depth[N],fa[N][20];
int e[N*2],ne[N*2],idx,h[N];
int n,m;
void add(int a,int b)
{
e[idx]=b;ne[idx]=h[a];h[a]=idx++;
}
void bfs(int root)
{
memset(depth,0x3f,sizeof depth);
depth[0]=0;//哨兵
depth[root]=1;
queue<int> q;
q.push(root);
while(q.size())
{
int t=q.front();
q.pop();
for(int i=h[t];~i;i=ne[i])
{
int j=e[i];
if(depth[j]>depth[t]+1)
{
depth[j]=depth[t]+1;
q.push(j);
fa[j][0]=t; //第j个点,跳0步是t
for(int k=1;k<=15;k++) fa[j][k]=fa[fa[j][k-1]][k-1];//第j个点,跳2^k步是谁
}
}
}
}
int lca(int x,int y)
{
if(depth[x]<depth[y]) swap(x,y);//把x和y跳到同一层中
for(int k=15;k>=0;k--)//从高往低跳
if(depth[fa[x][k]]>=depth[y]) x=fa[x][k];
if(x==y) return x;//x和y就是最近公共祖先
for(int k=15;k>=0;k--)//在同一层后,同时向上跳
{
if(fa[x][k]!=fa[y][k])//x和y还是没有跳到最近公共 祖先上
{
x=fa[x][k];
y=fa[y][k];
}
}
return fa[x][0];//这边返回fa[x][0],fa[y][0]都可以因为都在最近公共祖先上
}
int main()
{
memset(h,-1,sizeof h);
cin>>n;
int root=-1;
for(int i=1;i<=n;i++)
{
int a,b;
cin>>a>>b;
if(b==-1) root=a;
else add(a,b);add(b,a);
}
bfs(root); //预处理深度和fa[i][j]
cin>>m;
while(m--)
{
int a,b;
cin>>a>>b;
int p=lca(a,b);
if(p==a) cout<<1<<endl;
else if(p==b)cout<<2<<endl;
else cout<<0<<endl;
}
return 0;
}