01分数规划是指这样一类问题:有n个物品,第i个物品有
ai
的价值,价格是
bi
,然后取出若干个物品使得
∑ai∑bi
取到最值。
不妨设这个最大值为x(最小值同理),则有
∑ai−x∑bi<=0
,然后我们就能发现左边这个式子关于x是单调递减的,这样我们就能通过二分来确定这个最大值x。
回到本题,把式子化一下就可以得到
∑(ai−xbi)<=0
,二分一个x之后把所有的
ai−xbi
排序,取出最大的n-k个看和是不是大于0就可以了。
#include<cmath>
#include<cstdio>
#include<vector>
#include<queue>
#include<cstring>
#include<iomanip>
#include<stdlib.h>
#include<iostream>
#include<algorithm>
#define ll long long
#define inf 1000000
#define mod 1000000007
#define N 2000000
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
int a[N],b[N];
double y[N];
int n,k,i;
double l,mid,r;
bool cmp(const int &x,const int &y) {return x > y;}
bool judge(double x)
{
int i; double sum = 0;
fo(i,1,n) y[i] = a[i] - x * b[i];
sort(y+1,y+n+1,cmp);
fo(i,1,n-k) sum += y[i];
return sum >= 0;
}
int main()
{
scanf("%d%d",&n,&k);
while (!(!n && !k))
{
fo(i,1,n) scanf("%d",&a[i]);
fo(i,1,n) scanf("%d",&b[i]);
l = 0.0; r = inf*1.0;
while (r - l >= 1.0e-6)
{
mid = (l + r) / 2;
if (judge(mid)) l = mid; else r = mid;
}
int p = (int) (l * 100 + 0.5);
printf("%d\n",p);
scanf("%d%d",&n,&k);
}
return 0;
}