bzoj1787 [Ahoi2008]Meet 紧急集合 树上倍增
来源:互联网 发布:慈溪楼市每日成交数据 编辑:程序博客网 时间:2024/05/22 04:48
树上给出任意三个点,找出一个点使它到三个点的距离之和最小
结论:枚举任意两个点的lca,求出距离取最小值 伪证一发:
首先考虑x,y两个点的情况,显然lca(x,y)是最小值,加入第三个点z后,可以看成a,b,已将在lca(a,b)的一个点,向上的花费为2.而z的花费为1.为了是花费最小.考虑让z跳到lca(a,b),枚举lca(a,b)即可
PS:倍增难得1A的题(bzoj1977WA了一天了)
#include<iostream>#include<cstdio>#include<cmath>#include<cstring>#include<algorithm>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fod(i,a,b) for(int i=a;i>=b;i--)using namespace std;const int N=5e5+10;int dep[N],f[N][21],dis[N],n,m,len=0,last[N];struct Edge{int to,next,val;Edge(int to=0,int next=0,int val=0):to(to),next(next),val(val){}}e[N<<1];void add_edge(int u,int v,int w){e[++len]=Edge(v,last[u],w);last[u]=len;}void dfs(int u,int fa){ for(int i=1;i<=20;i++) { if(dep[u]<(1<<i))break; f[u][i]=f[f[u][i-1]][i-1]; } for(int i=last[u];i;i=e[i].next) { int id=e[i].to; if(id==fa)continue; dep[id]=dep[u]+1;dis[id]=dis[u]+e[i].val; f[id][0]=u; dfs(id,u); }}int lca(int x,int y){ if(dep[x]<dep[y])swap(x,y); int d=dep[x]-dep[y]; for(int i=0;i<=20;i++) if(d&(1<<i)) x=f[x][i]; for(int i=20;i>=0;i--) if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i]; if(x==y)return x; else return f[x][0];}int cal(int x,int y,int z){return dis[x]+dis[y]-dis[z]*2;}void solve(int x,int y,int z){ int lc1=lca(x,y),lc2=lca(x,z),lc3=lca(y,z); int dis1=cal(x,y,lc1)+cal(lc1,z,lca(lc1,z)),dis2=cal(x,z,lc2)+cal(y,lc2,lca(y,lc2)),dis3=cal(y,z,lc3)+cal(x,lc3,lca(x,lc3)); int mn=min(dis1,min(dis2,dis3)); if(dis1==mn)printf("%d %d\n",lc1,dis1); else if(dis2==mn)printf("%d %d\n",lc2,dis2); else if(dis3==mn)printf("%d %d\n",lc3,dis3); }int main(){ scanf("%d%d",&n,&m); for(int u,v,i=1;i<n;i++) { scanf("%d%d",&u,&v); add_edge(u,v,1);add_edge(v,u,1); }dfs(1,0); for(int x,y,z,i=1;i<=m;i++) { scanf("%d%d%d",&x,&y,&z); solve(x,y,z); } return 0;}
阅读全文
0 0
- bzoj1787 [Ahoi2008]Meet 紧急集合 树上倍增
- [bzoj1787][Ahoi2008]Meet 紧急集合 倍增LCA
- [BZOJ1787][Ahoi2008]Meet 紧急集合
- bzoj1787 [Ahoi2008]Meet 紧急集合
- BZOJ1787 [Ahoi2008]Meet 紧急集合
- bzoj1787 [Ahoi2008]Meet 紧急集合
- BZOJ1787: [Ahoi2008]Meet 紧急集合
- [BZOJ1787][Ahoi2008]Meet 紧急集合
- 【Ahoi2008】【lca】【bzoj1787】Meet 紧急集合
- [BZOJ1787][AHOI2008]Meet 紧急集合(LCA)
- BZOJ1787: [Ahoi2008]Meet 紧急集合(LCA)
- 11.30 bzoj1787 [Ahoi2008]Meet 紧急集合
- BZOJ1787 [Ahoi2008]Meet 紧急集合 【LCA】
- [BZOJ1787][Ahoi2008]Meet 紧急集合&&[BZOJ1832][AHOI2008]聚会
- [bzoj1787][Ahoi2008]Meet 紧急集合&&[bzoj1832][AHOI2008]聚会
- BZOJ1787: [Ahoi2008]Meet 紧急集合&1832: [AHOI2008]聚会
- bzoj1787: [Ahoi2008]Meet 紧急集合&&bzoj1832: [AHOI2008]聚会
- 【BZOJ1787】【Ahoi2008】Meet 紧急集合 LCA、双倍经验
- 基于领域的社会化推荐
- 网络干货,无论是运维还是开发都要知道的网络知识系列之(一)
- jsp取得绝对路径
- 哈哈
- bzoj3408 [Usaco2009 Oct]Heat Wave 热浪(Dijkstra裸题)
- bzoj1787 [Ahoi2008]Meet 紧急集合 树上倍增
- Access 与C++编程
- qt静态链接库心得
- 正则表达式
- Elasticsearch 架构
- BZOJ1801 Ahoi2009 chess 中国象棋
- 关联容器
- C. The Tag Game【dfs】
- 佳能iR2018一体机驱动下载 v14.2官方版