首页 > 代码库 > BZOJ1787: [Ahoi2008]Meet 紧急集合
BZOJ1787: [Ahoi2008]Meet 紧急集合
1787: [Ahoi2008]Meet 紧急集合
Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 1155 Solved: 488
[Submit][Status]
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
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
题解:
求3个点的LCA?样例过不了,那就猜吧,估计是其中某两个点的LCA,暴力枚举一下就行了。
有人说:若lca(x,y)=lca(y,z) 那么最后集合的地点就是 lca(x,z)?
不懂,求解释。
不过考场上暴力就可以,只是常数大一点。
唉,连个LCA都能写错。。。
代码:
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cmath> 4 #include<cstring> 5 #include<algorithm> 6 #include<iostream> 7 #include<vector> 8 #include<map> 9 #include<set>10 #include<queue>11 #include<string>12 #define inf 100000000013 #define maxn 500000+100014 #define maxm 500+10015 #define eps 1e-1016 #define ll long long17 #define pa pair<int,int>18 #define for0(i,n) for(int i=0;i<=(n);i++)19 #define for1(i,n) for(int i=1;i<=(n);i++)20 #define for2(i,x,y) for(int i=(x);i<=(y);i++)21 using namespace std;22 inline int read()23 {24 int x=0,f=1;char ch=getchar();25 while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}26 while(ch>=‘0‘&&ch<=‘9‘){x=10*x+ch-‘0‘;ch=getchar();}27 return x*f;28 }29 int n,m,tot,ans,pos,head[maxn],dep[maxn],f[maxn][22];30 struct edge{int next,go;}e[2*maxn];31 void insert(int x,int y)32 {33 e[++tot].go=y;e[tot].next=head[x];head[x]=tot;34 }35 void dfs(int x)36 {37 for1(i,20)38 if(dep[x]>=(1<<i))39 f[x][i]=f[f[x][i-1]][i-1];40 else break;41 for(int i=head[x],y;i;i=e[i].next)42 if(dep[y=e[i].go]==-1) 43 {44 dep[y]=dep[x]+1;45 f[y][0]=x;46 dfs(y);47 }48 }49 int lca(int x,int y)50 {51 if(dep[x]<dep[y])swap(x,y);52 int t=dep[x]-dep[y];53 for0(i,20)54 if(t&(1<<i))55 x=f[x][i];56 if(x==y)return x;57 for(int i=20;i>=0;i--)58 if(f[x][i]!=f[y][i])59 x=f[x][i],y=f[y][i];60 return f[x][0]; 61 }62 void update(int x,int y,int z)63 {64 int t=lca(x,y),tt=lca(t,z),ttt=dep[x]+dep[y]+dep[z]-dep[t]-2*dep[tt];65 if(ttt<ans)ans=ttt,pos=t;66 }67 int main()68 {69 freopen("input.txt","r",stdin);70 freopen("output.txt","w",stdout);71 n=read();m=read();72 int x,y,z;73 for1(i,n-1)x=read(),y=read(),insert(x,y),insert(y,x);74 memset(dep,-1,sizeof(dep));75 dep[1]=0;76 dfs(1);77 while(m--)78 {79 x=read();y=read();z=read();80 ans=inf;pos=0;81 update(x,y,z);82 update(y,z,x);83 update(z,x,y);84 printf("%d %d\n",pos,ans);85 }86 return 0;87 }
BZOJ1787: [Ahoi2008]Meet 紧急集合
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。