关于dp:聊聊动态规划

34次阅读

共计 5294 个字符,预计需要花费 14 分钟才能阅读完成。

文章首发于公众号👉 青梅主码,欢送关注

简介

动静布局(Dynamic programming,简称 DP)是美国数学家 Richard Bellman 在钻研 决策过程和控制系统实践 时创立的新办法。它在数学上属于运筹学的一个分支,在数学、管理科学、计算机科学、经济学和生物信息学中均有利用,外围是 通过把原问题合成为绝对简略的子问题的形式 来求解简单问题,次要利用是求解决策过程最优的数学方法。

简略来讲,动静布局是一种算法思维,其外围是把问题合成为子问题,通过求解子问题进而解决以后问题。理论中,并非所有问题都能够通过动静布局来求解,通过动静布局解决的问题,对问题和其合成对子问题都有肯定场景要求的,动静布局实用于有 重叠子问题 最优子结构 性质的问题,这类问题应用动静布局所耗时间往往比奢侈解法更少。

如下图,咱们对动静布局问题从 利用场景 题解思路 常见示例题目 三个方面来开展。

动静布局 Wiki

利用场景

重叠子问题

重叠子问题关注的是一个问题被合成为多个子问题时,子问题和子问题之间的关联。子问题会被反复计算,特地是咱们应用递归自上向下对问题进行求解时,那动静布局是如何解决和优化这些子问题呢?咱们以最简略的斐波那契数列为例来阐明。

斐波那契数列:
写一个函数,输出 n,求斐波那契(Fibonacci)数列的第 n 项(即 F(N))。斐波那契数列的定义如下:

F(0) = 0,F(1) = 1

F(N) = F(N – 1) + F(N – 2), 其中 N > 1.

  • 示例 1:

    • 输出:n = 2
    • 输入:1
  • 示例 2:

    • 输出:n = 5
    • 输入:5

自上而下

通过咱们能够通过自上而下的递归办法来解决该问题,代码如下

function fib(n: number): number {// n = 0, f(n) = 0;
    // n = 1, f(n) = 1;
    if (n < 2) return n

    return fib(n - 1) + fib(n - 2)
};

这样通过递归来解决,代码看起来很简洁,然而仔细分析下,当 n 取值较大时,这其中的反复计算可就太多了,如下,能够简略合成下

能够看到,f(n-2)f(n-3) 等均进行了反复计算,可想而知,当数据量较大时,反复计算会更多,所以,递归实现代码很简略,计算却尤为消耗性能,能够尝试计算 n 取值为 100, 看下计算工夫。
那么,应用动静布局能够怎么求解这个问题呢

自下而上

咱们能够从最简略的、规模最小的 f(0)f(1) 开始计算,直到求解 f(n),所以动静布局个别是不采纳递归的形式,而是经由循环来实现计算的,看下代码

function fib(n: number): number {let dp = new Array(n).fill(0)

    // n = 0, f(n) = 0;
    // n = 1, f(n) = 1;
    dp[0] = 0, dp[1] = 1

    for (let i = 2; i <= n; i++){dp[i] = dp[i-1]+dp[i-2]
    }

    return dp[n]
};

能够看出,这种自下而上的解决办法,每次计算时都会用到上次的计算结果,从而能够躲避大量的子问题反复计算

数组最大索引越界,针对本题,可能会超过数组的索引边界,所以须要额定解决下

function fib(n: number): number {let dp = new Array(n).fill(0), mod = 1000000007

    dp[0] = 0, dp[1] = 1

    for (let i = 2; i <= n; i++) {dp[i] = dp[i - 1] % mod + dp[i - 2] % mod
    }

    return dp[n] % mod
};

最优子结构

最优子结构关注的是当一个问题被合成为多个子问题时,子问题和原问题之间的关联。即问题的最优解所蕴含的子问题的解也是最优的,也就是说咱们能够通过求解每个子问题的最优解来最初决定该问题的最优解

无后效性

子问题的解一旦确定,就不受其后问题的影响。具体来讲,就是子问题的解一旦确定,就不再扭转

题解思路

动静布局问题的数学表达式能够简略的形象为如下: dp[y] = f(dp[x]),x<y,个别常见解题思路,次要蕴含以下四步,其中最外围的是确定 转移方程

  • 状态定义:定义动静布局过程中波及的变量状态
  • 转移方程:演绎和形象出子问题对应的数学方程
  • 初始值:初始化时变量的取值
  • 返回值:返回求解后所须要的值

咱们以如下题目为动手来理解动静布局的常见解题思路

常见题目

动静布局常见的几种类型题目,⭐ 代表题目难度

最大子序和 ⭐

👉 Leetcode 链接

给你一个整数数组 nums,请你找出一个具备最大和的间断子数组(子数组起码蕴含一个元素),返回其最大和。

子数组 是数组中的一个间断局部。

用例

  • 示例 1:

    • 输出: nums = [-2,1,-3,4,-1,2,1,-5,4]
    • 输入: 6
    • 解释: 间断子数组 [4,-1,2,1] 的和最大,为 6。
  • 示例 2:

    • 输出: nums = [1]
    • 输入: 1

思路剖析

  • 状态定义,本题外围是剖析什么是 最大和的间断子数组,这里咱们能够用 f(n) 示意第 n 个数结尾的最大和的间断子数组
  • 转移方程剖析,不言而喻,咱们的问题就是求 [0, n] 这个区间外面 f(n) 的最大值,利用动静布局的划分子问题的办法,对 f(n) 来讲,其最大值能够通过比拟 nums[n]nums[n] + f(n-1),取两者中较大值即可,从而,咱们能够得出其转移方程为:
    f(n) = max{(f(n-1) + nums[n], nums[n])}
  • 初始值剖析,咱们能够应用数组中第一个数来临时示意最大值 max,随后,通过循环数组中所有项,代入转移方程,不断更新 max 取值即可
  • 返回值即为 max

题解

function maxSubArray(nums: number[]): number {let max = nums[0], dp = [nums[0]], len = nums.length

    for (let n = 1; n < len; n++) {dp[n] = Math.max(0, dp[n - 1]) + nums[n]
        max = Math.max(max, dp[n])
    }

    return max
};

股票的最大利润 ⭐

👉 Leetcode 链接

假如把某股票的价格依照工夫先后顺序存储在数组中,请问交易该股票一次可能取得的最大利润是多少?

用例

  • 示例 1:

    • 输出: [7,1,5,3,6,4]
    • 输入: 5
    • 解释: 在第 2 天(股票价格 = 1)的时候买入,在第 5 天(股票价格 = 6)的时候卖出,最大利润 = 6-1 = 5。
      留神利润不能是 7-1 = 6, 因为卖出价格须要大于买入价格。
  • 示例 2:

    • 输出: [7,6,4,3,1]
    • 输入: 0
    • 解释: 在这种状况下, 没有交易实现, 所以最大利润为 0。

思路剖析

和上题相似,咱们首先思考🤔,如何获取 最大利润 ?依照惯例思路,天然是 最低点买入,最高点卖出

  • 状态定义,咱们能够用 f(n) 示意前 n 天的最大利润,min 示意前 n 天的最低价格
  • 转移方程剖析,能够预感的是,在 [0, n] 这个区间外面,f(n)的最大值就是第 n-1 天的利润,第 n 天价格和前 n 天的最低价格 min 的较大值,即 f(n-1)prices[n] - min 两者中的较大值,从而,咱们能够得出其转移方程为:
    $f(n) = max{(f(n-1), prices[n] – min)}$
  • 初始值剖析,最小值 min 间接应用 第一天的股票价格即可,max 赋值为 0,通过循环更新
  • 返回值剖析,返回值为 max

    function maxProfit(prices: number[]): number {let min = prices[0], max = 0
    
      for (let n = 0; n < prices.length; n++) {min = Math.min(min, prices[n])
          max = Math.max(max, prices[n] - min)
      }
    
      return max
    };

当然,本题也能够通过奢侈循环求解,这里就不剖析了,贴下代码

暴力法

function maxProfit(prices: number[]): number {
        let max = 0;
        for (let i = 0; i < prices.length - 1; i++) {for (let j = i + 1; j < prices.length; j++) {let profit = prices[j] - prices[i];
                if (profit > max) {max = profit;}
            }
        }
        return max;
    }
}

最长回升子序列 ⭐⭐

👉 Leetcode 链接

给你一个整数数组 nums,找到其中最长严格递增子序列的长度。

子序列是由数组派生而来的序列,删除(或不删除)数组中的元素而不扭转其余元素的程序。例如,[3,6,2,7] 是数组 [0,3,1,6,2,2,7] 的子序列。

用例

  • 输出:nums = [10,9,2,5,3,7,101,18]
  • 输入:4
  • 解释:最长递增子序列是 [2,3,7,101],因而长度为 4

思路剖析

做了两道题目之后,你对动静布局有没有些感觉了呢😜,咱们接着来看这道经典的动静布局题,如何获取 最长严格递增子序列 ?惯例思路剖析就是 数组中后一个数比后面每一个都要大,一旦呈现小的,连续性就不存在了

  • 状态定义,咱们能够用 f(n) 示意前 n 个回升子序列
  • 转移方程剖析,什么状况下咱们能够确定数组中后一个数比后面每一个都要大?这里能够预感的是,应用双重循环,[0, n] 这个区间为外层循环,[0, m] 这个区间为内层循环,且 m<n, 如果第 n 个数为回升子序列,则前 n 个数势必比任何一个 [f(0), f(m)]的取值要大,所以其转移方程为:
    $ f(n) = max{(f(n), f(m)+1)}, m<n$
  • 初始值剖析,最大值 max 赋值为 0,通过循环更新,dp 为长度为数组长度的新数组,其初始值均为 1
  • 返回值剖析,返回值为 max

    function lengthOfLIS(nums: number[]): number {let len = nums.length, dp = new Array(len).fill(1), max = 0
    
      if (len === 0) return 0
      
      for (let i = 0; i < len; i++){for (let j = 0; j < i; j++){if (nums[i] > nums[j]) {dp[i] = Math.max(dp[i], dp[j] + 1)
              }
          }
    
          max = Math.max(max, dp[i])
      }
    
      return max
    };

这道题工夫复杂度曾经是 $O(n^2)$ 了,思考下,有没有其余更优的解法呢?

最长递增子序列的个数 ⭐⭐⭐

👉 Leetcode 链接

给定一个未排序的整数数组,找到最长递增子序列的个数。

用例

  • 输出:[1,3,5,4,7]
  • 输入:2
  • 解释:有两个最长递增子序列,别离是 [1, 3, 4, 7] 和[1, 3, 5, 7]。

思路剖析

这道题目是上一题的变形,上一题求解了 最长递增子序列 ,那 最长递增子序列 是否只有一个呢?当然不是了,有几个,正是本题的求解目标

首先思考下,这道题和上道题不同的中央在哪里,就是要在 求解最长递增子序列的同时,记录其个数 ,那要怎么记录这个 个数 ,就是本题的要害了,对于上个题目,咱们给出的转移方程是 $ f(n) = max{(f(n), f(m)+1)}, m<n$,
那对于有多少个这样的最长子序列,咱们只须要将合乎 $f(n)=f(m) + 1$
的组合筛选进去即可,所以给予之前的做法,咱们只须要在对应条件内,再进行判断即可

  • 状态定义,咱们能够用 f(n) 示意前 n 个回升子序列
  • 转移方程剖析,什么状况下咱们能够确定数组中后一个数比后面每一个都要大?这里能够预感的是,应用双重循环,[0, n] 这个区间为外层循环,[0, m] 这个区间为内层循环,且 m< n, 如果第 n 个数为回升子序列,则前 n 个数势必比任何一个 [f(0), f(m)]的取值要大,所以其转移方程为:
    $ f(n) = max{(f(n), f(m)+1)}, m<n$,判断条件 $ f(n)=f(m) + 1$
  • 初始值剖析,最大值 max 赋值为 0,通过循环更新,dp 为长度为数组长度的新数组,其初始值均为 1,count 为长度为数组长度的新数组,其初始值均为 1
  • 返回值剖析,返回值为 count

    function findNumberOfLIS(nums: number[]): number {let len = nums.length, max = 0, res = 0, dp = new Array(len).fill(1), count = new Array(len).fill(1)
    
      for (let i = 0; i < len; i++){for (let j = 0; j < i; j++){if (nums[i] > nums[j]) {if (dp[i] < dp[j] + 1) {dp[i] = dp[j] + 1
                      count[i] = count[j]
                  } else if (dp[i] === dp[j] + 1) {count[i] += count[j]
                  }
              }
          }
    
          if (dp[i] > max) {max = dp[i]
              res = count[i]
          } else if(max === dp[i]) {res += count[i]
          }
      }
    
      return res
    };

    总结

    本期动静布局通过四个常见题目,简要介绍了动静布局问题的常见解题思路,从状态定义、转移方程、初始值和返回值四方面进行了剖析,其中 转移方程 是最外围的,而转移方程的确立来自对问题和子问题的梳理和合成,当你有了一个正确的转移方程,就曾经胜利了一大半,再思考好相干边界取值,根本就功败垂成了。

当然,动静布局可不止这些问题,比方:矩阵运算、前缀和、背包和门路压缩等问题,下一期咱们将介绍进阶版本的动静布局,敬请期待。

正文完
 0