【POJ1986】【模板】LCA

LCA 最近公共祖先

倍增法
f[i][j]表示i的2^j祖先，f[i][0]即表示i的父亲
首先预处理出f数组，一次dfs然后递推即可
接下来求LCA，思路就是让两个需要求LCA的点上升到同一层，然后一直向上找到他们LCA的子节点即可
详细说明代码里有

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define Maxn 500010
int ecnt,head[Maxn],deep[Maxn],f[Maxn][30],dis[Maxn],n,m;
struct node
{
    int next,to,dis;    
}e[4*Maxn+5];
void Add_edge(int from,int to,int dis)
{
    e[++ecnt].next=head[from];
    e[ecnt].to=to;
    e[ecnt].dis=dis;
    head[from]=ecnt;
}

void dfs(int u) //更新深度及dis值
{
    for(int i=head[u];i!=-1;i=e[i].next)
    {
        int v=e[i].to;
        int d=e[i].dis;
        if(!deep[v])
        {
            deep[v]=deep[u]+1;
            dis[v]=dis[u]+d;
            f[v][0]=u;
            dfs(v);
        }
    }
}

void init() //i的2^j祖先就是i的2^（j-1）祖先的 2^（j-1）祖先
{
    for(int j=1;(1<<j)<=n;j++)
        for(int i=1;i<=n;i++)
        {
            if(f[i][j-1])
                f[i][j]=f[ f[i][j-1] ][j-1];
        }   
}

int LCA(int a,int b)
{
    if(deep[a]<deep[b]) swap(a,b); //使a成为深度较深的那个点
    int i;
    for(i=0;(1<<i)<=deep[a];i++);
    i--;
    for(int j=i;j>=0;j--) //使a，b深度相同
    {
        if(deep[a]-(1<<j)>=deep[b])
            a=f[a][j];
    }
    if(a==b) return a; //同一层同一节点
    for(int j=i;j>=0;j--)//一直向上找，找到两个节点不同的最靠上的那一个
    {
        if(f[a][j]&&f[a][j]!=f[b][j])
        {
            a=f[a][j];
            b=f[b][j];
        }
    }
    return f[a][0];//返回他们的父亲
}
int main()
{
    int aa,bb,cc;
    int k;
    char d;
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d %c",&aa,&bb,&cc,&d);
        Add_edge(aa,bb,cc);
        Add_edge(bb,aa,cc);
    }
    deep[1]=1;
    f[1][0]=1;
    dfs(1);
    init();
    scanf("%d",&k);
    for(int i=1;i<=k;i++)
    {
        scanf("%d%d",&aa,&bb);
        printf("%d\n",dis[aa]+dis[bb]-2*dis[LCA(aa,bb)]);
    }
    return 0;
}