zl程序教程

您现在的位置是:首页 >  其他

当前栏目

【POJ2976】Dropping Tests(分数规划)

规划 分数 tests
2023-09-11 14:14:41 时间

【POJ2976】Dropping Tests(分数规划)

题面

Vjudge
翻译在\(Vjudge\)上有(而且很皮)

题解

简单的\(01\)分数规划
需要我们做的是最大化\(\frac{\sum a[i]}{\sum b[i]}\)
考虑二分答案
将最大化问题转换为判定问题
\(\sum{a[i]}-mid\sum{b[i]}\geq 0\)
因为所有选定的\(i\)是一样的
所以可以将权值化为\(a[i]-mid·b[i]\),这样只需要贪心的选择最大的那部分检查是否大于零就行了

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 1111
inline int read()
{
    RG int x=0,t=1;RG char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}
int a[MAX],b[MAX],n,K;
double c[MAX];
int main()
{
	while(233)
	{
		n=read();K=read();
		if(!n&&!K)break;
		for(int i=1;i<=n;++i)a[i]=read();
		for(int i=1;i<=n;++i)b[i]=read();
		double l=0,r=100;
		while(r-l>1e-5)
		{
			double mid=(l+r)/2;
			for(int i=1;i<=n;++i)c[i]=a[i]-mid*b[i];
			sort(&c[1],&c[n+1]);
			double t=0;
			for(int i=n;i>K;--i)t+=c[i];
			if(t>=0)l=mid;
			else r=mid;
		}
		printf("%.0f\n",l*100);
	}
	return 0;
}