POJ2976 Dropping tests (分数规划)

传送门

今天在看Amber神犇的《最小割模型在信息学竞赛中的应用》的时候看到了分数规划这玩意，然后就小水了一发（虽然这道题和最小割并没有什么关系）。

Dropping tests 是一道裸的01分数规划，大意就是给出n个分数，在里面选出n-k个要求这些分数 分子加分子 : 分母加分母最大。

关于分数规划的介绍，我觉得[Algorithm]01分数规划——Update：2012年7月27日这篇文章写得相当不错

上一下我的代码（采用二分做的，迭代不是很熟，以后再更新）

#include<cstdio>#include<algorithm>#define MAXN 1005#define eps 1e-7#define max(a, b) (a)>(b) ? (a):(b)int a[MAXN], b[MAXN];double l, r, ans, d[MAXN];int main(){    int n, k;    while(~scanf("%d%d", &n, &k) && n+k)    {        l = 0;r = 1;        for(int i = 1; i <= n; i ++)            scanf("%d", a + i);        for(int i = 1; i <= n; i ++)            scanf("%d", b + i);        double mid;        while(l + eps < r)        {            ans = 0;            mid = (l + r) / 2;            for(int i = 1; i <= n; i ++)                d[i] = a[i] - mid * b[i];            std::sort(d+1, d+n+1);            for(int i = k + 1; i <= n; i ++)                ans += d[i];            if(ans > 0) l = mid;            else r = mid;        }        printf("%.0lf\n", mid*100);    }    return 0;}