hdu 2586 How far away ？ lca求最短路

给定一棵n个结点的树，有m个询问，问a,b之间的最短路是多少？。

设dis[x]为树根到x节点的路程，所以只要求出a,b的lca，答案就是dis[a]+dis[b]-2*dis[lca(a,b)]

dis的记录在遍历树的时候可以求得。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <vector>
#include <map>
#include <cmath>
#include <set>
#include <queue>
using namespace std;

const int INF=1e9+10;
const double EPS = 1e-10;  
typedef long long ll;
const int mod=1e9+7;

const int N = 40010;  
const int M = 25;
int n,m;
int dir[N];
int dp[2*N][M];  //这个数组记得开到2*N，因为遍历后序列长度为2*n-1  
bool vis[N];  
struct edge  
{  
    int u,v,w,next;  
}e[2*N];  
int tot,head[N];  
inline void add(int u ,int v ,int w ,int &k)  
{  
    e[k].u = u; e[k].v = v; e[k].w = w;  
    e[k].next = head[u]; head[u] = k++;  
    u = u^v; v = u^v; u = u^v;  
    e[k].u = u; e[k].v = v; e[k].w = w;  
    e[k].next = head[u]; head[u] = k++;  
}  
int ver[2*N],R[2*N],first[N];  
//ver:节点编号 R：深度 first：点编号位置   
void dfs(int u ,int dep)  
{  
    vis[u] = true; ver[++tot] = u; first[u] = tot; R[tot] = dep;  
    for(int k=head[u]; k!=-1; k=e[k].next)  
        if( !vis[e[k].v] )  
        {  
            int v = e[k].v,w = e[k].w;  
            dir[v] = dir[u] + w;  
            dfs(v,dep+1);  
            ver[++tot] = u; R[tot] = dep;  
        }  
}
//返回的不是最小值，而是最小值时对应的下标
void ST(int n)  
{  
    for(int i=1;i<=n;i++)  
        dp[i][0] = i;  
    for(int j=1;(1<<j)<=n;j++)  
    {  
        for(int i=1;i+(1<<j)-1<=n;i++)  
        {  
            int a = dp[i][j-1] , b = dp[i+(1<<(j-1))][j-1];  
            dp[i][j] = R[a]<R[b]?a:b;  
        }  
    }  
}
//中间部分是交叉的。  
int RMQ(int l,int r)  
{  
    int k=0;  
    while((1<<(k+1))<=r-l+1)  
        k++;  
    int a = dp[l][k], b = dp[r-(1<<k)+1][k]; //保存的是编号  
    return R[a]<R[b]?a:b;  
}  
  
int LCA(int u ,int v)  
{  
    int x = first[u] , y = first[v];  
    if(x > y) swap(x,y);  
    int res = RMQ(x,y);  
    return ver[res];  
}  

void init1(){
	 memset(head,-1,sizeof(head));  
     memset(vis,false,sizeof(vis));
}
void init2(){
	tot=0;
	dfs(1,1);
	ST(2*N-1);
}
int main(){
	//freopen("out.txt","w",stdout);
	int t;
	scanf("%d",&t);
	while(t--){
		int num=0;
		init1();
		scanf("%d %d",&n,&m);
		for(int i=0;i<n-1;i++){
			int u,v,w;  
            scanf("%d%d%d",&u,&v,&w);  
            add(u,v,w,num);
		}
		dir[1]=0;
		init2();
		while(m--){
			int u,v;
			scanf("%d %d",&u,&v);
			int k=LCA(u,v);
			printf("%d\n",dir[u]+dir[v]-2*dir[k]);
		}
	}	
    return 0;
}

const int N = 40010;  
const int M = 25;
int dp[2*N][M];  //这个数组记得开到2*N，因为遍历后序列长度为2*n-1  
bool vis[N];  
struct edge  
{  
    int u,v,w,next;  
}e[2*N];  
int tot,head[N];  
inline void add(int u ,int v ,int w ,int &k)  
{  
    e[k].u = u; e[k].v = v; e[k].w = w;  
    e[k].next = head[u]; head[u] = k++;  
    u = u^v; v = u^v; u = u^v;  
    e[k].u = u; e[k].v = v; e[k].w = w;  
    e[k].next = head[u]; head[u] = k++;  
}  
int ver[2*N],R[2*N],first[N];  
//ver:节点编号 R：深度 first：点编号位置   
void dfs(int u ,int dep)  
{  
    vis[u] = true; ver[++tot] = u; first[u] = tot; R[tot] = dep;  
    for(int k=head[u]; k!=-1; k=e[k].next)  
        if( !vis[e[k].v] )  
        {  
            int v = e[k].v,w = e[k].w;  
            dfs(v,dep+1);  
            ver[++tot] = u; R[tot] = dep;  
        }  
}
//返回的不是最小值，而是最小值时对应的下标,n=2*N-1
void ST(int n)  
{  
    for(int i=1;i<=n;i++)  
        dp[i][0] = i;  
    for(int j=1;(1<<j)<=n;j++)  
    {  
        for(int i=1;i+(1<<j)-1<=n;i++)  
        {  
            int a = dp[i][j-1] , b = dp[i+(1<<(j-1))][j-1];  
            dp[i][j] = R[a]<R[b]?a:b;  
        }  
    }  
}
//中间部分是交叉的。  
int RMQ(int l,int r)  
{  
    int k=0;  
    while((1<<(k+1))<=r-l+1)  
        k++;  
    int a = dp[l][k], b = dp[r-(1<<k)+1][k]; //保存的是编号  
    return R[a]<R[b]?a:b;  
}  
  
int LCA(int u ,int v)  
{  
    int x = first[u] , y = first[v];  
    if(x > y) swap(x,y);  
    int res = RMQ(x,y);  
    return ver[res];  
}  

void init1(){
	 memset(head,-1,sizeof(head));  
     memset(vis,false,sizeof(vis));
}
void init2(){
	tot=0;
	dfs(1,1);
	ST(2*N-1);
}