完全背包问题

本文介绍了完全背包问题和零钱兑换问题的动态规划方法,包括状态转移方程、代码实现以及空间优化策略。完全背包问题中物品可以无限次选取,而零钱兑换问题关注的是最少硬币数量。
摘要由CSDN通过智能技术生成

'); --md-admonition-pythontutor-color: #eee;">

给定 \(n\) 个物品,第 \(i\) 个物品的重量为 \(wgt[i-1]\)、价值为 \(val[i-1]\) ,和一个容量为 \(cap\) 的背包。每个物品可以重复选取,问在限定背包容量下能放入物品的最大价值。示例如图 14-22 所示。

完全背包问题的示例数据

图 14-22   完全背包问题的示例数据

1.   动态规划思路

完全背包问题和 0-1 背包问题非常相似,区别仅在于不限制物品的选择次数

  • 在 0-1 背包问题中,每种物品只有一个,因此将物品 \(i\) 放入背包后,只能从前 \(i-1\) 个物品中选择。
  • 在完全背包问题中,每种物品的数量是无限的,因此将物品 \(i\) 放入背包后,仍可以从前 \(i\) 个物品中选择

在完全背包问题的规定下,状态 \([i, c]\) 的变化分为两种情况。

  • 不放入物品 \(i\) :与 0-1 背包问题相同,转移至 \([i-1, c]\) 。
  • 放入物品 \(i\) :与 0-1 背包问题不同,转移至 \([i, c-wgt[i-1]]\) 。

从而状态转移方程变为:

\[ dp[i, c] = \max(dp[i-1, c], dp[i, c - wgt[i-1]] + val[i-1]) \]

2.   代码实现

对比两道题目的代码,状态转移中有一处从 \(i-1\) 变为 \(i\) ,其余完全一致:

C++

unbounded_knapsack.cpp

/* 完全背包:动态规划 */
int unboundedKnapsackDP(vector<int> &wgt, vector<int> &val, int cap) {
    int n = wgt.size();
    // 初始化 dp 表
    vector<vector<int>> dp(n + 1, vector<int>(cap + 1, 0));
    // 状态转移
    for (int i = 1; i <= n; i++) {
        for (int c = 1; c <= cap; c++) {
            if (wgt[i - 1] > c) {
                // 若超过背包容量,则不选物品 i
                dp[i][c] = dp[i - 1][c];
            } else {
                // 不选和选物品 i 这两种方案的较大值
                dp[i][c] = max(dp[i - 1][c], dp[i][c - wgt[i - 1]] + val[i - 1]);
            }
        }
    }
    return dp[n][cap];
}

可视化运行

3.   空间优化

由于当前状态是从左边和上边的状态转移而来的,因此空间优化后应该对 \(dp\) 表中的每一行进行正序遍历

这个遍历顺序与 0-1 背包正好相反。请借助图 14-23 来理解两者的区别。

<1><2><3><4><5><6>

unbounded_knapsack_dp_comp_step3

图 14-23   完全背包问题在空间优化后的动态规划过程

代码实现比较简单,仅需将数组 dp 的第一维删除:C++

unbounded_knapsack.cpp

/* 完全背包:空间优化后的动态规划 */
int unboundedKnapsackDPComp(vector<int> &wgt, vector<int> &val, int cap) {
    int n = wgt.size();
    // 初始化 dp 表
    vector<int> dp(cap + 1, 0);
    // 状态转移
    for (int i = 1; i <= n; i++) {
        for (int c = 1; c <= cap; c++) {
            if (wgt[i - 1] > c) {
                // 若超过背包容量,则不选物品 i
                dp[c] = dp[c];
            } else {
                // 不选和选物品 i 这两种方案的较大值
                dp[c] = max(dp[c], dp[c - wgt[i - 1]] + val[i - 1]);
            }
        }
    }
    return dp[cap];
}

可视化运行

14.5.2   零钱兑换问题

背包问题是一大类动态规划问题的代表,其拥有很多变种,例如零钱兑换问题。

Question

给定 \(n\) 种硬币,第 \(i\) 种硬币的面值为 \(coins[i - 1]\) ,目标金额为 \(amt\) ,每种硬币可以重复选取,问能够凑出目标金额的最少硬币数量。如果无法凑出目标金额,则返回 \(-1\) 。示例如图 14-24 所示。

零钱兑换问题的示例数据

图 14-24   零钱兑换问题的示例数据

1.   动态规划思路

零钱兑换可以看作完全背包问题的一种特殊情况,两者具有以下联系与不同点。

  • 两道题可以相互转换,“物品”对应“硬币”、“物品重量”对应“硬币面值”、“背包容量”对应“目标金额”。
  • 优化目标相反,完全背包问题是要最大化物品价值,零钱兑换问题是要最小化硬币数量。
  • 完全背包问题是求“不超过”背包容量下的解,零钱兑换是求“恰好”凑到目标金额的解。

第一步:思考每轮的决策,定义状态,从而得到 \(dp\) 表

状态 \([i, a]\) 对应的子问题为:前 \(i\) 种硬币能够凑出金额 \(a\) 的最少硬币数量,记为 \(dp[i, a]\) 。

二维 \(dp\) 表的尺寸为 \((n+1) \times (amt+1)\) 。

第二步:找出最优子结构,进而推导出状态转移方程

本题与完全背包问题的状态转移方程存在以下两点差异。

  • 本题要求最小值,因此需将运算符 \(\max()\) 更改为 \(\min()\) 。
  • 优化主体是硬币数量而非商品价值,因此在选中硬币时执行 \(+1\) 即可。

\[ dp[i, a] = \min(dp[i-1, a], dp[i, a - coins[i-1]] + 1) \]

第三步:确定边界条件和状态转移顺序

当目标金额为 \(0\) 时,凑出它的最少硬币数量为 \(0\) ,即首列所有 \(dp[i, 0]\) 都等于 \(0\) 。

当无硬币时,无法凑出任意 \(> 0\) 的目标金额,即是无效解。为使状态转移方程中的 \(\min()\) 函数能够识别并过滤无效解,我们考虑使用 \(+ \infty\) 来表示它们,即令首行所有 \(dp[0, a]\) 都等于 \(+ \infty\) 。

2.   代码实现

大多数编程语言并未提供 \(+ \infty\) 变量,只能使用整型 int 的最大值来代替。而这又会导致大数越界:状态转移方程中的 \(+ 1\) 操作可能发生溢出。

为此,我们采用数字 \(amt + 1\) 来表示无效解,因为凑出 \(amt\) 的硬币数量最多为 \(amt\) 。最后返回前,判断 \(dp[n, amt]\) 是否等于 \(amt + 1\) ,若是则返回 \(-1\) ,代表无法凑出目标金额。代码如下所示:

PythonC++JavaC#GoSwiftJSTSDartRustCZig

coin_change.cpp

/* 零钱兑换:动态规划 */
int coinChangeDP(vector<int> &coins, int amt) {
    int n = coins.size();
    int MAX = amt + 1;
    // 初始化 dp 表
    vector<vector<int>> dp(n + 1, vector<int>(amt + 1, 0));
    // 状态转移:首行首列
    for (int a = 1; a <= amt; a++) {
        dp[0][a] = MAX;
    }
    // 状态转移:其余行和列
    for (int i = 1; i <= n; i++) {
        for (int a = 1; a <= amt; a++) {
            if (coins[i - 1] > a) {
                // 若超过目标金额,则不选硬币 i
                dp[i][a] = dp[i - 1][a];
            } else {
                // 不选和选硬币 i 这两种方案的较小值
                dp[i][a] = min(dp[i - 1][a], dp[i][a - coins[i - 1]] + 1);
            }
        }
    }
    return dp[n][amt] != MAX ? dp[n][amt] : -1;
}

可视化运行

图 14-25 展示了零钱兑换的动态规划过程,和完全背包问题非常相似。

<1><2><3><4><5><6><7><8><9><10><11><12><13><14><15>

coin_change_dp_step3

图 14-25   零钱兑换问题的动态规划过程

3.   空间优化

零钱兑换的空间优化的处理方式和完全背包问题一致:

PythonC++JavaC#GoSwiftJSTSDartRustCZig

coin_change.cpp

/* 零钱兑换:空间优化后的动态规划 */
int coinChangeDPComp(vector<int> &coins, int amt) {
    int n = coins.size();
    int MAX = amt + 1;
    // 初始化 dp 表
    vector<int> dp(amt + 1, MAX);
    dp[0] = 0;
    // 状态转移
    for (int i = 1; i <= n; i++) {
        for (int a = 1; a <= amt; a++) {
            if (coins[i - 1] > a) {
                // 若超过目标金额,则不选硬币 i
                dp[a] = dp[a];
            } else {
                // 不选和选硬币 i 这两种方案的较小值
                dp[a] = min(dp[a], dp[a - coins[i - 1]] + 1);
            }
        }
    }
    return dp[amt] != MAX ? dp[amt] : -1;
}

可视化运行

14.5.3   零钱兑换问题 II

Question

给定 \(n\) 种硬币,第 \(i\) 种硬币的面值为 \(coins[i - 1]\) ,目标金额为 \(amt\) ,每种硬币可以重复选取,问凑出目标金额的硬币组合数量。示例如图 14-26 所示。

零钱兑换问题 II 的示例数据

图 14-26   零钱兑换问题 II 的示例数据

1.   动态规划思路

相比于上一题,本题目标是求组合数量,因此子问题变为:前 \(i\) 种硬币能够凑出金额 \(a\) 的组合数量。而 \(dp\) 表仍然是尺寸为 \((n+1) \times (amt + 1)\) 的二维矩阵。

当前状态的组合数量等于不选当前硬币与选当前硬币这两种决策的组合数量之和。状态转移方程为:

\[ dp[i, a] = dp[i-1, a] + dp[i, a - coins[i-1]] \]

当目标金额为 \(0\) 时,无须选择任何硬币即可凑出目标金额,因此应将首列所有 \(dp[i, 0]\) 都初始化为 \(1\) 。当无硬币时,无法凑出任何 \(>0\) 的目标金额,因此首行所有 \(dp[0, a]\) 都等于 \(0\) 。

2.   代码实现

PythonC++JavaC#GoSwiftJSTSDartRustCZig

coin_change_ii.cpp

/* 零钱兑换 II:动态规划 */
int coinChangeIIDP(vector<int> &coins, int amt) {
    int n = coins.size();
    // 初始化 dp 表
    vector<vector<int>> dp(n + 1, vector<int>(amt + 1, 0));
    // 初始化首列
    for (int i = 0; i <= n; i++) {
        dp[i][0] = 1;
    }
    // 状态转移
    for (int i = 1; i <= n; i++) {
        for (int a = 1; a <= amt; a++) {
            if (coins[i - 1] > a) {
                // 若超过目标金额,则不选硬币 i
                dp[i][a] = dp[i - 1][a];
            } else {
                // 不选和选硬币 i 这两种方案之和
                dp[i][a] = dp[i - 1][a] + dp[i][a - coins[i - 1]];
            }
        }
    }
    return dp[n][amt];
}

可视化运行

3.   空间优化

空间优化处理方式相同,删除硬币维度即可:

PythonC++JavaC#GoSwiftJSTSDartRustCZig

coin_change_ii.cpp

/* 零钱兑换 II:空间优化后的动态规划 */
int coinChangeIIDPComp(vector<int> &coins, int amt) {
    int n = coins.size();
    // 初始化 dp 表
    vector<int> dp(amt + 1, 0);
    dp[0] = 1;
    // 状态转移
    for (int i = 1; i <= n; i++) {
        for (int a = 1; a <= amt; a++) {
            if (coins[i - 1] > a) {
                // 若超过目标金额,则不选硬币 i
                dp[a] = dp[a];
            } else {
                // 不选和选硬币 i 这两种方案之和
                dp[a] = dp[a] + dp[a - coins[i - 1]];
            }
        }
    }
    return dp[amt];
}

欢迎在评论区留下你的见解

  • 25
    点赞
  • 15
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值