最近公共祖先（LCA）

LCA

lca在处理一个树上问题非常的常用，接下来就是利用倍增来求LCA的介绍。
首先我们发现在往上跳的过程中，不需要一步一步的走，我们可以直接大步子往上跳直接跳到LCA上，这是我们就想到使用倍增来求LCA，我们先预处理出来大步跳到的点，然后我们就要二进制那样拆解的向上跳，但是我们要从大数开始跳。下面为代码。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
using namespace std;

int n,m,s,tot=0;
int head[1000010],nxt[1000010],to[1000001];
int d[500100],f[30][1000100];

int read()
{
    int x=0,flag=0;
    char ch=getchar();
    if(ch=='-') flag=1;
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9')x*=10,x+=ch-'0',ch=getchar();
    if(flag) return -x;
    return x;
}

void add(int x,int y)
{
	tot++;
	nxt[tot]=head[x];
	head[x]=tot;
	to[tot]=y;
	return;
}

void dfs(int u,int dep)
{
	d[u]=dep;
	for(int i=head[u];i!=-1;i=nxt[i])
	{
		int v=to[i];
		if(!d[v])dfs(v,dep+1),f[0][v]=u;
	}
}

int lca(int x,int y)
{
	int l=0;
	while((1<<l)<=n)l++;
	l--;
	if(d[x]<d[y])swap(x,y);
	for(int i=20;i>=0;i--)
	  if(d[y]<=d[x]-(1<<i))x=f[i][x];
	if(x==y)return x;
	for(int i=20;i>=0;i--)
	{
		if(f[i][x]!=f[i][y])
		{
			x=f[i][x];
			y=f[i][y];
		}
	}
	return f[0][x];
}

int main()
{
	memset(head,-1,sizeof(head));
	n=read(),m=read(),s=read();
	for(int i=1;i<n;i++)
	{
		int x,y;
		x=read(),y=read();
		add(x,y);
		add(y,x);
	}
	f[0][s]=s;
	dfs(s,1);
	for(int i=1;(i<<i)<=n;i++)
	  for(int j=1;j<=n;j++)
	    f[i][j]=f[i-1][f[i-1][j]];
	for(int i=1;i<=m;i++)
	{
		int l,r;
		l=read(),r=read();
		printf("%d\n",lca(l,r));
	}
	return 0;
}