自然数幂和 斯特林数

2018 UPD：
其实第二类斯特林数做自然数幂和更简单，这里简单写一下：
由一个基本式子出发

$$n^k=\sum_{i=0}^k\big\{^k_{\ i}\big\}[n]_i$$
考虑对 $n$求和
$$Ans=\sum_{i=0}^ni^k=\sum_{i=0}^n\sum_{j=0}^k\big\{^k_j\big\}[i]_j$$
$$=\sum_{j=0}^k\big\{^k_j\big\}\sum_{i=0}^n[i]_j$$
$$=\sum_{j=0}^k\big\{^k_j\big\}j!\sum_{i=0}^n\big(^i_j\big)$$
$$=\sum_{j=0}^k\big\{^k_j\big\}j!\big(^{n+1}_{j+1}\big)$$
然后就没了， $\big\{^k_j\big\}=\big\{^{k-1}_{j-1}\big\}+j\big\{^{k-1}_{\ \ j}\big\}$可以 $O(k^2)$递推，而且与下面第一类同理对模数没有要求

---

一年多前接触过这个东西，但是当时理解很浅，很快就忘了
今天Drin_E讲题的时候手推的，感觉很劲，浑然天成，逻辑性强，推理、转化自然，易理解
还是要mark一下的=w=

Description

求

$$S_{n,k}=\sum_{i=0}^ni^k$$

Analysis

构造一个多项式

$$P_x^k=x(x-1)...(x-k+1)$$
显然也可以写成如下形式，容易发现系数只与 $k,i$有关，设为 $Ss_{k,i}$
$$P_x^k=\sum_{i=1}^kx^iSs_{k,i}$$
考虑如何求出 $Ss_{k,i}$，为了使用归纳法，构造 $P_x^{k+1}$，观察其与 $P_x^k$各项系数的关系
$$P_x^{k+1}=\sum_{i=1}^{k+1}x^iSs_{k+1,i}=x(x-1)...(x-k)$$
$$=(x-k)P_x^k=(x-k)\sum_{i=1}^kx^iSs_{k,i}$$
$$=\sum_{i=1}^kx^{i+1}Ss_{k,i}-kx^iSs_{k,i}$$
通过比较各项系数，可以发现
$$Ss_{k+1,i}=Ss_{k,i-1}-kSs_{k,i}$$
于是 $Ss_{k,i}$可以 $O(n^2)$递推求出
实际上，这里的 $Ss$表示带符号 第一类斯特林数，有实际意义
上面铺垫结束，为了靠拢所求式，单独分出 $x^k$，移项
$$x^k=P_x^k-\sum_{i=1}^{k-1}x^iSs_{k,i}$$
靠拢所求式，两边求和
$$S_{n,k}=\sum_{x=0}^nx^k=k!\sum_{x=0}^nC_x^k-\sum_{x=0}^n\sum_{i=1}^{k-1}x^iSs_{k,i}$$
$$=k!C_{n+1}^{k+1}-\sum_{i=1}^{k-1}Ss_{k,i}\sum_{x=0}^nx^i$$
$$=\frac{P_{n+1}^{k+1}}{k+1}-\sum_{i=1}^{k-1}Ss_{k,i}S_{n,i}$$
减号前面部分暴力做，分子为k+1个连续自然数相乘，其中一定有k+1这个因数，所以可以避免除法
$Ss$可以预处理， $S$递推，$O(k^2)$
再贴一些别人的其他方法
伯努利数_from_PhilipsWeng
WerkeyTom_FTD