#include<bits/stdc++.h>
using namespace std;

typedef long long LL;

const int N=1000,M=20,mod=1000000007;

int add(int a,int b,int p=mod){return a+b>=p?a+b-p:a+b;}
int sub(int a,int b,int p=mod){return a-b<0?a-b+p:a-b;}
int mul(int a,int b,int p=mod){return (LL)a*b%p;}
void sadd(int &a,int b,int p=mod){a=add(a,b,p);}
void ssub(int &a,int b,int p=mod){a=sub(a,b,p);}
void smul(int &a,int b,int p=mod){a=mul(a,b,p);}

int n,m,a[N+9];

void into(){
  scanf("%d%d",&n,&m);
  for (int i=1;i<=m;++i)
	scanf("%d",&a[i]);
}

int inv[N+9],fac[N+9],ifac[N+9];

void Get_inv(){
  inv[1]=1;
  fac[0]=fac[1]=1;
  ifac[0]=ifac[1]=1;
  for (int i=2;i<=n;++i){
	inv[i]=mul(mod-mod/i,inv[mod%i]);
	fac[i]=mul(fac[i-1],i);
	ifac[i]=mul(ifac[i-1],inv[i]);
  }
}

int c[N+9][N+9];

void Get_c(){
  c[0][0]=1;
  for (int i=1;i<=n;++i){
	c[i][0]=c[i][i]=1;
	for (int j=1;j<i;++j) c[i][j]=add(c[i-1][j-1],c[i-1][j]);
  }
}

int dp[M+9][N+9];

void Get_dp(){
  dp[0][0]=1;
  for (int i=1;i<=m;++i)
	for (int j=0;j<=n;++j)
	  for (int k=0;k<=j&&k<=a[i];++k)
		sadd(dp[i][j],mul(mul(c[n-k][a[i]-k],ifac[k]),dp[i-1][j-k]));
}

void work(){
  Get_inv();
  Get_c();
  Get_dp();
}

void outo(){
  printf("%d\n",mul(dp[m][n],fac[n]));
}

int main(){
  into();
  work();
  outo();
  return 0;
}