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

Posted 友人A

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了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=50007;
int read(){
    int ans=0,f=1,c=getchar();
    while(c<0||c>9){if(c==-) f=-1; c=getchar();}
    while(c>=0&&c<=9){ans=ans*10+(c-0); 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-6;}
bool check(double mid){
    for(int i=1;i<=n;i++) e[i].v=1.0*w[i]-1.0*c[i]*mid,e[i].pos=i; 
    sort(e+1,e+1+n,cmp);
    double sum=0; sumx=0; sumy=0;
    for(int i=1;i<=k;i++){
        sum+=e[i].v;
        sumx+=w[e[i].pos];
        sumy+=c[e[i].pos];
    }
    return sum>=0;
}
int main()
{
    n=read(); k=read();
    for(int i=1;i<=n;i++) c[i]=read(),w[i]=read();
    double l=0,r=50000;
    while(r-l>1e-6){
        double mid=(l+r)/2;
        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 0;
}
View Code

 

以上是关于1257 背包问题 V3——分数规划的主要内容,如果未能解决你的问题,请参考以下文章

51nod 1257 01分数规划/二分

[51nod] 1257 背包问题 V3

51nod 1257 背包问题 V3(这不是背包问题是二分)

用动态规划求解分数背包问题

P4322 [JSOI2016]最佳团体(分数规划&树上背包)

01分数规划