CF993D Compute Power（二分+Dp）

一看到这种求\(min/max\left \{ \frac{\sum a_i}{\sum b_i} \right \}\)的题肯定是\(01\)分数规划，大概可以算作一种二分？

设\(ans\)为当前枚举到的答案，对于本题，如果想把答案往更小的方向更新，必须满足：

\[\sum a_i-ans*\sum b_i \leq 0\]

即：

\[\sum a_i-ans* b_i \leq 0\]

二分从来都是套板子的事，难点是如何写一个\(check\)函数

发现枚举第一个做的任务不现实（有后面跟着第二个任务和不跟第二个任务两种情况，不好讨论），考虑枚举第二个做的任务的最大总价值，用 \(总价值-第二个任务的价值\) 即可,

这个可以用\(Dp\)求：设\(f[i][j]\)表示枚举到在数列\(a\)中排名为\(i\)的任务，一共做了\(j\)个任务时的最大价值。转移方程为：

\[f[i][j+k]=max\left \{ f[i+1][j]+W[i][k](0\leq k \leq cnt[i]) \right \}\]

\(W[i][j]\)为把排名为\(i\)的任务中，价值\((a-ans*b)\) 最大的\(j\)个任务的价值和

注意的事：

1、向上取整

2 、用cout一定要输出long long！！

代码：

#include <bits/stdc++.h>
#define rd(x) scanf("%d",&x)
#define N 55
#define getrank(x) lower_bound(p+1,p+sz+1,x)-p
#define rnk(i) x[i].rnk
#define ll long long
using namespace std;

const double eps=1e-7;
struct qwq{
    int a,b,rnk;
    double v;
}x[N];
int n,cnt[N],pre[N],num[N],sz;
double f[N][N],p[N],W[N][N];

bool cmp(const qwq &x,const qwq &y){ return (x.a==y.a)?x.v>y.v:x.a>y.a; }

double check(double tmp){
    double tot=0,ans=0;int i,j,k;
    for(i=1;i<=n;++i){
        x[i].v=(double)x[i].a-tmp*x[i].b;
        tot+=x[i].v;
    }
    sort(x+1,x+n+1,cmp);memset(f,0xc0,sizeof(f));memset(num,0,sizeof(num));
    for(i=1;i<sz;++i) W[i][0]=0;
    for(i=1;i<=n;++i) W[rnk(i)][++num[rnk(i)]]=W[rnk(i)][num[rnk(i)]-1]+x[i].v;
    ans=f[0][0];
    if(sz==1) return tot-W[sz][num[1]];
    for(int j=0;j<=min(num[sz],num[sz-1]);++j) f[sz-1][j]=W[sz-1][j];
    for(i=sz-2;i>=1;--i){
        for(j=0;j<=pre[sz-1]-pre[i];++j)
            for(k=0;k<=num[i];++k)
                if(n-pre[i]-2*j>=k) f[i][j+k]=max(f[i+1][j]+W[i][k],f[i][j+k]);
    }
    for(i=0;i<=pre[sz-1];++i) ans=max(ans,f[1][i]);
    return tot-ans;
}

int main(){
    int i;
    rd(n);
    for(i=1;i<=n;++i) rd(x[i].a),p[i]=x[i].a;
    for(i=1;i<=n;++i) rd(x[i].b);
    sort(p+1,p+n+1);sz=unique(p+1,p+n+1)-p-1;
    for(i=1;i<=n;++i) cnt[x[i].rnk=getrank(x[i].a)]++;
    for(i=1;i<=sz;++i) pre[i]=pre[i-1]+cnt[i];
    if(sz==1){
        ll aa=0,bb=0;
        for(i=1;i<=n;++i) aa+=1ll*x[i].a,bb+=1ll*x[i].b;
        double qwq=(double)aa/bb*1000.0;
        printf("%I64d",(ll)ceil(qwq));return 0;
    }
    double l=0,r=1e8;
    while(r-l>eps){
        double mid=(r-l)*0.5+l;
        if(check(mid)<=0) r=mid;
        else l=mid;
    }
    l*=1000.0;
    printf("%I64d",(ll)ceil(l));
}