标签：modp 22 2020.6 ekx sum equiv 同余 Bn mathrm

今天粉兔同学问了一个问题：如何证明贝尔数的 Touchard’s Congruence 性质：
$B_{n+p} \equiv B_{n+1} + B_n \pmod p$Bn+p​≡Bn+1​+Bn​(modp)
其中 $p$p 是质数，$B_n$Bn​ 是贝尔数。

为了证明这个问题，我们首先证明一个引理：
引理 1：$\sum_k {p\brace k} x^k \equiv x + x^p \pmod p$∑k​{kp​}xk≡x+xp(modp)。
证明：我们考虑多项式 $x^p = \sum_k {p\brace k} x^{\underline k}$xp=∑k​{kp​}xk​，这是由下降幂转化得到的。而 $x^p \equiv x \pmod p$xp≡x(modp)，$x^{\underline p} \equiv 0 \pmod p$xp​≡0(modp)，因此我们可以知道对于 $\sum_{k<p} {p\brace k} x^{\underline k} \equiv x \pmod p$∑k<p​{kp​}xk​≡x(modp)，等式两边都是小于 $p$p 次的多项式，必然是对应相等的，而 ${p \brace p} = 1${pp​}=1 显然。故 $\sum_k {p\brace k} x^k \equiv x + x^p \pmod p$∑k​{kp​}xk≡x+xp(modp) 得证。

接下来考虑这样一件事，记贝尔数的 EGF 是 $B(x) = \exp (\mathrm e^x - 1)$B(x)=exp(ex−1)，那么 $B'(x) = B(x)\mathrm e^x$B′(x)=B(x)ex，我们考虑对 $B(x)$B(x) 连续求导 $n$n 次，我们知道求导完必然是一个 $B^{(n)}(x) = B(x)\sum_k a_{n,k} \mathrm e^{kx}$B(n)(x)=B(x)∑k​an,k​ekx 的形式，归纳则有

$\begin{aligned} B^{(n+1)}(x) & = \left(B^{(n)}(x)\right)'\\ B(x)\sum_k a_{n+1,k} \mathrm e^{kx} &= \left(B(x)\sum_k a_{n,k} \mathrm e^{kx}\right)'\\ &= \sum_k a_{n,k} \left(B(x)\mathrm e^{kx}\right)'\\ &= \sum_k a_{n,k} \left(kB(x)\mathrm e^{kx} + B(x)\mathrm e^{(k+1)x}\right)\\ \sum_k a_{n+1,k} \mathrm e^{kx} &= \sum_k a_{n,k} \left(k\mathrm e^{kx} + \mathrm e^{(k+1)x}\right)\\ a_{n+1,k} &= ka_{n,k} + a_{n,k-1} \end{aligned}$B(n+1)(x)B(x)k∑​an+1,k​ekxk∑​an+1,k​ekxan+1,k​​=(B(n)(x))′=(B(x)k∑​an,k​ekx)′=k∑​an,k​(B(x)ekx)′=k∑​an,k​(kB(x)ekx+B(x)e(k+1)x)=k∑​an,k​(kekx+e(k+1)x)=kan,k​+an,k−1​​

带入初值，我们容易得到
$B^{(n)}(x) = B(x) \sum_k {n\brace k} \mathrm e^{kx}$B(n)(x)=B(x)k∑​{kn​}ekx

令 $n=p$n=p，有
$\begin{aligned} B^{(p)}(x) &= B(x) \sum_k {p\brace k} \mathrm e^{kx}\\ & \equiv B(x) (\mathrm e^x + \mathrm e^{px})\\ & \equiv B(x) (\mathrm e^x + 1)\\ & \equiv B'(x) + B(x)\\ B_{n+p} &\equiv B_{n+1} + B_n \end{aligned}$B(p)(x)Bn+p​​=B(x)k∑​{kp​}ekx≡B(x)(ex+epx)≡B(x)(ex+1)≡B′(x)+B(x)≡Bn+1​+Bn​​

故原问题得证。

后面还提到了 $B_{n+p^m} \equiv mB_{n+1} + B_n \pmod p$Bn+pm​≡mBn+1​+Bn​(modp)，这个其实就不是很有新东西了，不过一个很爽的事情是我们可以用 Umbral Calculus 来快速得到：现在挪用上指标，我们有

$B^p = 1 + B$Bp=1+B

因此 $B^{p^m} = (1+B)^{p^{m-1}} = 1 + B^{p^{m-1}} = 2 + B^{p^{m-2}} = \cdots = m + B$Bpm=(1+B)pm−1=1+Bpm−1=2+Bpm−2=⋯=m+B

得证。

标签：modp,22,2020.6,ekx,sum,equiv,同余,Bn,mathrm
来源： https://blog.csdn.net/EI_Captain/article/details/106909220

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。