SDUT 2408 Pick apples

Pick apples

Time Limit: 1000ms   Memory limit: 165536K  有疑问？点这里^_^

题目描述

Once ago, there is a mystery yard which only produces three kinds of apples. The number of each kind is infinite. A girl carrying a big bag comes into the yard. She is so surprised because she has never seen so many apples before. Each kind of apple has a size and a price to be sold. Now the little girl wants to gain more profits, but she does not know how. So she asks you for help, and tell she the most profits she can gain.

输入

In the first line there is an integer T (T <= 50), indicates the number of test cases.
In each case, there are four lines. In the first three lines, there are two integers S and P in each line, which indicates the size (1 <= S <= 100) and the price (1 <= P <= 10000) of this kind of apple.

In the fourth line there is an integer V，(1 <= V <= 100,000,000)indicates the volume of the girl's bag.

输出

For each case, first output the case number then follow the most profits she can gain.

示例输入

1
1 1
2 1
3 1
6

示例输出

Case 1: 6

提示

 

来源

2012年"浪潮杯"山东省第三届ACM大学生程序设计竞赛

　　

　　一开始以为是完全背包问题，但是背包的体积V太大，直接背包果断不过。　

　　正解应该是 大范围贪心，小范围背包。

　　直接谈心不能保证充分利用背包的体积，从而不能保证找到最优解。

　　但是背包找到的最优解也肯定是尽可能将性价比高的物品放进去。

　　这样答题思路就出来了，先分出一部分空间 V1 来用来贪心，剩余部分 V2 = V - V1用来背包，贪心即选择性价比最高的物品放进去。

　　如果贪心后 V1 有剩余的空间 ， 设剩余空间为V3，则应将V3回收用来背包，即 V2 += V3。

　　这样得出的结构就是正确答案了。

　　但是这里还牵扯到一个问题，那就是V1和V2的划分问题。

　　这里有两种方法，一种就是直接分出较大的一部分空间来用来背包，显然此部分要在时间和空间上都能承受。

　　还有一种就是找给出的这三个体积的最小公倍数，至于为什么这样能行，我也不知道了，坐等大神来证明 o(∩_∩)o 

　　下面是在SDUT OJ上 AC 的代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>

using namespace std;

long long int gcd(long long int a,long long int b)//辗转相处求最大公约数，再用此数求最小共倍数
{
    long long int c = a%b;
    if(c)
    {
        return gcd(b,c);
    }
    return b;
}

long long int bag[1000010];

struct N
{
    long long int p,w;
    double pc;
} a[5];

bool cmp(N a1,N a2)//按性价比降序排序 为大范围贪心做准备
{
    if(a1.pc > a2.pc)
        return true;
    return false;
}

long long int Max;

int main()
{
    long long int t,js = 1;
    long long int v,i,j;
    cin>>t;
    while(t--)
    {
        for(i = 0; i < 3; i++)
        {
            cin>>a[i].w>>a[i].p;
            a[i].pc = a[i].p*1.0/a[i].w;
        }
        cin>>v;

        sort(a,a+3,cmp);



        long long int g1 = gcd(a[0].w,a[1].w);
        g1 = a[0].w/g1*a[1].w;
        long long int g = gcd(a[2].w,g1);
        g = a[2].w/g*g1;//寻找这个三个数的最小共倍数；

        g += 100;//一开始没有这一步在SDUST OJ上AC了  但是在SDUT上 死活A不了 

        if(g <= v)
        {
            long long int temp_v = v-g;
            Max = (long long int)temp_v / a[0].w * a[0].p;
            temp_v %= a[0].w;
            g += temp_v;
        }
        else
        {
            Max = 0;
            g = v;
        }

        memset(bag,0,(g+2)*sizeof(long long int));

        long long int Max2 = 0;

        for(i = 0;i < 3; i++)//中规中矩的完全背包
            for(j = a[i].w;j <= g; j++)
            {
                bag[j] = max(bag[j],bag[j-a[i].w]+a[i].p);
                if(bag[j] > Max2)
                    Max2 = bag[j];
            }

        cout<<"Case "<<js++<<": "<<Max+Max2<<endl;
    }
    return 0;
}