poj2976（0-1 分数规划）

poj2976（0-1 分数规划）

题意

给定n个二元组$(a_i,b_i)$除去n个，剩下的n-k组成集合S作$\frac{{\sum }_{i\in S}{a}_i}{{\sum }_{i\in S}{b}_i}$使其达到最大，原题还让将其取到最近整数！

0-1分数规划

这道题是一道很典型的0-1分数规划，网上的资源看了许多，但是好像没有人介绍在求解答案的时候二分的依据，也没有人解释为何要求解max(F(L))，在这里，我将给出解释。
首先还是对整个算法大概做个描述，具体细节一会解释
1.记$F(L)={\sum }_{i=1}^{i=n}(a_i-L{b}_i)x_i$其中${\sum }_{i=1}^{i=n}{x}_i=n-k$
2.由上面的表达式可以知道对于固定的L使得F(L)最大时即$max(F(L))$为将$a_i-L{b}_i$进行排序取前n-k大的元素，使其$x_i=1$其余为0，注意这里是在求F(L)的最值，而不是$\frac{{\sum }_{i\in S}{a}_i}{{\sum }_{i\in S}{b}_i}$的最值，
3.记$max（\frac{{\sum }_{i\in S}{a}_i}{{\sum }_{i\in S}{b}_i}）=ans$对L进行二分，
当maxF(L)>0的时候则有L<max,
当maxF(L)<0时，L>max,
于是根据上两条进行二分，直到找到使得max(F(L))的L，此时的L==ans

在这里我们着重解释第3步中的二分的想法是怎么来的，其余未明白部分见大佬文章

好啦下面就开始了

对第3条的分析

1.首先对，“当maxF(L)>0的时候则有L<max”,进行解释

$F(L)={\sum }_{i=1}^{i=n}(a_i{x}_i)-L{\sum }_{i=1}^{i=n}(b_i{x}_i)=C$
我们先忘掉上述的结论，重新给出一个结论：
对于某一L，若存在序列$x_i$满足其和未n-k，且使得F(L)>0的时候则有L<max，证明如下：
$F(L)={\sum }_{i=1}^{i=n}(a_i{x}_i)-L{\sum }_{i=1}^{i=n}(b_i{x}_i)=C>0$
${\sum }_{i=1}^{i=n}(a_i{x}_i)=C+L{\sum }_{i=1}^{i=n}(b_i{x}_i)$
$\frac{{\sum }_{i=1}^{i=n}(a_i{x}_i)}{{\sum }_{i=1}^{i=n}(b_i{x}_i)}=\frac{C}{{\sum }_{i=1}^{i=n}(b_i{x}_i)}+L=D$
根据上式，可以发现，存在xi使得表达式$\frac{{\sum }_{i=1}^{i=n}(a_i{x}_i)}{{\sum }_{i=1}^{i=n}(b_i{x}_i)}$的值D>L，于是L比目标表达式最值max要小，因为max(F(L))>F(L),若存在xi使得F(L)>0,则根据不等式传递性必有max(F(L))>F(L)>0， 并且使用max(F(L))在判断下面证明也能用到，节约了代码量！

2.接下来对"当maxF(L)<0的时候则有L>max"解释

F(L)<maxF(L)<0,则可以知道对所有$x_i$有$F(L)={\sum }_{i=1}^{i=n}(a_i{x}_i)-L{\sum }_{i=1}^{i=n}(b_i{x}_i)=C<0$,
对等式按照上面的方式进行变形
$\frac{{\sum }_{i=1}^{i=n}(a_i{x}_i)}{{\sum }_{i=1}^{i=n}(b_i{x}_i)}=\frac{C}{{\sum }_{i=1}^{i=n}(b_i{x}_i)}+L=D<L$对任意$x_i$成立，对于使得$\frac{{\sum }_{i=1}^{i=n}(a_i{x}_i)}{{\sum }_{i=1}^{i=n}(b_i{x}_i)}$取的最大值max得$x_i$也有max<L

有了1,2得判据,就可以根据二分法进而得出答案了.

代码如下:

#include<iostream>
#include<cmath>
#include<algorithm>
#define maxsize 1004
using namespace std;
int n,k;
double a[maxsize],b[maxsize];
bool justify(double r)
{
	double d[maxsize];
	for(int i=0;i<n;i++){
		d[i]=a[i]-r*b[i];
		// cout<<d[i]<<endl;
	}
	sort(d,d+n);
	double sum=0;
	for(int i=k;i<n;i++)sum+=d[i];
	// cout<<sum<<endl;
	return sum>0;
}
int main()
{
	// freopen("ini","r",stdin); 
	while(~scanf("%d%d",&n,&k)&&(n!=0||k!=0)){
		for(int i=0;i<n;i++)scanf("%lf",&a[i]);
		for(int i=0;i<n;i++)scanf("%lf",&b[i]);
		double l=0,r=1;
		while(fabs(r-l)>1e-6)
		{
			double mid=(l+r)/2;
			if(justify(mid))l=mid;
			else r=mid;
		}
		int ans=(int)((r+0.005)*100);
		cout<<ans<<endl;
	}
	return 0;
}