Unit2_1:动态规划DP

news2024/12/23 13:00:57

文章目录

  • 一、介绍
  • 二、0-1背包问题
    • 问题描述
    • 分析
    • 伪代码
    • 时间复杂度
  • 三、钢条切割问题
    • 问题描述
    • 分析
    • 伪代码
    • 过程
  • 四、矩阵链乘法
    • 背景
    • 性质
    • 分析
    • 案例
    • 伪代码

一、介绍

动态规划类似于分治法,它们都将一个问题划分为更小的子问题
最优子结构:问题的最优解包含子问题的最优解。DP适用的原因就在这
当子问题重叠时,即它们共享公共子问题时,可减小时间复杂度
DP通常用于优化问题,有许多解决方案的问题,我们想找到最好的一个
DP问题的求解思路一般就是
   先描述最优解的结构
   递归地定义最优解的值
   计算最优解的值(通常是自下而上)
   根据计算出的信息构造最优解(如果需要)

二、0-1背包问题

问题描述

n 个商品 , v i 表示第 i 个物品的价值 , w i 表示第 i 个物品的重量 一个能装入 W 重的背包 , 使用背包装下价值最多的物品 限制条件 : 我们不能取物品的一部分,我们取整个物品,或者什么都不取。 ( 这就是为什么它被称为 0 − 1 背包。 ) n个商品,v_i表示第i个物品的价值,w_i表示第i个物品的重量\\ 一个能装入W重的背包,使用背包装下价值最多的物品\\ 限制条件:我们不能取物品的一部分,我们取整个物品,或者什么都不取。\\(这就是为什么它被称为0-1背包。) n个商品,vi表示第i个物品的价值,wi表示第i个物品的重量一个能装入W重的背包,使用背包装下价值最多的物品限制条件:我们不能取物品的一部分,我们取整个物品,或者什么都不取。(这就是为什么它被称为01背包。)

分析

V [ i , w ] 表示重量为 w 背包 , 在前 i 种商品选择的最大价值 V[i,w]表示重量为w背包,在前i种商品选择的最大价值 V[i,w]表示重量为w背包,在前i种商品选择的最大价值
对于第i个物品,我们要么选取它,要么不选择,因此最大价值转移方程为

V [ i , w ] = m a x ( V [ i − 1 , w ] , v i + V [ i − 1 , w − w i ] ) V[i,w]=max(V[i-1,w],v_i+V[i-1,w-w_i]) V[i,w]=max(V[i1,w],vi+V[i1,wwi])

若使用递归重复计算很多值,时间复杂度为 T ( n ) = O ( 2 W ) T(n)=O(2^W) T(n)=O(2W),因此要重复利用最优子结构的性质.
在这里插入图片描述
初始化: V [ 0 , w ] = 0 V[0,w]=0 V[0,w]=0 f o r for for 0 ≤ w ≤ W 0 \leq w \leq W 0wW 此时没有商品,自然没有价值
接下来按顺序填表:
在这里插入图片描述
在这里插入图片描述

伪代码

Knapsack(v,w,n,W)
for w=0 to W do
	V[0,w]=0;
end
for i=1 to n do
	for w=0 to W
		if w[i]<W then
			V[i,w] = max{V[i-1,w],v[i]+V[i-1,w-w[i]]}
		else
			V[i,w] = V[i-1,w]
		end
	end
end
return V[n,W]

若是想要记录最优解的路径,需要维护一个 k e e p [ i ] [ w ] keep[i][w] keep[i][w],如果选择i作为 V [ i , w ] V[i,w] V[i,w],则 k e e p [ i ] [ w ] = 1 keep[i][w]=1 keep[i][w]=1
路径只需要
i f k e e p [ n , w ] = 1. 则选择 n 且继续从 k e p p [ n − 1 ] [ w − w n ] 开始 if keep[n,w]=1.则选择n且继续从kepp[n-1][w-w_n]开始 ifkeep[n,w]=1.则选择n且继续从kepp[n1][wwn]开始
i f k e e p [ n , w ] = 0. 则不选择 n 且继续从 k e p p [ n − 1 ] [ w ] 开始 if keep[n,w]=0.则不选择n且继续从kepp[n-1][w]开始 ifkeep[n,w]=0.则不选择n且继续从kepp[n1][w]开始
因此路径输出代码:

K ← W
for i ← n to 1 do
	if keep[i][K] is equal to 1 then
		Output i
		K ← K-w[i]
	end
end

在这里插入图片描述

时间复杂度

两层循环,时间复杂度 T ( n ) = O ( n W ) T(n)=O(nW) T(n)=O(nW)
本质上,DP问题是用空间换时间,将结果放在空间中而不用下次花费时间再去计算

三、钢条切割问题

问题描述

给定一个长度为 n 的棒材和一个价格表,其中 p i 为长度为 i 的棒材的价格 确定最大的收入 r n , 以及切割钢条的方案 给定一个长度为n的棒材和一个价格表,其中pi为长度为i的棒材的价格\\ 确定最大的收入r_n,以及切割钢条的方案 给定一个长度为n的棒材和一个价格表,其中pi为长度为i的棒材的价格确定最大的收入rn,以及切割钢条的方案

分析

此题暴力解法,即遍历长度.每个点有两种选择,切 o r or or不切,判断哪种选择最合适即可,时间复杂度 T ( n ) = O ( 2 n ) T(n)=O(2^n) T(n)=O(2n).
考虑到最优子结构,可利用较短的杆最优收益来确定较长的,此题的状态转移如下:
r n = m a x ( p n , r 1 + r n − 1 , r 2 + r n − 2 , . . . . . , r n − 1 + r 1 ) r_n=max(p_n,r_1+r_{n-1},r_2+r_{n-2},.....,r_{n-1}+r_1) rn=max(pn,r1+rn1,r2+rn2,.....,rn1+r1)
简化定义:
r n = m a x ( p i + r n − i ) r_n=max(p_i+r_{n-i}) rn=max(pi+rni)      1 ≤ i ≤ n 1 \leq i \leq n 1in

伪代码

r[0] ← 0
for j ← 0 to n do
	q  ← -∞
	for i  ← 1 to j do
		q  ← max(q,p[i]+r[j-i])
	end
	r[j]  ← q if j != 0
end
return r[n]

这种做法时间复杂度 T ( n ) = O ( n 2 ) T(n)=O(n^2) T(n)=O(n2)
若是需要保存切割的方案,则需要维护一个 s [ n ] s[n] s[n]数组. s [ n ] 保存前一次切割的长度 s[n]保存前一次切割的长度 s[n]保存前一次切割的长度:

r[0] ← 0
for j ← 0 to n do
	q  ← -∞
	for i  ← 1 to j do
		if q < p[i]+r[j-i] then
			q ← p[i]+r[j-i]
			s[j] ← i
		end
	end
	r[j]  ← q if j != 0
end
while n>0 do
	Output s[n]
	n ← n-s[n]
end

过程

在这里插入图片描述

四、矩阵链乘法

背景

p × q 矩阵 a 和 q × r 矩阵 B 的乘积 C = A B 是 p × r 矩阵 p × q矩阵a和q × r矩阵B的乘积C = AB是p × r矩阵 p×q矩阵aq×r矩阵B的乘积C=ABp×r矩阵

c [ i ] [ j ] = ∑ k = 1 q a [ i ] [ k ] b [ k ] [ j ] c[i][j]=\sum_{k=1}^{q}a[i][k]b[k][j] c[i][j]=k=1qa[i][k]b[k][j] f o r for for 1 ≤ i ≤ p 1 \leq i \leq p 1ip a n d and and 1 ≤ j ≤ r 1 \leq j \leq r 1jr
时间复杂度:注意 C C C p r pr pr个条目,每个条目需要 O ( q ) O(q) O(q)时间来计算,所以整个过程需要 O ( p q r ) O(pqr) O(pqr)时间

性质

矩阵乘法有结合律, A 1 A 2 A 3 = ( A 1 A 2 ) A 3 = A 1 ( A 2 A 3 ) A_1A_2A_3=(A_1A_2)A_3=A_1(A_2A_3) A1A2A3=(A1A2)A3=A1(A2A3)
因此当计算 A B C ABC ABC时,有两种选择, ( A B ) C = A ( B C ) (AB)C=A(BC) (AB)C=A(BC)
m u l t [ ( A B ) C ] = p q r + p r s mult[(AB)C]=pqr+prs mult[(AB)C]=pqr+prs
m u l t [ A ( B C ) ] = q r s + p q s mult[A(BC)]=qrs+pqs mult[A(BC)]=qrs+pqs
每种选择的时间复杂度不一样
因此矩阵链乘法需要解决的就是怎么样结合能使得计算的量最小

分析

令 A i . . j = A i A i + 1 . . . A j 令A_{i..j}=A_iA_{i+1}...A{j} Ai..j=AiAi+1...Aj,显然, A i . . j A_{i..j} Ai..j P i − 1 × P j P_{i-1}×P_j Pi1×Pj的矩阵
A i . . j A_{i..j} Ai..j可以表示为 A i . . j = ( A i . . . A k ) ( A k + 1 . . . A j ) = A i . . k A k + 1.. j A_{i..j}=(A_i...A_k)(A_{k+1}...A{j})=A_{i..k}A_{k+1..j} Ai..j=(Ai...Ak)(Ak+1...Aj)=Ai..kAk+1..j

1 ≤ i ≤ j ≤ n 1 \leq i \leq j \leq n 1ijn时,令 m [ i , j ] m[i,j] m[i,j]表示计算 A i . . j A_{i..j} Ai..j所需的最小乘法次数。最优成本可以用下面的递归定义来描述

m ( i , j ) = { 0 i f   i = j m i n i ≤ k ≤ j ( m ( i , k ) + m ( k + 1 , j ) + p i − 1 p k p j ) i f   n = 1 m(i,j)=\left\{ \begin{array}{ll} 0 & if \space i=j \\ min_{i \leq k \leq j} (m(i,k)+m(k+1,j)+p_{i-1}p_kp_j)& if \space n=1 \nonumber \end{array} \right. m(i,j)={0minikj(m(i,k)+m(k+1,j)+pi1pkpj)if i=jif n=1

注意计算并保存 m [ i , j ] m[i, j] m[i,j]的顺序是,当计算 m [ i , j ] m[i, j] m[i,j]时, m [ i , k ] m[i, k] m[i,k] m [ k + 1 , j ] m[k + 1, j] m[k+1,j]的值已经可用,因此按矩阵链长度的递增顺序计算它们:

m [ 1 , 2 ] , m [ 2 , 3 ] , m [ 3 , 4 ] , … , m [ n − 3 , n − 2 ] , m [ n − 2 , n − 1 ] , m [ n − 1 , n ] m[1,2], m[2,3], m[3,4],…, m[n-3,n-2], m[n-2,n-1], m[n-1,n] m[1,2],m[2,3],m[3,4],,m[n3,n2],m[n2,n1],m[n1,n]
m [ 1 , 3 ] , m [ 2 , 4 ] , m [ 3 , 5 ] , … , m [ n − 3 , n − 1 ] , m [ n − 2 , n ] m[1,3], m[2,4], m[3,5],…, m[n-3,n-1], m[n-2,n] m[1,3],m[2,4],m[3,5],,m[n3,n1],m[n2,n]
m [ 1 , 4 ] , m [ 2 , 5 ] , m [ 3 , 6 ] , … , m [ n − 3 , n ] m[1,4], m[2,5], m[3,6],…, m[n-3,n] m[1,4],m[2,5],m[3,6],,m[n3,n]
… …
m [ 1 , n − 1 ] , m [ 2 , n ] m[1,n-1], m[2,n] m[1,n1],m[2,n]
m [ 1 , n ] m[1,n] m[1,n]

若需要记录分隔括号路径,需要维护一个二维数组 s [ 1.. n , 1.. n ] s[1..n, 1..n] s[1..n,1..n],里面存储 A i . . j A_{i..j} Ai..j的最优分隔k
s [ 1 , n ] s[1,n] s[1,n]           ( A 1 . . A s [ 1 , n ] ) ( A s [ 1 , n ] + 1 . . . A n ) (A_1..A_{s[1,n]})(A_{s[1,n]+1}...A{n}) (A1..As[1,n])(As[1,n]+1...An)
s [ 1 , s [ 1. n ] ] s[1,s[1.n]] s[1,s[1.n]]           ( A 1 . . A s [ 1 , s [ 1. n ] ] ) ( A s [ 1 , s [ 1. n ] ] + 1 . . . A n ) (A_1..A_{s[1,s[1.n]]})(A_{s[1,s[1.n]]+1}...A{n}) (A1..As[1,s[1.n]])(As[1,s[1.n]]+1...An)

案例

在这里插入图片描述

伪代码

MatrixChain(p,n)
for i ← 1 to n do
	m[i,i] ← 0;
end
for l ← 2 to n do
	for i ← 1 to n-l+1 do
		j ← i+l-1;
		m[i,j] ← ∞
		for k ← i to j-1 do
			q ← m[i,k]+m[k+1,j]+p[i-1]*p[k]*p[j]
			if q<m[i,j] then
				m[i,j] ← q
				s[i,j] ← k
			end
		end
	end
end
return m[1,n] and s

三层循环,时间复杂度 T ( n ) = O ( n 3 ) T(n)=O(n^3) T(n)=O(n3)

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

一种单总线串口通信的调试方法

单总线的优点&#xff1a; 节省IO口&#xff0c;发送时可以将单片机的RXD设置为普通IO进行软件模拟发送&#xff0c;发送完设置为串口接收。避免通信干扰&#xff0c;由于是通过IO口对三极管/MOS管进行拉高拉低&#xff0c;外部信号不易对IO口进行干扰&#xff0c;EMI&#xf…

代码随想录算法训练营第四十六天 | LeetCode 139. 单词拆分、多重背包、背包总结

代码随想录算法训练营第四十六天 | LeetCode 139. 单词拆分、多重背包、背包总结 文章链接&#xff1a;单词拆分 多重背包 背包总结 视频链接&#xff1a;单词拆分 1. LeetCode 139. 单词拆分 1.1 思路 本题的那些单词就是物品&#xff0c;字符串就是背包&#xff0c;问用这些…

大数据学习之一文学会Spark【Spark知识点总结】

文章目录 什么是SparkSpark的特点Spark vs HadoopSparkHadoopSpark集群安装部署Spark集群安装部署StandaloneON YARN Spark的工作原理什么是RDDRDD的特点Spark架构相关进程Spark架构原理 Spark实战&#xff1a;单词统计Scala代码开发java代码开发任务提交 Transformation与Acti…

echarts 类目轴设置xAxis.interval不起效果

情景&#xff1a;在xAxis.typecategory&#xff1b;设置xAxis.interval不起效果&#xff1b; 解决方案&#xff1a;通过设置xAxis.axisLabel.interval属性

【H616_语言小美_控制安卓刷抖音项目 orangePi zero2 (已开源) 】.md uptada:23/11/07

文章目录 H616_语言小美_控制安卓刷抖音项目小美效果展示H616 ubuntu系统 安装adb智能公元 SU-03T 离线语音模组 固件制作配合串口实现 小美_控制安卓刷抖音 H616_语言小美_控制安卓刷抖音项目 注意&#xff1a;orangePi zero2 H616 安装系统为ubuntu 小美效果展示 语言小美 …

企业级私有化部署数字办公工具集合,解决企业全方位需求,快收藏

随着信息技术的日新月异&#xff0c;企业对IT建设的依赖程度越来越高。为了满足这一需求&#xff0c;软开企服凭借积累了多年的项目经验&#xff0c;推出了一款名为JVS的私有化部署企业信息化底座。JVS涵盖了多个产品领域&#xff0c;如协同办公、低代码开发、数据分析、风控决…

KaiOS APN配置文件apn.json调试验证方法(无需项目全编)

1、KaiOS 的应用就类似web应用&#xff0c;结合文件夹路径webapp字面意思理解。 2、KaiOS APN配置文件源代码在apn.json&#xff0c; &#xff08;1&#xff09;apn.json可以自定义路径&#xff0c;通过配置脚本实现拷贝APN在编译时动态选择路径在机器中生效。 &#xff08;…

linux复习笔记03(小滴课堂)

find命令&#xff1a; d查找目录&#xff1a; 按照文件权限查找&#xff1a; 600全部权限&#xff1a; -user根据所属主&#xff1a; 上面的例子是找出文件并打印有多少行。 我们也可以把我们查询到的结果复制到其它文件位置中去&#xff1a; 复制成功。 -mtime根据修改时间…

DOCTYPE是什么,有何作用、 使用方式、渲染模式、严格模式和怪异模式的区别?

前言 持续学习总结输出中&#xff0c;今天分享的是DOCTYPE是什么&#xff0c;有何作用、 使用方式、渲染模式、严格模式和怪异模式的区别。 DOCTYPE是什么&#xff0c;有何作用&#xff1f; DOCTYPE是HTML5的文档声明&#xff0c;通过它可以告诉浏览器&#xff0c;使用那个H…

MySQL的高阶语句

数据库的权限一般很小&#xff0c;工作中使用最多的场景就是查 排序、分组、子查询、视图、多表连接查询&#xff08;左连接、右连接、内连接&#xff09; create TABLE info ( id int(4) primary key, NAME varchar(5) not null, score decimal(5,2), address varchar(20)…

Django实战项目-学习任务系统-发送短信通知

接着上期代码内容&#xff0c;继续完善优化系统功能。 本次增加发送短信通知功能&#xff0c;学习任务系统发布的任务&#xff0c;为了更加及时通知到学生用户&#xff0c;再原有发送邮件通知基础上&#xff0c;再加上手机短信通知功能。 第一步&#xff1a;开通短信通知服务…

JSP通用材料收集归档系统eclipse定制开发mysql数据库BS模式java编程jdbc

一、源码特点 JSP 通用材料收集归档系统是一套完善的web设计系统&#xff0c;对理解JSP java编程开发语言有帮助&#xff0c;系统具有完整的源代码和数据库&#xff0c;系统主要采用B/S模式开发。开发环境为TOMCAT7.0,eclipse开发&#xff0c;数据库为Mysql5.0&#xff0c…

计算机服务器中了eking勒索病毒怎么办,勒索病毒解密,文件数据恢复

随着科技的不断发展&#xff0c;网络技术也得到了不断更新&#xff0c;极大地方便了人们的生产与生活&#xff0c;但随之而来的网络安全威胁也不断增加&#xff0c;近期&#xff0c;网络上的eking勒索病毒开始流行&#xff0c;它是phobos勒索家族中的一种病毒&#xff0c;具有较…

红黑树的性质和实现

红黑树 由于AVL树为了保持平衡需要经常旋转&#xff0c;开销是很大的。因此&#xff0c;红黑树比起AVL树放宽了平衡的限制。 概念 红黑树&#xff0c;在每个节点上增加一个存储位表示节点的颜色&#xff0c;可以是RED或BLACK。 通过对任何一条从根到叶子的路径上各个结点着…

EGpKa

利用数据增强的可解释图神经网络预测碳氢化合物的 pKa&#xff0c;原文 Explainable Graph Neural Networks with Data Augmentation for Predicting pKa of C–H Acids&#xff0c;代码在 DATA AND SOFTWARE AVAILABILITY。模型框架如下&#xff1a;

如何使用Linux DataEase数据可视化分析工具结合内网穿透实现远程办公

文章目录 前言1. 安装DataEase2. 本地访问测试3. 安装 cpolar内网穿透软件4. 配置DataEase公网访问地址5. 公网远程访问Data Ease6. 固定Data Ease公网地址 前言 DataEase 是开源的数据可视化分析工具&#xff0c;帮助用户快速分析数据并洞察业务趋势&#xff0c;从而实现业务…

git从入门到会用

git从入门到大牛 什么是git理解版本控制版本控制软件 远程仓库git安装windows 安装git 客户端安装下载客户端安装界面介绍 客户端操作仓库管理创建本地仓库切换本地仓库删除本地仓库 文件操作创建文件修改文件删除文件新增/修改/删除 标志 分支管理创建分支切换分支合并分支 标…

【嵌入式 C 常用算法 3 -- 非线性存储结构 二叉树】

文章目录 树概念介绍树相关术语结点子树结点的度结点的层次有序树和无序树森林空树&#xff08;简单了解即可&#xff09; 二叉树二叉树性质满二叉树完全二叉树 二叉树的深度遍历前序遍历中序遍历后续遍历 二叉树的层次遍历二叉树的顺序存储结构二叉树的顺序存储结构C代码实现 …

如何解决DC电源模块的电源噪声问题

BOSHIDA 如何解决DC电源模块的电源噪声问题 在电子设备的设计和制作过程中&#xff0c;电源噪声是一个非常重要的考虑因素。DC电源模块的电源噪声问题是电子设备中普遍存在的问题之一。它不仅会影响设备的性能&#xff0c;还会对设备的寿命和稳定性产生负面影响。因此&#xf…

第二证券:消费电子概念活跃,博硕科技“20cm”涨停,天龙股份斩获10连板

消费电子概念7日盘中再度拉升&#xff0c;到发稿&#xff0c;博硕科技“20cm”涨停&#xff0c;光大同创、波长光电涨超10%&#xff0c;易德龙、向阳科技、得润电子、天龙股份、同兴达等涨停。 博硕科技强势涨停&#xff0c;公司昨日在接受安排调研时表明&#xff0c;公司从上…