【POJ 2976】Dropping tests&01分数规划详解

01分数规划是指这样一类问题：有n个物品，第i个物品有$a_i$的价值，价格是$b_i$，然后取出若干个物品使得$\frac{\sum a_i}{\sum b_i}$取到最值。
不妨设这个最大值为x（最小值同理），则有$\sum a_i-x\sum b_i <= 0$，然后我们就能发现左边这个式子关于x是单调递减的，这样我们就能通过二分来确定这个最大值x。
回到本题，把式子化一下就可以得到$\sum (a_i-xb_i)<=0$，二分一个x之后把所有的$a_i-xb_i$排序，取出最大的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;
}