【SPOJ913】Query on a tree II

題意：

給定一棵n個點的樹，邊具有邊權。要求作以下操作：

DIST a b 詢問點a至點b路徑上的邊權之和

KTH a b k 詢問點a至點b有向路徑上的第k個點的編號

有多組測試數據，每組數據以DONE結尾。

 

解析：
       樹上倍增。

 

代碼：
 

#include <bits/stdc++.h>
using namespace std;

const int Max=20005;
int t,n,m,size,k;
int first[Max],dis[Max],dep[Max];
int f[Max][17];
struct shu{int to,next,len;};
shu edge[Max<<1];
char ch[6];

inline int get_int()
{
    int x=0,f=1;
    char c;
    for(c=getchar();(!isdigit(c))&&(c!='-');c=getchar());
    if(c=='-') f=-1,c=getchar();
    for(;isdigit(c);c=getchar()) x=(x<<3)+(x<<1)+c-'0';
    return x*f;
}
inline void print(int x)
{
    if(x<0) putchar('-'),x=-x;
    if(x>9) print(x/10);
    putchar('0'+x%10);
}

inline void clean()
{
    size=0;
    for(int i=1;i<=n;i++) first[i]=dis[i]=dep[i]=0;
    memset(f,0,sizeof(f));
}

inline void build(int x,int y,int z)
{
    edge[++size].next=first[x];
    first[x]=size;
    edge[size].to=y,edge[size].len=z;
}

inline void dfs(int point)
{
    for(int u=first[point];u;u=edge[u].next)
    {
      int to=edge[u].to;
      if(to==f[point][0]) continue;
      f[to][0]=point,dep[to]=dep[point]+1,dis[to]=dis[point]+edge[u].len;
      dfs(to);
    }
}

inline int LCA(int x,int y)
{
    if(dep[x]<dep[y]) swap(x,y);
    int len=dep[x]-dep[y];
    for(int i=15;i>=0;i--)
      if(len>=(1<<i)) len-=1<<i,x=f[x][i];
    if(x==y) return x;
    for(int i=15;i>=0;i--)
      if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
    return f[x][0];
}
inline int solve(int x,int y,int k)
{
    int fa=LCA(x,y),s=x;
    int len1=min(dep[x]-dep[fa],k-1),len2=dep[y]-dep[fa]-(k-1-(dep[x]-dep[fa]));
    for(int i=15;i>=0;i--)
      if(len1>=(1<<i))len1-=1<<i,x=f[x][i];
    if(dep[s]-dep[fa]>=k-1) return x;
    for(int i=15;i>=0;i--)
      if(len2>=(1<<i))len2-=1<<i,y=f[y][i];
    return y;
}

int main()
{
    t=get_int();
    while(t--)
    {
      n=get_int();
      clean();
      for(int i=1;i<n;i++)
      {
      	int x=get_int(),y=get_int(),z=get_int();
      	build(x,y,z),build(y,x,z);
      }
      dfs(1);
      for(int j=0;j<=15;j++)
        for(int i=1;i<=n;i++)
          if(f[i][j-1]) f[i][j]=f[f[i][j-1]][j-1];
      while(1)
      {
     	scanf("%s",ch);
     	if(ch[1]=='O') break;
     	int x=get_int(),y=get_int();
     	if(ch[1]=='I') print(dis[x]+dis[y]-2*dis[LCA(x,y)]),putchar('\n');
     	else k=get_int(),print(solve(x,y,k)),putchar('\n');
      }
    }
    return 0;
}