基准时间限制:1 秒 空间限制:131072 KB 分值: 40
难度:4级算法题
有N种物品,每种物品的数量为C1,C2......Cn。从中任选若干件放在容量为W的背包里,每种物品的体积为W1,W2......Wn(Wi为整数),与之相对应的价值为P1,P2......Pn(Pi为整数)。求背包能够容纳的最大价值。
Input
第1行,2个整数,N和W中间用空格隔开。N为物品的种类,W为背包的容量。(1 <= N <= 100,1 <= W <= 50000) 第2 - N + 1行,每行3个整数,Wi,Pi和Ci分别是物品体积、价值和数量。(1 <= Wi, Pi <= 10000, 1 <= Ci <= 200)
Output
输出可以容纳的最大价值。
Input示例
3 6 2 2 5 3 3 8 1 4 1
Output示例
9
李陶冶
(题目提供者)
典型的多重背包问题,传统做法是将其转化成01背包,但是时间复杂度还是没有降低
因此提供一种优化方式:二进制优化。
思路:
我们考虑把第
i
种物品换成若干件物品,使得原问题中第
i
种物品可取的每种策略
——
取
0..n[i]
件
——
均能等价于取若干件代换以后的物品。另外,取超过
n[i]
件的策略必不能出现。
方法是:将第i种物品分成若干件物品,其中每件物品有一个系数,这件物品的费用和价值均是原来的费用和价值乘以这个系数。使这些系数分别为 1,2,4,...,2^(k-1),n[i]-2^k+1,且k是满足n[i]-2^k+1>0的最大整数。例如,如果n[i]为13,就将这种物品分成系数分别为1,2,4,6的四件物品。
分成的这几件物品的系数和为n[i],表明不可能取多于n[i]件的第i种物品。另外这种方法也能保证对于0..n[i]间的每一个整数,均可以用若干个系数的和表示,这个证明可以分0..2^k-1和2^k..n[i]两段来分别讨论得出,并不难,希望你自己思考尝试一下。
这样就将第i种物品分成了O(log n[i])种物品,将原问题转化为了复杂度为O(V*∑log n[i])的01背包问题,是很大的改进。
分成的这几件物品的系数和为n[i],表明不可能取多于n[i]件的第i种物品。另外这种方法也能保证对于0..n[i]间的每一个整数,均可以用若干个系数的和表示,这个证明可以分0..2^k-1和2^k..n[i]两段来分别讨论得出,并不难,希望你自己思考尝试一下。
这样就将第i种物品分成了O(log n[i])种物品,将原问题转化为了复杂度为O(V*∑log n[i])的01背包问题,是很大的改进。
大牛链接:http://www.cnblogs.com/jbelial/articles/2116074.html
具体代码如下:
#include<stdio.h>
#include<algorithm>
using namespace std;
typedef long long ll;
ll dp[50005];
int n,m;
struct node
{
int u,v,w;
}a[200];
void work2(int x,int y)//01背包
{
int i;
for(i=m;i>=x;i--)
dp[i]=max(dp[i],dp[i-x]+y);
}
void work3(int x,int y)//完全背包
{
int i;
for(i=x;i<=m;i++)
dp[i]=max(dp[i],dp[i-x]+y);
}
void work1(int x,int y,int z)//多重背包
{
if(x*z>=m)
work3(x,y);
else
{
int t=1;
while(t<z)
{
work2(t*x,t*y);
z-=t;
t*=2;
}
if(z>0)
work2(z*x,z*y);
}
}
int main()
{
int i;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)
scanf("%d%d%d",&a[i].u,&a[i].v,&a[i].w);
for(i=1;i<=n;i++)
work1(a[i].u,a[i].v,a[i].w);
printf("%lld\n",dp[m]);
}