POJ Dropping tests

                                                                  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;
}