BZOJ 1787 Ahoi2008 Meet紧急集合

1787: [Ahoi2008]Meet 紧急集合

Time Limit: 20 Sec  Memory Limit: 162 MB
Submit: 3590  Solved: 1642
[Submit][Status][Discuss]

Description

Input

Output

Sample Input

6 4
1 2
2 3
2 4
4 5
5 6
4 5 6
6 3 1
2 4 4
6 6 6

Sample Output

5 2
2 5
4 1
6 0

HINT

Source

Day1

设自己在x点，小可可在y点，小可可的朋友在z点，那么设xy为x和y的lca，xz为x和z的lca，yz为y和z的lca

那么汇聚点要么在xy要么在yz要么在xz

三个进行比较即可

代码如下

#include <bits/stdc++.h>
using namespace std;
inline int read()
{
 int x=0;int f=1;char ch=getchar();
 while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
 while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
 return x*f;
}
const int MAXN=1e6+10;
int fa[MAXN][25],linkk[MAXN*2],depth[MAXN],len,n,m;
struct node
{
 int y,next;
}e[MAXN*2];
inline void insert(int xx,int yy)
{
 e[++len].y=yy;e[len].next=linkk[xx];linkk[xx]=len;
}
void dfs(int father,int st,int dep)
{
 fa[st][0]=father;depth[st]=dep;
 for(int i=linkk[st];i;i=e[i].next)
 {
  if(e[i].y!=father) dfs(st,e[i].y,dep+1);
 }
}
void getfather()
{
  for(int i=1;i<=20;i++)
  {
   for(int j=1;j<=n;j++)
   {
    if(fa[j][i-1]!=0) fa[j][i]=fa[fa[j][i-1]][i-1];
   }
  }
}
int LCA(int x,int y)
{
 if(x==y) return x;
 if(depth[x]<depth[y]) swap(x,y);
 for(int i=20;i>=0;i--)
 {
  if(depth[x]-(1<<i)>=depth[y]) x=fa[x][i];
 }
 if(x==y) return x;
 for(int i=20;i>=0;i--)
 {
  if(fa[x][i]!=fa[y][i]) {x=fa[x][i];y=fa[y][i];}
 }
 return fa[x][0];
}
void find()
{
 int xx=read();int yy=read();int zz=read();
 int xy=LCA(xx,yy);int xz=LCA(xx,zz);int yz=LCA(yy,zz);
 //cout<<xy<<' '<<xz<<' '<<yz<<endl;
 int ans1;int ans2=10000000;
 int xyz=LCA(xy,zz);
 if(depth[xx]+depth[yy]-depth[xy]-depth[xyz]*2+depth[zz]<ans2)
 {
  ans2=depth[xx]+depth[yy]-depth[xy]-depth[xyz]*2+depth[zz];
  ans1=xy;
 }
 xyz=LCA(xz,yy);
 if(depth[xx]+depth[zz]-depth[xz]-depth[xyz]*2+depth[yy]<ans2)
 {
  ans2=depth[xx]+depth[zz]-depth[xz]-depth[xyz]*2+depth[yy];
  ans1=xz;
 }
 xyz=LCA(yz,xx);
 if(depth[yy]+depth[zz]-depth[yz]-depth[xyz]*2+depth[xx]<ans2)
 {
  ans2=depth[yy]+depth[zz]-depth[yz]-depth[xyz]*2+depth[xx];
  ans1=yz;
 }
 printf("%d %d
",ans1,ans2);
}
void init()
{
 n=read();m=read();
 len=0;
 for(int i=1;i<n;i++)
 {
  int xx=read();
  int yy=read();
  insert(xx,yy);
  insert(yy,xx);
 }
 dfs(0,1,1);
 getfather();
 for(int i=1;i<=m;i++)
 {
  find();
 }
}
int main()
{
 init();
 return 0;
}