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的深度差，加起来就是结果。
代码：
1 #include<cstdio> 2 #include<iostream> 3 #include<algorithm> 4 #define N 500010 5 #define M 1000010 6 using namespace std; 7 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; 8 void add(int false_from,int false_to){ 9 next[++num]=head[false_from]; 10 to[num]=false_to; 11 head[false_from]=num; 12 } 13 void dfs1(int x){ 14 size[x]=1; 15 deep[x]=deep[father[x]]+1; 16 for(int i=head[x];i;i=next[i]) 17 if(father[x]!=to[i]){ 18 father[to[i]]=x; 19 dfs1(to[i]); 20 size[x]+=size[to[i]]; 21 } 22 } 23 void dfs2(int x){ 24 int mmax=0; 25 if(!top[x]) 26 top[x]=x; 27 for(int i=head[x];i;i=next[i]) 28 if(father[x]!=to[i]&&size[to[i]]>size[mmax]) 29 mmax=to[i]; 30 if(mmax){ 31 top[mmax]=top[x]; 32 dfs2(mmax); 33 } 34 for(int i=head[x];i;i=next[i]) 35 if(father[x]!=to[i]&&to[i]!=mmax) 36 dfs2(to[i]); 37 } 38 int lca(int x,int y){ 39 while(top[x]!=top[y]){ 40 if(deep[top[x]]<deep[top[y]]) 41 swap(x,y); 42 x=father[top[x]]; 43 } 44 if(deep[x]<deep[y])return x; 45 return y; 46 } 47 int main(){ 48 scanf("%d%d",&n,&m); 49 for(int i=1;i<n;++i){ 50 scanf("%d%d",&a,&b); 51 add(a,b); 52 add(b,a); 53 } 54 dfs1(1); 55 dfs2(1); 56 for(int i=1;i<=m;++i){ 57 scanf("%d%d%d",&a,&b,&c); 58 ans=0; 59 c1=lca(a,b); 60 c2=lca(a,c); 61 c3=lca(b,c); 62 if(c1==c2&&c2==c3){ 63 ans=abs(deep[c1]-deep[a])+abs(deep[c1]-deep[b])+abs(deep[c1]-deep[c]); 64 printf("%d %d\n",c1,ans); 65 } 66 else 67 if(c1==c2){ 68 ans+=deep[b]+deep[c]-deep[c3]*2; 69 ans+=deep[c3]+deep[a]-deep[c1]*2; 70 printf("%d %d\n",c3,ans); 71 } 72 else 73 if(c1==c3){ 74 ans+=deep[a]+deep[c]-deep[c2]*2; 75 ans+=deep[c2]+deep[b]-deep[c1]*2; 76 printf("%d %d\n",c2,ans); 77 } 78 else{ 79 ans+=deep[a]+deep[b]-deep[c1]*2; 80 ans+=deep[c1]+deep[c]-deep[c2]*2; 81 printf("%d %d\n",c1,ans); 82 } 83 } 84 return 0; 85 }
