BZOJ1787 [Ahoi2008]Meet 紧急集合（洛谷P4281、BZOJ1832）

LCA

BZOJ1787
BZOJ1832
洛谷P4281

划水

求出两两点对之间的LCA，显然其中有两个一定相等，而且这个点必然是深度更小的那个。

那么集结点就是深度大，也就是不一样的那个，化简后可得答案为这三个点的深度之和减去它们LCA的深度之和。

代码：

#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 500005
#define F inline
using namespace std;
struct edge{ int nxt,to; }ed[N<<1];
int n,m,k,h[N],dep[N],fa[N][20];
F char readc(){
    static char buf[100000],*l=buf,*r=buf;
    if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
    return l==r?EOF:*l++;
}
F int _read(){
    int x=0; char ch=readc();
    while (!isdigit(ch)) ch=readc();
    while (isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=readc();
    return x;
}
F void writec(int x){ if (x>9) writec(x/10); putchar(x%10+48); }
#define addedge(x,y) ed[++k]=(edge){h[x],y},h[x]=k
void dfs(int x){
    dep[x]=dep[fa[x][0]]+1;
    for (int i=h[x],v;i;i=ed[i].nxt)
        if ((v=ed[i].to)!=fa[x][0])
            fa[v][0]=x,dfs(v);
}
F void Make(){
    for (int j=1;j<20;j++)
        for (int i=1;i<=n;i++)
            fa[i][j]=fa[fa[i][j-1]][j-1];
}
F int lca(int x,int y){
    if (dep[x]<dep[y]) swap(x,y);
    for (int j=19;~j;j--)
        if (dep[fa[x][j]]>=dep[y]) x=fa[x][j];
    if (x==y) return x;
    for (int j=19;~j;j--)
        if (fa[x][j]!=fa[y][j])
            x=fa[x][j],y=fa[y][j];
    return fa[x][0];
}
int main(){
    n=_read(),m=_read();
    for (int i=1,x,y;i<n;i++)
        x=_read(),y=_read(),addedge(x,y),addedge(y,x);
    dfs(1),Make();
    while (m--){
        int x=_read(),y=_read(),z=_read();
        int a=lca(x,y),b=lca(y,z),c=lca(x,z),ans1,ans2=a^b^c;
        ans1=dep[x]+dep[y]+dep[z]-dep[a]-dep[b]-dep[c];
        writec(ans2),putchar(' '),writec(ans1),puts("");
    }
    return 0;
}