BZOJ——1787: [Ahoi2008]Meet 紧急集合

http://www.lydsy.com/JudgeOnline/problem.php?id=1787

题目描述

输入

输出

样例输入

6 4
1 2
2 3
2 4
4 5
5 6
4 5 6
6 3 1
2 4 4
6 6 6

样例输出

5 2
2 5
4 1
6 0

提示

易发现：三个点两两求出LCA，一定至少有两个LCA相等

　　　　若只有两个相等，那聚集点就是剩下的LCA

　　　　若三个相等，那聚集点就是LCA

#include <algorithm>
#include <cstdio>
#include <vector>

using namespace std;

const int N(1e5*5+15);
vector<int>vec[N];
int n,m,x,y,z;
int anspoint,anssum;
int deep[N],dad[N],size[N],top[N];

void DFS(int x)
{
    size[x]=1; deep[x]=deep[dad[x]]+1;
    for(int i=0;i<vec[x].size();i++)
        if(dad[x]!=vec[x][i])
        {
            dad[vec[x][i]]=x;
            DFS(vec[x][i]);
            size[x]+=size[vec[x][i]];
        }
}

void DFS_(int x)
{
    int t=0; if(!top[x]) top[x]=x;
    for(int i=0;i<vec[x].size();i++)
        if(vec[x][i]!=dad[x]&&size[t]<size[vec[x][i]]) t=vec[x][i];
    if(t) top[t]=top[x],DFS_(t);
    for(int i=0;i<vec[x].size();i++)
        if(vec[x][i]!=dad[x]&&t!=vec[x][i]) DFS_(vec[x][i]);
}

int LCA(int x,int y)
{
    while(top[x]!=top[y])
    {
        if(deep[top[x]]<deep[top[y]]) swap(x,y);
        x=dad[top[x]];
    }
    if(deep[x]>deep[y]) swap(x,y);
    return x;
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<n;i++)
    {
        scanf("%d%d",&x,&y);
        vec[x].push_back(y);
        vec[y].push_back(x);
    }
    DFS(1); DFS_(1);
    for(;m;m--)
    {
        scanf("%d%d%d",&x,&y,&z);
        anspoint=LCA(x,y)^LCA(x,z)^LCA(y,z);
        anssum =deep[x]+deep[anspoint]-(deep[LCA(x,anspoint)]<<1);
        anssum+=deep[y]+deep[anspoint]-(deep[LCA(y,anspoint)]<<1);
        anssum+=deep[z]+deep[anspoint]-(deep[LCA(z,anspoint)]<<1);
        printf("%d %d
",anspoint,anssum);
    }
    return 0;
}