[二分][DP]表白(love)

题目描述

鸡腿是CZYZ的著名DS，但是不想追妹子的DS不是好GFS，所以鸡腿想通过表白来达到他追到妹子的目的！虽然你对鸡腿很无语，但是故事的设定是你帮助鸡腿找到了妹子，所以现在你必须帮助鸡腿安排表白来实现故事的结局 ！
鸡腿想到了一个很高（sha）明（bi）的做法，那就是去找人来组成表白队伍来增强气势 ！鸡腿有很多好基友来帮忙，鸡腿数了数一共有N个人。但是鸡腿觉得大家排成两队来比较好看，而且鸡腿经过计算，第一队N1个人，第二队N2个人是最佳的队伍。问题来了…有些好基友们虽然很好心但是可能造成不好的影响（形象猥琐），所以鸡腿就给每个人打了分。Q1i表示第i个好基友排到第一队里时的好影响，C1i表示第i个好基友排到第一队里时的不良影响，Q2i表示第i个好基友排到第二队里时的好影响，C2i表示第i个好基友排到第二队里时的不良影响。请给鸡腿一种安排使得Q的和与C的和的比值最大，给出最大值。

Input
第一行给出三个整数N、N1、N2。
第2到N+1行，每行四个整数Q1,C1,Q2,C2。

Output
一行输出一个小数d表示最优化比例是d（保留6位小数）

Sample Input
5 2 2
12 5 8 3
9 4 9 4
7 3 16 6
11 5 7 5
18 10 6 3

Sample Output
2.444444

Data Constraint
对于50%的数据0 < N1 + N2 ≤ N ≤ 50；
对于100%的数据0 < N1 + N2 ≤ N ≤ 500，1 ≤ Q1, Q2 ≤ 2000，1 ≤ C1, C2 ≤ 50。

分析

我们可以二分答案，下界0.01，上界2000
然后对于判断二分的条件，我们先知道这些东西：
对于一个可行答案w，一定有
Σp1[i]+Σp2[i]
—————— ≥ w
Σc1[i]+Σc2[i]
∀i≠j

我们移一波项，得

Σp1[i]-Σc1[i]*w+Σp2[i]-Σc2[i]*w≥0
这个作为判断条件，但是如何选出较为优秀的值使得二分无错误呢？
容易想到最优解，DP
对于p,c比例，我们知道：
如果p1/c1>p2/c2那么这个人给1队是较为好的
然后这个比例我们已经搞成了不是比例的东西：p1[i]-c1[i]*w
那么对Σp1[i]-Σc1[i]*w+Σp2[i]-Σc2[i]*w≥0进行排序（带p1[i]-c1[i]*w和p2[i]-c2[i]*w）
可以有转移方程：
fi,j表示从前i个选择j个出来加入某队
f1[i][j]=max{f1[i-1][j],f1[i-1][j-1]+p1[i]-c1[i]*w
f2[i][j]=max{f2[i-1][j],f2[i-1][j-1]+p2[i]-c2[i]*w
然后由于排序是先优于1队后次于1队，所有f2要倒过来即
n-i+1
然后一波判断即可

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <memory.h>
#define rep(i,a,b) for (i=a;i<=b;i++)
const double eps=1e-9;
using namespace std;
int n,n1,n2;
struct peo {
    int p1,c1,p2,c2;
}p[501];
struct cnt {
    double a,b,c;
}a[501];
double f1[501][501],f2[501][501];

bool Cmp(cnt a,cnt b) {
    return a.c>b.c;
}

bool Judge(double w) {
    int i,j;
    rep(i,1,n) {
        a[i].a=p[i].p1-p[i].c1*w;
        a[i].b=p[i].p2-p[i].c2*w;
        a[i].c=a[i].a-a[i].b;
    }
    sort(a+1,a+n+1,Cmp);
    memset(f1,200,sizeof(f1));
    memset(f2,200,sizeof(f2));
    f1[0][0]=0;
    rep(i,1,n) {
        f1[i][0]=0;
        rep(j,1,min(n1,i))
        f1[i][j]=max(f1[i-1][j],f1[i-1][j-1]+a[i].a);
    }
    f2[0][0]=0;
    for (i=n;i>=1;i--) {
        f2[n-i+1][0]=0;
        rep(j,1,min(n2,n-i+1))
        f2[n-i+1][j]=max(f2[n-i][j],f2[n-i][j-1]+a[i].b);
    }
    rep(i,n1,n-n2)
    if (f1[i][n1]+f2[n-i][n2]>=0) return 1;
    return 0;
}

int main() {
    int i;
    double mid;
    scanf("%d%d%d",&n,&n1,&n2);
    rep(i,1,n)
    scanf("%d%d%d%d",&p[i].p1,&p[i].c1,&p[i].p2,&p[i].c2);
    double l=0.01,r=2000;
    while (l+eps<r) {
        mid=(l+r)/2;
        if (Judge(mid)) l=mid; else r=mid-eps;
    }
    printf("%.6lf",mid);
}