01分数规划

所谓的01分数规划问题就是说，给定两个等长数组a[i]，b[i]，求取第三个数组X[i]，其中X[i]中只含0和1两种元素，也就是一串和a[i],b[i]等长的0 1 序列，使得R=sigma(a[i]*x[i])/sigma(b[i]*x[i])，可以看出，X[i]的作用就是决定选或不选a,b数组中的第i号元素，使得R取得最值。

01分数规划问题主要包含一般的01分数规划、最优比率生成树问题、最优比率环问题等。这里我选择一般的01分数规划和大家谈谈自己所学。

我们要求得这样的数组X[i]使R最大，直接试遍各种0 1 组合算出R再求个最值的方法是我们的第一反应，但细细想来这方法可实现性不强，随便给个不长的数组，就要判断很多种不同的排列组合，时间/空间复杂度都是极其的高。所以我们要换一种思路进行下去：

首先我们对所要求的R式进行调整：

        R=sigma(a[i]*x[i])/sigma(b[i]*x[i])

得到

        sigma(a[i]*x[i])-R*sigma(b[i]*x[i]) = 0

这不难理解吧，但是却变复杂的除运算为减法了，我们只需求得表达式为零的R值，这样的求方程根的问题用二分法不难解决

在这就有很多人不解了，这样的R值岂不是有很多，怎样得到最优的X[i]呢？
所以说我们在求R值之前得先把最优的X[i]确定下来

大家知道，对于一组X[i]来说，上式 F(R) = sigma(a[i]*x[i])-R*sigma(b[i]*x[i]) 是随R值的增大而递减的函数，R大到一定程度，F（R）始终为负。所以我们二分求R的时候，利用R值算出F(R)的最大值（此时必然对应着一组X[i]），若最大值大于零，则R可以再大一点，即可以在二分的后半部分去找。否则在二分的前半部分去找。

那么问题来了，怎么由一个R值算出此时F(R)的最大值呢？我们这里再化简一下F(R)就可以得到 F(R) = sigma((a[i]-R*b[i])*x[i])，若再建一个数组的d[i] = a[i]-R*b[i]，那么F(R) = sigma(d[i]*x[i]) . 看，这不是求一个数组d哪几项的和最大的问题嘛，遍历一遍d[i]数组就行了。

大致体现思想的伪代码：

double  low = -1e6,high = 1e6;   //R的最值开始查找区间，
for(int i = 0;i<200;i++){       //循环二分200次大致得到比较精准的R
  double mid = ( low + high)/2; //二分（mid即使当前R值）
 double F_R = findMax(mid);    //找到该R值对应的函数F(R)的最大值
 if(F_R>0)  low = mid;
     else   high = mid;
}
double R = low;  //这里取low和high都行，二分200次low与high都很接近了

关于最大值那就要看具体题目的具体要要求了0_=,有的只是单纯求个最大，有的要求连续和的最大值，有的要求连续至少m个数的最大值……

下一次blog我将分享一道与之相关的题，展现0 1分数规划的具体代码实现
大家对于我说的不足之处欢迎补充，有哪个地方叙述不能理解也可以评论区给我留言0_0 谢谢大家！