文章详情

短信预约-IT技能 免费直播动态提醒

请输入下面的图形验证码

提交验证

短信预约提醒成功

Python之动态规划

2023-08-30 14:56

关注

序言

最近在学习python语言,语言有通用性,此文记录复习动态规划并练习python语言。

动态规划(Dynamic Programming)

动态规划是运筹学的一个分支,是求解决策过程最优化的过程。20世纪50年代初,美国数学家贝尔曼(R.Bellman)等人在研究多阶段决策过程的优化问题时,提出了著名的最优化原理,从而创立了动态规划。动态规划的应用极其广泛,包括工程技术、经济、工业生产、军事以及自动化控制等领域,并在背包问题、生产经营问题、资金管理问题、资源分配问题、最短路径问题和复杂系统可靠性问题等中取得了显著的效果。

基本思想:将待求解问题分解成若干个子问题,先求解子问题,然后从这些子问题的解得到原问题的解。与分治法不同的是,适合于用动态规划求解的问题,经分解得到子问题往往不是互相独立的。

斐波那契数列(Fibonacci sequence)

斐波那契数列,又称黄金分割数列,因数学家莱昂纳多·斐波那契(Leonardo Fibonacci)以兔子繁殖为例子而引入,故又称“兔子数列”,其数值为:1、1、2、3、5、8、13、21、34……在数学上,这一数列以如下递推的方法定义:F(0)=1,F(1)=1, F(n)=F(n - 1)+F(n - 2)(n ≥ 2,n ∈ N*)。

先以斐波那契数列为例,了解动态规划。

def fibonacci(num):    if num == 0:        return 1    if num == 1:        return 1    return fibonacci(num - 1) + fibonacci(num - 2)if __name__ == "__main__":    print(fibonacci(10))

在这里插入图片描述
上述是以递归的方式实现的,然而递归方式存在以下几个缺点:

结果(n-1)项(n-2)项
f(n)f(n-1)f(n-2)
f(5)f(4)f(3)
f(4)f(3)f(2)
f(3)f(2)f(1)
f(2)f(1)f(0)

以上述表格为例,可以看到在求下一个递归结果时,计算了之前已经计算出来的结果,存在重复计算项。

如果采用动态规划的方式,那么可以节省计算,采用数组暂存之前已经计算出来的结果。如下,

def fibonacci_dp(num):    # 定义一个数组暂存dp结果,数组初始值为-1    dp = [-1] * (num + 1)    dp[0] = 1    dp[1] = 1    for i in range(2, num + 1):        dp[i] = dp[i - 1] + dp[i - 2]    return dp[num]if __name__ == "__main__":    print(fibonacci_dp(10))

在这里插入图片描述

不同路径

上面的斐波那契数列是一维数组,较为简单,下面以二维数组为例。

题目描述

一个机器人位于一个 m x n 网格的左上角 (起始点在下图中标记为 “Start” )。

机器人每次只能向下或者向右移动一步。机器人试图达到网格的右下角(在下图中标记为 “Finish” )。

问总共有多少条不同的路径?在这里插入图片描述

示例1: 输入:m = 3, n = 7 输出:28

示例2: 输入:m = 3, n = 2 输出:3
解释: 从左上角开始,总共有 3 条路径可以到达右下角。
1.向右 -> 向下 -> 向下
2.向下 -> 向下 -> 向右
3.向下 -> 向右 -> 向下
示例3:
输入:m = 7, n = 3 输出:28
示例4:
输入:m = 3, n = 3 输出:6
1 <= m, n <= 100
题目数据保证答案小于等于 2 * 10^9

python代码

class UniquePaths(object):    def uniquePaths(self, m: int, n: int) -> int:        """        :type m: int        :type n: int        :rtype: int        """        # 初始化一个二维数组        dp = [[0] * n for _ in range(m)]        for i in range(m):            dp[i][0] = 1        for j in range(n):            dp[0][j] = 1        for i in range(1, m):            for j in range(1, n):                dp[i][j] = dp[i - 1][j] + dp[i][j - 1]        return dp[m - 1][n - 1]if __name__ == "__main__":    demo = UniquePaths()    print(demo.uniquePaths(7, 3))

在这里插入图片描述

最小路径和

题目描述

给定一个包含非负整数的 m x n 网格 grid ,请找出一条从左上角到右下角的路径,使得路径上的数字总和为最小。

说明:一个机器人每次只能向下或者向右移动一步。

在这里插入图片描述

示例1:
输入:grid = [[1,3,1],[1,5,1],[4,2,1]]
输出:7
解释:因为路径 1→3→1→1→1 的总和最小。
示例2:
输入:grid = [[1,2,3],[4,5,6]]
输出:12
提示:
m == grid.length
n == grid[i].length
1 <= m, n <= 200
0 <= grid[i][j] <= 100

python代码

class MinPathSum(object):    def minPathSum(self, grid):        """        :type grid: List[List[int]]        :rtype: int        """        row = len(grid)        column = len(grid[0])        # 定义dp[i][j]为到(i,j)处的最小路径和        dp = [[0] * column for _ in range(row)]        dp[0][0] = grid[0][0]        # 第0行j列        for j in range(1, column):            dp[0][j] = dp[0][j - 1] + grid[0][j]        # 第i行0列        for i in range(1, row):            dp[i][0] = dp[i - 1][0] + grid[i][0]        # 非第0行或第0列        for i in range(1, row):            for j in range(1, column):                dp[i][j] = min(dp[i - 1][j], dp[i][j - 1]) + grid[i][j]        return dp[row - 1][column - 1]if __name__ == "__main__":    demo = MinPathSum()    grid = [[1, 3, 1], [1, 5, 1], [4, 2, 1]]    print(demo.minPathSum(grid))

在这里插入图片描述

零钱兑换

题目描述

给你一个整数数组 coins ,表示不同面额的硬币;以及一个整数 amount ,表示总金额。

计算并返回可以凑成总金额所需的 最少的硬币个数 。如果没有任何一种硬币组合能组成总金额,返回 -1 。

你可以认为每种硬币的数量是无限的。

示例 1:
输入:coins = [1, 2, 5], amount = 11
输出:3
解释:11 = 5 + 5 + 1
示例 2:
输入:coins = [2], amount = 3
输出:-1
示例 3:
输入:coins = [1], amount = 0
输出:0
提示:
1 <= coins.length <= 12
1 <= coins[i] <= 231 - 1
0 <= amount <= 104

python代码

class CoinChange(object):    def coinChange(self, coins: list[int], amount: int) -> int:        """        :type coins: List[int]        :type amount: int        :rtype: int        """        # 状态转移方程dp(i) = min(dp(i-Cj)) + 1,Cj为货币面值        """        i<0  忽略        i==0 dp[0] = 0        i==1 dp[1] = min(dp[1-1], dp[1-2], dp[1-5]) + 1 = 1        i==2 dp[2] = min(dp[2-1], dp[2-2], dp[2-5]) + 1 = 1        i==3 dp[3] = min(dp[3-1], dp[3-2], dp[3-5]) + 1 = 2        i==4 dp[4] = min(dp[4-1], dp[4-2], dp[4-5]) + 1 = 2        ... ...        """        dp = [0] * (amount + 1)        dp[0] = 0        for i in range(1, amount + 1):            mini = int(1e9)            for coin in coins:                if i >= coin:                    res = dp[i - coin]                    if 0 <= res < mini:                        mini = res            dp[i] = mini + 1 if mini < int(1e9) else -1        if amount < 1:            return 0        return dp[amount]if __name__ == "__main__":    demo = CoinChange()    coins = [1, 2, 5]    amount = 11    print(demo.coinChange(coins, amount))

在这里插入图片描述

来源地址:https://blog.csdn.net/zkkzpp258/article/details/132549864

阅读原文内容投诉

免责声明:

① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的,并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据,供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。

② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341

软考中级精品资料免费领

  • 历年真题答案解析
  • 备考技巧名师总结
  • 高频考点精准押题
  • 2024年上半年信息系统项目管理师第二批次真题及答案解析(完整版)

    难度     807人已做
    查看
  • 【考后总结】2024年5月26日信息系统项目管理师第2批次考情分析

    难度     351人已做
    查看
  • 【考后总结】2024年5月25日信息系统项目管理师第1批次考情分析

    难度     314人已做
    查看
  • 2024年上半年软考高项第一、二批次真题考点汇总(完整版)

    难度     433人已做
    查看
  • 2024年上半年系统架构设计师考试综合知识真题

    难度     221人已做
    查看

相关文章

发现更多好内容

猜你喜欢

AI推送时光机
位置:首页-资讯-后端开发
咦!没有更多了?去看看其它编程学习网 内容吧
首页课程
资料下载
问答资讯