LeetCode 周赛 342（2023/04/23）容斥原理、

大家好，我是小彭。

前天刚举办 2023 年力扣杯个人 SOLO 赛，昨天周赛就出了一场 Easy - Easy - Medium - Medium 的水场，不得不说 LeetCode 是懂礼数的 😁。

接下来，请你跟着小彭的思路，一步步将问题做难，再将问题做简单。

往期回顾：LeetCode 单周赛 341 · 难度上来了，图论的问题好多啊！

LeetCode 周赛 342 概览

Q1. 计算列车到站时间（Easy）

简单模拟题，不多讲解。

Q2. 倍数求和（Easy）

题解 1：暴力解法 时间复杂度

题解 2：分析数据特征后发现数据存在等差数列性质，我们利用容斥原理和等差数列求和公式，可以把优化到 时间复杂度

Q3. 滑动子数组的美丽值（Medium）

题解 1：在滑动窗口的基础上，结合快速选择查找滑动窗口中第 x 小的元素，时间复杂度是

题解 2：分析数据特征后发现题目的值域非常小，我们可以用计数排序代替快速选择，时间复杂度为

Q4. 使数组所有元素变成 1 的最少操作次数（Medium）

在问题分析后我们将原问题抽象为 “寻找 GCB 为 1 的最短子数组”，关联相似的 “和为 k 的最短子数组” 问题，我们有从暴力 → 有序集合 → 单调性优化的解法：

题解 1：暴力 时间复杂度

题解 2：有序集合 时间复杂度

题解 3：单调性优化 时间复杂度

---

Q1. 计算列车到站时间（Easy）

题目地址

https://leetcode.cn/problems/calculate-delayed-arrival-time/

题目描述

给你一个正整数 arrivalTime 表示列车正点到站的时间（单位：小时），另给你一个正整数 delayedTime 表示列车延误的小时数。

返回列车实际到站的时间。

注意，该问题中的时间采用 24 小时制。

示例 1：

输入：arrivalTime = 15, delayedTime = 5
输出：20
解释：列车正点到站时间是 15:00 ，延误 5 小时，所以列车实际到站的时间是 15 + 5 = 20（20:00）。

示例 2：

输入：arrivalTime = 13, delayedTime = 11
输出：0
解释：列车正点到站时间是 13:00 ，延误 11 小时，所以列车实际到站的时间是 13 + 11 = 24（在 24 小时制中表示为 00:00 ，所以返回 0）。

提示：

• 1 <= arrivaltime < 24
• 1 <= delayedTime <= 24

题解（模拟）

简单模拟题，按题意实现即可。

class Solution {
    fun findDelayedArrivalTime(arrivalTime: Int, delayedTime: Int): Int {
        return (arrivalTime + delayedTime) % 24
    }
}

复杂度分析：

• 时间复杂度：
• 空间复杂度：

---

Q2. 倍数求和（Easy）

题目地址

https://leetcode.cn/problems/sum-multiples/

题目描述

给你一个正整数 n ，请你计算在 [1，n] 范围内能被 3、5、7 整除的所有整数之和。

返回一个整数，用于表示给定范围内所有满足约束条件的数字之和。

示例 1：

输入：n = 7
输出：21
解释：在[1, 7] 范围内能被 3、5、7 整除的所有整数分别是 3、5、6、7 。数字之和为21 。

示例 2：

输入：n = 10
输出：40
解释：在[1, 10] 范围内能被 3、5、7 整除的所有整数分别是 3、5、6、7、9、10 。数字之和为40 。

示例 3：

输入：n = 9
输出：30
解释：在[1, 9] 范围内能被 3、5、7 整除的所有整数分别是 3、5、6、7、9 。数字之和为30 。

提示：

• 1 <= n <= 103

预备知识 - 容斥原理

定义集合 A 表示能够被 3 整除的数，定义集合 B 表示能够被 5 整除的数，定义集合 C 表示能够被 7 整除的数。如果把这 3 个集合直接加起来，会多出来一些元素重复统计了，因此需要扣除 A ∩ B，A ∩ C 和 B ∩ C ，但是又有一小部分元素多扣了，反而再需要加上 A ∩ B ∩ C。这就是 容斥原理：

其中：

• A ∪ B ∪ C 表示能够被 3 或 5 或 7 整除的数，也就是原问题的解；
• A ∩ B 表示能够同时被 3 和 5 整除的数；
• A ∩ C 表示能够同时被 3 和 7 整除的数；
• B ∩ C 表示能够同时被 5 和 7 整除的数。

预备知识 - 等差数列求和

• 等差数列求和公式：(首项 + 尾项) * 项数 / 2

题解一（暴力）

先思考暴力解法：

算法：枚举每个数，检查该数字是否为 3 / 5 / 7 的倍数，是则计入结果中。

class Solution {
    fun sumOfMultiples(n: Int): Int {
        var ret = 0
        for (i in 1 .. n) {
            if(i % 3 == 0 || i % 5 == 0 || i % 7 == 0) ret += i
        }
        return ret
    }
}

复杂度分析：

• 时间复杂度： 其中 n 为 nums 数组的长度，每个元素最多访问 1 次；
• 空间复杂度：

题解二（容斥原理 + 等差数列求和公式）

暴力解法是否有优化空间呢，先分析重复计算：

• 要点 1：可以观察到 [1, n] 范围中的目标数是存在关联的，以 3 的倍数为例，3、6、9、12 是以 3 为等差的等差数列，而等差数列的和可以使用公式计算。数字 m 在 [1, n] 范围内中的倍数为 n / m 个，可以使用等差数列求和公式以 O(1) 算出这部分元素之和；
• 要点 2：结合容斥原理，可以在 O(1) 时间复杂度求出原问题。那么能够同时被 3 和 5 整除的等差数列如何计算呢？其实就是计算 15 的倍数。同理能够同时被 3 和 5 和 7 整除的等差数列就是 105 的倍数。

至此，结合容斥原理模拟即可：

class Solution {
    fun sumOfMultiples(n: Int): Int {
        return sum(n, 3) + sum(n, 5) + sum(n, 7) - sum(n, 15) - sum(n, 21) - sum(n, 35) + sum(n, 105)
    }

    private fun sum(n:Int, k:Int) :Int {
        // 等差数列求和公式：(首项 + 尾项) * 项数 / 2
        return (k + (n / k * k)) * (n / k) / 2
    }
}

复杂度分析：

• 时间复杂度：
• 空间复杂度：

---

Q3. 滑动子数组的美丽值（Medium）

题目地址

https://leetcode.cn/problems/sliding-subarray-beauty/description/

题目描述

给你一个长度为 n 的整数数组 nums ，请你求出每个长度为 k 的子数组的 美丽值 。

一个子数组的 美丽值 定义为：如果子数组中第 x 小整数 是 负数 ，那么美丽值为第 x 小的数，否则美丽值为 0 。

请你返回一个包含 n - k + 1 个整数的数组，依次** 表示数组中从第一个下标开始，每个长度为 k 的子数组的 美丽值 。

• 子数组指的是数组中一段连续 非空 的元素序列。

示例 1：

输入：nums = [1,-1,-3,-2,3], k = 3, x = 2
输出：[-1,-2,-2]
解释：总共有 3 个 k = 3 的子数组。
第一个子数组是[1, -1, -3] ，第二小的数是负数 -1 。
第二个子数组是[-1, -3, -2] ，第二小的数是负数 -2 。
第三个子数组是[-3, -2, 3] ，第二小的数是负数 -2 。

示例 2：

输入：nums = [-1,-2,-3,-4,-5], k = 2, x = 2
输出：[-1,-2,-3,-4]
解释：总共有 4 个 k = 2 的子数组。
[-1, -2] 中第二小的数是负数 -1 。[-2, -3] 中第二小的数是负数 -2 。[-3, -4] 中第二小的数是负数 -3 。[-4, -5] 中第二小的数是负数 -4 。

示例 3：

输入：nums = [-3,1,2,-3,0,-3], k = 2, x = 1
输出：[-3,0,-3,-3,-3]
解释：总共有 5 个 k = 2 的子数组。
[-3, 1] 中最小的数是负数 -3 。[1, 2] 中最小的数不是负数，所以美丽值为 0 。[2, -3] 中最小的数是负数 -3 。[-3, 0] 中最小的数是负数 -3 。[0, -3] 中最小的数是负数 -3 。

提示：

• n == nums.length
• 1 <= n <= 105
• 1 <= k <= n
• 1 <= x <= k
• 50 <= nums[i] <= 50

预备知识

求出每个长度为 k 的子数组的美丽值，容易想到可以用滑动窗口。

伪代码为：

// 伪代码
for (i in 0 until n) {
    // 进入窗口
    list.add(i)
    // 离开窗口
    if (i >= k) list.remove(i - k)
    if (i >= k - 1) {
        // 计算窗口答案
    }
}

题解一（滑动窗口 + 快速选择 · 超出时间限制）

在滑动窗口的基础上，使用快速选择查找窗口中第 x 小的数：

class Solution {

    private val random = Random(0)

    fun getSubarrayBeauty(nums: IntArray, k: Int, x: Int): IntArray {
        val n = nums.size
        val ret = LinkedList<Int>()
        val list = ArrayList<Int>()
        for (i in 0 until n) {
            // 进入窗口
            list.add(i)
            // 离开窗口
            if (i >= k) list.remove(i - k)
            if (i >= k - 1) {
                // 计算窗口答案
                quickSelect(nums, list, x)
                val num = nums[list[x - 1]]
                ret.add(if (num < 0) num else 0)
            }
        }
        return ret.toIntArray()
    }

    private fun quickSelect(nums: IntArray, list: ArrayList<Int>, x: Int) {
        val target = x - 1
        var left = 0
        var right = list.size - 1
        while (left < right) {
            val pivot = partition(nums, list, left, right)
            if (pivot == target) {
                return
            } else if (pivot < target) {
                left = pivot + 1
            } else {
                right = pivot - 1
            }
        }
    }

    private fun partition(nums: IntArray, list: ArrayList<Int>, left: Int, right: Int): Int {
        val random = random.nextInt(right - left + 1) + left
        list.swap(random, right)
        var p = left
        for (i in left until right) {
            if (nums[list[i]] < nums[list[right]]) list.swap(i, p++)
        }
        list.swap(p, right)
        return p
    }

    private fun ArrayList<Int>.swap(i: Int, j: Int) {
        val temp = this[i]
        this[i] = this[j]
        this[j] = temp
    }
}

复杂度分析：

• 时间复杂度： 其中 n 是 nums 数组的长度，单次窗口快速选择的时间复杂度是 ；
• 空间复杂度： 滑动窗口空间。

题解二（滑动窗口 + 计数排序）

注意到题目的值域非常小，能否利用起来：

我们可以用计数排序代替快速选择，用 cnts 计数数组计算窗口内每个元素的出现次数，再根据计数数组计算出第 x 小的数：

class Solution {

    private val random = Random(0)

    fun getSubarrayBeauty(nums: IntArray, k: Int, x: Int): IntArray {
        val n = nums.size
        val OFFSET = 50
        val cnts = IntArray(OFFSET * 2 + 1)
        val ret = IntArray(n - k + 1)
        outer@ for (i in 0 until n) {
            // 进入窗口
            cnts[OFFSET + nums[i]]++
            // 离开窗口
            if (i >= k) cnts[OFFSET + nums[i - k]]--
            if (i >= k - 1) {
                // 计算窗口美丽值
                var count = x
                // for (num in -OFFSET .. -1) {
                for (num in -OFFSET .. OFFSET) {
                    count -= cnts[num + 50]
                    if (count <= 0) {
                        // 找到第 x 小的数
                        // ret[i - k + 1] = num
                        ret[i - k + 1] = if(num < 0) num else 0
                        continue@outer
                    }
                }
            }
        }
        return ret
    }
}

另外，由于题目要求美丽值是负数，所以在计算窗口美丽值时，我们只需要枚举 [-50, -1] 的元素值。

复杂度分析：

• 时间复杂度： 其中 n 是 nums 数组的长度，U 是值域大小 101。每次滑动窗口求第 x 小的元素时间是 ；
• 空间复杂度： 计数数组空间。

---

Q4. 使数组所有元素变成 1 的最少操作次数（Medium）

题目地址

https://leetcode.cn/problems/minimum-number-of-operations-to-make-all-array-elements-equal-to-1/description/

题目描述

给你一个下标从 0 开始的 正 整数数组 nums 。你可以对数组执行以下操作 任意 次：

• 选择一个满足 0 <= i < n - 1 的下标 i ，将 nums[i] 或者 nums[i+1] 两者之一替换成它们的最大公约数。

请你返回使数组 nums 中所有元素都等于 1 的 最少 操作次数。如果无法让数组全部变成 1 ，请你返回 -1 。

两个正整数的最大公约数指的是能整除这两个数的最大正整数。

示例 1：

输入：nums = [2,6,3,4]
输出：4
解释：我们可以执行以下操作：
- 选择下标 i = 2 ，将 nums[2] 替换为 gcd(3,4) = 1 ，得到 nums = [2,6,1,4] 。
- 选择下标 i = 1 ，将 nums[1] 替换为 gcd(6,1) = 1 ，得到 nums = [2,1,1,4] 。
- 选择下标 i = 0 ，将 nums[0] 替换为 gcd(2,1) = 1 ，得到 nums = [1,1,1,4] 。
- 选择下标 i = 2 ，将 nums[3] 替换为 gcd(1,4) = 1 ，得到 nums = [1,1,1,1] 。

示例 2：

输入：nums = [2,10,6,14]
输出：-1
解释：无法将所有元素都变成 1 。

提示：

• 2 <= nums.length <= 50
• 1 <= nums[i] <= 106

问题分析

分析目标结果：

使得数组中所有元素都变成 1 的最少操作次数。

分析题目示例：

• 由于在每次操作中最多只能将一个数字修改为最大公约数，那么将 1 个元素操作为 “1” 的最小操作次数（如果可行）不会低于 1 次，将 n 个大于 1 的元素操作为 “1” 的最小次数不会低于 n 次，例如样例 [2,6,1,4]。
• 如果数组中至少存在 1 个 “1” 时，我们只需要将每个 “1” 与相邻的 “非 1” 元素组合操作，就能将所有元素，例如样例 [2,6,1,4]。这说明，问题的最小操作次数正好就是数组中不是 “1” 的个数。
• 如果数组中不存在 “1”，需要先操作出原始的 “1”：
  • 如果数组中所有元素的最大公约数大于 1，那么无论如何也无法操作出数字 1，例如样例 [2, 10, 6, 14]；
  • 否则，我们总可以操作 x 次获得原始 “1”，那么问题就等于 count + n - 1；

至此，程序整体框架确定。伪代码为：

if (所有元素的最大公约数 > 1) return -1
if (1 的个数 > 0) return n - (1 的个数)
操作 count 次得到原始的 “1”
return count + n - 1

接下来，我们需要思考如何计算出操作出原始 “1” 的最小次数：

回归到原问题操作，我们在每次操作中可以将一个数修改为最大公约数，那么对于连续的一段子数组（长度为 subSize），我们总可以用 subSize - 1 次操作将其中一个数变为整个子数组的最大公约数。如果这个最大公约数是 1，那么操作次数正好是 subSize - 1，反之无法操作出 1。

至此，可以想出暴力解法：

题解一（暴力枚举子数组）

在暴力解法中，我们枚举所有子数组，记录出所有子数组操作出原始 “1” 的最少操作次数。

class Solution {
    fun minOperations(nums: IntArray): Int {
        val n = nums.size
        // 1 的个数
        var cnt1 = 0
        var gcbAll = 0
        for (x in nums) {
            gcbAll = gcb(gcbAll, x)
            if (x == 1) cnt1++
        }
        // 所有数的最大公约数大于 1
        if (gcbAll > 1) return -1
        // 1 的个数大于 0
        if (cnt1 > 0) return n - cnt1

        // 操作出原始 “1” 的最小次数
        var minCount = n
        // 枚举子数组
        for (i in 0 until n) {
            var gcb = 0
            for (j in i until n) {
                gcb = gcb(gcb, nums[j])
                if (gcb == 1) {
                    minCount = Math.min(minCount, j - i /* 子数组长度 - 1 */)
                    break // 继续枚举 i 为起点的子数组不会得到更优解
                }
            }
        }
        return minCount + n - 1
    }

    // 求 x 和 y 的最大公约数（辗转相除法）
    private fun gcb(x: Int, y: Int): Int {
        var a = x
        var b = y
        while (b != 0) {
            val temp = a % b
            a = b
            b = temp
        }
        return a
    }
}

复杂度分析：

• 时间复杂度： 其中 n 是 nums 数组的长度，U 是数组元素的最大值。单次 GCD 计算的时间复杂度是 ，乍看起来算法整体的时间复杂度是 ，其实不对。因为在每层循环中，每次 GCD 计算并不是独立的，而是贯穿整个内层循环的，所以 GCD 的总时间取决于数据的最大值 U，在辗转相除中取余的次数也取决于 U。
• 空间复杂度： 不考虑结果数组，仅使用常量级别空间。

题解一的复杂度是平方级别的，如果放大题目数据量到 10^5 要怎么做？

问题抽象

在分析暴力解法的重复计算之前，我先向你抛出一个 “题外话”：

请你回答：“给定一个整数数组 nums 和目标和 k，如何求和为 k 的最短子数组？”

• 解法 1：暴力枚举所有子数组，记录出所有子数组和为 k 的最短子数组长度（这与题解一暴力枚举子数组求操作出原始 “1” 的最少操作次数类似）；
• 解法 2：我们从左向右线性遍历，并维护以 num[j] 为右端点的前缀和映射表 <preSum to index>。在此基础上，我们将当前位置 nums[i] 的前缀和与前缀和映射表中的每个元素取差值，就可以快速地获得以 num[i] 为右端点所有子数组的和。另外，由于我们是从左向右遍历的，所以前缀和映射表记录的索引正好是可以构造最短子数组的索引，子数组长度为 i - j + 1（当然，我们可以直接 O(1) 查询目标前缀和出现时的索引，而不需要真的用前缀和映射表的每个元素取差值）。

注：这个 “题外话” 与 LeetCode 560. 和为 K 的子数组 类似，如果你不熟悉可以先做做看。

那么，这个 “题外话” 与今天这道题有什么关系：

根据 GCB 运算的性质，当我们以 nums[i] 为左端点，不断向右扩展子数组的右端点时，我们的目标是求 “GCB 为 1 的子数组” 对吧。与 “求和为 k 的最短子数组” 类似，我们可以维护以 nums[j] 为左端点的 GCB 映射表 <gcb to 左端点 index>。在此基础上，我们将当前位置 nums[i] 与 GCB 映射表中的每个元素取 GCB，就可以快速的获得以 nums[i] 为右端点的所有子数组的 GCB。

那听起来这个算法依然是 O(n^2)？不对。

原因在题解一的时间复杂度分析中讲到了，因为每次 GCD 计算并不是独立的，而是贯穿整个循环的，GCB 映射表的大小取决于数据的最大值 U，而不是数据量，最多有 logU 种 GCB。因此优化后算法的时间复杂度是 O(n·lgU)，但增加了空间复杂度为 O(lgU)。

示意图

题解二（有序集合）

至此，在题解一的基础上修改 “枚举子数组计算操作出原始 “1” 的最小次数” 不分代码即可：

class Solution {
    fun minOperations(nums: IntArray): Int {
        // 略...

        // 计算操作出原始 “1” 的最小次数
        var minCount = n
        // gcb 散列表 <gcd to 左端点 index>
        var gcbMap = TreeMap<Int, Int>()
        // 枚举子数组
        for (i in 0 until n) {
            val newGcbMap = TreeMap<Int, Int>()
            // 枚举 gcb 映射表
            for ((gcb, index) in gcbMap) {
                newGcbMap[gcb(gcb, nums[i])] = index
            }
            newGcbMap[nums[i]] = i
            // 检查最小的 gcb 是否为 1
            val minEntry = newGcbMap.firstEntry()
            if (1 == minEntry.key) {
                minCount = Math.min(minCount, i - minEntry.value /* 子数组长度 - 1 */)
            }
            gcbMap = newGcbMap
        }
        return minCount + n - 1
    }

    // 求 x 和 y 的最大公约数
    private fun gcb(x: Int, y: Int): Int {
        // 略...
    }
}

复杂度分析：

• 时间复杂度： 由于使用了有序集合，所以每一轮迭代中要算上排序时间 ；
• 空间复杂度： GCB 映射表空间。

题解三（单调性优化）

思路参考：灵茶山艾府的题解

题解二的时间复杂度比我们分析的复杂度略要一些，如何寻找优化空间？

继续分析 GCB 的数据特征，可以发现：当我们从左向右遍历时，随着子数组的长度增大，子数组的 GCB 要么不变，要么变小，存在 单调性。 所以，我们并不需要维护有序集合，GCB 列表中最靠前的元素一定是最小的 GCB。

class Solution {
    fun minOperations(nums: IntArray): Int {
        // 略...

        // 计算操作出原始 “1” 的最小次数
        var minCount = n
        // gcb 列表 <gcd to 左端点 index>
        var gcbs = ArrayList<IntArray>()
        // 枚举子数组
        for (i in 0 until n) {
            val newGcbs = ArrayList<IntArray>()
            // 枚举 gcb 列表
            for (element in gcbs) {
                val gcb = gcb(element[0], nums[i])
                if (newGcbs.isEmpty() || newGcbs[newGcbs.size - 1][0] != gcb) {
                    // 增加 GCB
                    newGcbs.add(intArrayOf(gcb, element[1]))
                } else {
                    // 原地去重
                    newGcbs[newGcbs.size - 1][1] = element[1]
                }
            }
            newGcbs.add(intArrayOf(nums[i], i))
            // 检查最小的 gcb 是否为 1
            val minEntry = newGcbs[0]
            if (1 == minEntry[0]) {
                minCount = Math.min(minCount, i - minEntry[1] /* 子数组长度 - 1 */)
            }
            gcbs = newGcbs

        }
        return minCount + n - 1
    }

    // 求 x 和 y 的最大公约数
    private fun gcb(x: Int, y: Int): Int {
        // 略...
    }
}

复杂度分析：

• 时间复杂度：
• 空间复杂度：

相似题目：

• 560. 和为 K 的子数组
• 898. 子数组按位或操作
• 1521. 找到最接近目标值的函数值

---