简单斯特林反演

定义

定义$n^{\underline{k}}$为n的k次下降幂，则$n^{\underline{k}}=n\cdot (n-1)\cdot ...\cdot (n-k+1)=\frac{n!}{(n-k)!}=\frac{n!k!}{(n-k)!k!}={(}_k^n)k!$。

定义$n^{\overline{k}}$为n的k次上升幂，则$n^{\overline{k}}=n\cdot (n+1)\cdot ...\cdot (n+k-1)=\frac{(n+k-1)!}{(n-1)!}=\frac{(n+k-1)!k!}{(n-1)!k!}={(}_k^{n+k-1})k!$

斯特林反演

斯特林反演揭示了上升幂，下降幂与普通幂的关系。

斯特林反演公式

• 下降幂转普通幂
  $x^{\underline{n}}=\stackrel{n}{\sum _{k=0}}\left[\begin{array}{c}n\\ k\end{array}\right](-1{)}^{n-k}{x}^k$
• 上升幂转普通幂
  $x^{\overline{n}}=\stackrel{n}{\sum _{k=0}}\left[\begin{array}{c}n\\ k\end{array}\right]x^k$
• 普通幂转下降幂（常用）
  $x^n=\stackrel{n}{\sum _{k=0}}\left\{\begin{array}{c}n\\ k\end{array}\right\}{x}^{\underline{k}}=\stackrel{n}{\sum _{k=0}}\left\{\begin{array}{c}n\\ k\end{array}\right\}\left(\begin{array}{c}x\\ k\end{array}\right)k!$
• 普通幂转上升幂
  $x^n=\stackrel{n}{\sum _{k=0}}\left\{\begin{array}{c}n\\ k\end{array}\right\}(-1{)}^{n-k}{x}^{\overline{k}}=\stackrel{n}{\sum _{k=0}}(-1{)}^{n-k}\left\{\begin{array}{c}n\\ k\end{array}\right\}\left(\begin{array}{c}x+k-1\\ k\end{array}\right)k!$

事实上，必须说到，若有$n,k$不合法，则斯特林数为$0$，比如$n<k,n<0或k<0$。
（你可能已经意识到了${(}_k^x)k!=A_x^k$，但是为了方便处理，我们不这么写。）

启发

这两个式子一般在斯特林反演中不直接用到，但是可以给人启发：
组合恒等式：$k\left(\begin{array}{c}n\\ k\end{array}\right)=n\left(\begin{array}{c}n-1\\ k-1\end{array}\right)$
因而有：$k^{\underline{x}}\left(\begin{array}{c}n\\ k\end{array}\right)=n^{\underline{x}}\left(\begin{array}{c}n-x\\ k-x\end{array}\right)$
因此我们注意到，在普通幂转下降幂的过程中，要选择原有组合数底项的次幂，作为转下降幂的对象。

证明

我们证明一下常用公式：
$x^n=\stackrel{n}{\sum _{k=0}}\left\{\begin{array}{c}n\\ k\end{array}\right\}{x}^{\underline{k}}=\stackrel{n}{\sum _{k=0}}\left\{\begin{array}{c}n\\ k\end{array}\right\}\left(\begin{array}{c}x\\ k\end{array}\right)k!$

组合意义：从组合意义上说，这个式子的成立是显然而明确的。

董老师的组合解释：
左边表示n个不同球放在x个不同盒子里，方案数显然有$x^n$。
右边表示枚举非空盒子的个数k，从x个盒子中选出k个盒子（组合数），把n个不同球放进去（斯特林数），对盒子做全排列。

我的组合解释：
左式表示有n个盒子，要往每个盒子里放一个小球，小球共有x种颜色的方案数。
右式表示的是，假设n个盒子已经放完了球，把放了颜色相同的球的盒子分为一组，总共有k组，现在指定这k组盒子分别放了什么颜色的球，首先从x种颜色中取出k种颜色作为用到的颜色，让颜色和盒子一一配对。注意到颜色和盒子是无关的，所以再对颜色做全排列。
（大雾）

然后还可以用数学归纳法证明。

事实上有：

• 下降幂与普通幂的关系：
  
• 上升幂与普通幂的关系；
  

证明就是直接拆开。

我们归纳常用公式，其他公式同理：
假设对于$x^{n-1}$符合归纳。

套一下刚才的式子：

注意到左式在$i=0$时等于$0$，由于斯特林数为0。

QED.

其他归纳过程，有类似的证明。

例题

题目要求：

这个n实在是太高了，绝对不能留。

首先考虑把f带进去：

注意存在$k^i$与$\left(\begin{array}{c}n\\ k\end{array}\right)$，考虑对$k$进行普通幂转下降幂：

注意到最后两个组合数构成分离式：

越来越乱了，整理一下：

最后一个和式显然是二项式定理：

感觉可以做了？似乎。

注意到这个组合数不太好处理，当然用快速阶乘也不是处理不了…，继续拆一下式子：

这个形式看起来很好做了：

递推可以写出两项来对比一下。

然后就做完了。

#include<iostream>
using namespace std;
long long pow(long long a,long long b,long long c) {
	long long ans=1;
	while(b) {
		if(b&1) ans=ans*a%c;
		a=a*a%c;
		b>>=1;
	}
	return ans;
}
long long a[10005],J[10005],x[10005],f[10005];
int n,m;
long long M;
int s[10005][10005];
int main() {
	cin>>n>>x[1]>>M>>m;
	for(int i=0;i<=m;i++) cin>>a[i];
	x[0]=1;
	f[0]=pow(x[1]+1,n,M);
	for(int i=1;i<=m;i++)
		x[i]=x[i-1]*x[1]%M,
		f[i]=pow(x[1]+1,n-i,M);
	s[0][0]=1;
	for(int i=1;i<=m;i++)
		for(int j=1;j<=i;j++)
			s[i][j]=((long long)s[i-1][j-1]+(long long)j*s[i-1][j])%M;
	J[0]=1;
	for(int i=1;i<=m;i++)
		J[i]=J[i-1]*(n-i+1)%M;
	long long ans=0;
	for(int i=0;i<=m;i++) {
		long long sum=0;
		for(int j=0;j<=i;j++)
			(sum+=(long long)s[i][j]*J[j]%M*x[j]%M*f[j]%M)%=M;
		(ans+=a[i]*sum)%=M; 
	}
	cout<<ans;
	return 0;
}

习题
注意存在n<k的情况。

后记

于是皆大欢喜。