Jzoj P3454 表白___二分+思维+dp

题目大意：

有$2$个队伍分别选人，一共有$N$个人，$1$队选$N_1$个人，$2$队选$N_2$个人，每一个成员给出在$1$队的好影响$Q_1[i]$以及不良影响$C_1[i]$，在$2$队的好影响$Q_2[i]$以及不良影响$C_2[i]$，每个人最多只能进入一个队伍，问最后是2个队的$Q$总和与$C$总和的比例最大，即让 $\frac{Sum_Q}{Sum_C}$最大，输出最大比例保留$6$位小数。
$0 < N1 + N2 ≤ N ≤ 500，1 ≤ Q1, Q2 ≤ 2000，1 ≤ C1, C2 ≤ 50$

分析：

我们可以发现这题可以进行对答案的二分，
即二分区间$[10^{-6}..2000]$，
但是因为避免精度的误差，
我选择的是二分$[10^{-8}..2000]$
然后我们当前二分到的答案，设为$x$，
则，我们要判断$x$是否可行，
显然当$x$满足，
$\frac{ΣQ_1[i]+ΣQ_2[j]}{ΣC_1[i]+ΣC_2[j]}≥x$
那么，我们化简一下，
约掉分母，
$ΣQ_1[i]+ΣQ_2[j]≥(ΣC_1[i]+ΣC_2[j])x$
然后去括号，
$ΣQ_1[i]+ΣQ_2[j]≥ΣC_1[i]*x+ΣC_2[j]*x$
移项得到：
$(ΣQ_1[i]-ΣC_1[i]*x)+(ΣQ_2[j]-ΣC_2[j]*x)≥0$

那么我们知道，对于第$i$个人而言，如果要选择，
那么选择$Q_1[i]-C_1[i]*x$， $Q_2[i]-Q_2[i]*x$较大的那一个，是更优的
那么每次我们设
$a[i]$表示$Q_1[i]-C_1[i]*x$
$b[i]$表示$Q_2[i]-C_2[i]*x$
$c[i]$表示$a[i]-b[i]$
那么显然我们可以发现，当$c[i]$排序后，为一个递减序列的时候（即从大到小排列），
那么此时前面的选$1$队比选$2$队的更优，后面的选$2$队比选$1$队更优
此时我们设
$f[i][j]$表示前$i$个人选了$j$个人可以得到的最大$(ΣQ_1[]-ΣC_1[]*x)$
$g[i][j]$表示后面的人$[i..n]$选了$j$个人可以得到的最大$(ΣQ_2[]-ΣC_2[]*x)$
那么当存在任意$f[i][n1]+g[i+1][n2]$$≥0$，
则此时的 x x <script type="math/tex" id="MathJax-Element-7181">x</script>是有效的，继续向后二分，否则向前

代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define N 505

using namespace std;

struct Node {
    double x1, x2, cd;
}a[N];
double f[N][N], g[N][N], Q1[N], C1[N], Q2[N], C2[N];
int n, n1, n2;

bool cmp(Node aa, Node bb) {
    return aa.cd > bb.cd;
}

bool Check(double x) {
    for (int i = 1; i <= n; i++) {
         a[i].x1 = Q1[i] - C1[i] * x;
         a[i].x2 = Q2[i] - C2[i] * x;
         a[i].cd = a[i].x1 - a[i].x2;
    }   
    sort(a + 1, a + n + 1, cmp);
    for (int i =  0; i <= n+1; i++) {
         f[i][0] = g[i][0] = 0;
         for (int j = 1; j <= n1; j++) f[i][j] = -1e9;
         for (int j = 1; j <= n2; j++) g[i][j] = -1e9;
    }
    for (int i = 1; i <= n; i++) 
         for (int j = 1; j <= min(i, n1); j++) f[i][j] = max(f[i-1][j], f[i-1][j-1] + (a[i].x1));

    for (int i = n; i >= 1; i--) 
         for (int j = 1; j <= min(n-i+1, n2); j++) g[i][j] = max(g[i+1][j], g[i+1][j-1] + (a[i].x2));

    for (int i = n1; i <= n-n2; i++)
         if (f[i][n1] + g[i+1][n2] >= 0) return 1;
    return 0;
}

int main() {
    freopen("love.in","r",stdin);
    freopen("love.out","w",stdout);
    scanf("%d %d %d", &n, &n1, &n2);
    for (int i = 1; i <= n; i++) 
         scanf("%lf %lf %lf %lf", &Q1[i], &C1[i], &Q2[i], &C2[i]);
    double l = 1e-8, r = 2000, ans = 0; 
    while (l <= r) {
           double mid = (l + r) / 2;
           if (Check(mid)) ans = mid, l = mid + 1e-8;
                      else r = mid - 1e-8;
    }
    printf("%.6lf\n", ans);
    return 0;
}