14549：[Ctsc2016]萨菲克斯.阿瑞

#include<stdio.h>
const int maxn=505,mod=1000000007;
int n,m,M,ans;
int c[maxn],C[maxn],f[maxn][maxn],sum[maxn][maxn],fac[maxn],inv[maxn],nfac[maxn];
int main(){
	scanf("%d%d",&n,&M);
	fac[0]=nfac[0]=1;
	for(int i=1;i<=n;i++)
		fac[i]=1ll*fac[i-1]*i%mod,inv[i]=i==1? 1:(mod-1ll*(mod/i)*inv[mod%i]%mod),nfac[i]=1ll*nfac[i-1]*inv[i]%mod;
	for(int i=1;i<=M;i++){
		scanf("%d",&C[i]);
		if(C[i])
			c[++m]=C[i];
	}
	f[0][0]=1;
	for(int r=1;r<=m;r++){
		for(int i=0;i<=n;i++)
			for(int j=0;j<=i;j++)
				sum[i][j]=((i==0||j==0? 0:sum[i-1][j-1])+f[i][j])%mod,f[i][j]=0;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=i;j++){
				f[i][0]=(f[i][0]+1ll*nfac[j]*(sum[i-1][j-1]-(j<=c[r]? 0:sum[i-c[r]-1][j-c[r]-1])+mod))%mod;
				f[i][j]=(0ll+f[i][j]-(sum[i-1][j-1]-(j<c[r]? 0:sum[i-(c[r]-1)-1][j-(c[r]-1)-1]))+mod)%mod;
			}
		ans=(ans+f[n][0])%mod;
	}
	printf("%d\n",1ll*ans*fac[n]%mod);
	return 0;
}