hdu2586 lca+rmq

#include<bits/stdc++.h>using namespace std;const int N=40010;const int M=20;int s;struct node{    int v,w;    node() {}    node(int vv,int ww)    {        v=vv;        w=ww;    }};vector<node> q[N];int fa[N][M];int dis[N][M];int dep[N];int n,m;void dfs(int u){    for(int i=1; i<=s; i++)    {        fa[u][i]=fa[fa[u][i-1]][i-1];        dis[u][i]=dis[u][i-1]+dis[fa[u][i-1]][i-1];        if(!fa[u][i]) break;    }    for(int i=0; i<q[u].size(); i++)    {        int v=q[u][i].v,w=q[u][i].w;        if(v!=fa[u][0])        {            dep[v]=dep[u]+1;            fa[v][0]=u;            dis[v][0]=w;            dfs(v);        }    }}void init(){    s=floor(log(n+0.0)/log(2.0));    //最多能跳的2^i祖先    dep[0]=-1;    for(int i=1; i<=n; i++)q[i].clear();}int LCA(int a,int b){    if(dep[a]>dep[b])swap(a, b);    int ans=0;    for(int i=s; i>=0; i--)        if(dep[a]<dep[b]&&dep[fa[b][i]]>=dep[a])            ans+=dis[b][i],b=fa[b][i];    for(int i=s; i>=0; i--)        if(fa[a][i]!=fa[b][i])//a,b的2^i祖先不一样=>没跳到同样深度位置，接着跳            ans+=dis[a][i],ans+=dis[b][i],a=fa[a][i],b=fa[b][i];//一起往上跳    if(a!=b)ans+=dis[a][0],ans+=dis[b][0];    //没跳到一个地方，那么再往上一个结点就是它们的LCA    return ans;}int main(){    int t;    scanf("%d",&t);    while(t--)    {        scanf("%d%d",&n,&m);        init();        for(int i=1; i<n; i++)        {            int u,v,w;            scanf("%d%d%d",&u,&v,&w);            q[u].push_back(node(v,w));            q[v].push_back(node(u,w));        }        dfs(1);        for(int i=0; i<m; i++)        {            int u,v;            scanf("%d%d",&u,&v);            printf("%d\n",LCA(u,v));        }    }}