背包问题之01背包 全详解(最浅显易懂)

01背包是一种非常经典的动态规划问题,这里对01背包问题进行详细解读。


01背包问题题目描述

N N N 件物品和一个容量为 V V V 的背包。第 i i i件物品的体积是 c [ i ] c[i] c[i] ,价值是 w [ i ] w[i] w[i] ,求将哪些物品装入背包可使价值总和最大。

01背包问题解析

对于所有的动态规划问题,第一步都是确定状态。我们定义状态 d p [ i ] [ j ] dp[i][j] dp[i][j] 是表示目前正在枚举第 i i i 个物品,目前已取的总体积为 j j j,最大价值为 d p [ i ] [ j ] dp[i][j] dp[i][j]
第二步自然是找状态转移方程
首先我们注意一点:物品只有取和不取两种选择,这是符合我们日常生活的。状态转移方程就需要从这里为突破口来思考:
(1):假如我们不取这个物品,那么 d p [ i ] [ j ] dp[i][j] dp[i][j] 肯定是能从上一个物品,同样体积转移过来的,所以 d p [ i ] [ j ] = d p [ i − 1 ] [ j ] dp[i][j] = dp[i - 1][j] dp[i][j]=dp[i1][j]
(2):假如我们取这个物品,那么 d p [ i ] [ j ] dp[i][j] dp[i][j] 是从什么情况转移呢?思考一下,首先可以得出, d p [ i ] [ j ] dp[i][j] dp[i][j] 肯定可以从上一个物品转移过来,那可以从什么体积转移呢?我们注意到,对于 d p [ i ] [ j ] dp[i][j] dp[i][j] ,它的当前取到的体积为 j j j,由于我们取了这个物品,所以上一个物品的体积为 j − c [ i ] j - c[i] jc[i]。所以我们可以得出, d p [ i ] [ j ] dp[i][j] dp[i][j] 可以从 d p [ i − 1 ] [ j − c [ i ] ] dp[i - 1][j - c[i]] dp[i1][jc[i]] 转移过来,由于我们取了这个物品,所以还要加上这个物品的价值 w [ i ] w[i] w[i]
所以我们可以得出01背包的状态转移方程(很重要,尽量背下来!!)
j > = c [ i ] j >= c[i] j>=c[i] 时, d p [ i ] [ j ] = m a x ( d p [ i − 1 ] [ j ] , d p [ i − 1 ] [ j − c [ i ] ] + w [ i ] ) dp[i][j] = max(dp[i - 1][j], dp[i - 1][j - c[i]] + w[i]) dp[i][j]=max(dp[i1][j],dp[i1][jc[i]]+w[i])
j < c [ i ] j j<c[i] 时, d p [ i ] [ j ] = d p [ i − 1 ] [ j ] dp[i][j] = dp[i - 1][j] dp[i][j]=dp[i1][j]

上代码:

#include 

#include 

#include 

using namespace std;

int dp[21][1010];
int w[21], c[21];
int main() {
    int N, V;
    cin >> N >> V;
    for (int i = 1; i <= N; i++) {
        cin >> w[i] >> c[i];
    }
    for (int i = 1; i <= N; i++) {
        for (int j = 0; j <= V; j++) {
            if (j >= c[i]) {
                dp[i][j] = max(dp[i - 1][j - c[i]] + w[i], dp[i - 1][j]);
            } else {
                dp[i][j] = dp[i - 1][j];
            }
        }
    }
    cout << dp[N][V] << endl;
    return 0;
}

01背包滚动数组空间优化

我们分析一下空间复杂度: O ( N V ) O(NV) O(NV),显然较大,我们要优化一下,用什么优化呢?
我们可以用滚动数组!
观察转移方程,易得 d p [ i ] [ j ] dp[i][j] dp[i][j] 只从 d p [ i − 1 ] [ j ] dp[i - 1][j] dp[i1][j] d p [ i − 1 ] [ j − c [ i ] ] dp[i - 1][j - c[i]] dp[i1][jc[i]]转移,所以可以用一个 f l a g flag flag代替 i i i,用 1 − f l a g 1 - flag 1flag 代替 i − 1 i - 1 i1

这样只需要定义 d p [ 2 ] [ m a x n ] dp[2][maxn] dp[2][maxn] ,大大节省了空间。

上代码:

#include 
#include 
#include 
using namespace std;

int dp[2][1010];
int w[21], c[21];

int main() {
    int N, V;
    cin >> N >> V;
    for (int i = 1; i <= N; i++) {
        cin >> w[i] >> c[i];
    }
    int flag = 1;
    for (int i = 1; i <= N; i++) {
        for(int j = 0;j <= V; j++) {
            if(j >= c[i]) {
                dp[flag][j] = max(dp[1 - flag][j - c[i]] + w[i], dp[1 - flag][j]);
	    } else {
                dp[flag][j] = dp[1 - flag][j];
	    }
	}
        flag = 1 - flag;
    }
    cout << dp[1 - flag][V] << endl;
    return 0;
}

01背包空间优化

这里给大家介绍真正的空间优化。

如果我们将 d p dp dp 数组只用来表示体积,那么我们可以让内层循环的 j j j V V V 0 0 0 枚举,那么当前状态转移方程的 d p [ j ] dp[j] dp[j] d p [ j − c [ i ] ] dp[j - c[i]] dp[jc[i]] 由于我们没有更新,所以仍然是计算上一轮 i − 1 i - 1 i1 个物品的,就是二维状态下的 d p [ 1 − f l a g ] [ j ] dp[1 - flag][j] dp[1flag][j] d p [ i − f l a g ] [ j − c [ i ] ] dp[i - flag][j - c[i]] dp[iflag][jc[i]]。所以现在我们的转移方程是:
d p [ j ] = m a x ( d p [ j ] , d p [ j − c [ i ] ] + w [ i ] ) dp[j] = max(dp[j], dp[j - c[i]] + w[i]) dp[j]=max(dp[j],dp[jc[i]]+w[i])

上代码:

#include 
#include 
using namespace std;

int dp[1010];
int w[21], c[21];

int main() {
    int N, V;
    cin >> N >> V;
    for (int i = 1; i <= N; i++) {
        cin >> w[i] >> c[i];
    }
    for (int i = 1; i <= N; i++) {
		for (int j = V; j >= c[i]; j--) {
            dp[j] = max(dp[j - c[i]] + w[i], dp[j]);
		}
    }
    cout << dp[V] << endl;
    return 0;
}

课后习题

采药
开心的金明

你可能感兴趣的:(动态规划专题讲解,上课笔记)