第一类斯特林数求自然数幂和学习小记

目标

求$$\sum _{i=0}^n{i}^k$$

前置技能

第一类斯特林数

第一类斯特林数$s(n,m)$定义为有n个人，编号分别为1-n，排成m个循环排列的方案数。
递推式：$s(n,m)=s(n-1,m-1)+(n-1)\ast s(n-1,m)$
证明：若第n个人自己排成一个循环排列，则方案数为$s(n-1,m)$；若第n个人加入到原来的循环排列中，则可以放在任意一个人的左边，方案数为$(n-1)\ast s(n-1,m)$

另一种形式

$众所周知，C_n^k表示从n件物品中选出k件的方案数，$P_n^k$表示从n件物品中选出k个排列的方案数$
显然有$$C_n^k=\frac{P_n^k}{k!}$$
根据定义我们可以得到$$P_n^k=n\ast (n-1)\ast ...\ast (n-k+1)$$
展开后可得$$P_n^k=S(k,k)\ast n^k-S(k,k-1)\ast n^{k-1}...$$
所以有$$P_n^k=\sum _{i=0}^k(-1{)}^{k-i}\ast s(k,i)\ast n^i$$
这个才是第一类斯特林数的定义。
所以第一类斯特林数S就是排列数公式的展开式的系数，也是如上所述的那个东西。

第一类斯特林数的递推

边界条件：$s(p,0)=0(p>0)，s(p,p)=1(p>=0)$

一个小结论

$$\sum _{i=0}^n{C}_i^k=C_{n+1}^{k+1}$$
证明：因为$C_{n+1}^{k+1}=C_n^k+C_n^{k+1}$，等式两边可以同时把$C_n^k$消掉，然后就又得到了一个相同形式的式子。所以通过归纳法便可以证明这个结论。

求自然数幂和

我们设$$S_k(n)=\sum _{i=0}^n{i}^k$$
根据第一类斯特林数的定义，有$$C_n^k=\frac{P_n^k}{k!}=\frac{{\sum }_{i=0}^k(-1{)}^{k-i}\ast s(k,i)\ast n^i}{k!}$$
把式子中的$n^k$提出来，得$$n^k=C_n^k\ast k!-\sum _{i=0}^{k-1}(-1{)}^{k-i}\ast s(k,i)\ast n^i$$
那么$$S_k(n)=\sum _{i=0}^n{i}^k$$
$$=\sum _{i=0}^n(C_i^k\ast k!-\sum _{j=0}^{k-1}(-1{)}^{k-j}\ast s(k,j)\ast i^j)$$
$$=k!\ast \sum _{i=0}^n{C}_i^k-\sum _{i=0}^n\sum _{j=0}^{k-1}(-1{)}^{k-j}\ast s(k,j)\ast i^j$$
$$=k!\ast C_{n+1}^{k+1}-\sum _{j=0}^{k-1}(-1{)}^{k-j}\ast s(k,j)\ast S_j(n)$$
$$=\frac{P_{n+1}^{k+1}}{k+1}-\sum _{j=0}^{k-1}(-1{)}^{k-j}\ast s(k,j)\ast S_j(n)$$

于是就可以愉快地$O(k^2)$递推$S_k(n)$啦。边界条件$S_1(n)=\frac{n\ast (n+1)}{2}$

因为$P_{n+1}^{k+1}=(n+1)\ast n\ast (n-1)\ast ...\ast (n-k+1)$，而连续k+1个数中必然有一个是k+1的倍数，所以在用这种方法求自然数幂和时，由于不存在除法，所以并不用考虑模数是什么。

代码

int get_sum(int n)
{
	sum[1]=(LL)n*(n+1)/2%MOD;
	for (int i=2;i<=k;i++)
	{
		sum[i]=1;
		for (int j=0;j<i+1;j++)
			if ((LL)(n-j+1)%(i+1)==0) sum[i]=(LL)(n-j+1)/(i+1)*sum[i]%MOD;
			else sum[i]=(LL)(n-j+1)*sum[i]%MOD;
		for (int j=1;j<i;j++)
			if ((i-j)%2==0) (sum[i]-=(LL)s[i][j]*sum[j]%MOD)%=MOD;
			else (sum[i]+=(LL)s[i][j]*sum[j]%MOD)%=MOD;
	}
	return (sum[k]+MOD)%MOD;
}