原来这就是01分数规划,以前听左老师讲过做了类似的题,结果并不知道这玩意儿的名字叫01分数规划,涨姿势了
#include <iostream>
#include <cstdio>
#include <algorithm>
#define eps 0.00000001
#define N 1000050
using namespace std;
int n,k,a[N],b[N],e[N];
double v[N];
bool cmp(int p1,int p2) { return v[p1] > v[p2]; }
bool jud(double x) {
for (int i=1;i<=n;i++) e[i] = i;
for (int i=1;i<=n;i++) v[i] = a[i] - x*b[i];
sort(e+1,e+n+1,cmp);
double fz = 0 , fm = 0;
for (int i=1;i<=n-k;i++) fz += a[ e[i] ] , fm += b[ e[i] ];
return fz - x * fm >= 0;
}
int main() {
for (;;) {
scanf("%d%d",&n,&k); if (n == 0 && k == 0) break;
for (int i=1;i<=n;i++) scanf("%d",&a[i]);
for (int i=1;i<=n;i++) scanf("%d",&b[i]);
double l = 0 , r = 1000000000;
while (r-l > eps) {
double mid = (l + r) / (double)2;
if (jud(mid)) l = mid; else r = mid;
}
l = (double)100 * l;
printf("%d\n",(int)(l+0.5));
}
return 0;
}