斯特林数Stirling小结

最新推荐文章于 2021-05-16 21:34:15 发布

Liang-梁

最新推荐文章于 2021-05-16 21:34:15 发布

阅读量202

点赞数

分类专栏：斯特林数

本文链接：https://blog.csdn.net/qq_37555704/article/details/106788639

版权

斯特林数专栏收录该内容

2 篇文章 0 订阅

订阅专栏

放球问题即可

考虑现在有 $n$ 个元素， $x$ 种颜色，有 $k$ 个环，那么 $k$ 个环的方案数就是 $S_1(n,k)$
也可以一个一个元素的放，要么加入之前的位置 $i - 1$ ，要么独自成环 $x$
那么就是 $x * (x + 1) * (x + 2) * . . . * (x + n - 1)$

$x^{\underline n }=(-1)^n(-x)^n$
$x^{\overline n }=(-1)^n(-x)^{\underline n}$
利用这个转化即可

[省选联考 2020 A 卷] 组合数问题
在这里插入图片描述
$\sum_{i=0}^nx^if(i)\binom{n}{i}=\sum_{i=0}^ma_i\sum_{j=0}^nj^ix^j\binom{n}{j}=\sum_{i=0}^n a_i\sum_{k=0}^i\begin{Bmatrix}i\\k\end{Bmatrix}\sum_{j=0}^nj^{\underline k}x^j\binom{n}{j}$
考虑右边三个的组合意义，
选 $j$ 个位置，并给这 $j$ 个染色（颜色从 $1$ 到 $x$ ），再选出 $k$ 个排队
可以先选出 $k$ 个排队，再给这 $k$ 个染色，再从剩下的 $n - k$ 个选 $j - k$ 个再染色,于是有
$\sum_{i=0}^n a_i\sum_{k=0}^i\begin{Bmatrix}i\\k\end{Bmatrix}\sum_{j=0}^nj^{\underline k}x^j\binom{n}{j}=\sum_{i=0}^n a_i\sum_{k=0}^i\begin{Bmatrix}i\\k\end{Bmatrix}\sum_{j=0}^nn^{\underline k}x^kx^{j-k}\binom{n-k}{j-k}$
$=\sum_{i=0}^n a_i\sum_{k=0}^i\begin{Bmatrix}i\\k\end{Bmatrix}n^{\underline k}x^k\sum_{j=0}^nx^{j-k}\binom{n-k}{j-k}=\sum_{i=0}^n a_i\sum_{k=0}^i\begin{Bmatrix}i\\k\end{Bmatrix}n^{\underline k}x^k\sum_{j=0}^{n-k}x^{j}\binom{n-k}{j}$
$=\sum_{i=0}^n a_i\sum_{k=0}^i\begin{Bmatrix}i\\k\end{Bmatrix}n^{\underline k}x^k(x+1)^{n-k}$
$=\sum_{i=0}^n a_i\sum_{j=0}^i\begin{Bmatrix}i\\j\end{Bmatrix}n^{\underline j}x^j(x+1)^{n-j}$
然后就能 $O(k^2)$ 算了

int n,x,m;
int a[MAXN+5];
int S2[MAXN+5][MAXN+5],nd[MAXN+5];
void Init(){
	S2[0][0]=1;
	for(int i=1;i<=MAXN;i++)
		for(int j=1;j<=i;j++)
			S2[i][j]=Add(Mul(j,S2[i-1][j]),S2[i-1][j-1]);
	return ;
}
int main(){
	n=read(),x=read(),Mod=read(),m=read();
	Init();
	for(int i=0;i<=m;i++)
		a[i]=read();
	nd[0]=1;
	for(int i=1;i<=min(n,m);i++)
		nd[i]=Mul(nd[i-1],(n-i+1));
	int ans=0;
	for(int t=0;t<=m;t++){
		int tmp=0;
		for(int i=0;i<=min(n,t);i++)
			tmp=Add(tmp,Mul(Mul(S2[t][i],nd[i]),Mul(Pow(x+1,n-i),Pow(x,i))));
		ans=Add(ans,Mul(a[t],tmp));
	}
	printf("%d\n",ans);
	return 0;
}