3242: [Noi2013]快餐店
来源:互联网 发布:成都安吉斯主机编程 编辑:程序博客网 时间:2024/06/11 14:55
题意
就是给你一个环基树,然后要你找出一个点(可以在边上),使得这个点到所有图上的点距离最大值最小
题解
先是有一个结论:这个点肯定是在这“树”上的直径的中点,答案就是这个直径的长度的一半
所以其实我们就是要找一条最长的链
答案出现的情况有两种:
1.不经过环上的任意一条边,然后这个扫一下就可以了
2.要经过环,那么这个怎么算呢?
我们知道,对于一个方案,他肯定不会吧环遍历完的,也就是说至少有1条边不会被用到,我们就尝试枚举这一条边,看一下断掉这条边后的最大值取最小就是答案了。
为什么是最大值最小呢。对于最大值,我们是保证他是一条最长链,然后最小值是因为路径不会舍近求远
给环上的点标号1~k
设pre1[i]表示1~i的树上的所有点到1的最大值,suf1[i]表示i~k的树上的所有点到k的最大值
设pre2[i]表示1~i的任意两棵树之间组成的最长链的最大值,suf2[i]类似,表示的是i~k
这个正着和倒着分别扫一遍就可以求出来,具体看看代码吧。。
然后断开i~i+1的边就是
max(pre1[i]+suf1[i+1]+(1和k之间的边长),pre2[i],suf2[i+1])
然后最后和不经过环的方案取一个最大值
这里有一个细节要注意一下:
就是说在取ans的时候,顺序是有要求的
下面本该是这么写的:
for (LL u=2;u<=cnt;u++){ LL tt=0; tt=max(tt,pr2[u-1]); tt=max(tt,suf2[u]); tt=max(tt,pr1[u-1]+suf1[u]+cval[cnt]); ans=min(ans,tt);}for (int u=1;u<=n;u++) ans=max(ans,d[u]);
但是一开始我偷懒,写成了下面这种
ans=max(ans,d[1]);for (LL u=2;u<=cnt;u++){ LL tt=0; tt=max(tt,pr2[u-1]); tt=max(tt,suf2[u]); tt=max(tt,pr1[u-1]+suf1[u]+cval[cnt]); ans=min(ans,tt); ans=max(ans,d[u]);}
我调了很久都没有调出来,最后才发现这样是有区别的
如果你在里面取了最大值,可能他的最小值就不见了,就会影响我们的答案
看来我比较菜,搞成了min和max都有交换律。。
CODE:
#include<cstdio>#include<cstring>#include<algorithm>#include<iostream>#include<stack>using namespace std;typedef long long LL;const LL N=100005;LL n;struct qq{ LL x,y,c,last;}e[N*2];LL num,last[N];void init (LL x,LL y,LL c){ num++; e[num].x=x;e[num].y=y;e[num].c=c; e[num].last=last[x]; last[x]=num;}bool vis[N];//是不是在栈里面 bool in[N];//这个点是不是在环上面 LL cir[N],cnt;//环 LL cval[N],c1[N];//这个点在环上面与父亲边的权值stack<LL> sta;bool find_cir (LL x,LL fa,LL c)//当前的节点 父亲 与父亲边的值{ if (vis[x]==true)//环 { cir[++cnt]=x;in[x]=true;cval[cnt]=c; while (!sta.empty()&&sta.top()!=x) { LL x=sta.top();sta.pop(); cir[++cnt]=x; cval[cnt]=c1[x]; in[x]=true; } return true; } vis[x]=true; sta.push(x);c1[x]=c; for (LL u=last[x];u!=-1;u=e[u].last) { LL y=e[u].y; if (y==fa) continue; if (find_cir(y,x,e[u].c)) return true; } sta.pop(); return false;}LL dis[N],d[N];//这颗子树里面的最长链 子树里面的最长距离 LL dfs (LL x,LL fa){ LL re=0;//子树中两点的最远距离 for (LL u=last[x];u!=-1;u=e[u].last) { LL y=e[u].y; if (y==fa||in[y]==true)//我们的目的是遍历这个森林 continue; re=max(re,dfs(y,x)); LL tmp=dis[y]+e[u].c; re=max(re,dis[x]+tmp);//和之前的链并起来 dis[x]=max(dis[x],tmp); } return re;}LL pr1[N],pr2[N];LL suf1[N],suf2[N];int main(){ num=0;memset(last,-1,sizeof(last)); scanf("%lld",&n); for (LL u=1;u<=n;u++) { LL x,y,c; scanf("%lld%lld%lld",&x,&y,&c); init(x,y,c);init(y,x,c); } memset(in,false,sizeof(in)); memset(vis,false,sizeof(vis)); find_cir(1,0,0);//寻找环 cval[0]=cval[cnt];//因为是环嘛 for (LL u=1;u<=cnt;u++)//现在这棵数被环分成了很多个森林 d[u]=dfs(cir[u],0); LL t1=-cval[0],t2=-cval[0]; for (LL u=1;u<=cnt;u++) { LL nd=dis[cir[u]]; t1=t1+cval[u-1];t2=t2+cval[u-1]; pr1[u]=max(pr1[u-1],nd+t1); pr2[u]=max(pr2[u-1],nd+t2); t2=max(t2,nd); } t1=-cval[0],t2=-cval[0]; for (LL u=cnt;u>=1;u--) { LL nd=dis[cir[u]]; t1=t1+cval[u];t2=t2+cval[u]; suf1[u]=max(suf1[u+1],nd+t1); suf2[u]=max(suf2[u+1],nd+t2); t2=max(t2,nd); } LL ans=suf2[1]; for (LL u=2;u<=cnt;u++) { LL tt=0; tt=max(tt,pr2[u-1]); tt=max(tt,suf2[u]); tt=max(tt,pr1[u-1]+suf1[u]+cval[cnt]); ans=min(ans,tt); } for (int u=1;u<=n;u++) ans=max(ans,d[u]); printf("%.1lf\n",(double)ans/2.0); return 0;}
阅读全文
0 0
- 3242: [Noi2013]快餐店
- NOI2013 快餐店
- [Noi2013]快餐店
- bzoj 3242: [Noi2013]快餐店 dfs&递推
- BZOJ3242: [Noi2013]快餐店
- bzoj3242: [Noi2013]快餐店
- [线段树优化 DP] BZOJ 3242 [Noi2013]快餐店
- NOI2013快餐店【图上找环+线段树】
- 【树DP+基环树】[NOI2013]快餐店
- [BZOJ3242][Noi2013]快餐店 && 环套树+线段树
- bzoj3242: [Noi2013]快餐店 树形dp+线段树
- [BZOJ3242][NOI2013]快餐店-基环树-动态规划
- [BZOJ3242][Noi2013]快餐店(树形dp+线段树)
- Bzoj-3242 快餐店(环套树)
- 3066 快餐店
- NOI2013行纪
- NOI2013 书法家
- 上海快餐店一览
- session.cpp
- Vue中的$set的使用
- 支付宝第三方登录
- session.h
- Python网络爬虫与信息提取(一)
- 3242: [Noi2013]快餐店
- uva 225 Golygons
- 和研究生博士学长聊天有感
- docker容器 最简单的java演示实例
- 判断数独是否合法-LintCode
- C++继承[详细]
- 文件管理,系统管理,系统安全常用指令整理
- 数据结构——图的广度优先搜索模板
- some tips about python One