#include<bits/stdc++.h>
using namespace std;
// #define LOCAL
#ifdef LOCAL
    #define eprintf(...) fprintf(stderr, __VA_ARGS__)
#else
    #define eprintf(...) 42
#endif // LOCAL
int const nax = 3e5 + 10, mod = 1e9 + 7;
void plusle(int &a,int b){a+=b;if(a>=mod)a-=mod; return;}
void minun(int &a,int b){a-=b; if(a<0)a+=mod; return;}
int add(int a,int b) {a+=b;return a>=mod?a-mod:a;}
int sub(int a,int b) {a-=b;return a<0?a+mod:a;}
int mul(int a,int b) {return (int)((long long)a*b%mod);}
int power(int a,int b) {int res=1;while (b>0) {if (b&1) {res=mul(res,a);}a=mul(a,a);b>>=1;}return res;}

int const Mx = 3000;
int fac[Mx],invf[Mx];
void prefact(){
    fac[0]=1;
    for(int i=1;i< Mx;i++)
        fac[i]=mul(fac[i-1],i);
    invf[Mx-1]=power(fac[Mx -1],mod-2);
    for(int i=Mx -2 ;i>=0;i--){
        invf[i]=mul(invf[i+1],i+1);
    }

}
int ncr(int N,int R){
    if(N < 0) return 0;
    if(R > N || R < 0) return 0;
    if(N == 0 || N == R)
        return 1;
    return mul(fac[N],mul(invf[N-R],invf[R]));
}
int n, k, a[nax];
int f[22][1050];

int main(){
    prefact();
    scanf("%d %d", &n, &k);
    for(int i = 1 ;i <= k ;i ++) {
        scanf("%d", a  + i);
    }
    for (int i = 0 ; i <= n ; ++ i) {
        f[1][i] = mul(ncr(n , i), ncr(n - i, a[1] - i));
    }
    for (int j = 2 ; j <= k ; ++ j) {
        for (int i = 0 ; i <= n ; ++ i) {
            for (int from = 0 ; from <= n - i ; ++ from) {
                int toAdd = mul(f[j - 1][from], mul(ncr(n - from, i), ncr(n - i, a[j] - i)));
                plusle(f[j][from + i] , toAdd);
            }
        }
    }
    printf("%d\n", f[k][n]);
}
/*
    Good Luck
        -Lucina
*/
