最近公共祖先LCA倍增算法

题目一般是给你一颗树和树根节点的序号，然后再问你两个节点的最近公共祖先。

先说朴素的思想：先将这两个节点的深度统一，即把更深的节点提升到和更浅的节点一样的深度，然后如果两个节点不重合就将两个节点一起往上一个个找，直到两个点重合为止。当然，这样很容易就会T掉。

所以在这个基础上我们给出了倍增的算法：就是想办法缩减向上爬的次数。在这里有一个重要的公理，就是任何正数都可以由2的n次方加和得到。因此我们可以用一个fa[N][M]数组来存储某个节点n往上2的m次方的父节点。用一个dfs来记录每一个节点的上一个节点（也就是2的0次方），然后再用递推公式fa[j][i]=fa[fa[j][i-1]][i-1]来得出每个节点往上2的m次方的父节点。接下来跑一个lca问题就迎刃而解了。模板题传送门

模板题的代码：
 

#include<bits/stdc++.h>
using namespace std;
// #define int long long
const int N=5e5+5;
vector<int>e[N];
//bool vis[N];
const int M=20;
int fa[N][M];//表示N节点向上2的M次方的节点
int dep[N];//N节点的深度
void dfs(int x,int f=0)
{
    fa[x][0]=f;
    dep[x]=dep[f]+1;
    for(auto i:e[x])
    {
        if(i==f)//很重要，防止重复计算上一个节点
        continue;
        dfs(i,x);
    }
}
int LCA(int x,int y)
{
    if(dep[x]>dep[y])
    swap(x,y);
    int gap=0;
    while((1<<(gap+1))<=dep[y])
    ++gap;
    for(int i=gap;i>=0;i--)
    {
        if(dep[y]-dep[x]>=(1<<i))
        y=fa[y][i];//一步步将y同步到x的深度，倍增更快
    }
    if(x==y)
    return x;
    for(int i=gap;i>=0;i--)
    {
        if(fa[x][i]!=fa[y][i])
        {
            x=fa[x][i];//可以这样2次方爬，因为任何数字都可以拆解成2的n次方的和
            y=fa[y][i];
        }
    }
    return fa[x][0];
}
int main()
{
    int n,m,s;
    cin>>n>>m>>s;
    dep[s]=1;
    int t=n-1;
    while(t--)
    {
        int a,b;
        cin>>a>>b;
        e[a].push_back(b);
        e[b].push_back(a);       
    }
    dfs(s);
    for(int i=1;(1<<i)<=n;i++)
    {
        for(int j=1;j<=n;j++)
        fa[j][i]=fa[fa[j][i-1]][i-1];//倍增记录2的n次方父节点，没必要一个个往上推
    }
    while(m--)
    {
        int a,b;
        cin>>a>>b;
        int ans=LCA(a,b);
        cout<<ans<<endl;
    }
    return 0;
}

顺便记录一下vscode的快捷键，ctrl加z是撤回，加f是查找替换。(

天天开心！