自然数幂和（差分表法）

问题：

求$\sum_{i = 1}^n i^k$
(1<=k<=2000)
模数为大于n的质数。

题解：

这个东西可以证明k次的自然数幂和一定是个k+1次的多项式，对拉格朗日插值法有点帮助。

二次项展开：
$(n+1)^k-n^k$
$=\sum_{i=0}^{k-1}C_k^i*n^i$

$\sum_{i=1}^n ((i+1)^k - i^k) = (n + 1)^k-1$
$=\sum_{i=1}^n \sum_{j=0}^{k-1}C_k^j*i^j$
$=\sum_{i=0}^{k-1}C_k^j*\sum_{j=1}^nj^i$
设$s(j) =\sum_{i=1}^n i^j$,那么：
$=\sum_{i=0}^{k-1}C_k^j*s(i)$
$=\sum_{i=0}^{k-2}C_k^j*s(i)+k*s(k-1)$

得方程：
$\sum_{i=0}^{k-2}C_k^j*s(i)+k*s(k-1) =(n+1)^k-1$

移项得：
$s(k-1)={{(n+1)^k-1-\sum_{i=0}^{k-2}C_k^i*s(i)}\over k}$

整理得：
$s(k)={{(n+1)^{k+1}-1-\sum_{i=0}^{k-1}C_{k+1}^i*s(i)}\over {k+1}}$

组合数可以杨辉三角预处理。
除以k+1就必须要逆元了，所以模数必须是质数。