LeetCode 周赛上分之旅 #43 计算机科学本质上是数学吗?

leetcode,计算机科学,本质,数学 · 浏览次数 : 235

小编点评

**算法思想:** 1. **枚举所有字母的子序列:**使用循环枚举所有可能的子序列,并计算每个子序列的方案数。 2. **排序子序列的方案数:**根据子序列中字符的顺序和频率,对方案数进行排序。 3. **选择最大方案数:**从所有子序列中选择方案数最大的一种。 **代码:** ```python class Solution: def countKSubsequencesWithMaxBeauty(self, s: str, k: int) -> int: # MOD 定义 MOD = 1000000007 # 计数器 cnt = Counter(Counter(s).values()) # 排序子序列的方案数 ans = 1 for c, num in sorted(cnt.items(), reverse=True): if num >= k: return ans * pow(c, k, MOD) * comb(num, k) % MOD ans *= pow(c, num, MOD) k -= num # 返回结果 return ans ``` **分析:** * 时间复杂度:$O(n + m)$,其中 $n$ 是字符串长度,$m$ 是字符种类。 * 空间复杂度:$O(m)$。 * 特别情况下,如果字符串只包含一个字符,则方案数为 1。 * 算法利用了排序和乘法原理来高效地计算方案数。

正文

⭐️ 本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 [彭旭锐] 和 BaguTree Pro 知识星球提问。

学习数据结构与算法的关键在于掌握问题背后的算法思维框架,你的思考越抽象,它能覆盖的问题域就越广,理解难度也更复杂。在这个专栏里,小彭与你分享每场 LeetCode 周赛的解题报告,一起体会上分之旅。

本文是 LeetCode 上分之旅系列的第 43 篇文章,往期回顾请移步到文章末尾~

LeetCode 双周赛 112

T1. 判断通过操作能否让字符串相等 I(Easy)

  • 标签:模拟

T2. 判断通过操作能否让字符串相等 II(Medium)

  • 标签:模拟、计数、排序

T3. 几乎唯一子数组的最大和(Medium)

  • 标签:滑动窗口、计数

T4. 统计一个字符串的 k 子序列美丽值最大的数目(Hard)

  • 标签:枚举、贪心、排序、乘法原理、组合数


T1. 判断通过操作能否让字符串相等 I(Easy)

https://leetcode.cn/problems/check-if-strings-can-be-made-equal-with-operations-i/

题解(模拟)

因为只能交换距离偶数倍的位置,因此相当于比较两个字符串相同奇偶性下标上的元素是否相等。

  • 写法 1:基于散列表
class Solution {
    fun canBeEqual(s1: String, s2: String): Boolean {
        return setOf(s1[0], s1[2]) == setOf(s2[0], s2[2]) && setOf(s1[1], s1[3]) == setOf(s2[1], s2[3])
    }
}
  • 写法 2:基于字符串
class Solution:
    def checkStrings(self, s1: str, s2: str) -> bool:
        return sorted(s1[0::2]) == sorted(s2[0::2]) and sorted(s1[1::2]) == sorted(s2[1::2])

复杂度分析:

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

T2. 判断通过操作能否让字符串相等 II(Medium)

https://leetcode.cn/problems/check-if-strings-can-be-made-equal-with-operations-ii/

题解(模拟)

同上,分别统计奇偶下标上的元素个数是否相等。

写法 1:基于计数;

class Solution {
    fun checkStrings(s1: String, s2: String): Boolean {
        val U = 26
        val cnts = Array(2) { IntArray(U) }
        for ((i, e) in s1.withIndex()) {
            cnts[i % 2][e - 'a']++
        }
        for ((i, e) in s2.withIndex()) {
            cnts[i % 2][e - 'a']--
        }
        return cnts[0].all {it == 0} && cnts[1].all {it == 0}
    }
}

复杂度分析:

  • 时间复杂度:$O(n + U)$ 线性遍历时间与计数时间;
  • 空间复杂度:$O(U)$ 计数数组空间。

写法 2:基于字符串:

class Solution:
    def checkStrings(self, s: str, t: str) -> bool:
        return all(sorted(s[p::2]) == sorted(t[p::2]) for p in range(2))

复杂度分析:

  • 时间复杂度:$O(nlgn)$ 排序时间;
  • 空间复杂度:$O(n)$ 构造字符串空间。

T3. 几乎唯一子数组的最大和(Medium)

https://leetcode.cn/problems/maximum-sum-of-almost-unique-subarray/

题解(滑动窗口 + 计数)

滑动窗口模板题,维护窗口中不同元素的种类数和总和:

class Solution {
    fun maxSum(nums: List<Int>, m: Int, k: Int): Long {
        var cnts = HashMap<Int, Int>()
        var type = 0
        var sum = 0L
        var ret = 0L
        for (j in nums.indices) {
            // 滑入
            cnts[nums[j]] = cnts.getOrDefault(nums[j], 0) + 1
            if (1 == cnts[nums[j]]!!) type++
            sum += nums[j]
            // 滑出
            if (j >= k) {
                val i = j - k
                cnts[nums[i]] = cnts[nums[i]]!! - 1
                if (0 == cnts[nums[i]]) type --
                sum -= nums[i]
            }
            // 记录
            if (j >= k - 1 && type >= m) {
                ret = max(ret, sum)
            }
        }
        return ret
    }
}

复杂度分析:

  • 时间复杂度:$O(n)$ 线性遍历时间;
  • 空间复杂度:$O(n)$ 散列表空间。

T4. 统计一个字符串的 k 子序列美丽值最大的数目(Hard)

https://leetcode.cn/problems/count-k-subsequences-of-a-string-with-maximum-beauty/

问题分析

  • 问题目标: 求所有长为 $k$ 的子序列中美丽值是最大值的子序列数目;
  • 问题要件: 先计算长为 $k$ 的子序列的最大美丽值,再计算满足美丽值是最大值的子序列方案数;
  • 关键信息 1: 子序列要选择不重复的字母;
  • 关键信息 2: 同一个字符在原字符串中的不同位置可以构造不同子序列;
  • 核心变量: $f(c)$ 是 字符 $c$ 的出现次数,美丽值是子序列中字符的 $f(c)$ 之和;
  • 边界情况: 既然子序列要选择不重复的字母,那么存在边界情况,当 $k$ > 字符串的字符种类数:那么一定不能构造 $k$ 子序列,返回 $0$。

题解一(暴力枚举 + 乘法原理)

最简单的做法,我们可以枚举所有可能的 $k$ 子序列,并记录出现最大美丽值的方案数,怎么实现呢?

  • 方法 1 - 考虑到子序列需要保留原字符串的顺序,直接的想法是枚举字符串中的每个下标 $s[i]$ 选和不选,但是时间复杂度是 $O(2^n)$ 显然不成立;
  • 方法 2 - 事实上我们不需要从原字符串的角度枚举,而是可以从字符集的角度枚举,那样时间复杂度就可以用乘法原理来优化。比如说 a 的出现次数是 $2$,而 b 的出现次数是 $3$,那么所有 ab 可以构造的子序列方案数就是 2 * 3 = 6

那么,方法会不会超时呢,我们来简单分析下:

由于字符集的大小 $U$ 最多只有 $26$ 个,那么子序列的方案数最多有 $C_{26}^k$ 个,而由于 $k$ 大于 $U$ 的方案是不存在的,因此合法的方案数最多只有 $C_{U}^{\frac{U}{2}} = C_{26}^{13} = 10400600$ 约等于 $10^7$。只要我们保证求解每个子问题的时间复杂度是 $O(1)$ 的话是可以通过的。

枚举实现:

class Solution {
    fun countKSubsequencesWithMaxBeauty(s: String, k: Int): Int {
        val MOD = 1000000007
        // 计数
        val cnts = HashMap<Char, Int>()
        for (e in s) {
            cnts[e] = cnts.getOrDefault(e, 0) + 1
        }
        val m = cnts.size
        if (m < k) return 0 // 特判
        // 枚举子序列
        val keys = cnts.toList()
        var maxCount = 0L
        var maxF = 0
        // 回溯
        fun count(index: Int, size: Int, curF: Int, curCount: Long) {
            // 终止条件
            if (size == k) {                
                if (curF > maxF) {
                    maxF = curF
                    maxCount = curCount // 更新最大美丽值方案数
                } else if (curF == maxF) {
                    maxCount = (maxCount + curCount) % MOD // 增加方案数
                }
                return
            }
            if (size + m - index < k) return // 剪枝(长度不够)
            for (i in index until m) {
                val (c, cnt) = keys[i]
                count(i + 1, size + 1, curF + cnt, curCount * cnt % MOD /* 乘法原理 */) 
            }
        }
        count(0, 0, 0, 1)
        return maxCount.toInt()
    }
}

复杂度分析:

  • 时间复杂度:$O(C_m^k)$ 其中 $m$ 为字符种类 ;
  • 空间复杂度:$O(m)$ 散列表空间与递归栈空间。

题解二(排序 + 贪心 + 乘法原理)

考虑 $k = 1$ 的边界情况:

显然需要选择 $f(c)$ 值最大的 $1$ 个字母,如果存在 $m$ 个字母的 $f(c)$ 等于最大值,那么存在 $C_m^1 = m$ 种方案。这说明我们没必要枚举所有字母的子序列: 由于子序列中的字符是不重复的,因此 $k$ 子序列必然要选择 $f(c)$ 值最大的 $k$ 个字母,我们可以将字母按照 $f(c)$ 倒序排序,优先取 $f(c)$ 更大的字母。

具体实现上:

我们将字母按照 $f(c)$ 分桶排序,如果桶内字母数量 $K$ 小于等于 $k$,那么桶内元素都需要选择,否则还要计算桶内元素选择 $k$ 个的方案数:

  • 选择桶内所有元素,方案数为 $cnt^K$
  • 选择桶内部分元素,方案数为 $C_K^k · cnt^k$

其中涉及到幂运算,本质是倍增思想:

// 快速幂 x^n
private fun powM(a: Int, b: Int, mod: Int) : Long {
    var x = a.toLong()
    var n = b.toLong()
    var ret = 1L
    while (n > 0L) {
        if (n % 2 == 1L) ret = ret * x % mod
        x = x * x % mod
        n /= 2
    }
    return ret
}

其中涉及到 组合数

  • 计算式:
// 组合数计算公式 O(k)
private fun comb(n: Int, k: Int, mod: Int) : Int {
    var ret = 1L
    for (i in 1 .. k) {
        ret = ret * (n - i + 1) / i % mod
    }
    return ret.toInt()
}
  • 递推式(杨辉三角):
// 递归 O(n^2)
private fun comb(n: Int, k: Int, mod: Int) : Int {
    if (n == k) {
        return 1
    } else if (k == 1) {
        return n
    } else {
        return (comb(n - 1, k - 1, mod) + comb(n - 1, k, mod)) % mod
    }
}

// 迭代 O(n^2)
private fun comb(n: Int, k: Int, mod: Int) : Int {
    val c = Array(n + 1) { IntArray(n + 1) }
    for (i in 1 .. n) {
        c[i][0] = 1
        c[i][i] = 1
        for (j in 1 until i) {
            c[i][j] = (c[i-1][j] + c[i-1][j-1]) % mod
        }
    }
    return c[n][k]
}

// 组合数计算公式
private fun comb(n: Long, k: Long, mod: Int) : Int {
    var n = n
    var ret = 1L
    for (i in 1 .. k) {
        ret = ret * n-- / i % mod
    }
    return ret.toInt()
}

// 卢卡斯定理
fun Lucas(n: Long, k: Long, mod: Int) : Long {
    if (k == 0L) return 1L;
    return (comb(n % mod, k % mod, mod) * Lucas(n / mod, k / mod, mod)) % mod;
}

完整代码:

class Solution {
    fun countKSubsequencesWithMaxBeauty(s: String, k: Int): Int {
        val MOD = 1000000007
        // 计数
        val cnts = HashMap<Char, Int>()
        var maxCnt = 0
        for (e in s) {
            cnts[e] = cnts.getOrDefault(e, 0) + 1
            maxCnt = max(maxCnt, cnts[e]!!)
        }
        val m = cnts.size
        if (m < k) return 0 // 特判
        // 有序集合
        val map = TreeMap<Int, Int>() { c1, c2 -> 
            c2 - c1
        }
        // 二次频率
        for ((_, c) in cnts) {
            map[c] = map.getOrDefault(c, 0) + 1
        }
        val cntCnts = map.toList()
        // println(cntCnts.joinToString())
        // 构造方案
        var ret = 1L
        var leftK = k
        for ((cnt, K) in cntCnts) {
            if (K > leftK) {
                ret = ret * powM(cnt, leftK, MOD) * comb(K, leftK, MOD) % MOD
            } else {
                ret = ret * powM(cnt, K, MOD) % MOD
            }
            leftK -= K
            if (leftK <= 0) break
        }
        return ret.toInt()
    }

    // 组合数计算公式 C_n^k
    private fun comb(n: Int, k: Int, mod: Int) : Int {
        if (n == k) {
            return 1
        } else if (k == 1) {
            return n
        } else {
            return (comb(n - 1, k - 1, mod) + comb(n - 1, k, mod)) % mod
        }
    }

    // 快速幂 x^n
    private fun powM(x_: Int, n_: Int, mod: Int) : Long {
        var x = x_.toLong()
        var n = n_.toLong()
        var ret = 1L
        while (n > 0L) {
            if (n % 2 == 1L) ret = ret * x % mod
            x = x * x % mod
            n /= 2
        }
        return ret
    }
}

Python 中组合数和幂运算可以很方便地使用库函数:

class Solution:
    def countKSubsequencesWithMaxBeauty(self, s: str, k: int) -> int:
        MOD = 10 ** 9 + 7
        ans = 1
        cnt = Counter(Counter(s).values())
        for c, num in sorted(cnt.items(), reverse=True): # 二次计数
            if num >= k:
                return ans * pow(c, k, MOD) * comb(num, k) % MOD
            ans *= pow(c, num, MOD)
            k -= num
        return 0

复杂度分析:

  • 时间复杂度:$O(n + m)$ 主要时间在枚举字符串的环节;
  • 空间复杂度:$O(m)$ 散列表空间。

推荐阅读

LeetCode 上分之旅系列往期回顾:

⭐️ 永远相信美好的事情即将发生,欢迎加入小彭的 Android 交流社群~

与LeetCode 周赛上分之旅 #43 计算机科学本质上是数学吗?相似的内容:

LeetCode 周赛上分之旅 #43 计算机科学本质上是数学吗?

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 \[彭旭锐] 和 [BaguTree Pro](https://www.mdnice.com/writing/85b28c4e60354865a423728e668fc570) 知识星球提问。** > > 学习数据

LeetCode 周赛上分之旅 #33 摩尔投票派上用场

> ⭐️ **本文已收录到 [AndroidFamily](https://github.com/pengxurui/AndroidFamily),技术和职场问题,请关注公众号 [彭旭锐] 和 [BaguTree Pro] 知识星球提问。** > > 学习数据结构与算法的关键在于掌握问题背后的算法思

LeetCode 周赛上分之旅 #35 两题坐牢,菜鸡现出原形

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 [彭旭锐] 和 [BaguTree Pro] 知识星球提问。** > > 学习数据结构与算法的关键在于掌握问题背后的算法思维框架,你的思考越抽象,它能覆盖的问题域就越广,理解难度也更复杂。在这个专栏里,小彭与你分享

LeetCode 周赛上分之旅 #34 按部就班地解决动态规划问题

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 [彭旭锐] 和 [BaguTree Pro] 知识星球提问。** > > 学习数据结构与算法的关键在于掌握问题背后的算法思维框架,你的思考越抽象,它能覆盖的问题域就越广,理解难度也更复杂。在这个专栏里,小彭与你分享

LeetCode 周赛上分之旅 # 36 KMP 字符串匹配殊途同归

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 \[彭旭锐] 和 [BaguTree Pro](https://files.mdnice.com/user/3257/de950859-eb71-4821-a36b-bebe5cff500d.png) 知识星球提问

LeetCode 周赛上分之旅 # 37 多源 BFS 与连通性问题

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 \[彭旭锐] 和 [BaguTree Pro](https://files.mdnice.com/user/3257/de950859-eb71-4821-a36b-bebe5cff500d.png) 知识星球提问

LeetCode 周赛上分之旅 #38 结合排序不等式的动态规划

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 \[彭旭锐] 和 [BaguTree Pro](https://www.mdnice.com/writing/85b28c4e60354865a423728e668fc570) 知识星球提问。** > > 学习数据

LeetCode 周赛上分之旅 #39 结合中心扩展的单调栈贪心问题

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 \[彭旭锐] 和 [BaguTree Pro](https://www.mdnice.com/writing/85b28c4e60354865a423728e668fc570) 知识星球提问。** > > 学习数据

LeetCode 周赛上分之旅 #40 结合特征压缩的数位 DP 问题

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 \[彭旭锐] 和 [BaguTree Pro](https://www.mdnice.com/writing/85b28c4e60354865a423728e668fc570) 知识星球提问。** > > 学习数据

LeetCode 周赛上分之旅 #41 结合离散化的线性 DP 问题

> ⭐️ **本文已收录到 AndroidFamily,技术和职场问题,请关注公众号 \[彭旭锐] 和 [BaguTree Pro](https://www.mdnice.com/writing/85b28c4e60354865a423728e668fc570) 知识星球提问。** > > 学习数据