动态规划解题方法

魔幻的 2020 让我们怀疑人生是否存在最优解？我们某个时间的决策究竟是否正确？历史不能改变，但却会重演，我们究竟要从过去中学到什么呢？

让我们一起从动态规划中，来找寻这些问题的答案吧～

（咳咳，今天开始回归算法系列，来聊一聊之前的算法文章中没有讲到的内容。

什么是动态规划

动态规划（Dynamic Programic，简称 DP）是一种求解最优解的方法，它是一种特殊的分治思想，利用它可以实现时间复杂度的优化，有时也可以进行空间复杂度的优化，有时是需要更多的空间的（相比其他方法）。

dynamic 是动态的意思，也就是变化的，programing 可以理解为方程（我瞎说的），结合起来就是动态规划是用状态转移方程来求得最优解的算法。

在解释动态规划的时候，我们顺便理一理和它相关的两种思想——分治和贪心算法。

分治

分治是把大问题分解成若干个子问题，这样的能分解性质就是最优子结构的。

最简单的例子就是小明在解决问题 A 的时候，发现问题 A 是由问题 B 和 C 一起组成的，所以他想要解决问题 A，就需要把 B、C 一起去解决。

动态规划

动态规划是分治法的特例。

动态规划比分治法多了一种，就是重叠的子问题。

什么是重叠的子问题呢？举个例子来讲，可爱的小明遇到了一个可爱的问题，那就是问题 A，但是在前面需要解决一连串的问题，我们用 A1，A2，A3，A4 ... A 来表示，在解决 A1 之后会用它的解去解决类似的问题 A2 ，

然后再去解决 A3 ，最终再去解决 A，这就是重叠的子问题的典型代表。（下面的例题还会解释这个概念）

贪心

贪心比动态规划更加的特殊，它还需要问题满足另外一个性质——贪心选择性质，每次都可以把原问题分解成为一个子问题。

实际上再用动规的例子来说明贪心，在解决 A1，A2，A3，A4 ... A 的时候，他发现解决不光有一种重叠子问题的性质在里面，更有趣的是，解决 A1 需要一种特殊的规则。

例如小明现在在玩电脑游戏，而电脑游戏的最终目的是到达 A ，而他又发现，只要一直往右边走就能到达最终的目的地了。这就是一种贪心的算法，在每次往右边走，就是一种特殊的规则，而走到目的地 A 需要很多重复的子问题，也即每次活动一个单位。

入门

其实在很久之前我写的一篇文章中，以斐波那契数列这道基本题为例，详细阐述了从递归到 DP 的优化方法和思路，以及简单题的不简单的答法，大家不妨先去复习一下：

这才是面试官想听的：从递归到 DP 的优化

然后我们再来看看一般的动态规划解题思路。

解题思路

回到动态规划，这里有四个基本的概念：

• state（状态表示）
• function（转移方程）
• initial（初始化）
• final state（最终的状态）

在刚开始的时候，我们首先需要构建一个存储数据的表格，一般是数组或者矩阵，然后设定好每一个格子到下一个格子需要的转移方程。

然后去执行重复的步骤，从初始化的状态一直计算到最终需要的状态。

回到小明的例子，刚开始的时候小明需要确定一个 state（ A0 代表的是什么），然后找到 A0 与 A1 之间的关系，从初始化开始一直计算到最终的状态。

接下来，我们以 Leetcode 120 来详细的讲解这个算法。

题目描述

现在我们来分析一下这个题目，首先我们分析一下为什么它是一个动态规划的问题。

题目是要找到一种路径的和，这种路径和是要最小的，也就是求一个 最优解 。

因为这是路径，我们就是在每一层里面选择一个合适的数字，然后连成一个路径，在这道题目里面，最小的路径是 2-3-5-1 ，在第一层挑了 2 ，在第二层挑了 3 。也就是说总的问题拆分成了每一层的问题，而每一层之间都有一种依赖性在里面，例如第二层选择了 3 之后只能在 6,5 之中选择一个，不能跳到 7 ，这就是 重叠子问题 。

我们用 f[i][j] 表示从三角形顶部走到位置 [i][j] 的最小路径和。这里的位置 i, j 指的是三角形中第 i 行第 j 列的位置。

由于只能是从一个节点到相邻的两个节点（树），因此要想走到位置 [i][j] ，上一步就只能在位置 [i-1][j-1] 或者位置 [i-1][j] 。

我们在这两个位置中选择一个路径和较小的来进行转移，状态转移方程为：

f[i][j]=min(f[i−1][j−1],f[i−1][j])+c[i][j] ，

其中 c[i][j] 指的是 triangle[i][j] 的数值。

方法一

当设定完通项方程之后，我们还需要设定一些特殊的转化方程：

• 当靠近左边界时，也就是 j = 0 时，于是没有 f[i-1][j−1] 这一项 ，状态转移方程变为：

f[i][0]=f[i−1][0]+c[i][0]

• 当靠近右边界时，我们直接用上一层斜上角位置的数值进行计算：

f[i][i]=f[i−1][i−1]+c[i][i]

最终，我们只需要在 dp 三角形的最后一行找到最小值就可以了。

那么初始的状态是什么呢？

实际上就是刚开始的时候设定 dp 的第一个单位的数值为 cp[0][0] ，也即是 dp[0][0] = c[0][0] 。

状态转换图如下所示：

代码如下：

class Solution {
    public int minimumTotal(List<List<Integer>> triangle) {
       // 创建表格
        int[][] dp = new int[triangle.size()][triangle.size()];
        dp[0][0] = triangle.get(0).get(0);
       // 动态规划的方程式
        for (int i = 1; i < dp.length; i++) {
            for (int j = 0; j <= i; j++) {
                int curr = triangle.get(i).get(j);
                if (j == 0) {
                    dp[i][j] = dp[i-1][j] + curr;
                } else if (j == i) {
                    dp[i][j] = dp[i-1][j-1] + curr;
                } else {
                    dp[i][j] = Math.min(dp[i-1][j-1], dp[i-1][j]) + curr;
                }
            }
        }
        int res = dp[triangle.size()-1][0];
        for (int j = 1; j < triangle.size(); j++) {
            res = Math.min(res, dp[triangle.size()-1][j]);
        }
        return res;
    }
}

下面来分析这个问题的时间复杂度以及空间复杂度，一般来说空间复杂度是就是 DP 表格的大小。

在这道问题中是 O(n^2) ，而对于时间复杂度来说，就是整个 dp 的遍历次数，而在这个问题中我们只进行了一次遍历，也即一个矩阵的遍历，所以是 O(n^2) 。

而如果想要优化到 O(n) ，我们需要怎么做呢？

实际上这个就涉及到了一种 状态压缩 的方法，也即压缩这个状态表。

那么怎么去压缩呢？

这个问题比较简单，因为 dp[i][j] 仅仅与上一层的状态有关，所以说与前两层的是没有任何关系的，因此我们不必存储这些无关的状态。

实际上最简单的状态压缩就是保留好前两个状态即可，例如在计算第四行的时候，保留第三行以及第二行的状态表，然后交替的进行更新就可以啦。

class Solution {
    public int minimumTotal(List<List<Integer>> triangle) {
      // 只保留最近 2 行
        int[][] dp = new int[2][triangle.size()];
        dp[0][0] = triangle.get(0).get(0);
        for (int i = 1; i < triangle.size(); i++) {
            int row = i % 2;
            int prevRow = (i-1) % 2;
            for (int j = 0; j <= i; j++) {
                if (j == 0) {
                    dp[row][j] = dp[prevRow][j] + triangle.get(i).get(j);
                } else if (j == i) {
                    dp[row][j] = dp[prevRow][j-1] + triangle.get(i).get(j);
                } else {
                    dp[row][j] = Math.min(dp[prevRow][j-1], dp[prevRow][j]) + triangle.get(i).get(j);
                }
            }
        }

        int res = dp[(triangle.size() - 1) % 2][0];
        for (int j = 1; j < triangle.size(); j++) {
            res = Math.min(res, dp[(triangle.size() - 1) % 2][j]);
        }
        return res;
    }
}

这个的空间复杂度是 O(2n) ，能不能压缩成严格意义上的 O(n) 呢？

那么再往后是否还能够进行状态的压缩呢？

答案是可以的，我们可以再想一种方程然后达到最优的空间复杂度的目标。

当我们在计算位置 [i][j] 时， f[j+1] 到 f[i] 已经是第 i 行的值，而 f[0] 到 f[j] 仍然是第 i-1 行的值。

此时我们直接通过 f[j] = min(f[j−1], f[j]) + c[i][j] 来计算，但是这个时候我们需要 j 是倒着遍历的，因为这样才不会影响之前记录下的状态。

如果从 1 开始，那么计算 2 的时候就会用到新的 1 的数值而不是上一层 1 的数值。

代码如下：

class Solution {
    public int minimumTotal(List<List<Integer>> triangle) {
        int[] dp = new int[triangle.size()];
        dp[0] = triangle.get(0).get(0);

        for (int i = 1; i < triangle.size(); i++) {
            dp[i] = dp[i-1] + triangle.get(i).get(i);
            for (int j = i-1; j > 0; j--) {
                dp[j] = Math.min(dp[j-1], dp[j]) + triangle.get(i).get(j);
            }
            dp[0] += triangle.get(i).get(0);
        }

        int res = dp[0];
        for (int j = 1; j < triangle.size(); j++) {
            res = Math.min(res, dp[j]);
        }
        return res;
    }
}

方法 2

方法 1 有点绕，但如果自下向上来理解，就会变得很简单，这个方法也叫 bottom-up ，方法 1 则是 top-down 。

从结果出发，这个问题是一个 发散三角树 的问题，从最后一行出发，然后每一行每一行的进行递推，那么第一行就是最终的结果了。

举个最简单的例子：

如果从最底下往上出发，实际上找最小值方法的规律很容易找到，那就是在第二行 [1， 2] 里面选择一个就可以了，因为他们两个都要走到根节点。

也就是在下一行的两个数里面取个小的就行了，那么结果就是第一行的数值。

我们先用二维的转移矩阵来解释这个问题，用这种方法也不需要考虑方法 1 里面的边界条件了：

dp[i][j] = min(dp[i+1][j], dp[i+1][j+1]) + c[i][j]

状态转换图如下所示：

class Solution {
    public int minimumTotal(List<List<Integer>> triangle) {
        int[][] dp = new int[triangle.size()][triangle.size()];
        // 创建 DP 空间
        for (int i = 0; i < triangle.size(); i++) {
            dp[triangle.size() - 1][i] = triangle.get(triangle.size() - 1).get(i);
        }

        for (int i = triangle.size() - 2; i >= 0; i--) {
            for (int j = 0; j < triangle.get(i).size(); j++) {
                dp[i][j] = Math.min(dp[i+1][j], dp[i+1][j+1]) + triangle.get(i).get(j);
            }
        }

        return dp[0][0];
    }
}

那么在进行状态压缩的时候，我们该怎么去做呢？

实际上就是只用一个状态表来表示所有的。

因为只是和上一个状态相关，所以说可以表示成如下的形式：

dp[j] = min(dp[j], dp[j + 1]) + triangle[i][j] ，

我们只用 j 来代表当前的状态，然后最终输出 dp[0] 即可。

class Solution {
    public int minimumTotal(List<List<Integer>> triangle) {
        int[] dp = new int[triangle.size()];
        // 创建 DP 空间
        for (int i = 0; i < triangle.size(); i++) {
            dp[i] = triangle.get(triangle.size() - 1).get(i);
        }

        for (int i = triangle.size() - 2; i >= 0; i--) {
            for (int j = 0; j <= i; j++) {
                dp[j] = Math.min(dp[j], dp[j+1]) + triangle.get(i).get(j);
            }
        }

        return dp[0];
    }
}

总结

以上就是动态规划解题方法的粗浅介绍，总的来说，我们需要注意动态规划的这么几件事情。

1. 确定是否需要用动态规划；
2. 确定动态规划的四个部分；
3. 写代码。

实际上难点就是 转移方程 ，这个确实需要大量的积累才能够在面试的时候看穿，甚至有些题没见过的话就是想不出来的。

但是没见过就做不出来的题面试一般也不会考，所以大家也不用太担心，重点还是掌握方法，举一反三。

接下来我也会归纳总结一些动态规划的常见题型，和大家一起探索 最优解 。

更多算法文章，建议收藏这个链接：

• 齐姐聊算法

我是小齐，后端开发工程师，坐标纽约，曾在投行做 Quant，后加入科技公司 FLAG 中的一家。

业余时间打理公粽号【NYCSDE】，更新略少，干货很多，建议加个星标防止错过。

#感谢您访问本站#
#本文转载自互联网，若侵权，请联系删除，谢谢！657271#qq.com#