码迷,mamicode.com
首页 > 其他好文 > 详细

[模板] 斯特林数,性质以及求法

时间:2019-05-23 16:47:58      阅读:96      评论:0      收藏:0      [点我收藏+]

标签:line   math   rac   时间复杂度   复杂   排列   n+1   降幂   单点   

斯特林数

第一类斯特林数

定义

定义(无符号)第一类斯特林数, 表示 \(n\) 个有标号元素分为 \(k\) (\(k \le n\)) 个无标号环排列 (翻转算两种) 的方案数.
记做 \({n \brack k}\), 或者 \(c(n,k)\), \(|s(n,k)|\), \(\left|s_n^k\right|\).

枚举最后一个元素的位置, 容易得到递推式
\[\begin{cases} {n \brack k} = {n-1 \brack k-1} + (n-1){n-1 \brack k} & (n > 0, k < n)\{n \brack 0} = 0 & (n>0) \{n \brack n} = 1 & (n \ge 0)\\end{cases}\]

有符号的斯特林数可以表示为
\[ s(n,k) = (-1)^{n-k} {n\brack k}\]

它的定义是这样的:
首先定义 \(n\) 次下降幂和上升幂
\[x^{\underline n} = \frac {x!}{(x-n)!} = \prod_{i=x-n+1}^x i\]
\[x^{\overline n} = \frac {(x+n-1)!}{(x-1)!} = \prod_{i=x}^{x+n-1} i = (-1)^n (-x)^{\underline n}\]
容易发现这是一个常数项为 \(0\)\(n\) 次多项式. 那么可以定义
\[x^{\underline n} = \sum_{i=0}^n s(n,i) x^i\]
\[x^{\overline n} = \sum_{i=0}^n {n \brack i} x^i\]

求法... 待更

性质... 待更

第二类斯特林数

定义

定义第二类斯特林数, 表示 \(n\) 个有标号元素分为 \(k\) (\(k \le n\)) 个无标号非空集合的方案数.
记做 \({n \brace k}\), 或者 \(S(n,k)\).

同样的, 可以得到递推式:
\[ \begin{cases} {n \brace k} = {n-1 \brace k-1} + k{n-1 \brace k} & (n > 0, k < n)\{n \brace 0} = 0 & (n>0) \{n \brace n} = 1 & (n \ge 0)\\end{cases} \]

第二类斯特林数也有关于下降幂/上升幂的定义:
\[ x^{n}=\sum_{k=0}^{n}{n \brace k}x^{\underline k} \]
\[ x^{n}=\sum_{k=0}^{n}(-1)^{n-k}{n \brace k}x^{\overline k} \]

求法

我们显然有一个 \(O(n^2)\) 的递推求法.

但有时我们想求出单点的值, 或者求出 \(S(n,k), k \in \{0,1,\dotsc, n\}\) 的所有值, 那么就需要更快的算法.

\(S'(n,k)\) 表示 \(n\) 个有标号元素分为 \(k\)有标号非空集合的方案数, 显然有
\[ S'(n,k) = S(n,k) * k!\]
\(g(n,k)\) 表示 \(n\) 个有标号元素分为 \(k\)有标号可空集合的方案数, 还有
\[ g(n,k) = k^n \]

枚举非空集合的个数, 可以发现
\[ g(n,k) = \sum_{i=0}^k {k \choose i}S'(n,i) \]
二项式反演:
\[ S'(n,k) = \sum_{i=0}^k (-1)^{k-i} {k \choose i}g(n,i) \]
那么就可以求出 \(S(n,k)\):
\[ S(n,k) = \frac{S'(n,k)}{k!} \]
\[ S(n,k) = \sum_{i=0}^k (-1)^{k-i}\frac{1}{(k-i)!} \frac{i^n}{i!} \]

这样, 我们就可以在 \(O(n)\) 的时间复杂度内求出 \(S(n,k)\);
容易发现这是一个卷积, 直接FFT即可求出 \(S(n,k), k \in \{0,1,\dotsc, n\}\), 时间复杂度 \(O(n \log n)\).

[模板] 斯特林数,性质以及求法

标签:line   math   rac   时间复杂度   复杂   排列   n+1   降幂   单点   

原文地址:https://www.cnblogs.com/ubospica/p/10912502.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!