leetcode刷题题解之寻找两个正序数组的中位数

给定两个大小为 m 和 n 的正序(从小到大)数组 nums1 和 nums2

请你找出这两个正序数组的中位数,并且要求算法的时间复杂度为 O(log(m + n))。

你可以假设 nums1 和 nums2 不会同时为空。

 

示例 1:

nums1 = [1, 3]
nums2 = [2]

则中位数是 2.0

示例 2:

nums1 = [1, 2]
nums2 = [3, 4]

则中位数是 (2 + 3)/2 = 2.5

📺视频题解

code:

vid:

uuid:

requestId:

播放时间:

提示信息

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

          📖文字题解

          方法一:二分查找

          给定两个有序数组,要求找到两个有序数组的中位数,最直观的思路有以下两种:

          • 使用归并的方式,合并两个有序数组,得到一个大的有序数组。大的有序数组的中间位置的元素,即为中位数。

          • 不需要合并两个有序数组,只要找到中位数的位置即可。由于两个数组的长度已知,因此中位数对应的两个数组的下标之和也是已知的。维护两个指针,初始时分别指向两个数组的下标 00 的位置,每次将指向较小值的指针后移一位(如果一个指针已经到达数组末尾,则只需要移动另一个数组的指针),直到到达中位数的位置。

          假设两个有序数组的长度分别为 mmnn,上述两种思路的复杂度如何?

          第一种思路的时间复杂度是 O(m+n)O(m+n),空间复杂度是 O(m+n)O(m+n)。第二种思路虽然可以将空间复杂度降到 O(1)O(1),但是时间复杂度仍是 O(m+n)O(m+n)。题目要求时间复杂度是 O(log(m+n))O(\log(m+n)),因此上述两种思路都不满足题目要求的时间复杂度。

          如何把时间复杂度降低到 O(log(m+n))O(\log(m+n)) 呢?如果对时间复杂度的要求有 log\log,通常都需要用到二分查找,这道题也可以通过二分查找实现。

          根据中位数的定义,当 m+nm+n 是奇数时,中位数是两个有序数组中的第 (m+n)/2(m+n)/2 个元素,当 m+nm+n 是偶数时,中位数是两个有序数组中的第 (m+n)/2(m+n)/2 个元素和第 (m+n)/2+1(m+n)/2+1 个元素的平均值。因此,这道题可以转化成寻找两个有序数组中的第 kk 小的数,其中 kk(m+n)/2(m+n)/2(m+n)/2+1(m+n)/2+1

          假设两个有序数组分别是 A\text{A}B\text{B}。要找到第 kk 个元素,我们可以比较 A[k/21]\text{A}[k/2-1]B[k/21]\text{B}[k/2-1],其中 // 表示整数除法。由于 A[k/21]\text{A}[k/2-1]B[k/21]\text{B}[k/2-1] 的前面分别有 A[0..k/22]\text{A}[0\,..\,k/2-2]B[0..k/22]\text{B}[0\,..\,k/2-2],即 k/21k/2-1 个元素,对于 A[k/21]\text{A}[k/2-1]B[k/21]\text{B}[k/2-1] 中的较小值,最多只会有 (k/21)+(k/21)k/22(k/2-1)+(k/2-1) \leq k/2-2 个元素比它小,那么它就不能是第 kk 小的数了。

          因此我们可以归纳出三种情况:

          • 如果 A[k/21]<B[k/21]\text{A}[k/2-1] < \text{B}[k/2-1],则比 A[k/21]\text{A}[k/2-1] 小的数最多只有 A\text{A} 的前 k/21k/2-1 个数和 B\text{B} 的前 k/21k/2-1 个数,即比 A[k/21]\text{A}[k/2-1] 小的数最多只有 k2k-2 个,因此 A[k/21]\text{A}[k/2-1] 不可能是第 kk 个数,A[0]\text{A}[0]A[k/21]\text{A}[k/2-1] 也都不可能是第 kk 个数,可以全部排除。

          • 如果 A[k/21]>B[k/21]\text{A}[k/2-1] > \text{B}[k/2-1],则可以排除 B[0]\text{B}[0]B[k/21]\text{B}[k/2-1]

          • 如果 A[k/21]=B[k/21]\text{A}[k/2-1] = \text{B}[k/2-1],则可以归入第一种情况处理。

          fig1

          可以看到,比较 A[k/21]\text{A}[k/2-1]B[k/21]\text{B}[k/2-1] 之后,可以排除 k/2k/2 个不可能是第 kk 小的数,查找范围缩小了一半。同时,我们将在排除后的新数组上继续进行二分查找,并且根据我们排除数的个数,减少 kk 的值,这是因为我们排除的数都不大于第 kk 小的数。

          有以下三种情况需要特殊处理:

          • 如果 A[k/21]\text{A}[k/2-1] 或者 B[k/21]\text{B}[k/2-1] 越界,那么我们可以选取对应数组中的最后一个元素。在这种情况下,我们必须根据排除数的个数减少 kk 的值,而不能直接将 kk 减去 k/2k/2

          • 如果一个数组为空,说明该数组中的所有元素都被排除,我们可以直接返回另一个数组中第 kk 小的元素。

          • 如果 k=1k=1,我们只要返回两个数组首元素的最小值即可。

          用一个例子说明上述算法。假设两个有序数组如下:

          A: 1 3 4 9
          B: 1 2 3 4 5 6 7 8 9
          

          两个有序数组的长度分别是 4499,长度之和是 1313,中位数是两个有序数组中的第 77 个元素,因此需要找到第 k=7k=7 个元素。

          比较两个有序数组中下标为 k/21=2k/2-1=2 的数,即 A[2]\text{A}[2]B[2]\text{B}[2],如下面所示:

          A: 1 3 4 9
                 
          B: 1 2 3 4 5 6 7 8 9
                 
          

          由于 A[2]>B[2]\text{A}[2] > \text{B}[2],因此排除 B[0]\text{B}[0]B[2]\text{B}[2],即数组 B\text{B} 的下标偏移(offset)变为 33,同时更新 kk 的值:k=kk/2=4k=k-k/2=4

          下一步寻找,比较两个有序数组中下标为 k/21=1k/2-1=1 的数,即 A[1]\text{A}[1]B[4]\text{B}[4],如下面所示,其中方括号部分表示已经被排除的数。

          A: 1 3 4 9
               
          B: [1 2 3] 4 5 6 7 8 9
                       
          

          由于 A[1]<B[4]\text{A}[1] < \text{B}[4],因此排除 A[0]\text{A}[0]A[1]\text{A}[1],即数组 A\text{A} 的下标偏移变为 22,同时更新 kk 的值:k=kk/2=2k=k-k/2=2

          下一步寻找,比较两个有序数组中下标为 k/21=0k/2-1=0 的数,即比较 A[2]\text{A}[2]B[3]\text{B}[3],如下面所示,其中方括号部分表示已经被排除的数。

          A: [1 3] 4 9
                   
          B: [1 2 3] 4 5 6 7 8 9
                     
          

          由于 A[2]=B[3]\text{A}[2]=B[3],根据之前的规则,排除 A\text{A} 中的元素,因此排除 A[2]\text{A}[2],即数组 A\text{A} 的下标偏移变为 33,同时更新 kk 的值: k=kk/2=1k=k-k/2=1

          由于 kk 的值变成 11,因此比较两个有序数组中的未排除下标范围内的第一个数,其中较小的数即为第 kk 个数,由于 A[3]>B[3]\text{A}[3] > \text{B}[3],因此第 kk 个数是 B[3]=4\text{B}[3]=4

          A: [1 3 4] 9
                     
          B: [1 2 3] 4 5 6 7 8 9
                     
          
          class Solution {
              public double findMedianSortedArrays(int[] nums1, int[] nums2) {
                  int length1 = nums1.length, length2 = nums2.length;
                  int totalLength = length1 + length2;
                  if (totalLength % 2 == 1) {
                      int midIndex = totalLength / 2;
                      double median = getKthElement(nums1, nums2, midIndex + 1);
                      return median;
                  } else {
                      int midIndex1 = totalLength / 2 - 1, midIndex2 = totalLength / 2;
                      double median = (getKthElement(nums1, nums2, midIndex1 + 1) + getKthElement(nums1, nums2, midIndex2 + 1)) / 2.0;
                      return median;
                  }
              }
          
          <span class="hljs-function"><span class="hljs-keyword">public</span> <span class="hljs-keyword">int</span> <span class="hljs-title">getKthElement</span><span class="hljs-params">(<span class="hljs-keyword">int</span>[] nums1, <span class="hljs-keyword">int</span>[] nums2, <span class="hljs-keyword">int</span> k)</span> </span>{
              <span class="hljs-comment">/* 主要思路:要找到第 k (k&gt;1) 小的元素,那么就取 pivot1 = nums1[k/2-1] 和 pivot2 = nums2[k/2-1] 进行比较
               * 这里的 "/" 表示整除
               * nums1 中小于等于 pivot1 的元素有 nums1[0 .. k/2-2] 共计 k/2-1 个
               * nums2 中小于等于 pivot2 的元素有 nums2[0 .. k/2-2] 共计 k/2-1 个
               * 取 pivot = min(pivot1, pivot2),两个数组中小于等于 pivot 的元素共计不会超过 (k/2-1) + (k/2-1) &lt;= k-2 个
               * 这样 pivot 本身最大也只能是第 k-1 小的元素
               * 如果 pivot = pivot1,那么 nums1[0 .. k/2-1] 都不可能是第 k 小的元素。把这些元素全部 "删除",剩下的作为新的 nums1 数组
               * 如果 pivot = pivot2,那么 nums2[0 .. k/2-1] 都不可能是第 k 小的元素。把这些元素全部 "删除",剩下的作为新的 nums2 数组
               * 由于我们 "删除" 了一些元素(这些元素都比第 k 小的元素要小),因此需要修改 k 的值,减去删除的数的个数
               */</span>
          
              <span class="hljs-keyword">int</span> length1 = nums1.length, length2 = nums2.length;
              <span class="hljs-keyword">int</span> index1 = <span class="hljs-number">0</span>, index2 = <span class="hljs-number">0</span>;
              <span class="hljs-keyword">int</span> kthElement = <span class="hljs-number">0</span>;
          
              <span class="hljs-keyword">while</span> (<span class="hljs-keyword">true</span>) {
                  <span class="hljs-comment">// 边界情况</span>
                  <span class="hljs-keyword">if</span> (index1 == length1) {
                      <span class="hljs-keyword">return</span> nums2[index2 + k - <span class="hljs-number">1</span>];
                  }
                  <span class="hljs-keyword">if</span> (index2 == length2) {
                      <span class="hljs-keyword">return</span> nums1[index1 + k - <span class="hljs-number">1</span>];
                  }
                  <span class="hljs-keyword">if</span> (k == <span class="hljs-number">1</span>) {
                      <span class="hljs-keyword">return</span> Math.min(nums1[index1], nums2[index2]);
                  }
                  
                  <span class="hljs-comment">// 正常情况</span>
                  <span class="hljs-keyword">int</span> half = k / <span class="hljs-number">2</span>;
                  <span class="hljs-keyword">int</span> newIndex1 = Math.min(index1 + half, length1) - <span class="hljs-number">1</span>;
                  <span class="hljs-keyword">int</span> newIndex2 = Math.min(index2 + half, length2) - <span class="hljs-number">1</span>;
                  <span class="hljs-keyword">int</span> pivot1 = nums1[newIndex1], pivot2 = nums2[newIndex2];
                  <span class="hljs-keyword">if</span> (pivot1 &lt;= pivot2) {
                      k -= (newIndex1 - index1 + <span class="hljs-number">1</span>);
                      index1 = newIndex1 + <span class="hljs-number">1</span>;
                  } <span class="hljs-keyword">else</span> {
                      k -= (newIndex2 - index2 + <span class="hljs-number">1</span>);
                      index2 = newIndex2 + <span class="hljs-number">1</span>;
                  }
              }
          }
          

          }


          class Solution {
          public:
          int getKthElement(const vector<int>& nums1, const vector<int>& nums2, int k) {
          /* 主要思路:要找到第 k (k>1) 小的元素,那么就取 pivot1 = nums1[k/2-1] 和 pivot2 = nums2[k/2-1] 进行比较
          * 这里的 "/" 表示整除
          * nums1 中小于等于 pivot1 的元素有 nums1[0 .. k/2-2] 共计 k/2-1 个
          * nums2 中小于等于 pivot2 的元素有 nums2[0 .. k/2-2] 共计 k/2-1 个
          * 取 pivot = min(pivot1, pivot2),两个数组中小于等于 pivot 的元素共计不会超过 (k/2-1) + (k/2-1) <= k-2 个
          * 这样 pivot 本身最大也只能是第 k-1 小的元素
          * 如果 pivot = pivot1,那么 nums1[0 .. k/2-1] 都不可能是第 k 小的元素。把这些元素全部 "删除",剩下的作为新的 nums1 数组
          * 如果 pivot = pivot2,那么 nums2[0 .. k/2-1] 都不可能是第 k 小的元素。把这些元素全部 "删除",剩下的作为新的 nums2 数组
          * 由于我们 "删除" 了一些元素(这些元素都比第 k 小的元素要小),因此需要修改 k 的值,减去删除的数的个数
          */

              <span class="hljs-keyword">int</span> m = nums1.size();
              <span class="hljs-keyword">int</span> n = nums2.size();
              <span class="hljs-keyword">int</span> index1 = <span class="hljs-number">0</span>, index2 = <span class="hljs-number">0</span>;
          
              <span class="hljs-keyword">while</span> (<span class="hljs-literal">true</span>) {
                  <span class="hljs-comment">// 边界情况</span>
                  <span class="hljs-keyword">if</span> (index1 == m) {
                      <span class="hljs-keyword">return</span> nums2[index2 + k - <span class="hljs-number">1</span>];
                  }
                  <span class="hljs-keyword">if</span> (index2 == n) {
                      <span class="hljs-keyword">return</span> nums1[index1 + k - <span class="hljs-number">1</span>];
                  }
                  <span class="hljs-keyword">if</span> (k == <span class="hljs-number">1</span>) {
                      <span class="hljs-keyword">return</span> min(nums1[index1], nums2[index2]);
                  }
          
                  <span class="hljs-comment">// 正常情况</span>
                  <span class="hljs-keyword">int</span> newIndex1 = min(index1 + k / <span class="hljs-number">2</span> - <span class="hljs-number">1</span>, m - <span class="hljs-number">1</span>);
                  <span class="hljs-keyword">int</span> newIndex2 = min(index2 + k / <span class="hljs-number">2</span> - <span class="hljs-number">1</span>, n - <span class="hljs-number">1</span>);
                  <span class="hljs-keyword">int</span> pivot1 = nums1[newIndex1];
                  <span class="hljs-keyword">int</span> pivot2 = nums2[newIndex2];
                  <span class="hljs-keyword">if</span> (pivot1 &lt;= pivot2) {
                      k -= newIndex1 - index1 + <span class="hljs-number">1</span>;
                      index1 = newIndex1 + <span class="hljs-number">1</span>;
                  }
                  <span class="hljs-keyword">else</span> {
                      k -= newIndex2 - index2 + <span class="hljs-number">1</span>;
                      index2 = newIndex2 + <span class="hljs-number">1</span>;
                  }
              }
          }
          
          <span class="hljs-function"><span class="hljs-keyword">double</span> <span class="hljs-title">findMedianSortedArrays</span><span class="hljs-params">(<span class="hljs-built_in">vector</span>&lt;<span class="hljs-keyword">int</span>&gt;&amp; nums1, <span class="hljs-built_in">vector</span>&lt;<span class="hljs-keyword">int</span>&gt;&amp; nums2)</span> </span>{
              <span class="hljs-keyword">int</span> totalLength = nums1.size() + nums2.size();
              <span class="hljs-keyword">if</span> (totalLength % <span class="hljs-number">2</span> == <span class="hljs-number">1</span>) {
                  <span class="hljs-keyword">return</span> getKthElement(nums1, nums2, (totalLength + <span class="hljs-number">1</span>) / <span class="hljs-number">2</span>);
              }
              <span class="hljs-keyword">else</span> {
                  <span class="hljs-keyword">return</span> (getKthElement(nums1, nums2, totalLength / <span class="hljs-number">2</span>) + getKthElement(nums1, nums2, totalLength / <span class="hljs-number">2</span> + <span class="hljs-number">1</span>)) / <span class="hljs-number">2.0</span>;
              }
          }
          

          };


          class Solution:
          def findMedianSortedArrays(self, nums1: List[int], nums2: List[int]) -> float:
          def getKthElement(k):
          """
          - 主要思路:要找到第 k (k>1) 小的元素,那么就取 pivot1 = nums1[k/2-1] 和 pivot2 = nums2[k/2-1] 进行比较
          - 这里的 "/" 表示整除
          - nums1 中小于等于 pivot1 的元素有 nums1[0 .. k/2-2] 共计 k/2-1 个
          - nums2 中小于等于 pivot2 的元素有 nums2[0 .. k/2-2] 共计 k/2-1 个
          - 取 pivot = min(pivot1, pivot2),两个数组中小于等于 pivot 的元素共计不会超过 (k/2-1) + (k/2-1) <= k-2 个
          - 这样 pivot 本身最大也只能是第 k-1 小的元素
          - 如果 pivot = pivot1,那么 nums1[0 .. k/2-1] 都不可能是第 k 小的元素。把这些元素全部 "删除",剩下的作为新的 nums1 数组
          - 如果 pivot = pivot2,那么 nums2[0 .. k/2-1] 都不可能是第 k 小的元素。把这些元素全部 "删除",剩下的作为新的 nums2 数组
          - 由于我们 "删除" 了一些元素(这些元素都比第 k 小的元素要小),因此需要修改 k 的值,减去删除的数的个数
          """

                  index1, index2 = <span class="hljs-number">0</span>, <span class="hljs-number">0</span>
                  <span class="hljs-keyword">while</span> <span class="hljs-literal">True</span>:
                      <span class="hljs-comment"># 特殊情况</span>
                      <span class="hljs-keyword">if</span> index1 == m:
                          <span class="hljs-keyword">return</span> nums2[index2 + k - <span class="hljs-number">1</span>]
                      <span class="hljs-keyword">if</span> index2 == n:
                          <span class="hljs-keyword">return</span> nums1[index1 + k - <span class="hljs-number">1</span>]
                      <span class="hljs-keyword">if</span> k == <span class="hljs-number">1</span>:
                          <span class="hljs-keyword">return</span> min(nums1[index1], nums2[index2])
          
                      <span class="hljs-comment"># 正常情况</span>
                      newIndex1 = min(index1 + k // <span class="hljs-number">2</span> - <span class="hljs-number">1</span>, m - <span class="hljs-number">1</span>)
                      newIndex2 = min(index2 + k // <span class="hljs-number">2</span> - <span class="hljs-number">1</span>, n - <span class="hljs-number">1</span>)
                      pivot1, pivot2 = nums1[newIndex1], nums2[newIndex2]
                      <span class="hljs-keyword">if</span> pivot1 &lt;= pivot2:
                          k -= newIndex1 - index1 + <span class="hljs-number">1</span>
                          index1 = newIndex1 + <span class="hljs-number">1</span>
                      <span class="hljs-keyword">else</span>:
                          k -= newIndex2 - index2 + <span class="hljs-number">1</span>
                          index2 = newIndex2 + <span class="hljs-number">1</span>
              
              m, n = len(nums1), len(nums2)
              totalLength = m + n
              <span class="hljs-keyword">if</span> totalLength % <span class="hljs-number">2</span> == <span class="hljs-number">1</span>:
                  <span class="hljs-keyword">return</span> getKthElement((totalLength + <span class="hljs-number">1</span>) // <span class="hljs-number">2</span>)
              <span class="hljs-keyword">else</span>:
                  <span class="hljs-keyword">return</span> (getKthElement(totalLength // <span class="hljs-number">2</span>) + getKthElement(totalLength // <span class="hljs-number">2</span> + <span class="hljs-number">1</span>)) / <span class="hljs-number">2</span>
          


          func findMedianSortedArrays(nums1 []int, nums2 []int) float64 {
          totalLength := len(nums1) + len(nums2)
          if totalLength%2 == 1 {
          midIndex := totalLength/2
          return float64(getKthElement(nums1, nums2, midIndex + 1))
          } else {
          midIndex1, midIndex2 := totalLength/2 - 1, totalLength/2
          return float64(getKthElement(nums1, nums2, midIndex1 + 1) + getKthElement(nums1, nums2, midIndex2 + 1)) / 2.0
          }
          return 0
          }

          func getKthElement(nums1, nums2 []int, k int) int {
          index1, index2 := 0, 0
          for {
          if index1 == len(nums1) {
          return nums2[index2 + k - 1]
          }
          if index2 == len(nums2) {
          return nums1[index1 + k - 1]
          }
          if k == 1 {
          return min(nums1[index1], nums2[index2])
          }
          half := k/2
          newIndex1 := min(index1 + half, len(nums1)) - 1
          newIndex2 := min(index2 + half, len(nums2)) - 1
          pivot1, pivot2 := nums1[newIndex1], nums2[newIndex2]
          if pivot1 <= pivot2 {
          k -= (newIndex1 - index1 + 1)
          index1 = newIndex1 + 1
          } else {
          k -= (newIndex2 - index2 + 1)
          index2 = newIndex2 + 1
          }
          }
          return 0
          }

          func min(x, y int) int {
          if x < y {
          return x
          }
          return y
          }

          复杂度分析

          • 时间复杂度:O(log(m+n))O(\log(m+n)),其中 mmnn 分别是数组 nums1\text{nums1}nums2\text{nums2} 的长度。初始时有 k=(m+n)/2k=(m+n)/2k=(m+n)/2+1k=(m+n)/2+1,每一轮循环可以将查找范围减少一半,因此时间复杂度是 O(log(m+n))O(\log(m+n))

          • 空间复杂度:O(1)O(1)

          方法二:划分数组

          说明

          方法一的时间复杂度已经很优秀了,但本题存在时间复杂度更低的一种方法。这里给出推导过程,勇于挑战自己的读者可以进行尝试。

          思路与算法

          为了使用划分的方法解决这个问题,需要理解「中位数的作用是什么」。在统计中,中位数被用来:

          将一个集合划分为两个长度相等的子集,其中一个子集中的元素总是大于另一个子集中的元素。

          如果理解了中位数的划分作用,就很接近答案了。

          首先,在任意位置 iiA\text{A} 划分成两个部分:

                     left_A            |          right_A
              A[0], A[1], ..., A[i-1]  |  A[i], A[i+1], ..., A[m-1]
          

          由于 A\text{A} 中有 mm 个元素, 所以有 m+1m+1 种划分的方法(i[0,m]i \in [0, m])。

          len(left_A)=i,len(right_A)=mi\text{len}(\text{left\_A}) = i, \text{len}(\text{right\_A}) = m - i.

          注意:当 i=0i = 0 时,left_A\text{left\_A} 为空集, 而当 i=mi = m 时, right_A\text{right\_A} 为空集。

          采用同样的方式,在任意位置 jjB\text{B} 划分成两个部分:

                     left_B            |          right_B
              B[0], B[1], ..., B[j-1]  |  B[j], B[j+1], ..., B[n-1]
          

          left_A\text{left\_A}left_B\text{left\_B} 放入一个集合,并将 right_A\text{right\_A}right_B\text{right\_B} 放入另一个集合。 再把这两个新的集合分别命名为 left_part\text{left\_part}right_part\text{right\_part}

                    left_part          |         right_part
              A[0], A[1], ..., A[i-1]  |  A[i], A[i+1], ..., A[m-1]
              B[0], B[1], ..., B[j-1]  |  B[j], B[j+1], ..., B[n-1]
          

          A\text{A}B\text{B} 的总长度是偶数时,如果可以确认:

          • len(left_part)=len(right_part)\text{len}(\text{left\_part}) = \text{len}(\text{right\_part})
          • max(left_part)min(right_part)\max(\text{left\_part}) \leq \min(\text{right\_part})

          那么,{A,B}\{\text{A}, \text{B}\} 中的所有元素已经被划分为相同长度的两个部分,且前一部分中的元素总是小于或等于后一部分中的元素。中位数就是前一部分的最大值和后一部分的最小值的平均值:

          median=max(left_part)+min(right_part)2\text{median} = \frac{\text{max}(\text{left}\_\text{part}) + \text{min}(\text{right}\_\text{part})}{2}

          A\text{A}B\text{B} 的总长度是奇数时,如果可以确认:

          • len(left_part)=len(right_part)+1\text{len}(\text{left\_part}) = \text{len}(\text{right\_part})+1
          • max(left_part)min(right_part)\max(\text{left\_part}) \leq \min(\text{right\_part})

          那么,{A,B}\{\text{A}, \text{B}\} 中的所有元素已经被划分为两个部分,前一部分比后一部分多一个元素,且前一部分中的元素总是小于或等于后一部分中的元素。中位数就是前一部分的最大值:

          median=max(left_part)\text{median} = \text{max}(\text{left}\_\text{part})

          第一个条件对于总长度是偶数和奇数的情况有所不同,但是可以将两种情况合并。第二个条件对于总长度是偶数和奇数的情况是一样的。

          要确保这两个条件,只需要保证:

          • i+j=mi+nji + j = m - i + n - j(当 m+nm+n 为偶数)或 i+j=mi+nj+1i + j = m - i + n - j + 1(当 m+nm+n 为奇数)。等号左侧为前一部分的元素个数,等号右侧为后一部分的元素个数。将 iijj 全部移到等号左侧,我们就可以得到 i+j=m+n+12i+j = \frac{m + n + 1}{2}。这里的分数结果只保留整数部分。

          • 0im0 \leq i \leq m0jn0 \leq j \leq n。如果我们规定 A\text{A} 的长度小于等于 B\text{B} 的长度,即 mnm \leq n。这样对于任意的 i[0,m]i \in [0, m],都有 j=m+n+12i[0,n]j = \frac{m + n + 1}{2} - i \in [0, n],那么我们在 [0,m][0, m] 的范围内枚举 ii 并得到 jj,就不需要额外的性质了。

            • 如果 A\text{A} 的长度较大,那么我们只要交换 A\text{A}B\text{B} 即可。

            • 如果 m>nm > n ,那么得出的 jj 有可能是负数。

          • B[j1]A[i]\text{B}[j-1] \leq \text{A}[i] 以及 A[i1]B[j]\text{A}[i-1] \leq \text{B}[j],即前一部分的最大值小于等于后一部分的最小值。

          为了简化分析,假设 A[i1],B[j1],A[i],B[j]\text{A}[i-1], \text{B}[j-1], \text{A}[i], \text{B}[j] 总是存在。对于 i=0i=0i=mi=mj=0j=0j=nj=n 这样的临界条件,我们只需要规定 A[1]=B[1]=A[-1]=B[-1]=-\inftyA[m]=B[n]=A[m]=B[n]=\infty 即可。这也是比较直观的:当一个数组不出现在前一部分时,对应的值为负无穷,就不会对前一部分的最大值产生影响;当一个数组不出现在后一部分时,对应的值为正无穷,就不会对后一部分的最小值产生影响。

          所以我们需要做的是:

          [0,m][0, m] 中找到 ii,使得:

          B[j1]A[i]\qquad \text{B}[j-1] \leq \text{A}[i]A[i1]B[j]\text{A}[i-1] \leq \text{B}[j],其中 j=m+n+12ij = \frac{m + n + 1}{2} - i

          我们证明它等价于:

          [0,m][0, m] 中找到最大的 ii,使得:

          A[i1]B[j]\qquad \text{A}[i-1] \leq \text{B}[j],其中 j=m+n+12ij = \frac{m + n + 1}{2} - i

          这是因为:

          • ii0m0 \sim m 递增时,A[i1]\text{A}[i-1] 递增,B[j]\text{B}[j] 递减,所以一定存在一个最大的 ii 满足 A[i1]B[j]A[i-1] \leq B[j]

          • 如果 ii 是最大的,那么说明 i+1i+1 不满足。将 i+1i+1 带入可以得到 A[i]>B[j1]A[i] > B[j-1],也就是 B[j1]<A[i]B[j - 1] < A[i],就和我们进行等价变换前 ii 的性质一致了(甚至还要更强)。

          因此我们可以对 ii[0,m][0, m] 的区间上进行二分搜索,找到最大的满足 A[i1]B[j]A[i-1] \leq B[j]ii 值,就得到了划分的方法。此时,划分前一部分元素中的最大值,以及划分后一部分元素中的最小值,才可能作为就是这两个数组的中位数。

          class Solution {
              public double findMedianSortedArrays(int[] nums1, int[] nums2) {
                  if (nums1.length > nums2.length) {
                      return findMedianSortedArrays(nums2, nums1);
                  }
          
              <span class="hljs-keyword">int</span> m = nums1.length;
              <span class="hljs-keyword">int</span> n = nums2.length;
              <span class="hljs-keyword">int</span> left = <span class="hljs-number">0</span>, right = m, ansi = -<span class="hljs-number">1</span>;
              <span class="hljs-comment">// median1:前一部分的最大值</span>
              <span class="hljs-comment">// median2:后一部分的最小值</span>
              <span class="hljs-keyword">int</span> median1 = <span class="hljs-number">0</span>, median2 = <span class="hljs-number">0</span>;
          
              <span class="hljs-keyword">while</span> (left &lt;= right) {
                  <span class="hljs-comment">// 前一部分包含 nums1[0 .. i-1] 和 nums2[0 .. j-1]</span>
                  <span class="hljs-comment">// 后一部分包含 nums1[i .. m-1] 和 nums2[j .. n-1]</span>
                  <span class="hljs-keyword">int</span> i = (left + right) / <span class="hljs-number">2</span>;
                  <span class="hljs-keyword">int</span> j = (m + n + <span class="hljs-number">1</span>) / <span class="hljs-number">2</span> - i;
          
                  <span class="hljs-comment">// nums_im1, nums_i, nums_jm1, nums_j 分别表示 nums1[i-1], nums1[i], nums2[j-1], nums2[j]</span>
                  <span class="hljs-keyword">int</span> nums_im1 = (i == <span class="hljs-number">0</span> ? Integer.MIN_VALUE : nums1[i - <span class="hljs-number">1</span>]);
                  <span class="hljs-keyword">int</span> nums_i = (i == m ? Integer.MAX_VALUE : nums1[i]);
                  <span class="hljs-keyword">int</span> nums_jm1 = (j == <span class="hljs-number">0</span> ? Integer.MIN_VALUE : nums2[j - <span class="hljs-number">1</span>]);
                  <span class="hljs-keyword">int</span> nums_j = (j == n ? Integer.MAX_VALUE : nums2[j]);
          
                  <span class="hljs-keyword">if</span> (nums_im1 &lt;= nums_j) {
                      ansi = i;
                      median1 = Math.max(nums_im1, nums_jm1);
                      median2 = Math.min(nums_i, nums_j);
                      left = i + <span class="hljs-number">1</span>;
                  }
                  <span class="hljs-keyword">else</span> {
                      right = i - <span class="hljs-number">1</span>;
                  }
              }
          
              <span class="hljs-keyword">return</span> (m + n) % <span class="hljs-number">2</span> == <span class="hljs-number">0</span> ? (median1 + median2) / <span class="hljs-number">2.0</span> : median1;
          }
          

          }


          class Solution {
          public:
          double findMedianSortedArrays(vector<int>& nums1, vector<int>& nums2) {
          if (nums1.size() > nums2.size()) {
          return findMedianSortedArrays(nums2, nums1);
          }

              <span class="hljs-keyword">int</span> m = nums1.size();
              <span class="hljs-keyword">int</span> n = nums2.size();
              <span class="hljs-keyword">int</span> left = <span class="hljs-number">0</span>, right = m, ansi = <span class="hljs-number">-1</span>;
              <span class="hljs-comment">// median1:前一部分的最大值</span>
              <span class="hljs-comment">// median2:后一部分的最小值</span>
              <span class="hljs-keyword">int</span> median1 = <span class="hljs-number">0</span>, median2 = <span class="hljs-number">0</span>;
          
              <span class="hljs-keyword">while</span> (left &lt;= right) {
                  <span class="hljs-comment">// 前一部分包含 nums1[0 .. i-1] 和 nums2[0 .. j-1]</span>
                  <span class="hljs-comment">// 后一部分包含 nums1[i .. m-1] 和 nums2[j .. n-1]</span>
                  <span class="hljs-keyword">int</span> i = (left + right) / <span class="hljs-number">2</span>;
                  <span class="hljs-keyword">int</span> j = (m + n + <span class="hljs-number">1</span>) / <span class="hljs-number">2</span> - i;
          
                  <span class="hljs-comment">// nums_im1, nums_i, nums_jm1, nums_j 分别表示 nums1[i-1], nums1[i], nums2[j-1], nums2[j]</span>
                  <span class="hljs-keyword">int</span> nums_im1 = (i == <span class="hljs-number">0</span> ? INT_MIN : nums1[i - <span class="hljs-number">1</span>]);
                  <span class="hljs-keyword">int</span> nums_i = (i == m ? INT_MAX : nums1[i]);
                  <span class="hljs-keyword">int</span> nums_jm1 = (j == <span class="hljs-number">0</span> ? INT_MIN : nums2[j - <span class="hljs-number">1</span>]);
                  <span class="hljs-keyword">int</span> nums_j = (j == n ? INT_MAX : nums2[j]);
          
                  <span class="hljs-keyword">if</span> (nums_im1 &lt;= nums_j) {
                      ansi = i;
                      median1 = max(nums_im1, nums_jm1);
                      median2 = min(nums_i, nums_j);
                      left = i + <span class="hljs-number">1</span>;
                  }
                  <span class="hljs-keyword">else</span> {
                      right = i - <span class="hljs-number">1</span>;
                  }
              }
          
              <span class="hljs-keyword">return</span> (m + n) % <span class="hljs-number">2</span> == <span class="hljs-number">0</span> ? (median1 + median2) / <span class="hljs-number">2.0</span> : median1;
          }
          

          };


          class Solution:
          def findMedianSortedArrays(self, nums1: List[int], nums2: List[int]) -> float:
          if len(nums1) > len(nums2):
          return self.findMedianSortedArrays(nums2, nums1)

              infinty = <span class="hljs-number">2</span>**<span class="hljs-number">40</span>
              m, n = len(nums1), len(nums2)
              left, right, ansi = <span class="hljs-number">0</span>, m, <span class="hljs-number">-1</span>
              <span class="hljs-comment"># median1:前一部分的最大值</span>
              <span class="hljs-comment"># median2:后一部分的最小值</span>
              median1, median2 = <span class="hljs-number">0</span>, <span class="hljs-number">0</span>
          
              <span class="hljs-keyword">while</span> left &lt;= right:
                  <span class="hljs-comment"># 前一部分包含 nums1[0 .. i-1] 和 nums2[0 .. j-1]</span>
                  <span class="hljs-comment"># // 后一部分包含 nums1[i .. m-1] 和 nums2[j .. n-1]</span>
                  i = (left + right) // <span class="hljs-number">2</span>
                  j = (m + n + <span class="hljs-number">1</span>) // <span class="hljs-number">2</span> - i
          
                  <span class="hljs-comment"># nums_im1, nums_i, nums_jm1, nums_j 分别表示 nums1[i-1], nums1[i], nums2[j-1], nums2[j]</span>
                  nums_im1 = (-infinty <span class="hljs-keyword">if</span> i == <span class="hljs-number">0</span> <span class="hljs-keyword">else</span> nums1[i - <span class="hljs-number">1</span>])
                  nums_i = (infinty <span class="hljs-keyword">if</span> i == m <span class="hljs-keyword">else</span> nums1[i])
                  nums_jm1 = (-infinty <span class="hljs-keyword">if</span> j == <span class="hljs-number">0</span> <span class="hljs-keyword">else</span> nums2[j - <span class="hljs-number">1</span>])
                  nums_j = (infinty <span class="hljs-keyword">if</span> j == n <span class="hljs-keyword">else</span> nums2[j])
          
                  <span class="hljs-keyword">if</span> nums_im1 &lt;= nums_j:
                      ansi = i
                      median1, median2 = max(nums_im1, nums_jm1), min(nums_i, nums_j)
                      left = i + <span class="hljs-number">1</span>
                  <span class="hljs-keyword">else</span>:
                      right = i - <span class="hljs-number">1</span>
          
              <span class="hljs-keyword">return</span> (median1 + median2) / <span class="hljs-number">2</span> <span class="hljs-keyword">if</span> (m + n) % <span class="hljs-number">2</span> == <span class="hljs-number">0</span> <span class="hljs-keyword">else</span> median1
          


          func findMedianSortedArrays(nums1 []int, nums2 []int) float64 {
          if len(nums1) > len(nums2) {
          return findMedianSortedArrays(nums2, nums1)
          }
          m, n := len(nums1), len(nums2)
          left, right := 0, m
          median1, median2 := 0, 0
          for left <= right {
          i := (left + right) / 2
          j := (m + n + 1) / 2 - i
          nums_im1 := math.MinInt32
          if i != 0 {
          nums_im1 = nums1[i-1]
          }
          nums_i := math.MaxInt32
          if i != m {
          nums_i = nums1[i]
          }
          nums_jm1 := math.MinInt32
          if j != 0 {
          nums_jm1 = nums2[j-1]
          }
          nums_j := math.MaxInt32
          if j != n {
          nums_j = nums2[j]
          }
          if nums_im1 <= nums_j {
          median1 = max(nums_im1, nums_jm1)
          median2 = min(nums_i, nums_j)
          left = i + 1
          } else {
          right = i - 1
          }
          }
          if (m + n) % 2 == 0 {
          return float64(median1 + median2) / 2.0
          }
          return float64(median1)
          }

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

          func min(x, y int) int {
          if x < y {
          return x
          }
          return y
          }

          复杂度分析

          • 时间复杂度:O(logmin(m,n)))O(\log\min(m,n))),其中 mmnn 分别是数组 nums1\text{nums1}nums2\text{nums2} 的长度。查找的区间是 [0,m][0, m],而该区间的长度在每次循环之后都会减少为原来的一半。所以,只需要执行 logm\log m 次循环。由于每次循环中的操作次数是常数,所以时间复杂度为 O(logm)O(\log m)。由于我们可能需要交换 nums1\text{nums1}nums2\text{nums2} 使得 mnm \leq n,因此时间复杂度是 O(logmin(m,n)))O(\log\min(m,n)))

          • 空间复杂度:O(1)O(1)

          posted @ 2020-06-04 17:20  刷题之路1  阅读(270)  评论(0)    收藏  举报