bzoj 1787 Meet 紧急集合

Meet 紧急集合

这个题是在脖子oj（清北某奆佬给起的名字）八中oj（大视野在线评测）上的。

给出bzoj链接。

这个题还是求最近公共祖先的问题。

而该题不同于别的题，它是需要求三个点的最近公共祖先。

我们就需要求出三个点两两之间的LCA。

而这三个LCA之间，必有两个是相同的。

如果两个点相同，那另一点就是那三个点的LCA。

如果三个点都相同，那么该点就是那三个点的LCA。

最后还需要统计走过的边的长度。

三个点都相同的情况很好搞，就是计算三个点与那个LCA的深度差，加起来就是答案。

但是两个点就难一点。如果lca_a_b与lca_a_c相同，就需要求出lca_b_c到a,b,c的长度，就可以转化成b,c到它们的最近公共祖先lca_b_c的长度，和lca_b_c,a到它们的最近公共祖先lca_a_b(或lca_a_c)的长度。

那么就先求出lca_b_c与b和c的深度差，再求出lca_a_b（或lca_a_c）与lca_b_c和a的深度差，加起来就是结果。

代码：

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #define N 500010
 #define M 1000010
 using namespace std;
 int next[M],to[M],head[N],num,size[N],deep[N],father[N],top[N],n,m,a,b,c,c1,c2,c3,ans;
 void add(int false_from,int false_to){
     next[++num]=head[false_from];
     to[num]=false_to;
     head[false_from]=num;
 }
 void dfs1(int x){
     size[x]=;
     deep[x]=deep[father[x]]+;
     for(int i=head[x];i;i=next[i])
         if(father[x]!=to[i]){
             father[to[i]]=x;
             dfs1(to[i]);
             size[x]+=size[to[i]];
         }
 }
 void dfs2(int x){
     int mmax=;
     if(!top[x])
         top[x]=x;
     for(int i=head[x];i;i=next[i])
         if(father[x]!=to[i]&&size[to[i]]>size[mmax])
             mmax=to[i];
     if(mmax){
         top[mmax]=top[x];
         dfs2(mmax);
     }
     for(int i=head[x];i;i=next[i])
         if(father[x]!=to[i]&&to[i]!=mmax)
             dfs2(to[i]);
 }
 int lca(int x,int y){
     while(top[x]!=top[y]){
         if(deep[top[x]]<deep[top[y]])
             swap(x,y);
         x=father[top[x]];
     }
     if(deep[x]<deep[y])return x;
     return y;
 }
 int main(){
     scanf("%d%d",&n,&m);
     for(int i=;i<n;++i){
         scanf("%d%d",&a,&b);
         add(a,b);
         add(b,a);
     }
     dfs1();
     dfs2();
     for(int i=;i<=m;++i){
         scanf("%d%d%d",&a,&b,&c);
         ans=;
         c1=lca(a,b);
         c2=lca(a,c);
         c3=lca(b,c);
         if(c1==c2&&c2==c3){
             ans=abs(deep[c1]-deep[a])+abs(deep[c1]-deep[b])+abs(deep[c1]-deep[c]);
             printf("%d %d\n",c1,ans);
         }
         else
             if(c1==c2){
                 ans+=deep[b]+deep[c]-deep[c3]*;
                 ans+=deep[c3]+deep[a]-deep[c1]*;
                 printf("%d %d\n",c3,ans);
             }
         else
             if(c1==c3){
                 ans+=deep[a]+deep[c]-deep[c2]*;
                 ans+=deep[c2]+deep[b]-deep[c1]*;
                 printf("%d %d\n",c2,ans);
             }
         else{
             ans+=deep[a]+deep[b]-deep[c1]*;
             ans+=deep[c1]+deep[c]-deep[c2]*;
             printf("%d %d\n",c1,ans);
         }
     }
     return ;
 }