AcWing 1172. 祖孙询问 题解（最近公共祖先lca）

AcWing 1172. 祖孙询问
寻找最近公共祖先lca模板题

#include<bits/stdc++.h>

using namespace std;

const int N = 4e4 + 10, M = N * 2;

int e[M], ne[M], h[N], idx;
int n, m;
int depth[N];  //记录每个点的深度
int fa[N][16];  //fa[i][j]记录的是编号为i的点向上遍历2^i次方的点的编号是几
int q[N];
int root;

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;
	int hh = 0, tt = 0;
	q[0] = root;
	
	while(hh <= tt){
		int t = q[hh ++ ];
		
		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[ ++ tt] = j;
				fa[j][0] = t;
				for(int k = 1; k <= 15; k ++ ){
					fa[j][k] = fa[fa[j][k - 1]][k - 1];
				}
			}
		}
	}
}

int lca(int a, int b){
	if(depth[a] < depth[b]) swap(a, b);  //把a置为更靠下的节点 
	for(int k = 15; k >= 0; k -- ){  //不断向上遍历a的祖先，找到第一个深度不大于b的a的祖先节点 
		if(depth[fa[a][k]] >= depth[b])
			a = fa[a][k];
	}
	if(a == b) return a;
	for(int k = 15; k > 0; k -- ){
		if(fa[a][k] != fa[b][k]){  //如果两个节点目前的父节点不同 
			a = fa[a][k];
			b = fa[b][k];
		}
	}
	return fa[a][0];  //遍历最终找到的祖先节点，注意可能不是公共的 
}

int main()
{
	cin>>n;
	memset(h, -1, sizeof h); 
	
	for(int i = 0; i < n; i ++ ){
		int a, b;
		cin>>a>>b;
		if(b == -1) root = a;
		else add(a, b), add(b, a);
	}
	bfs(root);
	cin>>m;
	for(int i = 0; i < m; i ++ ){
		int a, b;
		cin>>a>>b;
		int t = lca(a, b);
		if(t == a) puts("1");
		else if(t == b) puts("2");
		else puts("0");
	}
	return 0;
}