BZOJ1787: [Ahoi2008]Meet 紧急集合&1832: [AHOI2008]聚会

以后根的dep要定1…不然LCA会炸…

因为好像没有这个类型只好随便分了一个类

树上的三个点找到一个点使得他们到这个点的距离和最小

A,B,C三点，对于A,B两点找一个点到他们的距离和最小，这个点一定在他们的路径上，对于A,B,C三点找一个点，这个点也一定在A,B的路径上，因为让一个人多走若干步一定比让两个人多走若干步优

求出两两的LCA，讨论一下就好了

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

inline void swap(int &x,int &y){x^=y;y^=x;x^=y;}
inline void read(int &x)
{
    char c;
    while(!((c=getchar())>='0'&&c<='9'));
    x=c-'0';
    while((c=getchar())>='0'&&c<='9') (x*=10)+=c-'0';
}
const int maxn = 510000;
const int maxm = 510000;
const int maxb = 21;

struct edge
{
    int y,nex;
    edge(){}
    edge(const int &_y,const int &_nex){y=_y;nex=_nex;}
}a[maxn<<1]; int len,fir[maxn];
int n,m;
inline void ins(const int &x,const int &y){a[++len]=edge(y,fir[x]);fir[x]=len;}

int fa[maxn][maxb],dep[maxn];

void dfs(const int &x)
{
    for(int i=1;fa[x][i-1]&&i<maxb;i++)
        fa[x][i]=fa[fa[x][i-1]][i-1];
    for(int k=fir[x];k;k=a[k].nex)
    {
        const int y=a[k].y;
        if(y!=fa[x][0])
        {
            fa[y][0]=x; dep[y]=dep[x]+1;
            dfs(y);
        }
    }
}
int LCA(int x,int y)
{
    if(dep[x]<dep[y]) swap(x,y);
    for(int i=maxb-1;i>=0&&dep[x]!=dep[y];i--)
        if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
    if(x==y) return x;
    for(int i=maxb-1;i>=0;i--)
        if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
    return fa[x][0];
}

int main()
{
    read(n); read(m);
    for(int i=1;i<n;i++)
    {
        int x,y; read(x); read(y);
        ins(x,y); ins(y,x);
    }
    dep[1]=1; dfs(1);
    for(int i=1;i<=m;i++)
    {
        int t1,t2,t3; read(t1); read(t2); read(t3);
        int tmp=LCA(t1,t2);
        int tmp2=LCA(t3,tmp);
        int c=dep[t1]+dep[t2]-2*dep[tmp];
        if(tmp2!=tmp) 
            printf("%d %d\n",tmp,dep[tmp]+dep[t3]-2*dep[tmp2]+c);
        else
        {
            int tmp3=LCA(t1,t3);
            int tmp4=LCA(t2,t3);
            if(tmp3==tmp)
                printf("%d %d\n",tmp4,c+dep[t3]-dep[tmp4]);
            else
                printf("%d %d\n",tmp3,c+dep[t3]-dep[tmp3]);
        }
    }

    return 0;
}