bzoj 5215: [Lydsy2017省队十连测]商店购物

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int ha=1000000007;
const int maxn=10000000;
int jc[maxn+5],ni[maxn+5];
int f[2][100005],n,m,K,tot;
int ans=0,a[305],now,pre;

inline int add(int x,int y){
	x+=y;
	return x>=ha?x-ha:x;
}

inline int ksm(int x,int y){
	int an=1;
	for(;y;y>>=1,x=x*(ll)x%ha) if(y&1) an=an*(ll)x%ha;
	return an;
}

inline void init(){
	jc[0]=1;
	for(int i=1;i<=maxn;i++) jc[i]=jc[i-1]*(ll)i%ha;
	ni[maxn]=ksm(jc[maxn],ha-2);
	for(int i=maxn;i;i--) ni[i-1]=ni[i]*(ll)i%ha;
}

inline int C(int x,int y){
	if(x<y) return 0;
	return jc[x]*(ll)ni[y]%ha*(ll)ni[x-y]%ha;
}

inline void dp(){
	f[0][0]=1;
	for(int i=1;i<=m;i++){
		tot+=a[i];
		pre=now,now^=1;
		for(int j=0,sum=0;j<=tot;j++){
			sum=add(sum,f[pre][j]);
			if(j>a[i]) sum=add(sum,ha-f[pre][j-a[i]-1]);
			f[now][j]=sum;
		}
	}
}

inline void calc(){
	if(!n){
		if(K<=tot) ans=f[now][K];
		else ans=0;
	}
	else for(int i=min(K,tot);i>=0;i--) ans=add(ans,f[now][i]*(ll)C(K-i+n-1,K-i)%ha);
}

int main(){
	init();
	scanf("%d%d%d",&n,&m,&K);
	for(int i=1;i<=m;i++) scanf("%d",a+i);
	n-=m;
	dp();
	calc();
	printf("%d\n",ans);
	return 0;
}