
本文针对在大数范围内(n git dp)技术,将问题复杂度降至o(logn),确保在大规模输入下也能快速准确地得出结果,并正确应用模运算。
1. 问题背景与挑战
我们需要解决一个计算问题:给定一个大整数 N (约束条件为 1 ≤ N
这个问题的核心挑战在于 N 的巨大范围。10^17 是一个非常大的数字,任何 O(N) 复杂度的算法都将导致严重的“时间限制超出”(Time Limit Exceed)问题。
2. 低效的暴力破解方案及其问题分析
最初的尝试通常是直接遍历从 1 到 N-1 的所有数字,并对每个数字计算其各位之和,判断奇偶性,然后累加。以下是这种暴力方法的示例代码:
MOD = 1000000007def getSum(number): """ 计算数字的各位之和。 注意:原始代码中对 number % MOD 的使用是错误的。 """ total = 0 # 原始代码中的错误:对 number % MOD 操作会改变数字本身, # 从而影响各位数字之和的计算。 # 正确的各位之和计算不应在此处引入取模。 while number > 0: total += number % 10 number //= 10 return (total % MOD) % 2 # 这里的 (total % MOD) 也是不必要的def sticker(number): """ 暴力遍历计算小于 number 且各位之和为奇数的数字之和。 """ stickerNeed = 0 for oneDigit in range(1, number): # 遍历范围过大 # 原始代码中的错误:getSum(oneDigit % MOD) 同理会改变数字 if (getSum(oneDigit) == 1): # 假设getSum已修正 stickerNeed += oneDigit return stickerNeed % MOD# number = int(input())# result = sticker(number)# print(result % MOD)
该方案存在的主要问题:
时间复杂度过高: 算法的运行时间与 N 成正比 (O(N))。当 N 达到 10^17 时,即使每秒处理 10^8 次操作,也需要 10^9 秒(约30年),这显然是不可接受的。模运算的错误使用:在 getSum(number) 函数中,对 number % MOD 的操作会改变 number 的值,导致计算出的各位数字之和不正确。例如,getSum(1000000008) 如果先取模会变成 getSum(1),结果显然错误。各位数字之和的计算应该在原始数字上进行,不涉及模运算。getSum 函数返回的是 (total % MOD) % 2。这里的 total % MOD 也是多余的,我们只需要判断 total 的奇偶性,直接 total % 2 即可。模运算 MOD 应该在对最终结果进行累加时使用,以防止中间结果溢出,而不是在计算数字本身的属性(如各位数字之和)时使用。
3. 高效解决方案:模式识别与数位动态规划 (Digit DP)
解决这类问题的关键在于避免逐一遍历,转而利用数字的结构性规律或更高级的算法。对于 N 达到 10^17 这种规模,数位动态规划 (Digit DP) 是最常用且高效的方法。
3.1 核心思想:分而治之与模式识别
尽管数位DP是最终方案,但理解其背后的模式识别思想很重要。
观察数字规律: 在任何一个完整的十进制数段中,数字的各位之和的奇偶性呈现出一定的规律。例如:在 1 到 9 中,1, 3, 5, 7, 9 的各位之和为奇数。在 10 到 19 中,10 (1+0=1), 12 (1+2=3), 14 (1+4=5), 16 (1+6=7), 18 (1+8=9) 的各位之和为奇数。在 20 到 29 中,21 (2+1=3), 23 (2+3=5), 25 (2+5=7), 27 (2+7=9), 29 (2+9=11) 的各位之和为奇数。这种规律表明,对于大范围的数字,我们可以通过数学公式或递推关系来快速计算满足条件的数字之和,而不是逐个检查。
3.2 数位动态规划 (Digit DP) 简介
数位DP是一种用于解决“在给定区间 [L, R] 内,有多少个/这些数的和是多少,满足某个与数字的各位数字相关的性质”的问题的通用技术。它通过递归和记忆化搜索来构建答案,避免重复计算。
数位DP 的基本思路:
问题转化: 通常将 [L, R] 区间的问题转化为 f(R) – f(L-1) 的形式,即计算 [1, X] 范围内的答案。递归函数定义: 定义一个递归函数 dp(index, tight, is_started, current_sum_parity):index: 当前正在考虑的数字位(从最高位开始)。tight: 布尔值,表示当前位是否受到 X 对应位的限制。如果为 True,则当前位只能取 0 到 X 对应位的数字;如果为 False,则可以取 0 到 9。is_started: 布尔值,表示是否已经开始放置非零数字。用于处理前导零。current_sum_parity: 到目前为止已构造数字的各位之和的奇偶性(0表示偶数,1表示奇数)。该函数通常返回一个元组,例如 (count, total_sum),表示在当前状态下,能构造出多少个满足条件的数字以及它们的总和。记忆化: 使用 memo 数组(或字典)存储 dp 函数的计算结果,避免重复计算相同状态。基线条件: 当 index 越界(所有位都已处理完毕)时,根据 current_sum_parity 返回 (1, 0)(如果和为奇数,表示找到一个数,其值为0)或 (0, 0)。状态转移: 遍历当前位可以放置的所有数字(digit),递归调用 dp 函数,并根据 digit 更新 current_sum_parity 和 is_started 等状态。需要特别注意如何计算 total_sum,因为它涉及到当前位的值以及后续位的贡献。
3.3 修正后的迭代求和函数 (用于处理尾部)
虽然数位DP适用于整个范围,但在某些情况下,如果 N 不是特别大,或者为了简化数位DP的实现,可以将问题分解为:一个大的、可以用公式或DP解决的部分,和一个小的、可以用修正后的迭代法解决的“尾部”。
以下是用于处理小范围(例如,一个不足以进行复杂DP计算的短
以上就是优化大数范围求和:避免超时,利用数位DP计算奇数数字和的详细内容,更多请关注创想鸟其它相关文章!
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 chuangxiangniao@163.com 举报,一经查实,本站将立刻删除。
发布者:程序猿,转转请注明出处:https://www.chuangxiangniao.com/p/1378519.html
微信扫一扫
支付宝扫一扫