#include "bits/stdc++.h"
using namespace std;
typedef long long ll;
const int INF = (1<<30);
const ll INFLL = (1ll<<60);
const ll MOD = (ll)(1e9+7);

#define l_ength size

void mul_mod(ll& a, ll b){
	a *= b;
	a %= MOD;
}

void add_mod(ll& a, ll b){
	a = (a<MOD)?a:(a-MOD);
	b = (b<MOD)?b:(b-MOD);
	a += b;
	a = (a<MOD)?a:(a-MOD);
}

ll a[32],fc[1024]={1ll},nv[1024]={1ll},dp[32][1024];

ll rwpw(ll a, ll p=MOD-2){
	ll ret = 1ll;
	while(p){
		if(p%2){
			mul_mod(ret,a);
		}
		mul_mod(a,a);
		p /= 2;
	}
	return ret;
}

int main(void){
	int n,d,i,j,k;
	ll tmp;
	cin >> n >> d;
	for(i=0; i<n; ++i){
		fc[i+1] = fc[i]*(i+1)%MOD;
		nv[i+1] = rwpw(fc[i+1]);
	}
	for(i=0; i<d; ++i){
		cin >> a[i];
	}
	dp[0][0] = 1ll;
	for(i=0; i<=d; ++i){
		for(j=0; j<=n; ++j){
			for(k=0; k<=a[i]; ++k){
				if(j+k>n){
					continue;
				}
				tmp = fc[n-k];
				mul_mod(tmp,nv[a[i]-k]);
				mul_mod(tmp,nv[n-a[i]]);
				mul_mod(tmp,nv[k]);
				mul_mod(tmp,dp[i][j]);
				add_mod(dp[i+1][j+k],tmp);
			}
		}
	}
	for(i=0; i<n; ++i){
		mul_mod(dp[d][n],(i+1));
	}
	cout << dp[d][n] << endl;
	return 0;
}
