#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int max_n=1005;
const int mod=1e9+7;
inline int Add(int x,int y){
    return (x+=y)>=mod?x-mod:x;
}
inline int Sub(int x,int y){
    return (x-=y)>=0?x:x+mod;
}
inline int Mul(int x,int y){
    return 1ll*x*y%mod;
}
inline int Pow(int x,int y){
    int res=1;
    while(y){
        if(y&1)res=1ll*res*x%mod;
        y>>=1;
        x=1ll*x*x%mod;
    }
    return res;
}
int n,k,a[max_n];
int f[max_n],inv[max_n];
int C(int n,int m){
    if(n<m||m<0)return 0;
    return 1ll*f[n]*inv[m]%mod*inv[n-m]%mod;
}
int dp[max_n];
int main(){
    scanf("%d%d",&n,&k);
    for(int i=1;i<=k;i++)scanf("%d",a+i);
    f[0]=1;
    for(int i=1;i<=n;i++)f[i]=Mul(f[i-1],i);
    inv[n]=Pow(f[n],mod-2);
    for(int i=n-1;i>=0;i--)inv[i]=Mul(inv[i+1],i+1);
    dp[0]=1;
    for(int i=1;i<=k;i++){
        for(int j=n;j>=0;j--){
            int res=0;
            for(int l=0;l<=min(j,a[i]);l++){
                int x=dp[j-l];
                x=Mul(x,inv[l]);
                x=Mul(x,C(n-l,a[i]-l));
                res=Add(res,x);
            }
            dp[j]=res;
        }
    }
    printf("%d\n",Mul(f[n],dp[n]));
    return 0;
}