求自然数幂和的各种方法（还有坑）

求

$$\sum_{i=1}^{n}{i^k}\mod p$$

高斯消元

一个定理，k次方和一定可以由0~k-1次方和表示出来，设方程组解出来就好了。
$O(k^3)$

倍增

我们设

$$f(n,k)=\sum_{i=1}^{n}i^k$$
怎么算呢？
我们采用分治思想。
如果n是奇数那么
$$f(n,k)=f(n-1,k)+n^k$$
否则，我们可以先求出n/2的f，然后对于n/2+1~n中的每个数都可以表示为n/2+i
那么
$$f(n,k)=f(n/2,k)+\sum_{i=1}^{n/2}(n/2+i)^k$$
二项式定理
$$f(n,k)=f(n/2,k)+\sum_{i=1}^{n/2}\sum_{j=0}^{k}C_k^j(n/2)^{k-j}i^j$$
调换一下
$$f(n,k)=f(n/2,k)+\sum_{j=0}^{k}C_j^kf(n/2,j)(n/2)^{k−j}$$
$O(k^2\log n \log k)$

第一类斯特林数

参考WerKeyTom_FTD的blog
设

$$S(n,k)=\sum_{i=1}^{n}i^k$$
引入 第一类斯特林数
$$P_m^n=\sum_{i=0}^{n}Ss(n,i)m^i$$
发现 $$j^k=Ss(k,k)j^k$$
于是 $$j^k=k!C_j^k-P_j^k+Ss(k,k)j^k$$
$$j^k=k!C_j^k-(P_j^k-Ss(k,k)j^k)$$
观察括号内的式子
$$P_j^k-Ss(k,k)j^k=\sum_{i=0}^{k}Ss(k,i)j^i-Ss(k,k)j^k$$
$$P_j^k-Ss(k,k)j^k=\sum_{i=0}^{k-1}Ss(k,i)j^i$$
于是 $$j^k=k!C_j^k-\sum_{i=0}^{k-1}Ss(k,i)j^i$$
$$S(n,k)=\sum_{j=1}^{n}(k!C_j^k-\sum_{i=0}^{k-1}Ss(k,i)j^i)$$
$$S(n,k)=k!\sum_{j=1}^{n}C_j^k-\sum_{i=0}^{k-1}\sum_{j=1}^{n}Ss(k,i)j^i$$
$$S(n,k)=k!C_{n+1}^{k+1}-\sum_{i=0}^{k-1}Ss(k,i)\sum_{j=1}^{n}j^i$$
$$S(n,k)={P_{n+1}^{k+1}\over {k+1}}-\sum_{i=0}^{k-1}Ss(k,i)S(n,i)$$
减号前面暴力做就好了。因为它是k+1个连续自然数相乘，必然会有k+1这个因数。预处理Ss，递推S。 $O(k^2)$