惯性聚合 高效追踪和阅读你感兴趣的博客、新闻、科技资讯
阅读原文 在惯性聚合中打开

推荐订阅源

Martin Fowler
Martin Fowler
D
DataBreaches.Net
F
Fortinet All Blogs
阮一峰的网络日志
阮一峰的网络日志
博客园_首页
Apple Machine Learning Research
Apple Machine Learning Research
H
Help Net Security
M
MIT News - Artificial intelligence
美团技术团队
人人都是产品经理
人人都是产品经理
freeCodeCamp Programming Tutorials: Python, JavaScript, Git & More
Cyber Security Advisories - MS-ISAC
Cyber Security Advisories - MS-ISAC
The Cloudflare Blog
有赞技术团队
有赞技术团队
L
LangChain Blog
博客园 - Franky
钛媒体:引领未来商业与生活新知
钛媒体:引领未来商业与生活新知
博客园 - 【当耐特】
S
SegmentFault 最新的问题
V
Visual Studio Blog
Blog — PlanetScale
Blog — PlanetScale
Hugging Face - Blog
Hugging Face - Blog
B
Blog
I
InfoQ

博客园_首页

Linux实操--组管理、权限管理和定时任务 Java + EasyExcel 实现单个接口导出多个Excel Mem0 源码解析系列(二):提示词工程的深度剖析 Openclaw TaskFlow究竟是什么?和普通Skill技能有什么区别 博文阅读密码验证 - 博客园 嘉立创开源:应该是全网MicroPython教程最多的开发板 Hermes Agent 集成实践:从协议到生产 2026年AI编程工具横评:Cursor、Codex、Claude Code、Zed、Windsurf Java程序员必看的RAG入门教程 2026 AI效率神器:Superpowers + Claude Code 保姆级教程 本地大模型部署全攻略:从 0 到 1 玩转 Ollama 【从0到1构建一个ClaudeAgent】内存管理-上下文压缩 .NET 高级开发 | 设计、实现一个事件总线框架 电子小白入门之NE555 3. WorkBuddy:隐藏玩法,一键召唤专家,让 AI 以"专家身份"给你干活 和AI一起搞事情#3:Claude Teammate 游戏开发翻车实录 【OpenClaw】通过 Nanobot 源码学习架构---(7)Memory C# .NET 周刊|2026年3月3期 我在 Debian 11 上把 K8s 单机搭起来了,过程没你想的那么顺(/opt 目录版) 深度学习进阶(七)Data-efficient Image Transformer CLI+Skill搭建浏览器AI自动化框架,告别一切重复枯燥任务 告别Token账单无底洞:OpenClaw本地部署,重塑企业数据主权的唯一解 FastAPI+Vue:文件分片上传+秒传+断点续传,这坑我帮你踩平了! SBTI 爆火后,我做了个程序员版的 CBTI。。已开源 + 附开发过程 多模态检索开始进入工程期:用 Sentence Transformers 搭建可落地的 Multimodal RAG 100多行代码实现一个最简单的Agent(用ReAct) Claude Code 通关手册(八):推荐 5 个 Hooks,代码质量提升 3 倍 老板:“有人截图了!”。安全部门:“收到,马上查暗水印!” - why技术 技术之外,皆是人间 C#/.NET/.NET Core技术前沿周刊 | 第 69 期(2026年4.01-4.12)
第二类斯特林数学习笔记
liduoduo2021 · 2026-05-13 · via 博客园_首页

定义与递推式

\(\left\{ \begin{matrix} n \\ k \end{matrix} \right\}\) ,也可记作 \(S(n,k)\) ,表示将 \(n\) 个有标号的数分成 \(k\) 个互不区分的组的方案数。

根据组合意义,能推出递推式为:

\[S(n,k)=S(n-1,k-1)+S(n-1,k)\cdot k \]

方幂转下降幂

\[n^m=\sum_{k=0}^{min(n,m)}S(m,k)\binom {n}{k}k! \]

左式组合意义是从 \(n\) 个数中选 \(m\) 个数,可重且有标号的方案数。右式是枚举本质不同的点的个数 \(k\),再将 \(m\) 个数 划分到 \(k\) 组 ,且组之间区分带上 \(k!\) 的系数,最后乘上 \(n\) 个点中选 \(k\) 个点的方案数。

观察式子两边与 \(n\) 关联的项,左边是方幂形式,右边是组合数,也就是下降幂形式。下降幂的形式能够保证与 \(n\) 有关的次数不会超过 \(n\) ,有时会起到优化作用。

应用

给定一棵树,且树有一个代价 \(cost(S)\),且 \(cost(S)\le|S|\) ,且 \(f\) 可以同子树合并得到,现在要求对于所有非空节点集合 \(S\) 的生成最小连通块的 \(cost(S)^M\) 之和。

一个平凡的思路是直接记录 \(f_{u,i}\) 表示以 \(u\) 为根的子树的所有方案中,\(cost(S)^i\) 的代价和。合并用二项式定理展开,假设两部分所有方案中的 \(cost\) 的序列分别是 \(L_0,L_1,…\)\(R_0,R_1,…\)

\[ \begin{aligned} f_{u,k}&=\sum_i\sum_j (L_i+R_j)^k\\ &=\sum_i\sum_j\sum_{t=0}^k\binom {k}{t}{L_i}^t{R_i}^{k-t}\\ &=\sum_{t=0}^k \binom {k}{t} (\sum_i{L_i}^t)(\sum_j{R_j}^{k-t})\\ &=\sum_{t=0}^k \binom {k}{t} f_{v1,t}\cdot f_{v2,k-t} \end{aligned} \]

这样做的时间复杂度是 \(O(nM^2)\)

考虑优化,记录 \(g_{u,i}\) 表示以 \(u\) 为根的子树的所有方案中,\(\binom {cost(S)}{i}\) 的代价和,通过方幂转下降幂的公式能够在 \(O(nM)\) 的时间复杂度复原答案。通过范德蒙德恒等式合并,有:

\[\begin{aligned} g_{u,k}&=\sum_i\sum_j\binom{L_i+R_j}{k}\\ &=\sum_i\sum_j\sum_{t=0}^k\binom{L_i}{t}\binom{R_j}{k-t}\\ &=\sum_{t=0}^k(\sum_i\binom{L_i}{t})(\sum_j\binom{R_j}{k-t})\\ &=\sum_{t=0}^kg_{v1,t}\cdot g_{v2,k-t} \end{aligned} \]

这样看似还是平方的,但是 \(g\) 的第二维不会超过子树大小,所以时间复杂度就是 \(O(nM)\)

本质上就是从维护 \(n^k\) 变成 \(\binom{n}{k}\) 枚举量就从 \(m\) 变成 \(min(n,m)\),最后通过公式将 \(\binom{n}{k}\) 复原成 \(n^k\)

显式公式

\[S(n,k)=\frac{1}{k!}\sum_{i=0}^k(-1)^{k-i}\binom{k}{i}i^n \]

考虑组合意义,可以先求出区分组的方案数,最后再带上 \(\frac{1}{k!}\) 。区分组的方案本质就是将 \(n\) 个数染色,出现颜色数为 \(k\) 的方案数,这个可以通过将 \(n\) 个数染成 \(1\)\(k\) 的颜色二项式反演得到。

生成函数

固定 \(k\) 的EGF:

\[\sum_{n=k}^\infty \frac{x^n}{n!}S(n,k)=\frac{(e^x-1)^k}{k!} \]

固定 \(k\) 的OGF:

\[\sum_{n=k}^\infty x^nS(n,k)=\sum_{i=0}^k(-1)^{k-i}\binom{k}{i}\frac{1}{1-ix}\\ \sum_{n=k}^\infty x^nS(n,k)=\frac{x^k}{\prod_{i=0}^{k}1-ix} \]