算法——动态规划:0/1 背包问题

文章目录

  • 一、问题描述
  • 二、解决方案
    • 1. DP 状态的设计
    • 2. 状态转移方程
    • 3. 算法复杂度
    • 4. 举例
    • 5. 实现
    • 6. 滚动数组
      • 6.1 两行实现
      • 6.2 单行实现
      • 6.3 优缺点
  • 三、总结


一、问题描述

问题的抽象:给定 n n n 种物品和一个背包,第 i i i 个物品的体积为 c i c_i ci,价值为 w i w_i wi,背包的总容量为 C C C。把物品装入背包时,第 i i i 种物品只有两种选择:装入背包或不装入背包。如何选择装入背包的物品,使装入背包中的物品的总价值最大?

具体的问题可以看这道洛谷题:P1048 [NOIP2005 普及组] 采药,将 物品 换成了 草药,将 容量 换成了 时间,将 背包的容量 换成了 规定的时间

二、解决方案

0/1 背包问题 是一道 经典 的使用 动态规划 思想的题目,同时也不是很难,掌握它之后就正式跨入 动态规划 的大门了。

1. DP 状态的设计

引入一个 ( N + 1 ) × ( C + 1 ) (N + 1) \times (C + 1) (N+1)×(C+1) 的二维数组 d p [ ] [ ] dp[][] dp[][],称为 DP 状态。其中, d p [ i ] [ j ] dp[i][j] dp[i][j] 表示把前 i i i 个(从第 1 1 1 个到第 i i i 个)物品装入容量为 j j j 的背包中获得的最大价值。

每个 d p [ i ] [ j ] dp[i][j] dp[i][j] 都是一个 0/1 背包问题:将前 i i i 个物品装入容量为 j j j 的背包。所以 d p [ N ] [ C ] dp[N][C] dp[N][C] 就代表将前 N N N 个物品装入容量为 C C C 的背包。

2. 状态转移方程

状态转移方程指的是 从一个状态转变到另一个状态的递推公式

一般使用 自底向上 的 递推 来解决背包问题,假设已经递推到 d p [ i ] [ j ] dp[i][j] dp[i][j],分两种情况:

  • i i i 个物品的体积比容量 j j j 还大,不能装进容量为 j j j 的背包。此时直接继承前 i − 1 i - 1 i1 个物品装进容量 j j j 的背包的情况即可,即 d p [ i ] [ j ] = d p [ i − 1 ] [ j ] dp[i][j] = dp[i - 1][j] dp[i][j]=dp[i1][j]
  • i i i 个物品的体积比容量 j j j 小,能装进背包。此时就体现 0/1 了( 0 0 0 表示不装, 1 1 1 表示装):
    • 装第 i i i 个物品。先在 前 i − 1 i - 1 i1 个物品装进容量为 j j j 的背包中,给第 i i i 个物品空出 c i c_i ci 的空间,从而得到背包 d p [ i − 1 ] [ j − c [ i ] ] dp[i - 1][j - c[i]] dp[i1][jc[i]]。然后将第 i i i 个物品放入这个背包,背包的价值增加 w i w_i wi,从而有当前背包 d p [ i ] [ j ] = d p [ i − 1 ] [ j − c [ i ] ] + w [ i ] dp[i][j] = dp[i - 1][j - c[i]] + w[i] dp[i][j]=dp[i1][jc[i]]+w[i]
    • 不装第 i i i 个物品。直接继承前 i − 1 i - 1 i1 个物品装进容量 j j j 的背包的情况即可,即 d p [ i ] [ j ] = d p [ i − 1 ] [ j ] dp[i][j] = dp[i - 1][j] dp[i][j]=dp[i1][j]
    • 此时取这两者中的较大值作为当前背包的价值,状态转移方程为: 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])

所以 0/1 背包问题 的特征如下:

  • 重叠子问题是 d p [ i ] [ j ] dp[i][j] dp[i][j]
  • 最优子结构是 d p [ i ] [ j ] dp[i][j] dp[i][j] 的状态转移方程。

3. 算法复杂度

  • 时间复杂度:在计算时,需要计算二维矩阵 d p [ ] [ ] dp[][] dp[][] 中的每一个值,矩阵的大小是 N C NC NC,每次计算的时间为 O ( 1 ) O(1) O(1),所以时间复杂度为 O ( N C ) O(NC) O(NC)
  • 空间复杂度:使用了大小为 N C NC NC 的二维数组,所以空间复杂度为 O ( N C ) O(NC) O(NC)

说明:如果 N , C N, C N,C 很大,则 N + 1 , C + 1 N + 1, C + 1 N+1,C+1 可以近似地看作 N , C N, C N,C,为了简化,所以复杂度不是 O ( ( N + 1 ) × ( C + 1 ) ) O((N + 1) \times (C + 1)) O((N+1)×(C+1))

4. 举例

初学者可能很难理解,这时举一个实际案例就懂了。对于 P1048 题,有一个测试用例:背包的容量为 70,物品的数量为 3,物品的体积数组为 [71, 69, 1],物品的价值数组为 [100, 1, 2]

定义一个 4 * 71 的二维数组,如下所示(由于长度限制,中间的数全部用 ... 代替,... 代表的数 和 ... 两边的数相同):
alt text
先填充索引为 1 的行(第一个物品的体积为 71,价值为 100):
alt text
然后填充索引为 2 的行(第二个物品的体积为 69,价值为 1):
alt text
接着填充索引为 3 的行(第三个物品的体积为 1,价值为 2):
alt text
最终,数组的 dp[3][70] 的位置存储着题目的结果——将前 3 个物品放入容量为 70 的背包的最大价值。

5. 实现

// 虽然这些代码看起来是 C 语言的代码,但如果选择 C 语言,则会编译失败。建议选择 C++14
#include <stdio.h>const int MAX_N = 105; // 最大的物品数量,与题目有关
const int MAX_C = 1005; // 最大的背包容量,与题目有关
int N, C; // N 是物品的数量,C 是背包的容量
int c[MAX_N], w[MAX_N]; // c[i], w[i] 分别表示第 i 个物品的 体积 和 价值
int dp[MAX_N][MAX_C]; // dp[i][j] 表示将前 i 个物品装入容量为 j 的背包中// 取 x 和 y 中的较大值
int max(int x, int y) {return x > y ? x : y;
}// 进行动态规划
int programming() {for (int i = 1; i <= N; i++) {for (int j = 1; j <= C; j++) {if (j < c[i]) { // 装不下第 i 个物品dp[i][j] = dp[i - 1][j];} else { // 能装下第 i 个物品dp[i][j] = max(dp[i - 1][j - c[i]] + w[i], // 装第 i 个物品dp[i - 1][j] // 不装第 i 个物品);}}}return dp[N][C];
}int main() {// 读取数据scanf("%d %d", &C, &N);for (int i = 1; i <= N; i++) {scanf("%d %d", &c[i], &w[i]);}// 进行动态规划printf("%d", programming());return 0;
}

6. 滚动数组

滚动数组动态规划 最常用的 空间优化技术。动态规划的状态方程通常是二维即以上,占用了太多空间,用滚动数组可以 极大程度上 减少空间占用,它能把 二维 状态方程的空间复杂度优化到 一维,更高维的数组也可以通过优化,减少一个维度。

从状态转移方程 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]) 中可以看出, d p [ i ] [ ] dp[i][] dp[i][] 只与 d p [ i − 1 ] [ ] dp[i - 1][] dp[i1][] 有关,与之前的 d p [ i − 2 ] [ ] , d p [ i − 3 ] [ ] , ⋯ dp[i - 2][], dp[i - 3][], \cdots dp[i2][],dp[i3][], 都无关。由于它们都没用了,索性就 复用 它们占用的空间,用新的一行覆盖已经无用的一行,只需要两行就够了。这就是“滚动”的含义。

滚动数组根据使用的行数不同,分为两种实现方式:

6.1 两行实现

仍然使用二维数组,不过不是 ( N + 1 ) × ( C + 1 ) (N + 1) \times (C + 1) (N+1)×(C+1) 的二维数组,而是 2 × ( C + 1 ) 2 \times (C + 1) 2×(C+1) 的二维数组,这就是所谓的“两行”:计算 本行 时用 上一行 的结果,然后将 本行 和 上一行 互换,计算 新的本行(新的本行 使用了 原来上一行 的内存空间)时使用 新的上一行(新的上一行 使用了 原来本行 的内存空间)。

#include <stdio.h>const int MAX_N = 105; // 最大的物品数量,与题目有关
const int MAX_C = 1005; // 最大的背包容量,与题目有关
int N, C; // N 是物品的数量,C 是背包的容量
int c[MAX_N], w[MAX_N]; // c[i], w[i] 分别表示第 i 个物品的 体积 和 价值
int dp[2][MAX_C]; // 使用两行滚动数组,dp[i][j] 成为不断更新的变量,没有固定的含义// 取 x 和 y 中的较大值
int max(int x, int y) {return x > y ? x : y;
}// 交换 i 和 j 的值
// 注意 &,它表示引用传递,使用传入的实参,而不是重新创建一个新变量来代表实参,这是 C++ 的特性
void swap(int &i, int &j) {int temp = i;i = j;j = temp;
}// 进行动态规划
int programming() {/*在状态转移方程 dp[i][j] = max(dp[i - 1][j], dp[i - 1][j - c[i]] + w[i]) 中,位于 dp 的行索引的 i 用 curr 来代替,i - 1 用 prev 来代替prev, curr 分别是 上一行 和 本行 索引*/int prev = 0, curr = 1;for (int i = 1; i <= N; i++) {swap(prev, curr); // 将本行和上一行互换for (int j = 1; j <= C; j++) {if (j < c[i]) { // 装不下第 i 个物品dp[curr][j] = dp[prev][j];} else { // 能装下第 i 个物品dp[curr][j] = max(dp[prev][j - c[i]] + w[i], // 装第 i 个物品dp[prev][j] // 不装第 i 个物品);}}}return dp[curr][C]; // 返回当前行的最后一个值
}int main() {// 读取数据scanf("%d %d", &C, &N);for (int i = 1; i <= N; i++) {scanf("%d %d", &c[i], &w[i]);}// 进行动态规划printf("%d", programming());return 0;
}

6.2 单行实现

实际上还能省,发现计算 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]] d p [ i − 1 ] [ j ] dp[i - 1][j] dp[i1][j],所以还可以复用 上一行 中的部分内存空间。

由于 j − c [ i ] j - c[i] jc[i] 不是一个确定的值,取值范围为 [ 0 , j ) [0, j) [0,j),所以 无法 在计算出 d p [ i ] [ j ] dp[i][j] dp[i][j] 的结果前 复用 上一行 中 j j j 之前的内存空间。既然如此,就考虑 复用 上一行 中 j j j 之后的内存空间。此时能想到:如果 从后向前计算 d p [ i ] [ j ] dp[i][j] dp[i][j],就不会影响 j j j 之前的值。进而就有 单行 的滚动数组实现:

#include <stdio.h>const int MAX_N = 105; // 最大的物品数量,与题目有关
const int MAX_C = 1005; // 最大的背包容量,与题目有关
int N, C; // N 是物品的数量,C 是背包的容量
int c[MAX_N], w[MAX_N]; // c[i], w[i] 分别表示第 i 个物品的 体积 和 价值
int dp[MAX_C]; // 使用单行滚动数组,dp[j] 成为不断更新的变量,没有固定的含义// 取 x 和 y 中的较大值
int max(int x, int y) {return x > y ? x : y;
}// 进行动态规划
int programming() {for (int i = 1; i <= N; i++) {for (int j = C; j >= 1; j--) { // 切记要从后往前计算if (j < c[i]) { // 装不下第 i 个物品dp[j] = dp[j];} else { // 能装下第 i 个物品dp[j] = max(dp[j - c[i]] + w[i], // 装第 i 个物品dp[j] // 不装第 i 个物品);}}}return dp[C]; // 返回单行的最后一个值
}int main() {// 读取数据scanf("%d %d", &C, &N);for (int i = 1; i <= N; i++) {scanf("%d %d", &c[i], &w[i]);}// 进行动态规划printf("%d", programming());return 0;
}

6.3 优缺点

  • 优点:优化了空间复杂度,很大程度上减少了内存的使用。
  • 缺点
    • 由于复用了空间,从而覆盖了之前的计算结果,数组中的值没有固定的实际含义,难以理解,通常都是先写出常规的动态规划,再使用滚动数组进行优化。
    • 还是因为复用了空间,导致只留下最后一次计算的状态,无法从数组中看出具体的方案。

三、总结

0/1 背包问题0/1 在于每种物品只有两种状态——放进背包 和 不放进背包,针对这一问题,使用了 动态规划 的解决方案,由于空间复杂度比较高,所以还使用 滚动数组 的思想进行优化,从而将占用的空间 降维。初学者适合用 两行 的实现,不容易犯错。熟练掌握动态规划后,可以使用 单行 的实现。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/399330.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

k8s分布式存储-ceph

文章目录 Cephdeploy-ceph部署1.系统环境初始化1.1 修改主机名&#xff0c;DNS解析1.2 时间同步1.3 配置apt基础源与ceph源1.4关闭selinux与防火墙1.5 **创建** ceph **集群部署用户** cephadmin1.6分发密钥 2. ceph部署2.1 **安装** ceph 部署工具2.2 **初始化** mon **节点**…

子串 前缀和 | Java | (hot100) 力扣560. 和为K的子数组

560. 和为K的子数组 暴力法&#xff08;连暴力法都没想出来……&#xff09; class Solution {public int subarraySum(int[] nums, int k) {int count0;int len nums.length;for(int i0; i<len; i) {int sum0;for(int ji; j<len; j) {sumnums[j];if(sum k) {count;}…

mysql注入-字符编码技巧

一、环境搭建 创建数据表 CREATE TABLE mysql_Bian_Man (id int(10) unsigned NOT NULL AUTO_INCREMENT,username varchar(255) COLLATE latin1_general_ci NOT NULL,password varchar(255) COLLATE latin1_general_ci NOT NULL,PRIMARY KEY (id) ) ENGINEMyISAM AUTO_INCREME…

Redis 缓存预热、雪崩、穿透、击穿

缓存预热 缓存预热是什么 缓存预热就是系统上线后&#xff0c;提前将相关的缓存数据直接加载到缓存系统。避免在用户请求的时候&#xff0c;先查询数据库&#xff0c;然后再将数据缓存的问题&#xff01;用户直接查询事先被预热的缓存数据&#xff01;解决方案 使用 PostConstr…

LeetCode旋转图像

题目描述&#xff1a; 给定一个 n n 的二维矩阵 matrix 表示一个图像。请你将图像顺时针旋转 90 度。 你必须在 原地 旋转图像&#xff0c;这意味着你需要直接修改输入的二维矩阵。请不要使用另一个矩阵来旋转图像。 示例 1&#xff1a; 输入&#xff1a;matrix [[1,2,3]…

opencv实战项目七:帧差法获取运动车辆蒙版

文章目录 一、简介二、实现方案三、算法实现步骤3.1 取出视频中间帧3.2 帧差法形成运动蒙版&#xff1a; 四、代码整体实现五&#xff1a;效果 一、简介 在智能交通系统中&#xff0c;车辆检测是一项重要的技术&#xff0c;而通常情况下一张图片中无用信息过多会带来额外的计算…

Linux--C语言之循环结构

文章目录 一、循环结构&#xff08;一&#xff09;循环的概念&#xff08;二&#xff09;循环的类型&#xff08;三&#xff09;循环的构成&#xff08;四&#xff09;当型循环的实现while死循环 &#xff08;五&#xff09;for...总结死循环 &#xff08;七&#xff09;循环实…

机器学习——逻辑回归(学习笔记)

目录 一、认识逻辑回归 二、二元逻辑回归&#xff08;LogisticRegression&#xff09; 1. 损失函数 2. 正则化 3. 梯度下降 4. 二元回归与多元回归 三、sklearn中的逻辑回归&#xff08;自查&#xff09; 1. 分类 2. 参数列表 3. 属性列表 4. 接口列表 四、逻辑回归…

大厂面试题分享第一期

大厂面试题分享第一期 Redis持久化方式AOF优缺点RDB优缺点 如何保证Redis和Myql的一致性索引下推输入url到浏览器发生了什么ReentranLock底层原理SpringBoot 的启动流程 Redis持久化方式 Redis提供了两种主要的持久化机制&#xff0c;分别是AOF&#xff08;Append-Only File&a…

Python 数据可视化,怎么选出合适数据的图表

数据可视化最佳实践 1. 引言&#xff1a;为什么数据可视化最佳实践很重要 数据可视化是数据分析和决策过程中不可或缺的一部分。通过有效的可视化&#xff0c;复杂的数据可以转化为易于理解的信息&#xff0c;从而帮助观众快速做出正确的判断。然而&#xff0c;糟糕的可视化可…

单片机IO灌入5V电压导致其他IO电压测量到大于供电电压问题

最近用GD32F103RCT6做项目&#xff0c;用了3个485收发器&#xff0c;都是直接接在单片机IO上的。 485收发器是5V供电的&#xff0c;这个时候就出现5V电平和3.3V电平兼容的问题了。 一开始只用了PA10、PC11这两个串口&#xff0c;他俩是兼容5V的&#xff0c;从手册可以看出IO最…

企业源代码也需要加密!十款好用的源代码加密软件排行榜

在当今竞争激烈的商业环境中&#xff0c;企业的源代码是其核心资产之一。为了保护这些宝贵的知识产权不被泄露&#xff0c;源代码加密成为了众多企业的重要举措。2024 年&#xff0c;市面上出现了众多功能强大的源代码加密软件。接下来&#xff0c;就让我们一同来探索十款备受好…

DockerCompose编排Nginx+Mysql并实现Nginx配置Mysql(TCP协议)负载均衡

场景 Nginx配置实例-负载均衡实例&#xff1a;平均访问多台服务器&#xff1a; Nginx配置实例-负载均衡实例&#xff1a;平均访问多台服务器_我想访问五个服务器的信息用nginx怎么做-CSDN博客 以上实现Nginx的http协议的负载均衡&#xff0c;如果使用Nginx实现TCP协议的负载…

Java的JVM中的概念之——新生代和老年代

JVM新生代和老年代是JVM中非常重要的概念&#xff0c;那么他们在JVM中扮演者什么样的角色和含义呢&#xff1f; 在Java虚拟机&#xff08;JVM&#xff09;的垃圾回收&#xff08;GC&#xff09;中&#xff0c;内存被分为不同的区域&#xff0c;其中两个主要区域是新生代&#…

PHP餐厅点餐系统小程序源码

&#x1f37d;️【餐厅点餐新纪元&#xff0c;点餐系统让用餐更便捷&#xff01;】&#x1f4f1; &#x1f50d; 一键浏览&#xff0c;菜单尽在掌握 &#x1f4f1; 走进餐厅&#xff0c;无需再担心找不到服务员或菜单被抢光&#xff01;餐厅点餐系统让你轻松扫描桌上的二维码…

HarmonyOS Developer之图片帧动画播放器

创建image-animator组件 在pages/index目录下的hml文件中创建一个image-animator组件&#xff0c;css文件中编写组件样式&#xff0c;js文件中引用图片。 设置image-animator组件属性 添加iteration&#xff08;播放次数&#xff09;、reverse&#xff08;播放顺序&#xf…

LVS原理及实例

目录 LVS原理 LVS概念 lvs集群的类型 lvs-nat 解释 传输过程 lvs-dr 解释 传输过程 特点 lvs-tun LVS&#xff08;Linux Virtual Server&#xff09;常见的调度算法 防火墙标记&#xff08;Firewall Marking&#xff09;结合轮询调度 实战案例 lvs的nat模式配置 …

使用Linux实现FTP云盘1

关于FTP服务器 FTP&#xff08;文件传输协议&#xff09;服务器是在互联网上提供文件存储和访问服务的计算机&#xff0c;它们依照FTP 协议提供服务。 FTP是File Transfer Protocol(文件传输协议)。 程序运行&#xff0c;服务端不断接收客户端指令&#xff0c;服务 端可同时处…

提取含有特定字符的行和列grep函数(含有替换)

目录 ①grep提取含有特定字符的列 ②grep提取含有特定字符的行 R语言进行字符的替换和删减gsub&#xff0c;substr函数R语言进行字符的替换和删减gsub&#xff0c;substr函数_r语言数据框字符替换-CSDN博客 ①grep提取含有特定字符的列 在一个dataframe中&#xff0c;需要提…

Element学习(axios异步加载数据、案例操作)(5)

1、这次学习的是上次还未完成好的恶element案例&#xff0c;对列表数据的异步加载&#xff0c;并渲染展示。 ——>axios来发送异步请求 &#xff08;1&#xff09; &#xff08;2&#xff09;在vue当中安装axios &#xff08;注意在当前的项目目录&#xff0c;并且安装完之后…