背包问题是一类经典的优化问题,它通常涉及到的是在限定容量下如何填充一定数量的物品,使得价值最大。
先看例题:
一个旅行者有一个最多能装M 公斤的背包,现在有 n 件物品,求旅行者能获得最大总价值。
贪心代码如下
#include <all>
using namespace std;
long long M,N;
struct wu
{
long long w;//重量
long long v;//价值
};
vector<wu>z;
stack<wu> s;
bool mb(wu w1,wu w2){
return w1.v-w1.w<w2.v-w2.w;
}
int main(){
cin>>M>>N;
z.reserve(N);
for(long long i=0;i<N;i++){
cin>>z[i].w>>z[i].v;//重量,价值
}
sort(z.begin(),z.end(),mb);
for(long long i=0;i<N;i++){
s.push(z[i]);
}
long long m=0,n=0;
while(!s.empty()){
if(s.top().w+m<=M){
m+=s.top().w;
n+=s.top().v;
}
s.pop();
}
cout<<n;
return 0;
}
为什么只得了40分?举个反例:输入
10 3
10 12
2 5
1 6
贪心算法结果为11,正确结果应为12,因为对于第一个物品"利润"为2,第二个物品"利润"为3,第三个物品"利润"为5,按照贪心算法程序会优先选择第三个和第二个物品,于是就装不下第一个物品了,因此会输出5+6,但是此时背包还有7个空空间——没有装满
对于此类题型应使用动态规划的01背包模型
参考代码:
#include <all>
using namespace std;
//f[i][j]=max(f[i-1][j],f[i-1][j-w[i]]+v[i])
long long n, m; // 种类,容量
long long w[35]; // 重量
long long v[35]; // 价值
long long f[35][210]; // 前i件物品,容量为j的最大价值
int main()
{
cin >> m >> n;
for (long long i = 1; i <= n; i++)
{
cin >> w[i] >> v[i];
}
for(long long i=1;i<=n;i++){
cin>>w[i]>>v[i];
}
for(long long i=1;i<=n;i++){
for(long long j=1;j<=m;j++){
if(j>=w[i]){//如果能装下
f[i][j]=max(f[i-1][j],f[i-1][j-w[i]]+v[i]);
}
else f[i][j]=f[i-1][j];
}
}
cout<<f[n][m];
return 0;
}
解析:
状态:最大价值
求解顺序:啊吧啊吧啊吧(先求出相同i的不同j的最大价值)
转移方程:f[i][j]的状态可以由两种情况得到,一种是不选第i个物[j]即为不选i的前i-1个物品可以获取的最大价值;另一种是选第i个物品,则f[i][j]的值为f[i-1][j-w[i]]+v[i],其中w[i]表示第i个物品的重量,v[i]表示第i个物品的价值。
边界条件:若使f[i][j]表示前i件物品背包容量为j时价值最大值,则i=0或j=0时取值为0(前0件物品或背包大小为0时价值为0)
求解目标:f[n][m]
空间优化
因为f的遍历顺序为从第1行开始一行一行遍历,且在遍历第i行时候不会用到第i-2行数据,也就是i-2行及以前的数据没有用了,可以清除。同时,第i-1行的数据每个只会用到一次。遍历每一行时候只用到当前容量j和j-w[i]的数据,也就是第 i 次遍历只需要 第 i-1 次遍历中容量小于等于 j 的数据 。
但如果如果我们顺序枚举,我们的f[i]可能被前面的状态影响.
所以我们考虑倒叙枚举,这样我们的f[i]不会被i以前的状态影响,而我们更新的话也不会影响其他位置的状态.
所以可以改为
#include <all>
using namespace std;
//f[i][j]=max(f[i-1][j],f[i-1][j-w[i]]+v[i])
long long n, m; // 种类,容量
long long w[1010]; // 重量
long long v[1010]; // 价值
long long f[1010]; // 前i件物品,容量为j的最大价值
int main()
{
cin >> m >> n;
for(long long i=1;i<=n;i++){
cin>>w[i]>>v[i];
}
for(long long i=1;i<=n;i++){
for(long long j=m;j>=w[i];j--){//倒序枚举
f[j]=max(f[j],f[j-w[i]]+v[i]);
}
}
cout<<f[m];
return 0;
}