1257 背包问题 V3——分数规划

N个物品的体积为W1，W2......Wn（Wi为整数），与之相对应的价值为P1,P2......Pn（Pi为整数），从中选出K件物品（K <= N)，使得单位体积的价值最大。

Input

第1行：包括2个数N, K(1 <= K <= N <= 50000)
第2 - N + 1行：每行2个数Wi, Pi（1 <= Wi, Pi <= 50000)

Output

输出单位体积的价值（用约分后的分数表示）。

Input示例

3 2

2 2

5 3

2 1

Output示例

3/4

————————————————————————————

第一眼看题目以为是贪心QAQ 后来发现不行

因为如果你现在已有的价值/体积是最佳

而现在有两个价值很小的物品 a b a价值比b大

但是a的体积远大于b的话 此时b肯定是更优的

所以正解应该是二分答案 判断是否合法就好辣

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
using namespace std;
const int M=;
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
int n,k;
int c[M],w[M];
LL sumx,sumy,ansx,ansy;
LL gcd(LL x,LL y){
    while(y){
        LL p=x%y;
        x=y;
        y=p;
    }
    return x;
}
struct node{double v; int pos;}e[M];
bool cmp(node a,node b){return a.v-b.v>1e-;}
bool check(double mid){
    for(int i=;i<=n;i++) e[i].v=1.0*w[i]-1.0*c[i]*mid,e[i].pos=i;
    sort(e+,e++n,cmp);
    double sum=; sumx=; sumy=;
    for(int i=;i<=k;i++){
        sum+=e[i].v;
        sumx+=w[e[i].pos];
        sumy+=c[e[i].pos];
    }
    return sum>=;
}
int main()
{
    n=read(); k=read();
    for(int i=;i<=n;i++) c[i]=read(),w[i]=read();
    double l=,r=;
    while(r-l>1e-){
        double mid=(l+r)/;
        if(check(mid)) l=mid,ansx=sumx,ansy=sumy;
        else r=mid;
    }
    LL d=gcd(ansx,ansy);
    printf("%lld/%lld\n",ansx/d,ansy/d);
    return ;
}