把多阶段过程转化为一系列单阶段问题,利用各阶段之间的关系,逐个求解。
动态规划过程:每次决策依赖于当前状态,又随即引起状态的转移。一个决策序列就是在变化的状态中产生出来的,所以,这种多阶段最优化决策解决问题的过程就称为动态规划。
假设问题是由交叠的子问题所构成,我们就能够用动态规划技术来解决它。一般来说,这种子问题出自对给定问题求解的递推关系中,这个递推关系包括了同样问题的更小子问题的解。动态规划法建议,与其对交叠子问题一次重新的求解,不如把每一个较小子问题仅仅求解一次并把结果记录在表中(动态规划也是空间换时间的)。这样就能够从表中得到原始问题的解。
基本思想与分治法类似,也是将待求解的问题分解为若干个子问题(阶段),按顺序求解子阶段,前一子问题的解,为后一子问题的求解提供了实用的信息。
在求解任一子问题时,列出各种可能的局部解,通过决策保留那些有可能达到最优的局部解,丢弃其它局部解。依次解决各子问题,最后一个子问题就是初始问题的解。
因为动态规划解决的问题多数有重叠子问题这个特点。为降低反复计算。对每个子问题仅仅解一次,将其不同阶段的不同状态保存在一个二维数组中。
与分治法最大的区别是:适合于用动态规划法求解的问题,经分解后得到的子问题往往不是互相独立的(即下一个子阶段的求解是建立在上一个子阶段的解的基础上,进行进一步的求解)。
能采用动态规划求解的问题的一般要具有3个性质:
(1)最优化原理:假设问题的最优解所包括的子问题的解也是最优的,就称该问题具有最优子结构,即满足最优化原理。
(2)无后效性:即某阶段状态一旦确定。就不受这个状态以后决策的影响。也就是说,某状态以后的过程不会影响曾经的状态。仅仅与当前状态有关;
(3)有重叠子问题:即子问题之间是不独立的,一个子问题在下一阶段决策中可能被多次使用到(该性质并非动态规划适用的必要条件,可是假设没有这条性质,动态规划算法同其它算法相比就不具备优势)。
动态规划所处理的问题是一个多阶段决策问题,一般由初始状态开始,通过对中间阶段决策的选择,达到结束状态。这些决策形成了一个决策序列,同时确定了完成整个过程的一条活动路线(通常是求最优的活动路线)。
(1)划分阶段:按照问题的时间或空间特征,把问题分为若干个阶段。在划分阶段时,注意划分后的阶段一定要是有序的或者是可排序的,否则问题就无法求解
(2)确定状态和状态变量:将问题发展到各个阶段时所处于的各种客观情况用不同的状态表示出来。当然,状态的选择要满足无后效性。
(3)确定决策并写出状态转移方程:状态转移就是根据上一阶段的状态和决策来导出本阶段的状态。如果确定了决策,状态转移方程也就可写出。但事实上常常是反过来做,根据相邻两个阶段的状态之间的关系来确定决策方法和状态转移方程。
(4)寻找边界条件:给出的状态转移方程是一个递推式,需要一个递推的终止条件或边界条件。
只要解决问题的阶段、状态和状态转移决策确定了,就可以写出状态转移方程(包括边界条件)。
1、找零钱问题(求最值动态规划)
有面值为1元、2元和7元的硬币若干枚,如何用最少的硬币凑够27元?
class Solution:
def coinChange(self, coins: List[int], amount: int) -> int:
# 开数组 27 个位置的数组
f = [float("inf")] * (amount+1)
# 初始条件:0元0种方法
f[0] = 0
n = len(coins)
#f[1],f[2]....f[27]
#每个数组空间装当前零钱最优的拼凑结果
for i in range(1,amount+1):
#对每种零钱进行遍历
#转移方程f[X] = min{f[X-2]+1 , f[X-5]+1 , f[x-7]+1}
for j in range(0,n):
if(i>=coins[j]):
f[i] = min(f[i-coins[j]]+1,f[i])
#如果凑不出目标数,返回-1
if f[-1] == float("inf"):
return -1
#返回次数
else:return f[-1]
2、走方格问题(计数型动态规划)
给定m行和n列网格,有一个机器人从左上角(0,0)出发,每一步可以向下或向右走一步,问多少种不同的方式走到右下角
1.确定状态
2.转移方程:对任意一个格子:f[i][j] = f[i-1][j]+f[i][j-1]
3.初始条件和边界情况:
4.计算顺序
求f[m-1][n-1](时间复杂度:O(MN) 空间复杂度(数组大小):O(MN))
5.代码实现
import numpy as np
def uniquePaths(m,n): #m:行 ,n:列 5 4
#如果用([0]*m)*n会涉及到拷贝问题
f = np.zeros((m,n),dtype=np.int)
for j in range(n):
for i in range(m):
#边界为0的情况下均为1步
if(i==0 or j==0):
f[i][j] = 1
#转移方程
else:
f[i][j] = f[i-1][j]+f[i][j-1]
return f[i][j]
uniquePaths(5,4)
3、青蛙过桥(存在型动态规划)
有n块石头分别在x轴的0,1,2,… ,n-1位置,一只青蛙在石头0处,如果青蛙在第i块石头上,它最多往右跳ai个石头。
问青蛙是否能够调到石头n-1处
1.确定状态
2.转移方程:
设状态为dp[j] = 青蛙能不能跳到石头j,有:
3.初始条件和边界情况
初始条件:dp[0] = True,青蛙一开始就在石头0,无边界情况。
4.计算顺序
由初始条件dp[0] = True,
计算dp[1],dp[2],…,dp[n-1],结果:dp[n-1]
时间复杂度:O(n^2)(用贪心算法复杂度为O(n)),空间复杂度:O(n)
5.代码实现