stirling（斯特林）数相关

一.两类stirling数的定义及递推公式.

第一类stirling数：定义第一类stirling数$S_1(n,m)$表示$n$个有标号的球分成$m$个无标号的非空轮换的方案数.

这里的轮换定义为旋转同构但翻转不同构，也就是说$1,2,3,4$和$2,3,4,1$是同一个轮换，但$1,2,3,4$与$4,3,2,1$不是同一个轮换.

第一类stirling数的递推公式分为两种情况转移：
1.把最后一个球单独放在一个轮换中，方案数为$S_1(n-1,m-1)$.
2.把最后一个球放在某个轮换的某个球之前，方案数为$(n-1)S_1(n-1,m)$.

那么得到第一类stirling数递推公式：
$$S_1(n,m)=S_1(n-1,m-1)+(n-1)S_1(n-1,m)$$

第二类stirling数：定义第二类stirling数$S_2(n,m)$表示$n$个有标号的球分成$m$个无标号的非空集合的方案数.

与第一类stirling数类似的，第二类stirling数的递推公式同样分为两种情况转移：
1.把最后一个球单独放在一个集合中，方案数为$S_2(n-1,m-1)$.
2.把最后一个球放在某个集合中，方案数为$m{S}_2(n-1,m)$.

那么得到第二类stirling数递推公式：
$$S_2(n,m)=S_2(n-1,m-1)+m{S}_2(n-1,m)$$

二.球与盒子的问题.

球与盒子的问题指的是将$n$个球放入$m$个盒子的计数问题，主要根据球有无标号，盒有无标号，是否可以有空盒分为$8$种情况.

接下来我们一一讨论这$8$种情况的解决方案.

情况1：球无标号，盒有标号，不允许空盒.

插板法，答案为$\left(\genfrac{}{}{0ex}{}{n-1}{m-1}\right)$.

情况2：球无标号，盒有标号，允许空盒.

插板法，答案为$\left(\genfrac{}{}{0ex}{}{n+m-1}{n}\right)$.

情况3：球无标号，盒无标号，不允许空盒.

设$f_{i,j}$表示$i$个球放入$j$个盒子的方案数，分两种情况讨论：
1.每个盒子中放入一个球，方案数为$f_{i-j,j}$.
2.将一个球放入一个盒子，方案数为$f_{i-1,j-1}$.

综合起来就是：
$$f_{i,j}=f_{i-j,j}+f_{i-1,j-1}$$

情况4：球无标号，盒无标号，允许空盒.

与情况3类似的有：
$$f_{i,j}=f_{i-j,j}+f_{i,j-1}$$

情况5：球有标号，盒有标号，允许空盒.

每个球可以放入有$m$种选择，答案为$m^n$.

情况6：球有标号，盒无标号，不允许空盒.

答案为第二类stirling数，即$S_2(n,m)$.

情况7：球有标号，盒无标号，允许空盒.

对第二类stirling数求和，即：
$$\sum _{i=0}^m{S}_2(n,i)$$

情况8：球有标号，盒有标号，不允许空盒.

一个办法是对第二类stirling数乘上$m$的全排列，即$m!S_2(n,m)$.

另一种办法是利用二项式反演，设$f_i$表示$n$个有标号的球在不允许空盒的情况下放入$i$个有标号的盒子的方案数，那么有：
$$m^n=\sum _{i=0}^m\left(\genfrac{}{}{0ex}{}{m}{i}\right)f_i\iff f_m=\sum _{i=0}^m(-1{)}^{m-i}(\genfrac{}{}{0ex}{}{m}{i})i^n$$

三.第二类stirling数的通项公式.

将$f_i=i!S_2(n,i)$代入第二个式子中我们可以得到：
$$\begin{gathered} m!S_2(n,m)=\sum _{i=0}^m(-1{)}^{m-i}(\genfrac{}{}{0ex}{}{m}{i})i^n \\ {S}_2(n,m)=\sum _{i=0}^m\frac{(-1{)}^i}{i!}\frac{(m-i{)}^n}{(m-i)!} \end{gathered}$$

这是第二类stirling数通项公式.

容易看出，第二类stirling数的通项公式是一个卷积形式，利用NTT，我们可以在$O(m(\log n+\log m))$的时间复杂度内求出一行第二类stirling数.


四.自然数幂和与第二类stirling数.

将$f_i=i!S_2(n,i)$代入第一个式子中我们可以得到：
$$m^n=\sum _{i=0}^m\left(\genfrac{}{}{0ex}{}{m}{i}\right)i!S_2(n,i)=\sum _{i=0}^n{S}_2(n,i)m^{\underline{i}}$$

将这个结论带入自然数幂和：
$$\begin{gathered} \sum _{i=1}^n{i}^k \\ =\sum _{i=1}^n\sum _{j=0}^i\left(\genfrac{}{}{0ex}{}{i}{j}\right)j!S_2(k,j) \\ =\sum _{i=0}^{n}i!S_2(k,i)\sum _{j=i}^n\left(\genfrac{}{}{0ex}{}{j}{i}\right) \\ =\sum _{i=0}^{n}i!S_2(k,i)\left(\genfrac{}{}{0ex}{}{n+1}{i+1}\right) \\ =\sum _{i=0}^{n}i!\frac{(n+1{)}^{\underline{i+1}}}{(i+1)!}{S}_2(k,i) \\ =\sum _{i=0}^n\frac{(n+1{)}^{\underline{i+1}}}{i+1}{S}_2(k,i) \\ =\sum _{i=0}^k\frac{(n+1{)}^{\underline{i+1}}}{i+1}{S}_2(k,i) \end{gathered}$$

这样就能做到$O(k\log k)$处理自然数幂和了.


五.Bell数.

对第二类stirling数有了更加深入的探讨后，我们再引入它的一行和Bell数.

Bell数：Bell数$B_n$定义为：
$$B_n=\sum _{i=0}^n{S}_2(n,i)$$

其组合意义为将$n$个有标号的球划分为若干个无标号的集合的方案数.

我们已经可以做到$O(n\log n)$求出单个Bell数，但是这个做法用到了常数大且较为难写的NTT，并不优秀.

来试着将其写成更好计算的形式：
$$\begin{gathered} B_n=\sum _{i=0}^n{S}_2(n,i) \\ =\sum _{i=0}^n\sum _{j=0}^i\frac{(-1{)}^j}{j!}\frac{(i-j{)}^n}{(i-j)!} \\ =\sum _{i=0}^n\frac{(-1{)}^i}{i!}\sum _{j=i}^n\frac{(j-i{)}^n}{(j-i)!} \\ =\sum _{i=0}^n\frac{(-1{)}^i}{i!}\sum _{j=0}^{n-i}\frac{j^n}{j!} \end{gathered}$$

现在只需要对$\frac{j^n}{j!}$这个式子求前缀和，就可以直接做到$O(n\log n)$求解单个Bell数了.


六.Bell数的生成函数.

这个方法可以做到$O(n\log n)$求解所有$i\in [1,n]$的Bell数$B_i$.

首先根据Bell数的组合意义，我们枚举最后一个数所在的集合大小，可以得到如下递推式：
$$B_{n+1}=\sum _{i=0}^n\left(\genfrac{}{}{0ex}{}{n}{i}\right)B_{n-i}$$

设$B_n$的指数型生成函数为$B(x)$，根据上式容易得到：
$$\begin{gathered} B^{\prime }(x)=B(x)e^x \\ \frac{B^{\prime }(x)}{B(x)}=e^x \\ \int \frac{B^{\prime }(x)}{B(x)}\mathrm{d}x=e^x-1 \\ \ln B(x)=e^x-1 \\ B(x)=e^{e^x-1} \end{gathered}$$

利用多项式$\exp$即可做到$O(n\log n)$求解.


七.两类stirling数与三种幂.

首先，我们可以尝试得到上升幂和下降幂之间的关系：
$$\begin{gathered} x^{\underline{n}}=\prod _{i=0}^{n-1}(x-i)=(-1{)}^n\prod _{i=0}^{n-1}(-x+i)=(-1{)}^n(-x{)}^{\overline{n}} \\ {x}^{\overline{n}}=\prod _{i=0}^{n-1}(x+i)=(-1{)}^n\prod _{i=0}^{n-1}(-x-i)=(-1{)}^n(-x{)}^{\underline{n}} \end{gathered}$$

然后，我们在上面已经得到了幂和下降幂之间的关系：
$$x^n=\sum _{i=0}^n{S}_2(n,i)x^{\underline{i}}$$

将$x$换成$-x$，我们可以得到幂和上升幂之间的关系：
$$\begin{gathered} (-x{)}^n=\sum _{i=0}^n{S}_2(n,i)(-x{)}^{\underline{i}} \\ (-x{)}^n=\sum _{i=0}^n(-1{)}^i{S}_2(n,i)x^{\overline{i}} \\ {x}^n=\sum _{i=0}^n(-1{)}^{n-i}{S}_2(n,i)x^{\overline{i}} \end{gathered}$$

我们发现只有第二类stirling数为系数的关系式，于是猜想一个第一类stirling数为系数的关系式：
$$x^{\overline{n}}=\sum _{i=0}^n{S}_1(n,i)x^i$$

归纳证明后发现它是对的，于是再次将$x$换$-x$可以得到：
$$x^{\underline{n}}=\sum _{i=0}^n(-1{)}^{n-i}{S}_1(n,i)x^i$$

于是我们得到了在幂和上升幂、下降幂之间转化的公式.


八.stirling反演.

我们将下降幂与幂的关系式代入幂与下降幂的关系式可以得到：
$$\begin{gathered} x^n=\sum _{i=0}^n{S}_2(n,i)x^{\underline{i}} \\ =\sum _{i=0}^n{S}_2(n,i)\sum _{j=0}^i(-1{)}^{i-j}{S}_1(i,j)x^j \\ =\sum _{i=0}^n\sum _{j=0}^i(-1{)}^{i-j}{S}_2(n,i)S_1(i,j)x^j \\ =\sum _{i=0}^n\sum _{j=i}^n(-1{)}^{j-i}{S}_2(n,j)S_1(j,i)x^i \end{gathered}$$

取出其中第$i$项：
$$[n=i]x^i=\sum _{j=i}^n(-1{)}^{j-i}{S}_2(n,j)S_1(j,i)x^i$$

同时除以$x^i$得到反转公式：
$$[n=i]=\sum _{j=i}^n(-1{)}^{j-i}{S}_2(n,j)S_1(j,i)$$

于是我们猜想stirling反演：
$$f_n=\sum _{i=0}^n{S}_2(n,i)g_i\iff g_n=\sum _{i=0}^n(-1{)}^{n-i}{S}_1(n,i)f_i$$

代入验证：
$$\begin{gathered} f_n=\sum _{i=0}^n{S}_2(n,i)g_i \\ =\sum _{i=0}^n{S}_2(n,i)\sum _{j=0}^i(-1{)}^{i-j}{S}_1(i,j)f_j \\ =\sum _{i=0}^n\sum _{j=0}^i(-1{)}^{i-j}{S}_2(n,i)S_1(i,j)f_j \\ =\sum _{i=0}^n\sum _{j=i}^n(-1{)}^{j-i}{S}_2(n,j)S_1(j,i)f_j \\ =\sum _{i=0}^n[n=i]f_i \\ =f_n \end{gathered}$$

同理我们可以得到另一个反转公式和stirling反演：
$$\begin{gathered} \sum _{j=i}^n(-1{)}^{j-i}{S}_1(n,j)S_2(j,i)=[n=i] \\ {f}_n=\sum _{i=0}^n{S}_1(n,i)g_i\iff g_n=\sum _{i=0}^n(-1{)}^{n-i}{S}_2(n,i)f_i \end{gathered}$$