zcmu 1111 松哥的困惑II（多重背包）

【题目】

1111: 松哥的困惑II

Time Limit: 5 Sec  Memory Limit: 128 MB
Submit: 121  Solved: 51
[Submit][Status][Web Board]

Description

松哥历尽千辛万苦终于找到女朋友了,但是一星期后又回到了单身的日子,松哥很生气后果很严重,所以松哥决定大吃一顿,但是由于松哥很胖,他要吃m千克的东西才会饱,松哥喜欢吃n样东西,每份东西有ai千克,需要bi元钱,只有ci份.松哥希望知道最少需要花多少钱才能使他吃饱,你能告诉他嘛.东西只能一份一份买,不能只买半份.

Input

多组测试数据.

每组测试数据的第一行包含两个正整数n,m(n<=100,m<=10000).

第二行有n个正整数a1,a2,a3,a4…an-1.(ai<=10000)

第三行有n个正整数b1,b2,b3,b4…bn-1.(bi<=100)

第四行有n个正整数c1,c2,c3,c4…cn-1.(ci<=100)

Output

对于每组数据输出能使松哥吃饱所需要花费最少的钱.如果松哥吃完了所有的东西还不能吃饱,请输出“impossible”.

Sample Input

3 1
1 1 2
1 2 1
1 1 1
3 100
1 1 2
1 2 1
1 1 1
2 2
2 1000
2 10
2 2
2 5
2 1000
2 10
2 2

Sample Output

1
Impossible
2
10

【题解】

【多重背包】

给你n种不同的物品，每个物品有自己的重量w[i]，和价值v[i]，如果每个物品最多只能拿c[i]个，给你容量为m的背包，怎样才能取得最大价值？

状态转移方程：dp[j]=MAX{dp[j],dp[j-k*w[i]]+k*v[i]};

基本操作：

for(i=1;i<=n;i++)
      for(j=m;j>=0;j--)
           for(k=0;k<=c[i];k++) 
           {
                if(j-k*w[i]<0)   break;
                dp[j]=max(dp[j],dp[j-k*w[i]]+k*v[i]);
          }

本题是要求对于每组数据输出能使松哥吃饱所需要花费最少的钱，所以状态转移方程为dp[j+k*a[i]]=MIN(dp[j+k*a[i]],dp[j]+k*b[i]);

dp[j]存储的是吃j千克食物的花费最小值，需要注意的是当j+k*a[i]超过m时存储在dp[m]中。

【代码】

#include<bits/stdc++.h>
using namespace std;
int main()
{
    int n,m,i,j,k;
    int a[105],b[105],c[105],dp[10005];
    while(~scanf("%d%d",&n,&m))
    {
        int sum=0,num=0;
        for(i=0;i<n;i++) scanf("%d",&a[i]);
        for(i=0;i<n;i++) scanf("%d",&b[i]);
        for(i=0;i<n;i++) scanf("%d",&c[i]);
        for(i=0;i<n;i++) num+=b[i]*c[i],sum+=a[i]*c[i];
        if(sum<m) printf("Impossible\n");
        else
        {
            for(i=1;i<=m;i++) dp[i]=num;
            dp[0]=0;
            for(i=0;i<n;i++)
                for(j=m;j>=0;j--)
                    if(dp[j]!=num) //可以达到吃j千克食物
                       for(k=0;k<=c[i];k++)
                       {
                           if(j+k*a[i]>=m)
                              dp[m]=min(dp[m],dp[j]+k*b[i]);
                           else
                              dp[j+k*a[i]]=min(dp[j+k*a[i]],dp[j]+k*b[i]);
                       }
            printf("%d\n",dp[m]);
        }
    }
    return 0;
}