高效求解非递减整数划分方案数:动态规划与数学优化实战

千墨同学_5339

千墨同学_5339

2026-05-06

628人浏览

原创

高效求解非递减整数划分方案数:动态规划与数学优化实战

本文详解如何将“分糖果”问题建模为带约束的整数划分问题,通过数学转化(零允许→正整数划分)、递推关系推导及空间优化dp实现,将时间复杂度稳定控制在 o(c×d),轻松应对 1000 组大规模输入。

本文详解如何将“分糖果”问题建模为带约束的整数划分问题,通过数学转化(零允许→正整数划分)、递推关系推导及空间优化dp实现,将时间复杂度稳定控制在 o(c×d),轻松应对 1000 组大规模输入。

该问题本质是:给定 D 个孩子和 C 颗糖果,要求分配方案满足 非递减顺序(即第 i 名的孩子所得糖果数 ≤ 第 i+1 名),等价于统计所有长度为 D、元素和为 C、且满足 $x_1 \le x_2 \le \cdots \le x_D$ 的非负整数序列个数。

关键洞察在于:这类序列与 整数 C 划分为至多 D 个部分(允许 0)的非递减划分 一一对应。但标准整数划分通常定义为正整数之和,因此我们采用经典平移技巧——令
$$ y_i = x_i + 1 \quad (i = 1,2,\dots,D) $$
则 $y_i \ge 1$ 且严格满足 $y_1 \le y_2 \le \cdots \le y_D$,其和为 $C + D$。此时,$(y_1,\dots,y_D)$ 正是整数 $C+D$ 的一个 恰好 D 部分的(正整数)非递减划分,即「D-分拆数」$p(C+D,\, D)$。

该分拆数满足如下经典递推式(可由 Ferrers 图共轭或生成函数导出):

$$ \begin{aligned} p(n, k) &= \begin{cases} 1 & n = k = 0 \ 0 & n \le 0 \text{ 或 } k \le 0 \ p(n - k,\, k) + p(n - 1,\, k - 1) & \text{其他} \end{cases} \end{aligned} $$

其中:

  • $p(n - k, k)$:表示所有部分 ≥ 2 的 D-分拆(每部分减 1 后仍为 D-分拆);
  • $p(n - 1, k - 1)$:表示至少有一个部分为 1 的 D-分拆(删去该 1 后变为 (D−1)-分拆)。

因此,原问题答案即为 p(C + D, D)。

✅ 推荐实现:空间优化的二维 DP(Bottom-up)

为避免递归开销与栈溢出,并支持 Q ≤ 1000 批量查询,我们预处理一个全局 DP 表 dp[d][c],其中 dp[d][c] = p(c, d),范围覆盖最大 D=100、最大 C=5000+100=5100:

MOD = (1 = d else 0
        term2 = dp[d - 1][c - 1] if c >= 1 and d >= 1 else 0
        dp[d][c] = (term1 + term2) % MOD

# 处理输入
import sys
input = sys.stdin.read().split()
q = int(input[0])
idx = 1
out_lines = []
for _ in range(q):
    d = int(input[idx])
    c = int(input[idx + 1])
    idx += 2
    # answer = p(c + d, d)
    res = dp[d][c + d] if (c + d) <h3>⚠️ 注意事项与优化要点</h3>
  • 边界安全:c + d 最大为 5000 + 100 = 5100,在 MAX_C = 5100 下完全覆盖,无需运行时判断越界。
  • 空间压缩(进阶):若内存受限,可仅保留两行 dp_prev[c] 和 dp_curr[c],将空间从 O(D×C) 降至 O(C),但需注意 dp_curr[c] 依赖 dp_curr[c-d](同层前向),故推荐使用完整二维表以保逻辑清晰。
  • 勿混淆模型:原始错误代码中尝试暴力枚举所有组合并事后排序校验,时间复杂度高达 $O(C^D)$,完全不可行;而早期 DP 版本中的状态定义有误(如 dp[i][j] += dp[i][j-i] 缺失语义支撑),导致结果偏差。
  • 模运算位置:务必在每次加法后取模,防止整数溢出(Python 虽支持大整数,但题目明确要求模 2^30−1)。

✅ 总结

本题是算法优化的经典范例:
? 建模转换 是提速前提——将“非递减非负分配”转化为“正整数 D-分拆”;
? 数学递推 提供最优子结构——p(n,k) = p(n−k,k) + p(n−1,k−1) 是核心加速引擎;
? 预处理 + 查表 是应对多查询的黄金策略——用 O(D×C) 空间换 O(1) 单次查询,总耗时稳定在毫秒级。

最终代码可在 1000 组输入下稳定通过时限,兼具正确性、效率与可维护性。

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1731

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

4424

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1709

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

25877

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3147

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3167

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1183

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

596

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2463

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
热门推荐
/
最新课程
phpStudy极速入门视频教程
phpStudy极速入门视频教程

共6课时 | 54.7万人学习

独孤九贱(4)_PHP视频教程
独孤九贱(4)_PHP视频教程

共89课时 | 133.7万人学习