最近公共祖先（三种方法）

1.向上标记法O(n)

如果求的是 x , y 两点的最近公共祖先，则任选一点，假设选的是 x 点，将从 x 开始（包括x）向父节点走直到根节点的路径 都做标记，再从 y 开始向上走，遇到的第一个有标记的点就是它们的最近公共祖先

2.倍增法(nlogn)

令 f [ i , j ] 表示从 i 开始，向上走 2 ^ j 步所能走到的节点  0 <= j <= logn

令depth [ i ] 表示深度  

设置哨兵：如果从 i 开始跳，走2^j 会跳过根节点，则f [ i , j ] = 0

跳法ti：如果需要跳 2 ^ j 步，先跳 2 ^ j-1步，再跳 2 ^ j-1 步

步骤：1.先让两个点跳到同一层

           2.让两个点同时往上跳，直到跳到它们的最近公共祖先的下一层

图示：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=40010,M=N*2;
int n,m;
int h[N],e[M],ne[M],w[M],idx;
int depth[N],fa[N][16];
int q[N];
void add(int a,int b)
{
	e[idx]=b;
	ne[idx]=h[a];
	h[a]=idx++;
}
void bfs(int root)//同时处理depth和fa数组
{
	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);
	for(int k=15;k>=0;k--)
	{
		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()
{
	scanf("%d",&n);
	int root=0;
	memset(h,-1,sizeof h);
	for(int i=0;i<n;i++)
	{
		int a,b;
		scanf("%d%d",&a,&b);
		if(b==-1)
		root=a;
		else
		add(a,b),add(b,a);
	}
	bfs(root);
	scanf("%d",&m);
	while(m--)
	{
		int a,b;
		scanf("%d%d",&a,&b);
		int p=lca(a,b);
		if(p==a)
		puts("1");
		else if(p==b)
		puts("2");
		else
		puts("0"); 
	}
	return 0;
}

在枚举怎么跳时，用二进制拼凑步数，如要跳11步，最大的小于11的数是跳8，剩3，再跳2……

Tarjan离线求法O(n+m)

在线是指进行一次询问，输出一次答案

离线是指将所有询问读进来之后统一处理

主要思想：

距离

给出 n 个点的一棵树，多次询问两点之间的最短距离。

注意：

• 边是无向的。
• 所有节点的编号是 1,2,…,n。

输入格式

第一行为两个整数 n 和 m。n 表示点数，m 表示询问次数；

下来 n−1 行，每行三个整数 x,y,k，表示点 x 和点 y 之间存在一条边长度为 k；

再接下来 m 行，每行两个整数 x,y，表示询问点 x 到点 y 的最短距离。

树中结点编号从 1 到 n。

输出格式

共 m 行，对于每次询问，输出一行询问结果。

数据范围

42 ≤n ≤ 104,
1 ≤ m ≤ 2×104,
0 < k ≤ 100,
1 ≤x,y≤ n

输入样例1：

2 2 
1 2 100 
1 2 
2 1

输出样例1：

100
100

输入样例2：

3 2
1 2 10
3 1 15
1 2
3 2

输出样例2：

10
25

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>

using namespace std;

typedef pair<int,int> PII;
const int N=20010,M=N*2;

int n,m;
int h[N],e[M],ne[M],w[M],idx;
int dist[N];
int p[N];
int res[N];
int st[N];
vector<PII> query[N];

void add(int a,int b,int c)
{
	e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

void dfs(int u,int fa)//每个点和一号点的距离 
{
	for(int i=h[u];~i;i=ne[i])
	{
		int j=e[i];
		if(j==fa)
		continue;
		dist[j]=dist[u]+w[i];
		dfs(j,u);
	}
}

int find(int x)
{
	if(p[x]!=x)
	p[x]=find(p[x]);
	return p[x];
}

void tarjan(int u)
{
	st[u]=1;//当前正在遍历类 
	for(int i=h[u];~i;i=ne[i])
	{
		int j=e[i];
		if(!st[j])
		{
			tarjan(j);
			p[j]=u;
		}
	}
	for(auto item:query[u])//所有和u相关的查询 
	{
		int y=item.first,id=item.second;
		if(st[y]==2)
		{
			int anc=find(y);//最近公共祖先 
			res[id]=dist[u]+dist[y]-dist[anc]*2;
		}
	}
	st[u]=2;
}

int main()
{
	scanf("%d%d",&n,&m);
	memset(h,-1,sizeof h);
	for(int i=0;i<n-1;i++)
	{
		int a,b,c;
		scanf("%d%d%d",&a,&b,&c);
		add(a,b,c),add(b,a,c);
	} 
	for(int i=0;i<m;i++)
	{
		int a,b;
		scanf("%d%d",&a,&b);
		if(a!=b)
		{
			query[a].push_back({b,i});
			query[b].push_back({a,i});
		}
	}
	for(int i=1;i<=n;i++)
	{
		p[i]=i;
	}
	dfs(1,-1);
	tarjan(1);
	for(int i=0;i<m;i++)
	printf("%d\n",res[i]); 
}