leetcode题解之53. 最大子序和

给定一个整数数组 nums ,找到一个具有最大和的连续子数组(子数组最少包含一个元素),返回其最大和。

示例:

输入: [-2,1,-3,4,-1,2,1,-5,4],
输出: 6
解释: 连续子数组 [4,-1,2,1] 的和最大,为 6。

进阶:

如果你已经实现复杂度为 O(n) 的解法,尝试使用更为精妙的分治法求解。

📺视频题解

code:

vid:

uuid:

requestId:

播放时间:

提示信息

字幕
    倍速
      清晰度
        音轨
          倍速正常
          字幕Off
          音轨
          清晰度

          📖文字题解

          方法一:动态规划

          思路和算法

          假设 nums 数组的长度是 nn,下标从 00n1n - 1

          我们用 aia_i 代表 nums[i],用 f(i)f(i) 代表以第 ii 个数结尾的「连续子数组的最大和」,那么很显然我们要求的答案就是:

          max0in1{f(i)}\max_{0 \leq i \leq n - 1} \{ f(i) \}

          因此我们只需要求出每个位置的 f(i)f(i),然后返回 f 数组中的最大值即可。那么我们如何求 f(i)f(i) 呢?我们可以考虑 aia_i 单独成为一段还是加入 f(i1)f(i - 1) 对应的那一段,这取决于 aia_if(i1)+aif(i - 1) + a_i 的大小,我们希望获得一个比较大的,于是可以写出这样的动态规划转移方程:

          f(i)=max{f(i1)+ai,ai}f(i) = \max \{ f(i - 1) + a_i, a_i \}

          不难给出一个时间复杂度 O(n)O(n)、空间复杂度 O(n)O(n) 的实现,即用一个 f 数组来保存 f(i)f(i) 的值,用一个循环求出所有 f(i)f(i)。考虑到 f(i)f(i) 只和 f(i1)f(i - 1) 相关,于是我们可以只用一个变量 pre 来维护对于当前 f(i)f(i)f(i1)f(i - 1) 的值是多少,从而让空间复杂度降低到 O(1)O(1),这有点类似「滚动数组」的思想。

          代码

          class Solution {
          public:
              int maxSubArray(vector<int>& nums) {
                  int pre = 0, maxAns = nums[0];
                  for (const auto &x: nums) {
                      pre = max(pre + x, x);
                      maxAns = max(maxAns, pre);
                  }
                  return maxAns;
              }
          };
          
          public class Solution {
              public int MaxSubArray(int[] nums) {
                  int pre = 0, maxAns = nums[0];
                  foreach (int x in nums) {
                      pre = Math.Max(pre + x, x);
                      maxAns = Math.Max(maxAns, pre);
                  }
                  return maxAns;
              }
          }
          
          var maxSubArray = function(nums) {
              let pre = 0, maxAns = nums[0];
              nums.forEach((x) => {
                  pre = Math.max(pre + x, x);
                  maxAns = Math.max(maxAns, pre);
              });
              return maxAns;
          };
          
          func maxSubArray(nums []int) int {
              max := nums[0]
              for i := 1; i < len(nums); i++ {
                  if nums[i] + nums[i-1] > nums[i] {
                      nums[i] += nums[i-1]
                  }
                  if nums[i] > max {
                      max = nums[i]
                  }
              }
              return max
          }
          

          复杂度

          • 时间复杂度:O(n)O(n),其中 nnnums 数组的长度。我们只需要遍历一遍数组即可求得答案。
          • 空间复杂度:O(1)O(1)。我们只需要常数空间存放若干变量。

          方法二:分治

          思路和算法

          这个分治方法类似于「线段树求解 LCIS 问题」的 pushUp 操作。 也许读者还没有接触过线段树,没有关系,方法二的内容假设你没有任何线段树的基础。当然,如果读者有兴趣的话,推荐看一看线段树区间合并法解决 多次询问 的「区间最长连续上升序列问题」和「区间最大子段和问题」,还是非常有趣的。

          我们定义一个操作 get(a, l, r) 表示查询 aa 序列 [l,r][l, r] 区间内的最大子段和,那么最终我们要求的答案就是 get(nums, 0, nums.size() - 1)。如何分治实现这个操作呢?对于一个区间 [l,r][l, r],我们取 m=l+r2m = \lfloor \frac{l + r}{2} \rfloor,对区间 [l,m][l, m][m+1,r][m + 1, r] 分治求解。当递归逐层深入直到区间长度缩小为 11 的时候,递归「开始回升」。这个时候我们考虑如何通过 [l,m][l, m] 区间的信息和 [m+1,r][m + 1, r] 区间的信息合并成区间 [l,r][l, r] 的信息。最关键的两个问题是:

          • 我们要维护区间的哪些信息呢?
          • 我们如何合并这些信息呢?

          对于一个区间 [l,r][l, r],我们可以维护四个量:

          • lSum 表示 [l,r][l, r] 内以 ll 为左端点的最大子段和
          • rSum 表示 [l,r][l, r] 内以 rr 为右端点的最大子段和
          • mSum 表示 [l,r][l, r] 内的最大子段和
          • iSum 表示 [l,r][l, r] 的区间和

          以下简称 [l,m][l, m][l,r][l, r] 的「左子区间」,[m+1,r][m + 1, r][l,r][l, r] 的「右子区间」。我们考虑如何维护这些量呢(如何通过左右子区间的信息合并得到 [l,r][l, r] 的信息)?对于长度为 11 的区间 [i,i][i, i],四个量的值都和 aia_i 相等。对于长度大于 11 的区间:

          • 首先最好维护的是 iSum,区间 [l,r][l, r]iSum 就等于「左子区间」的 iSum 加上「右子区间」的 iSum
          • 对于 [l,r][l, r]lSum,存在两种可能,它要么等于「左子区间」的 lSum,要么等于「左子区间」的 iSum 加上「右子区间」的 lSum,二者取大。
          • 对于 [l,r][l, r]rSum,同理,它要么等于「右子区间」的 rSum,要么等于「右子区间」的 iSum 加上「左子区间」的 rSum,二者取大。
          • 当计算好上面的三个量之后,就很好计算 [l,r][l, r]mSum 了。我们可以考虑 [l,r][l, r]mSum 对应的区间是否跨越 mm——它可能不跨越 mm,也就是说 [l,r][l, r]mSum 可能是「左子区间」的 mSum 和 「右子区间」的 mSum 中的一个;它也可能跨越 mm,可能是「左子区间」的 rSum 和 「右子区间」的 lSum 求和。三者取大。

          这样问题就得到了解决。

          代码

          class Solution {
          public:
              struct Status {
                  int lSum, rSum, mSum, iSum;
              };
          
          <span class="hljs-function">Status <span class="hljs-title">pushUp</span><span class="hljs-params">(Status l, Status r)</span> </span>{
              <span class="hljs-keyword">int</span> iSum = l.iSum + r.iSum;
              <span class="hljs-keyword">int</span> lSum = max(l.lSum, l.iSum + r.lSum);
              <span class="hljs-keyword">int</span> rSum = max(r.rSum, r.iSum + l.rSum);
              <span class="hljs-keyword">int</span> mSum = max(max(l.mSum, r.mSum), l.rSum + r.lSum);
              <span class="hljs-keyword">return</span> (Status) {lSum, rSum, mSum, iSum};
          };
          
          <span class="hljs-function">Status <span class="hljs-title">get</span><span class="hljs-params">(<span class="hljs-built_in">vector</span>&lt;<span class="hljs-keyword">int</span>&gt; &amp;a, <span class="hljs-keyword">int</span> l, <span class="hljs-keyword">int</span> r)</span> </span>{
              <span class="hljs-keyword">if</span> (l == r) <span class="hljs-keyword">return</span> (Status) {a[l], a[l], a[l], a[l]};
              <span class="hljs-keyword">int</span> m = (l + r) &gt;&gt; <span class="hljs-number">1</span>;
              Status lSub = get(a, l, m);
              Status rSub = get(a, m + <span class="hljs-number">1</span>, r);
              <span class="hljs-keyword">return</span> pushUp(lSub, rSub);
          }
          
          <span class="hljs-function"><span class="hljs-keyword">int</span> <span class="hljs-title">maxSubArray</span><span class="hljs-params">(<span class="hljs-built_in">vector</span>&lt;<span class="hljs-keyword">int</span>&gt;&amp; nums)</span> </span>{
              <span class="hljs-keyword">return</span> get(nums, <span class="hljs-number">0</span>, nums.size() - <span class="hljs-number">1</span>).mSum;
          }
          

          };


          public class Solution {
          public class Status {
          public int lSum, rSum, mSum, iSum;
          public Status(int lSum_, int rSum_, int mSum_, int iSum_) {
          lSum = lSum_; rSum = rSum_; mSum = mSum_; iSum = iSum_;
          }
          }

          <span class="hljs-function"><span class="hljs-keyword">public</span> Status <span class="hljs-title">pushUp</span>(<span class="hljs-params">Status l, Status r</span>)</span> {
              <span class="hljs-keyword">int</span> iSum = l.iSum + r.iSum;
              <span class="hljs-keyword">int</span> lSum = Math.Max(l.lSum, l.iSum + r.lSum);
              <span class="hljs-keyword">int</span> rSum = Math.Max(r.rSum, r.iSum + l.rSum);
              <span class="hljs-keyword">int</span> mSum = Math.Max(Math.Max(l.mSum, r.mSum), l.rSum + r.lSum);
              <span class="hljs-keyword">return</span> <span class="hljs-keyword">new</span> Status(lSum, rSum, mSum, iSum);
          }
          
          <span class="hljs-function"><span class="hljs-keyword">public</span> Status <span class="hljs-title">getInfo</span>(<span class="hljs-params"><span class="hljs-keyword">int</span>[] a, <span class="hljs-keyword">int</span> l, <span class="hljs-keyword">int</span> r</span>)</span> {
              <span class="hljs-keyword">if</span> (l == r) <span class="hljs-keyword">return</span> <span class="hljs-keyword">new</span> Status(a[l], a[l], a[l], a[l]);
              <span class="hljs-keyword">int</span> m = (l + r) &gt;&gt; <span class="hljs-number">1</span>;
              Status lSub = getInfo(a, l, m);
              Status rSub = getInfo(a, m + <span class="hljs-number">1</span>, r);
              <span class="hljs-keyword">return</span> pushUp(lSub, rSub);
          }
          
          <span class="hljs-function"><span class="hljs-keyword">public</span> <span class="hljs-keyword">int</span> <span class="hljs-title">MaxSubArray</span>(<span class="hljs-params"><span class="hljs-keyword">int</span>[] nums</span>)</span> {
              <span class="hljs-keyword">return</span> getInfo(nums, <span class="hljs-number">0</span>, nums.Length - <span class="hljs-number">1</span>).mSum;
          }
          

          }


          function Status(l, r, m, i) {
          this.lSum = l;
          this.rSum = r;
          this.mSum = m;
          this.iSum = i;
          }

          const pushUp = (l, r) => {
          const iSum = l.iSum + r.iSum;
          const lSum = Math.max(l.lSum, l.iSum + r.lSum);
          const rSum = Math.max(r.rSum, r.iSum + l.rSum);
          const mSum = Math.max(Math.max(l.mSum, r.mSum), l.rSum + r.lSum);
          return new Status(lSum, rSum, mSum, iSum);
          }

          const getInfo = (a, l, r) => {
          if (l === r) return new Status(a[l], a[l], a[l], a[l]);
          const m = (l + r) >> 1;
          const lSub = getInfo(a, l, m);
          const rSub = getInfo(a, m + 1, r);
          return pushUp(lSub, rSub);
          }

          var maxSubArray = function(nums) {
          return getInfo(nums, 0, nums.length - 1).mSum;
          };


          func maxSubArray(nums []int) int {
          return get(nums, 0, len(nums) - 1).mSum;
          }

          func pushUp(l, r Status) Status {
          iSum := l.iSum + r.iSum
          lSum := max(l.lSum, l.iSum + r.lSum)
          rSum := max(r.rSum, r.iSum + l.rSum)
          mSum := max(max(l.mSum, r.mSum), l.rSum + r.lSum)
          return Status{lSum, rSum, mSum, iSum}
          }

          func get(nums []int, l, r int) Status {
          if (l == r) {
          return Status{nums[l], nums[l], nums[l], nums[l]}
          }
          m := (l + r) >> 1
          lSub := get(nums, l, m)
          rSub := get(nums, m + 1, r)
          return pushUp(lSub, rSub)
          }

          func max(x, y int) int {
          if x > y {
          return x
          }
          return y
          }

          type Status struct {
          lSum, rSum, mSum, iSum int
          }

          复杂度分析

          假设序列 aa 的长度为 nn

          • 时间复杂度:假设我们把递归的过程看作是一颗二叉树的先序遍历,那么这颗二叉树的深度的渐进上界为 O(logn)O(\log n),这里的总时间相当于遍历这颗二叉树的所有节点,故总时间的渐进上界是 O(i=1logn2i1)=O(n)O(\sum_{i = 1}^{\log n} 2^{i - 1}) = O(n),故渐进时间复杂度为 O(n)O(n)
          • 空间复杂度:递归会使用 O(logn)O(\log n) 的栈空间,故渐进空间复杂度为 O(logn)O(\log n)

          题外话

          「方法二」相较于「方法一」来说,时间复杂度相同,但是因为使用了递归,并且维护了四个信息的结构体,运行的时间略长,空间复杂度也不如方法一优秀,而且难以理解。那么这种方法存在的意义是什么呢?

          对于这道题而言,确实是如此的。但是仔细观察「方法二」,它不仅可以解决区间 [0,n1][0, n - 1],还可以用于解决任意的子区间 [l,r][l, r] 的问题。如果我们把 [0,n1][0, n - 1] 分治下去出现的所有子区间的信息都用堆式存储的方式记忆化下来,即建成一颗真正的树之后,我们就可以在 O(logn)O(\log n) 的时间内求到任意区间内的答案,我们甚至可以修改序列中的值,做一些简单的维护,之后仍然可以在 O(logn)O(\log n) 的时间内求到任意区间内的答案,对于大规模查询的情况下,这种方法的优势便体现了出来。这棵树就是上文提及的一种神奇的数据结构——线段树。

          posted @ 2021-02-27 16:12  乐天的java  阅读(113)  评论(0)    收藏  举报