【块状树】【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<<],first[maxn<<],next[maxn<<],w[maxn<<],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< ? (-x) : x;}
 int Res,Num;char C,CH[];
 inline int R()
 {
     Res=;C='*';
     while(C<''||C>'')C=getchar();
     while(C>=''&&C<=''){Res=Res*+(C-'');C=getchar();}
     return Res;
 }
 inline void P(long long x)
 {
     Num=;if(!x){putchar('');return;}
     while(x>)CH[++Num]=x%,x/=;
     while(Num)putchar(CH[Num--]+);
 }
 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]+;
           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=;i<n;i++)
       {
         x=R();y=R();
         G.AddEdge(x,y);
         G.AddEdge(y,x);
       }
     for(int i=;i<=n;i++) {top[i]=i; siz[i]=;}
     sz=sqrt((double)n*5.0); makeblock();
     for(;m>;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 ;
 }