LCA之ST算法模板 poj-1986

n个点，m条边，k次询问。

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>
struct node{
    int to,next,w;
}a[100010];
int top,jl[40010],p[40010],deep[40010],begin[40010],dp[100010][20],first[40010];
void add(int x,int y,int w){
    a[++top].to=y;
    a[top].w=w;
    a[top].next=begin[x];
    begin[x]=top;
    return;
}
void dfs(int d){
    dp[++top][0]=d;
    first[d]=top;
    int h=begin[d];
    while(h){
        if(!p[a[h].to]){
            p[a[h].to]=1;
            deep[a[h].to]=deep[d]+1;
            jl[a[h].to]=jl[d]+a[h].w;
            dfs(a[h].to);
            dp[++top][0]=d;
        }

        h=a[h].next;
    }
    return;
}
void st(){
    int i,j,k;
    k=log(top) / log(2);
    for(i=1;i<=k;i++)
        for(j=1;j<=top;j++){
            if(j+(1<<i)>top)
                break;
            if(deep[dp[j][i-1] ]<deep[dp[j+(1<<i-1)][i-1] ])
                dp[j][i]=dp[j][i-1];
            else
                dp[j][i]=dp[j+(1<<i-1)][i-1];
        }
    return;
}
int rmq(int x1,int y1){
    int k;
    int x=first[x1];
    int y=first[y1];
    if(x>y){
        int xx=y;
        y=x;
        x=xx;
    }
    k = log(y-x+1) / log(2);
    if(deep[dp[x][k]]<deep[dp[y-(1<<k)+1][k]])
        return jl[x1] + jl[y1] - 2*jl[dp[x][k]];
    else 
        return jl[x1] + jl[y1] - 2*jl[dp[y-(1<<k)+1][k]];

}
int main(){
    int i,j,k,n,m,x,y,w;
    char c[1010];
    scanf("%d%d",&n,&m);
    for(i=1;i<=m;i++){
        scanf("%d%d%d",&x,&y,&w);
        add(x,y,w);
        add(y,x,w);
        scanf("%s",c);
    }
    top=0;
    p[1]=1;
    dfs(1);
    st();
    scanf("%d",&k);
    for(i=1;i<=k;i++){
        scanf("%d%d",&x,&y);
        printf("%d\n",rmq(x,y));
    }
    return 0;
}