01分数规划算法 信息学竞赛 OI ACM 二分 Dinkelbach 最优比率生成树 最优比率环
01分数规划
张天翔
blog.csdn.net/hzoi_ztxztx97@QQ.com
前置技能
二分思想最短路算法一些数学脑细胞？
问题模型1
基本01分数规划问题
给定n个二元组(valuei,costi)，valuei是选择此二元组获得的价值（非负），costi是选择此二元组付出的代价（非负），设xi(xi∈{0,1})代表第i个二元组的选与不选，最大(小)化下式 maximize(or minimize) r=∑valuei⋅xi∑costi⋅xi
解决方法
二分法
设r最大值为r∗， r∗=∑valuei⋅xi∑costi⋅xi
∑valuei⋅xi−r∗⋅∑costi⋅xi=0
设一个函数，自变量为r值， f(r)=∑valuei⋅xi−r∗⋅∑costi⋅xi
观察这个函数，假如{xi}固定，则这个函数就是坐标系中一条直线(y=B−A⋅x)，每一组{xi}对应着一条直线，这些直线斜率非正(因为−A=−∑costi⋅xi≤0)，纵截距非负(因为B=∑valuei⋅xi≥0 )，如图1。 
对于每一条直线，当f(r)=0时，横截距就是这一组的r，那么r∗就是每条直线横截距的最大值（每组{xi}对应r的最大值）如图2。 
在图中上任取一条垂直x轴的竖线， 如果存在直线与这条竖线的交点纵坐标为正，那么最优值一定在当前竖线的右侧； 如果所有直线与这条竖线交点纵坐标为负，那么最优值一定在当前竖线的左侧； 如果所有直线与这条竖线交点纵坐标非正且存在直线与这条竖线交点纵坐标为0，那么当前竖线横坐标即为最优值r∗。 
按照这个思想，可以二分答案r，那么二分时如何进行判断呢？
选择一个r时需要判断所有f(r)的最大值是否为0，如果max{f(r)}>0则r<r∗；如果max{f(r)}<0 则 r>r∗。 怎样求max{f(r)}? f(r)=∑valuei⋅xi−r⋅∑costi⋅xi=∑(valuei−r⋅costi)⋅xi
二分一个r时，每个二元组的valuei−r⋅costi 都可以求出，设其为weighti，现在的目标就是找到一组{xi}使得∑wighti⋅xi最大（即求max{f(r)}）。怎么找到这一组{xi}，或者直接求得max{f(r)}呢？具体问题具体分析，经常借助最短路算法判断是否存在负环。下面会有几道例题。
01分数规划还会与其他问题结合，如网络流等。
Dinkelbach算法
这个算法我是在写这篇文章时才知道的。
思考上述二分算法的思路，设二分过程中某一个二分值为r，二分时的判断条件是max{f(r)}的正负性，而这个r除了让L右移或者R左移就没有用了。现在思考某一过程中r与max{f(r)}能否再被利用。 二分时，假如max{f(r)}>0这说明最优解在当前r的右侧，于是让L=r，但是，如果将L移动到max{f(r)}对应直线的横截距呢？显然，算法会变得更快。这个思想就是Dinkelbach算法的内涵。 Dinkelbach实质上是一种迭代算法，基于这样的思想：不去二分答案，而是先随便给定一个答案，然后根据更优的解（max{f(r)}对应直线的横截距）不断移动答案，逼近最优解。理论上它比二分快些。 在这个算法中，一般将r初始化为0。
两种算法的比较
Dinkelbach算法的弊端就是需要保存解。这两个算法解决统一问题实际上都有可能快些。 我觉着我一般还是用二分。。。。
例题
只先写了文章，还差相关题目的分析，代码等。得等个一两天再。2017-02-06留。
Poj2976Dropping tests
问题模型2
最优比率生成树
带权无向图G， 对于图中每条边ei， 都有valuei和costi，现在求一棵生成树T，最大（小）化∑valuei∑costi,ei∈T
解决方法
套用01分数规划模型，如果ei∈T则xi=1否则xi=0。
二分法
二分答案r，边赋值weighti=valuei−r⋅costi，因为是生成树，边的数量确定，那么max{f(r)}需要选取前|G|−1大的weighti，也就是求最大生成树，按最大生成树权值的正负性就可以二分了。
Dinkelbach算法
当前答案r，边赋值weighti=valuei−r⋅costi，同样求最大生成树，找到max{f(r)}对应的边集{xi}，也就是最大生成树的边集。对这个边集找横截距当做下一次答案。横截距是啥呢？ f(r)=B−A⋅rrr=0=B/A=∑valuei⋅xi∑costi⋅xi
例题
poj2728Desert King(最优比率生成树)
问题模型3
最优比率环
给定有点权和边权的图，求一个环，使得环的点权和与边权和的比值最大。
解决方法
套用01分数规划模型，点权为valuei，边权为costi，一个环为C 问题要求最大化∑valuei∑costi,(i∈C) 边数和点数是相同的，但上述式子表述不是很正确，意会即可。 若答案为r∗，那么任意一个环 ∑valuei∑costi∑valueir∗⋅∑costi−∑valuei≤r∗≤r∗⋅∑costi≥0
二分法
设当前答案r， r<r∗，至少存在一个环，r⋅∑costi−∑valuei<0，即存在负权回路（将边权设为r⋅costi−valuei，不是提前算出，而是在更新路径的时候从哪个点访问到这条边的就将这条边设为相应点权与边权的对应值）； r≥r∗，则不存在负环。
求负环可以用Bellman-Ford，但是比较慢，一般用spfa算法求负环 具体判断方法为，一个点不能入队n次，否则有负环；一条最短路径长度不能到n，否则有负环。两个判断方法可以同时使用。
Dinkelbach算法
如果用这个算法需要记录下来一个负环，实现还是能实现的，但是没有二分+spfa好写。
例题
poj 3621(最优比率环)
问题模型4
最大密度子图
这个问题会写在网络流总结中。
参考资料
【1】KirisameMarisa - NYIST 914Yougth的最大化【二分搜索/Dinkelbach算法】 【2】PerSeAwe - ［Algorithm］01分数规划
