luogu CF932E Team Work

背景：

斯特林数学习中。

题目传送门：

https://www.luogu.org/problem/CF932E

题意：

求${\sum }_{i=1}^n{C}_n^i{i}^k$。

思路：

大力化化简式子。
$$\sum _{i=1}^n{C}_n^i{i}^k$$

考虑到：$i^k={\sum }_{j=0}^k{S}_{i,j}{C}_{i,j}j!$，因此有：
$$=\sum _{i=1}^n{C}_n^i\sum _{j=0}^k{S}_{k,j}{C}_{i,j}j!$$

$$=\sum _{i=1}^n\frac{n!}{i!(n-i)!}\sum _{j=0}^k{S}_{k,j}\frac{i!}{j!(i-j)!}j!$$

$$=\sum _{i=1}^n\frac{n!}{(n-i)!}\sum _{j=0}^k{S}_{k,j}\frac{1}{(i-j)!}$$

$$=\sum _{j=0}^k{S}_{k,j}\sum _{i=1}^n\frac{n!}{(n-i)!(i-j)!}$$

上下同时乘上$(n-j)!$，有：
$$=\sum _{j=0}^k{S}_{k,j}\sum _{i=1}^n\frac{(n-j)!}{(n-i)!(i-j)!}\frac{n!}{(n-j)!}$$

$$=\sum _{j=0}^k{S}_{k,j}\frac{n!}{(n-j)!}\sum _{i=1}^n{C}_{n-j}^{n-i}$$

考虑到${\sum }_{i=1}^n{C}_n^i=2^i$，二项式定理可证（即使$n-i>n-j$也没有关系，因为$C_{n-j}^{n-i}=0$）因此有：

$$=\sum _{j=0}^k{S}_{k,j}\frac{n!}{(n-j)!}{2}^{n-j}$$

$\Theta (n^2)$处理$S_{k,j}$即可，然后带入求解，即可通过。

能更优秀吗？
考虑到$S_{k,j}={\sum }_{i=0}^j\frac{(-1{)}^i}{i!}\frac{(j-i{)}^k}{(j-i)!}$，因此有：

$$=\sum _{j=0}^k\sum _{i=0}^j\frac{(-1{)}^i}{i!}\frac{(j-i{)}^k}{(j-i)!}\frac{n!}{(n-j)!}{2}^{n-j}$$

$$=\sum _{j=0}^k\frac{n!}{(n-j)!}{2}^{n-j}\sum _{i=0}^j\frac{(-1{)}^i}{i!}\frac{(j-i{)}^k}{(j-i)!}$$

后面的一个$\sum$是一个卷积的形式，因此$\Theta (k\log k)$求解即可。
然而这道题的模数不和谐，因此我才懒得去打。

显然$\frac{n!}{(n-j)!}$不能直接求，怎么办呢？
$\frac{n!}{(n-j)!}=(n-j+1)\ast (n-j+2)\ast (n-j+3)...\ast n$，即可。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define mod 1000000007
using namespace std;
	int n,k,ans=0;
	int s[2][5010];
void init()
{
	s[1][1]=1;
	for(int i=2;i<=k;i++)
		for(int j=1;j<=i;j++)
			s[i&1][j]=((LL)s[(i&1)^1][j-1]+(LL)j*s[(i&1)^1][j]%mod)%mod;
}
int ksm(int x,int k)
{
	int tot=1;
	for(;k;k>>=1)
	{
		if(k&1) tot=(LL)tot*x%mod;
		x=(LL)x*x%mod;
	}
	return tot;
}
int main()
{
	scanf("%d %d",&n,&k);
	init();
	int tmp=n;
	for(int i=1;i<=min(n,k);i++)
	{
		ans=((LL)ans+(LL)s[k&1][i]*tmp%mod*ksm(2,n-i)%mod)%mod;
		tmp=(LL)tmp*(n-i)%mod;
	}
	printf("%d",ans);
}