bzoj1602 [Usaco2008 Oct]牧场行走

题目

LCA模板题。

在树上求两个点之间的距离，可以用类似差分的方法。设dis[i]表示结点i到根的距离，那么s到t的距离为：

dis[s]+dis[t]-2*dis[lca(s,t)];

lca（s，t）表示它们的最近公共祖先。

#include<bits/stdc++.h>#define N 2000using namespace std;int n,q,x,y,z;int first[N+5],nxt[N+5],to[N+5],val[N+5],siz;int f[N+5][25],dis[N+5],dep[N+5];inline char nc(){    static char buf[100000],*p1=buf,*p2=buf;    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}inline int read(){    int x=0,b=1;    char c=nc();    for(;!(c<='9'&&c>='0');c=nc())if(c=='-')b=-1;    for(;c<='9'&&c>='0';c=nc())x=x*10+c-'0';    return x*b;}inline void write(int x){    if(x==0)putchar('0');    else    {        char buf[15];        int len=0;        if(x<0)putchar('-'),x=-x;        while(x)buf[++len]=x%10+'0',x/=10;        for(int i=len;i>=1;i--)putchar(buf[i]);    }    putchar('\n');}inline void add(int x,int y,int z){    nxt[siz]=first[x];    first[x]=siz;    to[siz]=y;    val[siz]=z;    siz++;}inline void init(int x,int fa){    f[x][0]=fa;    for(int i=1;i<=16;i++)        if(f[x][i-1])f[x][i]=f[f[x][i-1]][i-1];else break;    for(int i=first[x];i!=-1;i=nxt[i])    {        int u=to[i];        if(u==fa)continue;        dep[u]=dep[x]+1;        dis[u]=dis[x]+val[i];        init(u,x);    }}inline int LCA(int x,int y){    if(dep[x]<dep[y])swap(x,y);    int fet=dep[x]-dep[y];    for(int i=16;i>=0;i--)if(fet>=(1<<i))x=f[x][i],fet-=(1<<i);    if(x==y)return x;    for(int i=16;i>=0;i--)if(f[x][i]!=f[y][i])        x=f[x][i],y=f[y][i];    return f[x][0];}int main(){    freopen("in.txt","r",stdin);    n=read(),q=read();    memset(first,-1,sizeof(first));    for(int i=1;i<n;i++)x=read(),y=read(),z=read(),add(x,y,z),add(y,x,z);    init(1,0);    for(int i=1;i<=q;i++)        x=read(),y=read(),write(dis[x]+dis[y]-2*dis[LCA(x,y)]);    return 0;}