有 N 种物品和一个容量是 V 的背包,
(不同之处)每种物品都有无限件可用。
第 i 种物品的体积是 vi,价值是 wi。
求解将哪些物品装入背包,可使这些物品的总体积不超过背包容量,且总价值最大。
输出最大价值。
输入格式
第一行两个整数,N,V,用空格隔开,分别表示物品种数和背包容积。
接下来有 N 行,每行两个整数 vi,wi,用空格隔开,分别表示第 i 种物品的体积和价值。
输出格式
输出一个整数,表示最大价值。
数据范围
0<N,V≤1000
0<vi,wi≤1000
输入样例
4 5
1 2
2 4
3 4
4 5
输出样例:
10
f[i-1] [j-k*v[i]] 表示第i个物品选了k个,需要再加上第i-1个物品的最优解,第i-1个物品的最优解是f[i-1][j-k*v[i]]
这是二维朴素做法,时间复杂度O(n^3) =10^9 超时
#include<iostream>
using namespace std;
int n,m;
const int N=1010;
int f[N][N],v[N],w[N];
int main()
{
cin>>n>>m;
for(int i=1;i<=n;i++)
cin>>v[i]>>w[i];
for(int i=1;i<=n;i++)
for(int j=0;j<=m;j++)
{
for(int k=0;k*v[i]<=j;k++)//表示第i个物品选几个(从0个~k个 体积大于为止)
f[i][j]=max(f[i][j],f[i-1][j-k*v[i]]+k*w[i]);
}
cout<<f[n][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-1][j-v]=max( f[i-1][j-v], f[i-1][j-v]+w f[i-1][]j-2v]+2w)
所以推得
f[i][j]=max(f[i][j-v]+w, f[i-1][j])
有了上述关系,k循环其实就可以不要了,核心代码优化:
for(int j=1;i<=n;j++)
for(int j=0;j<=m;j++)
{
f[i][j]=f[i-1][j];
if(j-v[i]>=0)
f[i][j]=max(f[i][j].f[i][j-v[i]]+w[i]);
}
进一步优化 j 的初始化
for(int j=1;i<=n;j++)
for(int j=v[i];j<=m;j++)
f[i][j]=max(f[i][j].f[i][j-v[i]]+w[i]);
写法和01背包极其像 (仅下标不同)
对比:
f[i][j] = max(f[i][j],f[i-1][j-v[i]]+w[i]);//01背包
f[i][j] = max(f[i][j],f[i] [j-v[i]]+w[i]);//完全背包问题
所以
进一步优化
for(int j=1;i<=n;j++)
for(int j=v[i];j<=m;j++)
f[j]=max(f[j].f[j-v[i]]+w[i]);
//f[i][j]=max(f[i][j].f[i][j-v[i]]+w[i]);
因为所有的式子都是从右边开始算,j-v[i]<j,所以先算j-v[i],不会影响左边的计算
所以最终代码为
#include<iostream>
using namespace std;
const int N = 1010;
int f[N];
int v[N],w[N];
int main()
{
int n,m;
cin>>n>>m;
for(int i=1;i<=n;i++)
cin>>v[i]>>w[i];
for(int i=1;i<=n;i++)
for(int j=v[i];j<=m;j++)
f[j]=max(f[j],f[j-v[i]]+w[i]);
cout<<f[m]<<endl;
return 0;
}