51nod 1257 01分数规划/二分

Posted zzq

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了51nod 1257 01分数规划/二分相关的知识,希望对你有一定的参考价值。

http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1257

1257 背包问题 V3技术分享

基准时间限制:3 秒 空间限制:131072 KB 分值: 80 难度:5级算法题
技术分享 收藏
技术分享 关注
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
第一次写分数规划,感觉就是数学真神奇- -
假设我们已知选了k件,他们的单位价值就是 x = ∑(pi) / ∑(wi) , 分解一下得到 ∑(pi) - x*∑(wi) = 0 ,对于所有合法的x这个式子肯定成立,所以对于max {x}也是,
我们有 ∑(pi) - x1*∑(wi) = 0 , 如果 x1==max(x) ,显然x1就是最大值 ; 如果x1<max(x)则他不是最优解 ; 如果x1>max(x) 这个式子<0说明不合法。
只要二分一下x的值就好了,为了使得x尽量合法显然我们将所有的 pi-x*wi排序之后优先选择较大的即可。
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 using namespace std;
 6 #define LL long long 
 7 int W[50005], P[50005];
 8 struct node { double d;int u; }D[50005];
 9 bool cmp(node A, node B) { return A.d < B.d; }
10 int gcd(int a, int b) { return b == 0 ? a : gcd(b, a%b); }
11 bool ok(double x, int N, int K, int &a, int &b)
12 {
13     for (int i = 1;i <= N;++i) {
14         D[i].d = 1.0*P[i] - x*W[i];
15         D[i].u = i;
16     }
17     sort(D + 1, D + 1 + N,cmp);
18     double s = 0; int w1 = 0, p1 = 0;
19     for (int i = N;i > N - K;i--) {
20         s += D[i].d;
21         w1 += W[D[i].u];
22         p1 += P[D[i].u];
23     }
24     if (s >= 0) {
25         int y = gcd(w1, p1);
26         a = w1 / y;
27         b = p1 / y;
28     }
29     return s >= 0;
30 }
31 int main()
32 {
33     int N, K, i, j, k;
34     cin >> N >> K;
35     for (i = 1;i <= N;++i)
36     {
37         scanf("%d%d", W + i, P + i);
38     }
39     double l = 0, r = 50000;
40     int w,p;
41     while (r - l >= 1e-8) {
42         double mid = r - (r - l) / 2;
43         if (ok(mid,N,K,w,p)) {
44             l = mid;
45         }
46         else {
47             r = mid;
48         }
49     }//cout << l << endl;
50     printf("%d/%d\n",p,w);
51     return 0;
52 }

 









以上是关于51nod 1257 01分数规划/二分的主要内容,如果未能解决你的问题,请参考以下文章

01分数规划

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

[51nod] 1257 背包问题 V3

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

[POJ2976][POJ2728]01分数规划问题的二分答案解法

[转]01分数规划算法 ACM 二分 Dinkelbach 最优比率生成树 最优比率环