LC 480. 滑动窗口中位数

题目描述

这是 LeetCode 上的 480. 滑动窗口中位数 ,难度为 困难

中位数是有序序列最中间的那个数。如果序列的长度是偶数,则没有最中间的数;此时中位数是最中间的两个数的平均数。

例如:

  • $[2,3,4]$,中位数是 $3$
  • $[2,3]$,中位数是 $(2 + 3) / 2 = 2.5$

给你一个数组 $nums$,有一个长度为 $k$ 的窗口从最左端滑动到最右端。

窗口中有 $k$ 个数,每次窗口向右移动 $1$ 位。

你的任务是找出每次窗口移动后得到的新窗口中元素的中位数,并输出由它们组成的数组。

示例:

1
2
3
4
5
6
7
8
9
10
11
12
给出 nums = [1,3,-1,-3,5,3,6,7],以及 k = 3。

窗口位置 中位数
--------------- -----
[1 3 -1] -3 5 3 6 7 1
1 [3 -1 -3] 5 3 6 7 -1
1 3 [-1 -3 5] 3 6 7 -1
1 3 -1 [-3 5 3] 6 7 3
1 3 -1 -3 [5 3 6] 7 5
1 3 -1 -3 5 [3 6 7] 6

 因此,返回该滑动窗口的中位数数组 [1,-1,-1,3,5,6]。

提示:

  • 你可以假设 $k$ 始终有效,即:$k$ 始终小于等于输入的非空数组的元素个数。
  • 与真实值误差在 $10 ^ {-5}$ 以内的答案将被视作正确答案。

朴素解法

一个直观的做法是:对每个滑动窗口的数进行排序,获取排序好的数组中的第 $\frac{k}{2}$ 和 $\frac{k - 1}{2}$ 个数(避免奇偶数讨论),计算中位数。

我们大概分析就知道这个做法至少 $O(n k)$ 的,算上排序的话应该是 $O(n (k + k\log{k}))$。

比较无奈的是,这道题的中文说明中没有给出数据范围。我们无法根据判断这样的做法会不会超时。

朴素做法通常是优化的开始,所以还是提供一下朴素做法的代码。

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
class Solution {
public double[] medianSlidingWindow(int[] nums, int k) {
int n = nums.length;
int cnt = n - k + 1;
double[] ans = new double[cnt];
int[] t = new int[k];
for (int l = 0, r = l + k - 1; r < n; l++, r++) {
for (int i = l; i <= r; i++) t[i - l] = nums[i];
Arrays.sort(t);
ans[l] = (t[k / 2] / 2.0) + (t[(k - 1) / 2] / 2.0);
}
return ans;
}
}

  • 时间复杂度:最多有 n 个窗口需要滑动计算。每个窗口,需要先插入数据,复杂度为 $O(k)$,插入后需要排序,复杂度为 $O(k\log{k})$。整体复杂度为 $O(n * (k + k\log{k}))$
  • 空间复杂度:使用了长度为 k 的临时数组。复杂度为 $O(k)$

优先队列(堆)

从朴素解法中我们可以发现,其实我们需要的就是滑动窗口中的第 $\frac{k}{2}$ 小的值和第 $\frac{k - 1}{2}$ 小的值。

我们知道滑动窗口求最值的问题,可以使用优先队列来做。

但这里我们求的是第 $k$ 小的数,而且是需要两个值。还能不能使用优先队列来做呢?

我们可以维护两个堆:

  • 一个大根堆维护着滑动窗口中一半较小的值(此时堆顶元素为滑动窗口中的第 $\frac{k - 1}{2}$ 小的值)
  • 一个小根堆维护着滑动窗口中一半较大的值(此时堆顶元素为滑动窗口中的第 $\frac{k}{2}$ 小的值)

滑动窗口的中位数就是两个堆的堆顶元素的平均值。

实现细节:

  1. 初始化时,先让 k 个元素直接入 right,再从 right 中倒出 $\frac{k}{2}$ 个到 left 中。这时候可以根据 leftright 得到第一个滑动窗口的中位值。

  2. 开始滑动窗口,每次滑动都有一个待添加和待移除的数:

    2.1 根据与右堆的堆顶元素比较,决定是插入哪个堆和从哪个堆移除

    2.2 之后调整两堆的大小(确保只会出现 left.size() == right.size()right.size() - left.size() == 1,对应了窗口长度为偶数或者奇数的情况)

    2.3 根据 left 堆 和 right 堆得到当前滑动窗口的中位值

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
class Solution {
public double[] medianSlidingWindow(int[] nums, int k) {
int n = nums.length;
int cnt = n - k + 1;
double[] ans = new double[cnt];
// 如果是奇数滑动窗口,让 right 的数量比 left 多一个
PriorityQueue<Integer> left = new PriorityQueue<>((a,b)->Integer.compare(b,a)); // 滑动窗口的左半部分
PriorityQueue<Integer> right = new PriorityQueue<>((a,b)->Integer.compare(a,b)); // 滑动窗口的右半部分
for (int i = 0; i < k; i++) right.add(nums[i]);
for (int i = 0; i < k / 2; i++) left.add(right.poll());
ans[0] = getMid(left, right);
for (int i = k; i < n; i++) {
// 人为确保了 right 会比 left 多,因此,删除和添加都与 right 比较(left 可能为空)
int add = nums[i], del = nums[i - k];
if (add >= right.peek()) {
right.add(add);
} else {
left.add(add);
}
if (del >= right.peek()) {
right.remove(del);
} else {
left.remove(del);
}
adjust(left, right);
ans[i - k + 1] = getMid(left, right);
}
return ans;
}
void adjust(PriorityQueue<Integer> left, PriorityQueue<Integer> right) {
while (left.size() > right.size()) right.add(left.poll());
while (right.size() - left.size() > 1) left.add(right.poll());
}
double getMid(PriorityQueue<Integer> left, PriorityQueue<Integer> right) {
if (left.size() == right.size()) {
return (left.peek() / 2.0) + (right.peek() / 2.0);
} else {
return right.peek() * 1.0;
}
}
}

  • 时间复杂度:调整过程中堆大小最大为 k,堆操作中的指定元素删除复杂度为 $O(k)$;窗口数量最多为 n。整体复杂度为 $O(n * k)$
  • 空间复杂度:最多有 n 个元素在堆内。复杂度为 $O(n)$

答疑

以下是针对一些具有代表性的问题进行的集中答疑:

  • 为什么 new PriorityQueue<>((x,y)->(y-x)) 的写法会有某些案例无法通过?和 new PriorityQueue<>((x,y)->Integer.compare(y,x)) 写法有何区别?

    (x,y)->(y-x) 的写法逻辑没有错,AC 不了是因为 int 溢出。

    在 Java 中 Integer.compare 的实现是 (x < y) ? -1 : ((x == y) ? 0 : 1)。只是单纯的比较,不涉及运算,所以不存在溢出风险。

    而直接使用 y - x,当 y = Integer.MAX_VALUE, x = Integer.MIN_VALUE 时,到导致溢出,返回的是 负数 ,而不是逻辑期望的 正数

同样具有溢出问题的还有计算第 $\frac{k}{2}$ 小的数和第 $\frac{k - 1}{2}$ 小的数的平均值时。

因此题解中使用的是 (a / 2.0) + (b / 2.0) 的形式,而不是采用 (a + b) / 2.0 的形式。后者有相加溢出的风险。


复杂度说明

JDK 中 PriorityQueueremove(Object o) 实现是先调用 indexOf(Object o) 方法进行线性扫描找到下标(复杂度为 $O(n)$),之后再调用 removeAt(int i) 进行删除(复杂度为 $O(\log{n})$)。

对于本题而言,如果需要实现 $O(\log{n})$ 的 remove(Object o), 只能通过引入其他数据结构(如哈希表)来实现快速查找元素在对堆数组中的下标。

对于本题,可以使用元素在原数组中的下标作为 key,在堆数组中的真实下标作为 val。

通过哈希表可以 $O(1)$ 的复杂度找到下标,之后的删除只需要算堆调整的复杂度即可(最多 down 一遍,up 一遍,复杂度为 $O(\log{n})$)。

至于 JDK 没有这样做的原因,猜测是因为基本类型的包装类型存在小数缓存机制,导致无法很好的使用哈希表来对应一个插入元素的下标。

举个🌰,我们调用三次 add(10), 会有 $3$ 个 $10$ 在堆内,但是由于小数(默认范围为 $[-128,127]$)包装类型存在缓存机制,使用哈希表继续记录的话,只会有 { Integer.valueOf(10) : 移动过程中最后一次访问的数组下标 } 这样一条记录(add 进去的 $10$ 均为同一对象)。这时候删除一个 $10$ 之后,哈希表无法正确指导我们找到下一个 $10$ 的位置。


最后

这是我们「刷穿 LeetCode」系列文章的第 No.480 篇,系列开始于 2021/01/01,截止于起始日 LeetCode 上共有 1916 道题目,部分是有锁题,我们将先把所有不带锁的题目刷完。

在这个系列文章里面,除了讲解解题思路以外,还会尽可能给出最为简洁的代码。如果涉及通解还会相应的代码模板。

为了方便各位同学能够电脑上进行调试和提交代码,我建立了相关的仓库:https://github.com/SharingSource/LogicStack-LeetCode

在仓库地址里,你可以看到系列文章的题解链接、系列文章的相应代码、LeetCode 原题链接和其他优选题解。


本博客所有文章除特别声明外,均采用 CC BY-SA 4.0 协议 ,转载请注明出处!