POJ 1986 Distance Queries (在线LCA转RMQ)

题目地址：POJ 1986

纯模板题。输入的最后一个字母是多余的，完全不用管。还有注意询问的时候有相同点的情况。

代码如下：

#include <iostream>
#include <string.h>
#include <math.h>
#include <queue>
#include <algorithm>
#include <stdlib.h>
#include <map>
#include <set>
#include <stdio.h>
using namespace std;
#define LL __int64
#define pi acos(-1.0)
const int mod=1e9+7;
const int INF=0x3f3f3f3f;
const double eqs=1e-9;
const int MAXN=50000+10;
int F[MAXN<<1], fir[MAXN], rmq[MAXN<<1], dis[MAXN];
int head[MAXN], cnt, tot;
struct node
{
        int u, v, w, next;
}edge[1000000];
void add(int u, int v, int w)
{
        edge[cnt].v=v;
        edge[cnt].w=w;
        edge[cnt].next=head[u];
        head[u]=cnt++;
}
void dfs(int u, int fa, int dep)
{
        F[++tot]=u;
        fir[u]=tot;
        rmq[tot]=dep;
        for(int i=head[u];i!=-1;i=edge[i].next){
                int v=edge[i].v;
                if(v==fa) continue ;
                dis[v]=dis[u]+edge[i].w;
                dfs(v,u,dep+1);
                F[++tot]=u;
                rmq[tot]=dep;
        }
}
struct ST
{
        int dp[MAXN<<1][30], i, j;
        void init()
        {
                for(i=1;i<=tot;i++){
                        dp[i][0]=i;
                }
                for(j=1;(1<<j)<=tot;j++){
                        for(i=1;i<=tot-(1<<j-1)+1;i++){
                                dp[i][j]=rmq[dp[i][j-1]]<rmq[dp[i+(1<<j-1)][j-1]]?dp[i][j-1]:dp[i+(1<<j-1)][j-1];
                        }
                }
        }
        int Query(int l, int r)
        {
                if(r<l) swap(l,r);
                int k=0;
                while((1<<k+1)<=r-l+1) k++;
                return rmq[dp[l][k]]<rmq[dp[r+1-(1<<k)][k]]?dp[l][k]:dp[r+1-(1<<k)][k];
        }
}st;
void init()
{
        memset(head,-1,sizeof(head));
        memset(dis,0,sizeof(dis));
        cnt=tot=0;
}
int main()
{
        int n, m, i, u, v, k, tmp, w;
        char c;
        scanf("%d%d",&n,&m);
        init();
        while(m--){
                scanf("%d%d%d %c",&u,&v,&w,&c);
                add(u,v,w);
                add(v,u,w);
        }
        dfs(1,-1,0);
        st.init();
        scanf("%d",&k);
        while(k--){
                scanf("%d%d",&u,&v);
                if(u==v){
                        printf("0\n");
                        continue ;
                }
                tmp=F[st.Query(fir[u],fir[v])];
                printf("%d\n",dis[u]+dis[v]-2*dis[tmp]);
        }
        return 0;
}