$B_n$ 是集合大小为 $n$ 的集合的划分数目，注意这里划分成多少个集合都是可以的。根据贝尔数的定义我们就可以得到贝尔数和第二类斯特林数之间的关系：

$$
\sum\limits_{i=0}^nS(n,i)=B_n
$$

递推公式

贝尔数有递推公式：

$$
B_n=\sum\limits_{i=0}^{n-1}\dbinom{n-1}{i}B_i
$$

考虑如何证明这个式子。我们考虑第 $n$ 个元素，如果它与 $k$ 个元素划分到一个集合中去，那么我们要做的就是把剩下的 $n-k-1$ 个元素进行划分，方案数就是 $\binom{n-1}{k}B_{n-k-1}=\binom{n}{n-k-1}B_{n-k-1}$ 由此不难得出原式。

性质

• $B_{n+m}=\sum_{j=0}^m\sum_{i=0}nj^{n-i}\left{ m\atop j \right}\binom{n}{i}B_i$

证明：我们考虑这个式子的组合意义，这 $m$ 个元素怎么分配，枚举 $j$ 表示分配到 $j$ 个集合里，然后考虑这 $n$ 个元素如何分配，我们枚举 $i$ 个单独分配，分配方案就是 $B_i$，然后把剩下的放在那 $j$ 个集合里，方案数为 $j^{n-i}$。所以这个式子成立。

• Touchard 同余：$B_{p+n}\equiv B_n+B_{n+1}\bmod p$

根据第二类斯特林数的递推公式 ：

$$
\left{ n\atop m \right} = \frac{1}{m!}\sum\limits_{i=0}^m(-1)i\binom{m}{i}(m-i)^n
$$

我们可以得到，当 $n$ 为质数时，$1<m<n$ 时有 $\left{ n\atop m \right}\equiv 0\bmod n$ 成立。具体证明在博主的另一篇博客《斯特林数与斯特林反演中》提到，这里不再赘述。

我们用第一条性质把 $B_{n+p}$ 展开得到：

$$
\begin{aligned}
B_{n+p}=\sum\limits_{i=0}^{n\sum\limits_{j=0}}pj^{n-i}\left{ p\atop j \right}\binom{n}iB_i\
=\sum\limits_{j=0}^p\left{ p\atop j \right}\sum\limits_{i=0}^nj\binom{n}iB_i\
\end{aligned}
$$

根据第二类斯特林数上面所述的性质，在 $\bmod\ p$ 意义下，当 $1<j<p$ 时为 $0$，而当 $j=0$ 时根据第二类斯特林数的定义，也为 $0$，所以上面这个式子可以化简如下：

$$
\begin{aligned}
&\equiv \left{ p\atop 1 \right}\sum\limits_{i=0}^n\binom{n}{i}B_i+\left{ p\atop p\right}\sum\limits_{i=0}^np\binom{n}{i}B_i\
&\equiv B_{n+1}+B_n
\end{aligned}
$$

后面那个式子之所以成立是因为当且仅当 $i=n$ 的时候后面那个式子才不会被 $\bmod\ p$ 给模掉。

• $\sum_{r=0}^k(-1)r\binom{k}{r} B_{n-r+1}=\sum_{r=0}^{n-k}\binom{n-k}{r}B_{n-r}$

注意到 $k=0$ 时的情况是平凡的。我们考虑对 $k$ 归纳证明，我们分别让 $F(n,k)$ 等于等号两边的式子，通过对两边都进行证明 $F(n,k)-F(n-1,k)=F(n,k+1)$ 就可以证明。

这里只证明等号左边满足这个式子，等号右边同理。

设 $F(n,k)=\sum_{r=0}^k(-1)r\binom kr B_{n-r+1}$，那么我们可以得到：

$$
\begin{aligned}
&F(n-1,k)+F(n,k+1)\&=\sum\limits_{r=0}^k(-1)r\binom kr B_{n-r}+\sum\limits_{r=0}^{k+1}(-1)r\binom { k+1}{r}G_{n-r+1}\
&=G_{n+1}+\sum\limits_{r=1}^{k+1}(-1)r\binom kr G_{n-r+1}=F(n,k)
\end{aligned}
$$

• $B_{n+p^m}\equiv mB_n+B_{n+1}\bmod p$

这个式子可以看作是上面那个式子的一个扩展，我们同样利用第一条性质把这个东西展开得到：

$$
B_{n+p^{m}=\sum_{j=0}}\sum_{i=0}^{pm}j^{pm-i}\left{ n\atop j \right}\binom{p^m}{i}B_i
$$

根据卢卡斯定理，我们可以得到 $\binom{p^m}{r}\equiv 0\bmod p(0<r<p^m)$，所以上面这个式子在 $\bmod \ p$ 意义下就可以写成：

$$
\begin{aligned}
B_{n+p^{m}&\equiv\sum\limits_{j=0}}n \left{ n\atop j \right}(j^{pm}B_0+B_{p^m})\
&\equiv \sum\limits_{j=0}^n\left{ n\atop j \right}(j+B_{p^m})
\end{aligned}

$$

通过让性质 $1$ 的 $n=1$ 我们可以得到一个式子：

$$
B_{1+m}=\sum\limits_{j=0}^m(j+1)\left{ m\atop j \right}
$$

所以我们可以得到：

$$
\begin{aligned}
B_{n+p^m}&\equiv (B_{p^{n}-1)\sum\limits_{j=0}}n\left{ n\atop j \right}+\sum\limits_{j=0}^n(j+1)\left{ n\atop j \right}\
&=(B_{p^{m}-1)\sum\limits_{j=0}}n\left{ n\atop j \right}+B_{n+1}\
&=(B_{p^m}-1)B_n+B_{n+1}
\end{aligned}
$$

所以现在我们只需要证明下面这个结论：

$$
B_{p^m}\equiv m+1 \bmod p
$$

我们考虑如何证明，这里会用得到一点抽象代数中的定义，不过基本还是初等方法。

这里 $\Bbb{Z}/p^m\Bbb{Z}$ 指的是 $p^m$ 的完全剩余系。对于 $x,y \in \Bbb{Z}/p^m\Bbb{Z}$，定义 $f_y(x)=x+y\bmod p^m$，设 $P$ 是这个剩余系的一个划分，那么 $f_y(P)$ 表示的就是对这个划分中的每一个集合的每一个元素进行操作所得到的划分，于是 $f_y(P)$ 就可以理解成一个类似于置换的东西。设 $X={f_y(x)|y\in \Bbb{Z} }$，定义一个划分 $P$ 为不动的指的是 $\forall f_y(x)\in X,P=f_y(P)$ .其实这个东西等价于 $P=f_y(P)$。同样我们也可以定义不动的划分，有以下性质：任何非不动的划分必然属于一个大小为 $P$ 的正整数幂的轨道。

容易发现这个轨道的大小就是最小的 $k$ 满足 $f_k(P)=P$，我们已经知道有 $f_{p^m}=P$，这其实就是一个不断执行 $p^m$ 次 $f_1$ 的过程，我们可以形象的把这个过程看做一个环，那么我们在这个环上走 $p^m$ 次一定能够走回源点，那么这也就是说，环的长度一定能够乘除 $p^m$，而在非不动的意义下，我们也可以得到 $f_1(P)\not=P$，所以我们最终可以得到以上性质，这告诉我们，所有非不动的划分的数量可以被 $p$ 整除，这也就是说，总的划分数量，即 $B_{p^m}$，和不动的划分的数量在模 $p$ 的意义下相等。

假设一个划分 $P$ 是不动的，$a,b$ 分别在这个划分的某个集合中 $A,B$ 中，明显 $f_{b-a}(a)=b$ 又因为这个划分是不动的，所以我们有 $f_{b-a}(A)=B$，因此在一个不动的划分中，所有划分出的集合大小都相等，不妨设为 $p^j$，我们考虑证明这个划分有以下性质：在一个集合中的所有元素模 $p^{m-j}$ 都相同。

我们用反证法来证明这个性质，设 $P$ 是一个不动的划分，$A\in P$ 且 $a,b\in A$，并且 $a\not\equiv b\bmod p^{m-j}$ 所以 $f_{b-a}(a)=b$，所以 $f_{b-a}$ 是一个把 $A$ 映射到 $A$ 的映射，因此对于任意的正整数 $r$，我们有 $f_{r(b-a)}(a)$ 为 $A$ 中的某个元素，显然 $f_{r(b-a)}(a)\equiv f_{s(b-a)}(a)\bmod p^{m}$，当且仅当 $a+r(b-a)\equiv a+s(b-a)\bmod p^m$，有假设可知 $r\equiv s$，因此对于 $1\le r\le p^{j+1}$ $f_{r(b-a)}(a)$ 都是不同的元素，这表明 $|A|>p^j$ 矛盾。

于是这告诉我们，对于划分中的集合大小都为 $p^j$ 的不动的划分，它是唯一的。于是这样的集合个数就为 $m+1$ 个，于是我们可以得到 $B_{p^m}\equiv m+1$。由此上面的性质得以证明。

设 $B^n=B_n$ 也就是说我们规定了一个符号。

我们考虑证明一个更强的性质：$B_{kp^s+n}\equiv B^n(B+s)k$

我们考虑证明当 $s$ 固定的时候，对于所有的 $k$ 和 $n$ 上面这个式子成立，我们考虑对 $k$ 进行归纳，当 $k=1$ 的时候我们已经证明这个是对的，则：

$$
B_{(k+1)p^s+n}\equiv B_{kp^s+ps+n}\equiv B^{ps+n}(B+s)^k\equiv B^n(B+s)
$$

由此我们可以得到这个性质成立。

其实如果我们不知道 $k=1$ 的时候式子成立也可以证明，我们考虑这个式子在 $s=0$ 的时候显然成立，我们考虑取 $k=1$ 对 $s$ 进行归纳也可以证明这个式子的正确性。

由此我们可以在把这个同余性质加强一下，贝尔数满足：

$$
G_{\sum k_rp^r}\equiv \prod(G+r)^{k_r}
$$

证明这个式子，只需要把上面那个更强的式子进行迭代使用。

由此我们可以得到这样的一个式子：

$$
G_{\frac{p^p-1}{p-1}+n}\equiv G_n\bmod p
$$

证明：

$$
G_{\frac{p^p-1}{p-1}+n}\equiv G_{p^{p-1}+p+...1+n}\equiv G^{n+1}(G+p-1)...(G+1)
$$

不难发现这就是一个下降幂的形式，根据第二类斯特林数的性质，我们可以得到：

$$
G^{n+1}(G-1)=G^{n+p}-G=G_n
$$

所以我们可以得到以下的性质：

• 贝尔数模素数 $p$ 的周期为 $\frac{p^p-1}{p-1}$

贝尔三角形

我们可以通过构造贝尔三角形来得到贝尔数，具体方法如下：$b_{1,1}=1$，$b_{k,1}=b_{k-1,k-1}$，$b_{x,y}=b_{x,y-1}+b_{x-1,y-1}$。

这样构造，每一行的第一个数就是贝尔数，即 $b_{k,1}=B_k$。

这样构造的原因如下：

我们知道：$\Delta B_m=B_{m+1}-B_m$，不难发现并证明：

$$
\Delta^n B_m=\sum\limits_{k=0}^n(-1)\binom{n}{k}B_{m+k}
$$

我们令 $m=1$ 可以得到：

$$
\begin{aligned}
\Delta^{nB_1&=\sum\limits_{k=0}}n(-1)^{n-k}\binom{n}{k}B_{k+1}\
&=\sum\limits_{k=0}^n(-1)\binom{n}{k}\sum\limits_{j=0}^k\binom{k}{j}B_j\
&=\sum\limits_{j=0}^{n\sum\limits_{k=j}}n(-1)^{n-k}\binom{n}{k}\binom{k}{j}B_j\
&=\sum\limits_{j=0}^{nB_j\binom{n}{j}\sum\limits_{k=j}}n(-1)^{n-k}\binom{n-k}{k-j}\
&=\sum\limits_{j=0}^nB_j\binom nj \sum\limits_{j=0}^{n-k}(-1)j\binom{n-k}j\
&=\sum\limits_{j=0}^nB_j\binom nj [k=n]=B_n
\end{aligned}
$$

那么我们就可以利用差分公式计算贝尔数。

这个公式提示我们，如果我们按照上面的方式去构造贝尔数的话，不难发现前几项的规律，把没一行的最后一个数往前拆，就会发现这个东西其实就是 $\Delta^n B_1$。其实就是不断往下做差分的过程。我们看这个三角形的局部：

  1
  1 2
  2 3 5
  5 7 10 15
  15 20 27 37 52

就会发现下一行的差分是这一行。

那么这个三角形怎么理解呢？我们这样来理解，设 $a_{i,j}$ 表示第 $i$ 行的第 $j$ 个数，我们首先来看第一列，不难发现，这一列数全部是 $\Delta^n B_1$，然后第二列数，不难发现，这一列数全是 $\Delta^n B_2$，后面依次类推，那么它是怎么保持这种性质的呢？原因在于：

$$
B_n=\Delta B_{n-1}+B_{n-1}\Rightarrow \Delta^{n} B_n=\Delta^{n+1} B_{n-1}+\Delta^n B_{n-1}
$$

注意我们消除而加上 $\Delta$ 算子必须要满足的条件是我们应用的式子没有依靠于这个数列的特殊性质。

由上面这个式子可以知道递推式子成立，再由 $\Delta^n B_1=B_n$ 就可以知道我们把每一行的最后一个数抄下来当做第一行的第一个数的原因了。

生成函数

设贝尔数的 EGF 为 $f(x)$，由此我们可以得到：

$$
f(x)=\sum\limits_{i\ge0}\frac{B_i}{i!}x^i
$$

我们对上面这个东西进行求导可以得到：

$$
\begin{aligned}
f'(x)&=\sum\limits_{i\ge 0}\frac{B_{i+1}}{i!}x^i\
&=\sum\limits_{i\ge 0}\frac{x^i}{i!}\sum\limits_{i\ge j}\binom{i}{j}B_j\
&=\sum\limits_{j\ge 0}\frac{B_jx^j}{j!}\sum\limits_{i\ge j}\frac{x^{i-j}}{(i-j)!}\
\end{aligned}
$$

众所周知，我们把 $e^x$ 在原点处的泰勒展开式写出来可以得到：

$$
e^x=\sum\limits_{i\ge 0}\frac{x^i}{i!}
$$

把 $e^{ix}$ 在原点处的泰勒展开式写出来可以得到：

$$
e^{ix}=\sum\limits_{j\ge 0}\frac{x^jij}{j!}
$$

所以原式我们可以继续写得到：

$$
LHS=\sum\limits_{j\ge 0}\frac{B_jx_j}{j!}e^x=f(x)ex
$$

由此我们可以得到 $f(x)$ 的封闭形式为 $f(x)=e^{ex+c}$。由 $f(0)=1$ 我们可以得到 $c=-1$，所以 $f(x)=e^{ex-1}$。

至此我们得到贝尔数 EGF 的封闭形式。

设 $g(x)=e^{ex}$，我们把指数上的 $e^x$ 看做一个整体，然后我们可以得到 $e^{ex}=\sum_{i\ge 0}\frac{e^{xi}}{i!}$

所以我们可以得到 $f(x)=\frac 1e g(x)=\frac 1e \sum_{i\ge 0}\frac{e^{xi}}{i!}$。

我们继续推可以得到：

$$
\begin{aligned}
LHS&=\frac 1e \sum\limits_{i\ge 0}\frac{e^{xi}}{i!}\
&=\frac 1e \sum\limits_{i\ge 0}\frac{1}{i!}\sum\limits_{j\ge 0}\frac{i^jxj}{j!}\
&=\frac 1e \sum\limits_{j\ge 0}\frac{x^j}{j!}\sum\limits_{i\ge 0}\frac{i^j}{i!}
\end{aligned}
$$

我们去式子两边 $\frac{x^n}{n!}$ 的系数，就可以得到：$B_n=\frac 1e\sum_{i\ge 0}\frac{i^n}{i!}$（Dobinski公式）。

引用

• 知乎回答

• 高等组合学笔记(七): Bell数,第一,二类Stirling数及其递推关系_zorchp-CSDN博客

• 论文《AN ELEMENTARY (NUMBER THEORY) PROOF OF TOUCHARD’S CONGRUENCE》

• 论文《numbers generated by the function $e^{ex-1}$》