7. 混合背包问题
有 N 种物品和一个容量是 V 的背包。
物品一共有三类:
第一类物品只能用1次(01背包);
第二类物品可以用无限次(完全背包);
第三类物品最多只能用 si 次(多重背包);
每种体积是 vi,价值是 wi。
求解将哪些物品装入背包,可使物品体积总和不超过背包容量,且价值总和最大。
输出最大价值。
输入格式
第一行两个整数,N,V,用空格隔开,分别表示物品种数和背包容积。
接下来有 N 行,每行三个整数 vi,wi,si,用空格隔开,分别表示第 i 种物品的体积、价值和数量。
si=−1 表示第 i 种物品只能用1次;
si=0 表示第 i 种物品可以用无限次;
si>0 表示第 i 种物品可以使用 si 次;
输出格式
输出一个整数,表示最大价值。
数据范围
0<N,V≤1000
0<vi,wi≤1000
−1≤si≤1000
输入样例
4 5
1 2 -1
2 4 1
3 4 0
4 5 2
输出样例:
8
本题思路
分别 01背包,完全背包,多重背包(需要二进制优化,否则会TLE)
#include <bits/stdc++.h>
using namespace std;
const int N=1010;
int n,m;
int f[N];
int main()
{
cin>>n>>m;
for(int i=0;i<n;i++)
{
int v,w,s;
cin>>v>>w>>s;
if(!s)//完全背包
{
for(int j=v;j<=m;j++)
f[j]=max(f[j],f[j-v]+w);
}
else
{
if(s==-1)s=1;//01背包转物品数为1的多重背包
for(int k=1;k<=s;k*=2)//多重背包二进制优化,优化选择次数
{
for(int j=m;j>=k*v;j--)
f[j]=max(f[j],f[j-k*v]+k*w);
s-=k;
}
if(s>0)
{
for(int j=m;j>=s*v;j--)
{
f[j]=max(f[j],f[j-s*v]+s*w);
}
}
}
}
cout<<f[m]<<endl;
return 0;
}
这里给出多重背包二进制优化的相关思路
- 为什么多重背包不能像二重背包一样优化??
答:
完全背包物品个数不限,递推公式会有无限长度,(即没有最后一项)多重背包物品个数会被限制,能够写出最后一项
完全背包优化:
f[i , j ] = max( f[i-1,j] , f[i-1,j-v]+w , f[i-1,j-2v]+2w , f[i-1,j-3v]+3w , .....)
f[i , j-v]= max( f[i-1,j-v] , f[i-1,j-2v] + w , f[i-1,j-2v]+2*w , .....)
由上两式,可得出如下递推关系:
f[i][j]=max(f[i,j-v]+w , f[i-1][j])
- 为什么多重背包可以二进制优化??
答:正常情况下,如果有11个物品,我们需要枚举12次,(0~11)
如果我们将这11个物品分为:{1,2,4,4}(1,2,4,4,可以组成0~11任何一个数字)四组”新的商品“只需要枚举四次
11=0111(B)+(11−0111(B))=0111(B)+0100(B)
如果仍然不是很能理解的话,取这样一个例子:要求在一堆苹果选出n个苹果。我们传统的思维是一个一个地去选,选够n个苹果就停止。这样选择的次数就是n次
二进制优化思维就是:现在给出一堆苹果和10个箱子,选出n个苹果。将这一堆苹果分别按照1,2,4,8,16,…512分到10个箱子里,那么由于任何一个数字x ∈[1,1024]
都可以从这10个箱子里的苹果数量表示出来,但是这样选择的次数就是 ≤10次 。
这样利用二进制优化,时间复杂度就从O(n3)降到O(n2logS),从4*109降到了2*107。
- 如何划分?
//设有s个商品,也就是将s划分
for(int k = 1 ; k <= s ;k*=2)
{
s-=k;
goods.push_back({v*k,w*k});
}
if(s>0)
goods.push_back({v*s,w*s});
这个划分思路很像货币系统,我们的货币可以组成任意个数字。
多重背包二进制优化AC代码:
#include <bits/stdc++.h>
using namespace std;
const int N=2010;
int f[N],n,m;
struct good
{
int w,v;
};
int main()
{
cin>>n>>m;
vector<good>Good;
good tmp;
for(int i=1;i<=n;i++)
{
int v,w,s;
cin>>v>>w>>s;
for(int k=1;k<=s;k*=2)
{
s-=k;
Good.push_back({k*w,k*v});
}
if(s>0)Good.push_back({s*w,s*v});
}
for(auto t:Good)
for(int j=m;j>=t.v;j--)
f[j]=max(f[j],f[j-t.v]+t.w);
cout<<f[m]<<endl;
return 0;
}
本文参考:https://www.acwing.com/solution/content/53633/
https://www.acwing.com/solution/content/20115/