给你正整数 low
,high
和 k
。
如果一个数满足以下两个条件,那么它是 美丽的 :
- 偶数数位的数目与奇数数位的数目相同。
- 这个整数可以被
k
整除。
请你返回范围 [low, high]
中美丽整数的数目。
示例 1:
**输入:** low = 10, high = 20, k = 3
**输出:** 2
**解释:** 给定范围中有 2 个美丽数字:[12,18]
- 12 是美丽整数,因为它有 1 个奇数数位和 1 个偶数数位,而且可以被 k = 3 整除。
- 18 是美丽整数,因为它有 1 个奇数数位和 1 个偶数数位,而且可以被 k = 3 整除。
以下是一些不是美丽整数的例子:
- 16 不是美丽整数,因为它不能被 k = 3 整除。
- 15 不是美丽整数,因为它的奇数数位和偶数数位的数目不相等。
给定范围内总共有 2 个美丽整数。
示例 2:
**输入:** low = 1, high = 10, k = 1
**输出:** 1
**解释:** 给定范围中有 1 个美丽数字:[10]
- 10 是美丽整数,因为它有 1 个奇数数位和 1 个偶数数位,而且可以被 k = 1 整除。
给定范围内总共有 1 个美丽整数。
示例 3:
**输入:** low = 5, high = 5, k = 2
**输出:** 0
**解释:** 给定范围中有 0 个美丽数字。
- 5 不是美丽整数,因为它的奇数数位和偶数数位的数目不相等。
提示:
0 < low <= high <= 109
0 < k <= 20
下午两点【b站@灵茶山艾府】 直播讲题,欢迎关注!
本文接着 数位 DP 通用模板 继续讲,因为核心逻辑是一样的。
转换成计算「不超过 high 的美丽整数的数目」减去「不超过 low}-1 的美丽整数的数目」。
对于不超过 high 美丽整数的数目,先将 high 转换成字符串 s,设 s 的长度为 n。
本题我们定义 f(i,\textit{val}, \textit{diff}, \textit{isLimit},\textit{isNum}) 表示构造第 i 位及其之后数位的合法方案数,其中:
- 已经构造的数位,模 k 等于 val。根据取模的原则(请看本文末尾),当我们在第 i 位填入数字 d 后,将 val 更新为 (\textit{val}\cdot 10 + d)\bmod k。在递归终点,判断是否满足 val}=0。
- diff 表示奇数数位的数目与偶数数位的数目的差。在递归终点,判断是否满足 diff}=0。注意我们无需对奇数数位的数目和偶数数位的数目分别各用一个参数表示,那样效率更低。假设填入的数字是 d,那么将 diff 加上 (d\bmod 2) \cdot 2 - 1,对应着奇数加一,偶数减一。
代码实现时,如果你用的不是 Python,需要注意 diff 会出现负数,可以将递归入口处的 diff 设置为 n,在递归终点时判断是否满足 diff}=n。
[sol-Python3]1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
| class Solution: def numberOfBeautifulIntegers(self, low: int, high: int, k: int) -> int: def calc(high: int) -> int: s = str(high) @cache def dfs(i: int, val: int, diff: int, is_limit: bool, is_num: bool) -> int: if i == len(s): return int(is_num and val == 0 and diff == 0) res = 0 if not is_num: res = dfs(i + 1, val, diff, False, False) d0 = 0 if is_num else 1 up = int(s[i]) if is_limit else 9 for d in range(d0, up + 1): res += dfs(i + 1, (val * 10 + d) % k, diff + d % 2 * 2 - 1, is_limit and d == up, True) return res return dfs(0, 0, 0, True, False) return calc(high) - calc(low - 1)
|
[sol-Java]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
| class Solution { public int numberOfBeautifulIntegers(int low, int high, int k) { return calc(high, k) - calc(low - 1, k); }
private int calc(int high, int k) { var s = Integer.toString(high).toCharArray(); int n = s.length; var memo = new int[n][k][n * 2 + 1]; for (int i = 0; i < n; i++) for (int j = 0; j < k; j++) Arrays.fill(memo[i][j], -1); return dfs(0, 0, n, true, false, k, s, memo); }
private int dfs(int i, int val, int diff, boolean isLimit, boolean isNum, int k, char[] s, int[][][] memo) { if (i == s.length) return isNum && val == 0 && diff == s.length ? 1 : 0; if (!isLimit && isNum && memo[i][val][diff] != -1) return memo[i][val][diff]; int res = 0; if (!isNum) res = dfs(i + 1, val, diff, false, false, k, s, memo); int up = isLimit ? s[i] - '0' : 9; for (int d = isNum ? 0 : 1; d <= up; d++) res += dfs(i + 1, (val * 10 + d) % k, diff + d % 2 * 2 - 1, isLimit && d == up, true, k, s, memo); if (!isLimit && isNum) memo[i][val][diff] = res; return res; } }
|
[sol-C++]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
| class Solution { int calc(int high, int k) { auto s = to_string(high); int n = s.length(), memo[n][k + 1][n * 2 + 1]; memset(memo, -1, sizeof(memo)); function<int(int, int, int, bool, bool)> dfs; dfs = [&](int i, int val, int diff, bool is_limit, bool is_num) -> int { if (i == n) return is_num && val == 0 && diff == n; if (!is_limit && is_num && memo[i][val][diff] != -1) return memo[i][val][diff]; int res = 0; if (!is_num) res = dfs(i + 1, val, diff, false, false); int up = is_limit ? s[i] - '0' : 9; for (int d = 1 - is_num; d <= up; d++) res += dfs(i + 1, (val * 10 + d) % k, diff + d % 2 * 2 - 1, is_limit && d == up, true); if (!is_limit && is_num) memo[i][val][diff] = res; return res; }; return dfs(0, 0, n, true, false); }
public: int numberOfBeautifulIntegers(int low, int high, int k) { return calc(high, k) - calc(low - 1, k); } };
|
[sol-Go]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 42 43 44 45 46 47 48 49
| func numberOfBeautifulIntegers(low, high, k int) int { calc := func(high int) int { s := strconv.Itoa(high) n := len(s) memo := make([][][]int, n) for i := range memo { memo[i] = make([][]int, k+1) for j := range memo[i] { memo[i][j] = make([]int, n*2+1) for k := range memo[i][j] { memo[i][j][k] = -1 } } } var dfs func(int, int, int, bool, bool) int dfs = func(i, val, diff int, isLimit, isNum bool) (res int) { if i == n { if isNum && val == 0 && diff == n { return 1 } return 0 } if !isLimit && isNum { p := &memo[i][val][diff] if *p >= 0 { return *p } defer func() { *p = res }() } if !isNum { res += dfs(i+1, val, diff, false, false) } up := 9 if isLimit { up = int(s[i] - '0') } d := 0 if !isNum { d = 1 } for ; d <= up; d++ { res += dfs(i+1, (val*10+d)%k, diff+d%2*2-1, isLimit && d == up, true) } return } return dfs(0, 0, n, true, false) } return calc(high) - calc(low-1) }
|
复杂度分析
- 时间复杂度:\mathcal{O}(n^2kD),其中 n=\mathcal{O}(\log\textit{high}),也就是 high 的十进制表示的长度;D=10。由于每个状态只会计算一次,动态规划的时间复杂度 = 状态个数 \times 单个状态的计算时间。本题中状态个数等于 \mathcal{O}(n^2k),单个状态的计算时间为 \mathcal{O}(D),所以动态规划的时间复杂度为 \mathcal{O}(n^2kD)。
- 空间复杂度:\mathcal{O}(n^2k)。即状态个数。
算法小课堂:模运算
如果让你计算 1234\cdot 6789 的个位数,你会如何计算?
由于只有个位数会影响到乘积的个位数,那么 4\cdot 9=36 的个位数 6 就是答案。
对于 1234+6789 的个位数,同理,4+9=13 的个位数 3 就是答案。
你能把这个结论抽象成数学等式吗?
一般地,涉及到取模的题目,通常会用到如下等式(上面计算的是 m=10):
(a+b)\bmod m = ((a\bmod m) + (b\bmod m)) \bmod m
(a\cdot b) \bmod m=((a\bmod m)\cdot (b\bmod m)) \bmod m
证明:根据带余除法,任意整数 a 都可以表示为 a=km+r,这里 r 相当于 a\bmod m。那么设 a=k_1m+r_1,\ b=k_2m+r_2。
第一个等式:
\begin{aligned}
&\ (a+b) \bmod m\
=&\ ((k_1+k_2) m+r_1+r_2)\bmod m\
=&\ (r_1+r_2)\bmod m\
=&\ ((a\bmod m) + (b\bmod m)) \bmod m
\end{aligned}
第二个等式:
\begin{aligned}
&\ (a\cdot b) \bmod m\
=&\ (k_1k_2m^2+(k_1r_2+k_2r_1)m+r_1r_2)\bmod m\
=&\ (r_1r_2)\bmod m\
=&\ ((a\bmod m)\cdot (b\bmod m)) \bmod m
\end{aligned}
根据这两个恒等式,可以随意地对代码中的加法和乘法的结果取模。