bzoj 5281 [Usaco2018 Open]Talent Show 二分+01分数规划

现在好像流行题目传送门。。。

这道题目，一开始还以为就是一个贪心的排一下序的奶牛题(按体重与才艺的比值排序)
当然这是有反例的，比如说：
3 10
1 50
8 1
9 1
这个例子中已经把才艺与体重的比值由小到大排好了。
按照我原来贪心的写法就会取#1#2#3三头奶牛，但是呢，其实取#1#3更优。
不扯那个错得离谱的贪心了，来讲讲正解。
求$∑wi>=W,∑ti/∑wi$最大。
分数规划，二分答案，
$∑ti/∑wi>=ans$
$∑ti−wi∗ans>=0$
$∑1000∗ti−wi∗ans>=0$
那么这道题目就解完了。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 255
#define maxw 1005
using namespace std;
const int INF=1<<30;
int n,allw;
double ans,f[maxw];
struct chj{
    int w,p;
}a[maxn];
inline int read(){
    int ret=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9') ret=ret*10+ch-'0',ch=getchar();
    return ret*f;
}
inline bool check(double x){
    for (int i=1;i<=allw;i++) f[i]=-INF;
    for (int i=1;i<=n;i++){
        double val=(double)a[i].p-x*a[i].w;
        for (int j=allw;j>=0;j--){
            int k=min(j+a[i].w,allw);
            if (f[j]==-INF) continue;
            f[k]=max(f[k],f[j]+val);
        }
    }
    return f[allw]>=1e-5;
}
void fnd(){
    double L=0,R=1e6;
    while (R-L>1e-5){
        double mid=(L+R)/2;
        if (check(mid)) ans=mid,L=mid;
        else R=mid;
    }
}
int main(){
    freopen("5281.in","r",stdin);
    freopen("5281.out","w",stdout);
    n=read();allw=read();
    for (int i=1;i<=n;i++) a[i]=(chj){read(),read()};
    fnd();
    printf("%d\n",(int)floor(ans*1000));
    return 0;
}