多重背包问题
有 N 种物品和一个容量是 V 的背包。
第 i 种物品最多有 si 件,每件体积是 vi,价值是 wi。
求解将哪些物品装入背包,可使物品体积总和不超过背包容量,且价值总和最大。
输出最大价值。
输入格式
第一行两个整数,N,V,用空格隔开,分别表示物品种数和背包容积。
接下来有 N 行,每行三个整数 vi,wi,si,用空格隔开,分别表示第 i 种物品的体积、价值和数量。
输出格式
输出一个整数,表示最大价值。
数据范围
0<N≤1000
0<V≤2000
0<vi,wi,si≤2000
与完全背包问题的区别无非多了一个每类物品的最大上限,从这一点入手可以很简单的得到
朴素解法如下:
1.#include<bits/stdc++.h>
using namespace std;
const int N=1005;
int f[N][N],w[N],v[N],s[N];
int main()
{
int n,V;
cin>>n>>V;
for(int i=1;i<=n;++i)
cin>>v[i]>>w[i]>>s[i];
for(int i=1;i<=n;++i)
for(int j=v[i];j<=V;++j)
{
for(int k=0;k<=s[i]&&k*v[i]<=j;++k)
f[i][j]=max(f[i][j],f[i-1][j-k*v[i]]+k*w[i]);
}
cout<<f[n][V];
}
2.#include<bits/stdc++.h>
using namespace std;
const int N=1005;
int f[N][N],w[N],v[N],s[N];
int main()
{
int n,V;
cin>>n>>V;
for(int i=1;i<=n;++i)
cin>>v[i]>>w[i]>>s[i];
for(int i=1;i<=n;++i)
for(int j=v[i];j<=V;++j)
{
f[i][j]=f[i-1][j];
for(int k=1;k<=s[i]&&k*v[i]<=j;++k)
f[i][j]=max(f[i][j],f[i-1][j-k*v[i]]+k*w[i]);
}
cout<<f[n][V];
}
两个解法几乎一样,但循环的思路有所区别,第一个是第二个的一般化把,k=0就是f[i-1][j];因此2是可行的
但这两个方法的时间复杂度是大O的(nv^2)计算次数也超过了
那么来试一下能不能用前面的优化方法优化这一个
观察f[i][j-v[i]]=max(f[i-1][j-v[i],f[i-1][j-v[i]*2]+w[i]+f[i-1][j-v[i]3]+2w[i]+…+f[i-1][j-s[i]*v[i]]+(n-1)w[i]+s[i+1]v[i]+nw[i])
而f[i][j]=max(f[i-1][j],f[i-1][j-v[i]]+w[i],f[i-1][j-v[i]2]+2w[i]+f[i-1][j-v[i]3]+3w[i]+…+f[i-1][j-nv[i]]+(n)*w[i];
区别在于这里f[i][j-v[i]]比完全背包中的多了一项,且 没有办法确定这一项是不是最终决定了f[i][j-v[i]]的属性
完全背包中后面有几项由这个j来独立决定而这里s【i】也参与了决定,完全背包中必有f[i][j]比f[i][j-v[i]]多一项f[i-1][j]也即是他们
起始不同但结尾同,就是j不够了
而在这里如果j相对于v[i]足够大会使s[i]成为决定性因素,此时不再有上述性质,反而这里f[i][j],f[i][j-v[i]]都是有s【i】项
所以前面的优化行不通
下面采用常用的二进制优化方法将s次循环改变为logs次
数据2000时间复杂度降低100倍大概是1e7可行
所谓二进制优化就是,起始只需要算出logs个状态就可以表示所有s个状态了,从二进制可以表示所有整数可以简单看出
那么试着把s打包乘几组,然后组合起来一定能够达到不漏!!
这就变成了01背包问题
把N开成1000*log2000=11000多点
ac:
#include<bits/stdc++.h>
using namespace std;
const int N=15000;
int w[N],v[N],f[N],cnt=1;
int main()
{
int n,V;
cin>>n>>V;
for(int i=1;i<=n;++i)
{
int tv,tw,ts;
cin>>tv>>tw>>ts;
for(int k=1;2*k<=ts;k*=2){
v[cnt]=k*tv;
w[cnt++]=k*tw;
ts -=k ;
}
if(ts){
v[cnt]=ts*tv,w[cnt++]=ts*tw;
}
}
for(int i=1;i<=cnt;++i)
for(int j=V;j>=v[i];--j)
f[j]=max(f[j],f[j-v[i]]+w[i]);
cout<<f[V];
}