【块状树】【LCA】bzoj1787 [Ahoi2008]Meet 紧急集合

分块LCA什么的，意外地快呢……

就是对询问的3个点两两求LCA，若其中两组LCA相等，则答案为第三者。

然后用深度减一减什么的就求出距离了。

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define maxn 500001
struct Graph
{int v[maxn<<1],first[maxn<<1],next[maxn<<1],w[maxn<<1],en;
void AddEdge(const int &a,const int &b)
{v[++en]=b;next[en]=first[a];first[a]=en;}}G;
int dep[maxn],fa[maxn],top[maxn],siz[maxn],x,y,z,sz,n,m;
int Abs(const int &x){return x<0 ? (-x) : x;}
int Res,Num;char C,CH[12];
inline int R()
{
    Res=0;C='*';
    while(C<'0'||C>'9')C=getchar();
    while(C>='0'&&C<='9'){Res=Res*10+(C-'0');C=getchar();}
    return Res;
}
inline void P(long long x)
{
    Num=0;if(!x){putchar('0');return;}
    while(x>0)CH[++Num]=x%10,x/=10;
    while(Num)putchar(CH[Num--]+48);
}
void makeblock(int cur)
{
    for(int i=G.first[cur];i;i=G.next[i])
      if(G.v[i]!=fa[cur])
        {
          dep[G.v[i]]=dep[cur]+1;
          fa[G.v[i]]=cur;
          if(siz[top[cur]]<sz)
            {
              siz[top[cur]]++;
              top[G.v[i]]=top[cur];
            }
          makeblock(G.v[i]);
        }
}
inline int QLCA(int u,int v)
{
    while(u!=v)
      {
        if(top[u]==top[v])
          {
            if(dep[u]<dep[v]) swap(u,v);
            u=fa[u];
          }
        else
          {
            if(dep[top[u]]<dep[top[v]]) swap(u,v);
            u=fa[top[u]];
          }
      }
    return u;
}
int main()
{
    n=R();m=R();
    for(int i=1;i<n;i++)
      {
        x=R();y=R();
        G.AddEdge(x,y);
        G.AddEdge(y,x);
      }
    for(int i=1;i<=n;i++) {top[i]=i; siz[i]=1;}
    sz=sqrt((double)n*5.0); makeblock(1);
    for(;m>0;m--)
      {
        x=R();y=R();z=R();
        int f1=QLCA(x,y),f2=QLCA(x,z),f3=QLCA(y,z);
        if(f1==f2)
          {
            int f4=QLCA(x,f3); P(f3); putchar(' ');
            P(Abs(dep[f3]-dep[y])+Abs(dep[f3]-dep[z])
            +Abs(dep[f4]-dep[x])+Abs(dep[f4]-dep[f3])); puts("");
          }
        else if(f1==f3)
          {
            int f4=QLCA(y,f2); P(f2); putchar(' ');
            P(Abs(dep[f2]-dep[x])+Abs(dep[f2]-dep[z])
            +Abs(dep[f4]-dep[y])+Abs(dep[f4]-dep[f2])); puts("");
          }
        else
          {
            int f4=QLCA(z,f1); P(f1); putchar(' ');
            P(Abs(dep[f1]-dep[x])+Abs(dep[f1]-dep[y])
            +Abs(dep[f4]-dep[z])+Abs(dep[f4]-dep[f1])); puts("");
          }
      }
    return 0;
}