【JZOJ4230】淬炼神体【数论，数学】【二分】

题目大意：

给出$n,a[]$和$b[]$，求在满足$sum^{n}_{i=1}x_i=m$时$max{frac{sum^{n}_{i=1}a_i imes x_i}{sum^{n}_{i=1}b_i imes x_i}}$。

---

思路：

$0/1$分数规划裸体。
推荐参考《算法竞赛进阶指南》$P181$。
随便取一个数$L$，判断是否有一组$x[]$满足$sum^{n}_{i=1}(a_i-L imes b_i)geq L$。
若有，那么变形得$(sum^{n}_{i=1}a_i imes x_i)-L imes(sum^{n}_{i=1}b_i imes x_igeq 0)$。则$frac{sum^{n}_{i=1}a_i imes x_i}{sum^{n}_{i=1}b_i imes x_i}geq L$。
所以答案就不小于$L$，否则答案就小于$L$。
那么就二分$L$，知道精度小于$10^{-3}$为止。

---

代码：

#include <cstdio>
#include <algorithm>
using namespace std;

const int N=100010;
int n,m;
double l,r,mid;

struct node
{
	int a,b;
	double s;
}a[N];

bool cmp(node x,node y)
{
	return x.s>y.s;
}

bool check()  //判断最大的m个之和是否大于L
{
	double sum=0.0;
	for (int i=1;i<=m;i++)
		sum+=a[i].s;
	return sum>=0;
}

int main()
{
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++)
		scanf("%d",&a[i].a);
	for (int i=1;i<=n;i++)
		scanf("%d",&a[i].b);
	l=0.0001;
	r=20000.0;
	while (r-l>=0.0001)  //二分L
	{
		mid=(l+r)/2.0;
		for (int i=1;i<=n;i++)
			a[i].s=(double)a[i].a/mid-(double)a[i].b;
		sort(a+1,a+1+n,cmp);
		if (check()) l=mid;
			else r=mid;
	}
	printf("%0.3lf
",l);
	return 0;
}