码迷,mamicode.com
首页 > 编程语言 > 详细

[转]01分数规划算法 ACM 二分 Dinkelbach 最优比率生成树 最优比率环

时间:2019-01-24 18:55:24      阅读:322      评论:0      收藏:0      [点我收藏+]

标签:二分法   test   最短路径   二分答案   长度   解决   com   初始   最短路   

01分数规划

 

前置技能

二分思想
最短路算法
一些数学脑细胞?
问题模型1

基本01分数规划问题

给定nn个二元组(valuei,costi)(valuei,costi),valueivaluei是选择此二元组获得的价值(非负),costicosti是选择此二元组付出的代价(非负),设xi(xi∈{0,1})xi(xi∈{0,1})代表第ii个二元组的选与不选,最大(小)化下式
maximize(or minimize)   r=∑valuei?xi∑costi?xi
maximize(or minimize)   r=∑valuei?xi∑costi?xi
下面先说最大化

解决方法

二分法

设rr最大值为r?r?,
r?=∑valuei?xi∑costi?xi
r?=∑valuei?xi∑costi?xi
∑valuei?xi?r??∑costi?xi=0
∑valuei?xi?r??∑costi?xi=0
设一个函数,自变量为rr值,
f(r)=∑valuei?xi?r??∑costi?xi
f(r)=∑valuei?xi?r??∑costi?xi
观察这个函数,假如{xi}{xi}固定,则这个函数就是坐标系中一条直线(y=B?A?xy=B?A?x),每一组{xi}{xi}对应着一条直线,这些直线斜率非正(因为?A=?∑costi?xi≤0?A=?∑costi?xi≤0),纵截距非负(因为B=∑valuei?xi≥0B=∑valuei?xi≥0 ),如图1。

技术分享图片

 

对于每一条直线,当f(r)=0f(r)=0时,横截距就是这一组的rr,那么r?r?就是每条直线横截距的最大值(每组{xi}{xi}对应rr的最大值)如图2。

技术分享图片

 

 

在图中上任取一条垂直xx轴的竖线,
如果存在直线与这条竖线的交点纵坐标为正,那么最优值一定在当前竖线的右侧;
如果所有直线与这条竖线交点纵坐标为负,那么最优值一定在当前竖线的左侧;
如果所有直线与这条竖线交点纵坐标非正且存在直线与这条竖线交点纵坐标为0,那么当前竖线横坐标即为最优值r?r?。

技术分享图片

 

按照这个思想,可以二分答案rr,那么二分时如何进行判断呢?

选择一个rr时需要判断所有f(r)f(r)的最大值是否为0,如果max{f(r)}>0max{f(r)}>0则r<r?r<r?;如果max{f(r)}<0max{f(r)}<0 则 r>r?r>r?。
怎样求max{f(r)}max{f(r)}?
f(r)=∑valuei?xi?r?∑costi?xi=∑(valuei?r?costi)?xi
f(r)=∑valuei?xi?r?∑costi?xi=∑(valuei?r?costi)?xi
二分一个rr时,每个二元组的valuei?r?costivaluei?r?costi 都可以求出,设其为weightiweighti,现在的目标就是找到一组{xi}{xi}使得∑wighti?xi∑wighti?xi最大(即求max{f(r)}max{f(r)})。怎么找到这一组{xi}{xi},或者直接求得max{f(r)}max{f(r)}呢?具体问题具体分析,经常借助最短路算法判断是否存在负环。下面会有几道例题。

01分数规划还会与其他问题结合,如网络流等。

DinkelbachDinkelbach算法

这个算法我是在写这篇文章时才知道的。

思考上述二分算法的思路,设二分过程中某一个二分值为rr,二分时的判断条件是max{f(r)}max{f(r)}的正负性,而这个rr除了让LL右移或者RR左移就没有用了。现在思考某一过程中rr与max{f(r)}max{f(r)}能否再被利用。
二分时,假如max{f(r)}>0max{f(r)}>0这说明最优解在当前rr的右侧,于是让L=rL=r,但是,如果将LL移动到max{f(r)}max{f(r)}对应直线的横截距呢?显然,算法会变得更快。这个思想就是DinkelbachDinkelbach算法的内涵。
DinkelbachDinkelbach实质上是一种迭代算法,基于这样的思想:不去二分答案,而是先随便给定一个答案,然后根据更优的解(max{f(r)}max{f(r)}对应直线的横截距)不断移动答案,逼近最优解。理论上它比二分快些。
在这个算法中,一般将rr初始化为00。

再说最小化

看上面的图,也很好理解,就是最左边的rr为r?r?,当前的rr确定时需要用到min{f(r)}min{f(r)}。
如果min{f(r)}>0min{f(r)}>0,那么r<r?r<r?;
如果min{f(r)}=0min{f(r)}=0,那么r=r?r=r?;
如果min{f(r)}<0min{f(r)}<0,那么r>r?r>r?。

做题时认清哪个是valueivaluei,哪个是costicosti,再记住上面几张图,基本不会出错了。

两种算法的比较

DinkelbachDinkelbach算法的弊端就是需要保存解。这两个算法解决统一问题实际上都有可能快些。
我觉着我一般还是用二分。。。。

例题

[POJ2976]Dropping tests

问题模型2

最优比率生成树

带权无向图GG, 对于图中每条边eiei, 都有valueivaluei和costicosti,现在求一棵生成树TT,最大(小)化∑valuei∑costi,ei∈T∑valuei∑costi,ei∈T

解决方法

套用01分数规划模型,如果ei∈Tei∈T则xi=1xi=1否则xi=0xi=0。

二分法

二分答案rr,边赋值weighti=valuei?r?costiweighti=valuei?r?costi,因为是生成树,边的数量确定,那么max{f(r)}max{f(r)}需要选取前|G|?1|G|?1大的weightiweighti,也就是求最大生成树,按最大生成树权值的正负性就可以二分了。最小化就求最小生成树。

DinkelbachDinkelbach算法

当前答案rr,边赋值weighti=valuei?r?costiweighti=valuei?r?costi,同样求最大生成树,找到max{f(r)}max{f(r)}对应的边集{xi}{xi},也就是最大生成树的边集。对这个边集找横截距当做下一次答案。横截距是啥呢?
f(r)=B?A?rrr=0=B/A=∑valuei?xi∑costi?xi
f(r)=B?A?r=0r=B/Ar=∑valuei?xi∑costi?xi

最小化就求最小生成树。
例题

[POJ2728]Desert King

问题模型3

最优比率环

给定有点权和边权的图,求一个环,使得环的点权和与边权和的比值最大。

解决方法

套用01分数规划模型,点权为valueivaluei,边权为costicosti,一个环为CC
问题要求最大化∑valuei∑costi,(i∈C)∑valuei∑costi,(i∈C)
边数和点数是相同的,但上述式子表述不是很正确,意会即可。
若答案为r?r?,那么任意一个环
∑valuei∑costi∑valueir??∑costi?∑valuei≤r?≤r??∑costi≥0
∑valuei∑costi≤r?∑valuei≤r??∑costir??∑costi?∑valuei≥0
最小化时
∑valuei∑costi∑valuei∑valuei?r??∑costi≥r?≥r??∑costi≥0
∑valuei∑costi≥r?∑valuei≥r??∑costi∑valuei?r??∑costi≥0
二分法

设当前答案rr,
r<r?r<r?,至少存在一个环,r?∑costi?∑valuei<0r?∑costi?∑valuei<0,即存在负权回路(将边权设为r?costi?valueir?costi?valuei,不是提前算出,而是在更新路径的时候从哪个点访问到这条边的就将这条边设为相应点权与边权的对应值);
r≥r?r≥r?,则不存在负环。

求负环可以用Bellman-Ford,但是比较慢,一般用spfa算法求负环
具体判断方法为,一个点不能入队nn次,否则有负环;一条最短路径长度不能到nn,否则有负环。两个判断方法可以同时使用。

最小化时边权设为 ∑valuei?r?∑costi∑valuei?r?∑costi即可,同样也是更新时算出此值。

以上具体实现看例题。

DinkelbachDinkelbach算法

如果用这个算法需要记录下来一个负环,实现还是能实现的,但是没有二分+spfa好写。

例题

[POJ3621]Sightseeing Cows
---------------------
作者:tianxiang971016
来源:CSDN
原文:https://blog.csdn.net/hzoi_ztx/article/details/54898323
版权声明:本文为博主原创文章,转载请附上博文链接!

[转]01分数规划算法 ACM 二分 Dinkelbach 最优比率生成树 最优比率环

标签:二分法   test   最短路径   二分答案   长度   解决   com   初始   最短路   

原文地址:https://www.cnblogs.com/lyp-Bird/p/10315898.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!