[AHOI2008]紧急集合 / 聚会 题解

题目

欢乐岛上有个非常好玩的游戏，叫做“紧急集合”。在岛上分散有N个等待点，有N-1条道路连接着它们，每一条道路都连接某两个等待点，且通过这些道路可以走遍所有的等待点，通过道路从一个点到另一个点要花费一个游戏币。参加游戏的人三人一组，开始的时候，所有人员均任意分散在各个等待点上（每个点同时允许多个人等待），每个人均带有足够多的游戏币（用于支付使用道路的花费）、地图（标明等待点之间道路连接的情况）以及对话机（用于和同组的成员联系）。当集合号吹响后，每组成员之间迅速联系，了解到自己组所有成员所在的等待点后，迅速在N个等待点中确定一个集结点，组内所有成员将在该集合点集合，集合所用花费最少的组将是游戏的赢家。小可可和他的朋友邀请你一起参加这个游戏，由你来选择集合点，聪明的你能够完成这个任务，帮助小可可赢得游戏吗？
输入格式：
第一行两个正整数N和M（N<=500000，M<=500000），之间用一个空格隔开。分别表示等待点的个数（等待点也从1到N进行编号）和获奖所需要完成集合的次数。 随后有N-1行，每行用两个正整数A和B，之间用一个空格隔开，表示编号为A和编号为B的等待点之间有一条路。 接着还有M行，每行用三个正整数表示某次集合前小可可、小可可的朋友以及你所在等待点的编号。
输出格式
一共有M行，每行两个数P,C，用一个空格隔开。其中第i行表示第i次集合点选择在编号为P的等待点，集合总共的花费是C个游戏币。
输入输出样例
输入样例#1：
6 4
1 2
2 3
2 4
4 5
5 6
4 5 6
6 3 1
2 4 4
6 6 6
输出样例#1：
5 2
2 5
4 1
6 0

分析

一道比较水的一道图论题，我们先计算只有两个点的情况，$dis=d[u]+d[v]-2*d[lca(u,v)]$，那么来说有三个点的情况，我们来说三个点的情况，可以先求出两点的lca，再求lca和另外一个点的lca，然而第一步有三个lca，略有麻烦，通过手玩样例，我们发现d最大的点就是答案，所以判断那个点最深就可以做了。
上代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
int ver[1000010],nxt[1000010],hed[500010],tot,f[500010][25],d[500010],t,n,T;
queue<int> q;
inline void add(int x,int y){
    ver[++tot]=y,nxt[tot]=hed[x],hed[x]=tot;
}
inline void bfs(){
    q.push(1);
    d[1]=1;
    while(q.size()){
        int x=q.front();
        q.pop();
        for(int i=hed[x];i;i=nxt[i]){
            int y=ver[i];
            if(d[y])
                continue;
            d[y]=d[x]+1;
            f[y][0]=x;
            for(int j=1;j<=t;j++)
                f[y][j]=f[f[y][j-1]][j-1];
            q.push(y);
        }
    }
}
inline int lca(int x,int y){
    if(d[x]>d[y])
        swap(x,y);
    for(int i=t;i>=0;i--){
        if(d[f[y][i]]>=d[x])
            y=f[y][i];
    }
    if(x==y)
        return x;
    for(int i=t;i>=0;i--)
        if(f[x][i]!=f[y][i])
            x=f[x][i],y=f[y][i];
    return f[x][0];
}
inline int readn(){
    int x=0;
    char ch=getchar();
    while(!isdigit(ch))
        ch=getchar();
    while(isdigit(ch))
    {
        x=(x<<1)+(x<<3)+(ch^'0');
        ch=getchar();
    }
    return x;
}
inline bool pd(int x,int y,int z,int fa){
    if((d[x]>max(d[y],d[z]))||(d[x]<d[fa]))
        return 0;
    if(d[x]>=d[y]){
        if(x==y)
            return 1;
        for(int i=t;i>=0;i--)
            if(d[x]<=d[f[z][i]])
                z=f[z][i];
        if(z==x)
            return 1;
        return 0;
    }
    if(d[x]>=d[z]){
        if(x==z)
            return 1;
        for(int i=t;i>=0;i--)
            if(d[x]<=d[f[y][i]])
                y=f[y][i];
        if(x==y)
            return 1;
        return 0;
    }
    for(int i=t;i>=0;i--)
        if(d[x]<=d[f[y][i]])
            y=f[y][i];
    for(int i=t;i>=0;i--)
        if(d[x]<=d[f[z][i]])
            z=f[z][i];
    if((x==y)||(x==z))
        return 1;
    return 0;
}
int main(){
    n=readn(),T=readn();
    t=(int)(log(n)/log(2))+1;
    for(int i=1,x,y;i<n;i++){
        x=readn(),y=readn();
        add(x,y),add(y,x);
    }
    bfs();
    while(T--){
        int a,b,c;
        a=readn(),b=readn(),c=readn();
        int qaq=lca(a,b),qbq=lca(b,c),qcq=lca(a,c);
        if(d[qaq]>=max(d[qbq],d[qcq])){
            printf("%d %d\n",qaq,d[a]+d[b]-d[qaq]+d[c]-2*d[lca(qaq,c)]);
            continue;
        }
        if(d[qbq]>=d[qcq])
            printf("%d %d\n",qbq,d[a]+d[b]+d[c]-d[qbq]-2*d[lca(qbq,a)]);
        else
            printf("%d %d\n",qcq,d[a]+d[b]+d[c]-d[qcq]-2*d[lca(qcq,b)]);
    }
    return 0;
}