Dropping tests
题目分析:
给你一个形如 . 的公式。给你n个ai和bi让你求出最大的值。
思路分析:
一开始是想到了二分最大值,但是判断函数一时没想出来,一开始是认为可以贪心确认。然后,各种的试都不对。最后看了别人说是0-1分数规划问题。才知道了解题思路。
其实,这题涉及的0-1分数规划并不是很难。先介绍一下0-1分数规划吧。
y = sigma(a)*xi / sigma(b)*xi (xi = 0|1)
则存在y的函数t(y) 。
当t(y) = 0 的时候表明sigma(a) == sigma(b)*y
t(y) > 0 的时候表明sigma(a) > sigma(b)*y
t(y) < 0 的时候表明sigma(a) < sigma(b)*y
由以上不等式我们转换可知t(y) = sigma(a) - sigma(b)*y。
结合这道题目,我们为了有解,这希望去除k个最小的sigma(a) - sigma(b)*y这样得到的结果肯定更优。
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
const double INF = 10000000000000000;
const int MAXN = 1000 + 10;
int n,k;
double a[MAXN],b[MAXN],c[MAXN];
bool C(double mid){
double sum = 0;
for(int i = 0;i < n;++i){
c[i] = 100*a[i] - mid*b[i];
sum += c[i];
}
sort(c,c+n);
for(int i = 0;i < k;++i)
sum -= c[i];
return sum >= 0;
}
void solve(){
double lb = 0,ub = INF;
for(int k = 0;k < 100;++k){
double mid = (lb + ub) / 2;
if(C(mid)) lb = mid;
else ub = mid;
}
printf("%.0lf\n",lb);
}
int main()
{
while(scanf("%d%d",&n,&k),(n||k)){
for(int i = 0;i < n;++i)
scanf("%lf",&a[i]);
for(int i = 0;i < n;++i)
scanf("%lf",&b[i]);
solve();
}
return 0;
}