/* cerberus97 - Hanit Banga */

#include <iostream>
#include <iomanip>
#include <cassert>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <vector>
#include <algorithm>

using namespace std;

#define pb push_back
#define fast_cin() ios_base::sync_with_stdio(false); cin.tie(NULL)

typedef long long ll;
typedef long double ld;
typedef pair <int, int> pii;
typedef pair <ll, ll> pll;

const int N = 1e3 + 10, K = 30, mod = 1e9 + 7;

ll add(ll x, ll y) {
	x += y;
	if (x >= mod) {
		x -= mod;
	}
	return x;
}

ll sub(ll x, ll y) {
	x -= y;
	if (x < 0) {
		x += mod;
	}
	return x;
}

ll mul(ll x, ll y) {
	x *= y;
	return x % mod;
}

ll pwr(ll x, ll e) {
	ll res = 1;
	while (e) {
		if (e & 1) {
			res = mul(res, x);
		}
		x = mul(x, x);
		e /= 2;
	}
	return res;
}

int a[K];
ll choose[N][N], dp[N][K];

int main() {
	fast_cin();
	int n, k;
	cin >> n >> k;
	for (int i = 1; i <= k; ++i) {
		cin >> a[i];
	}
	for (int i = 0; i <= n; ++i) {
		choose[i][i] = choose[i][0] = 1;
		for (int j = 1; j < i; ++j) {
			choose[i][j] = add(choose[i - 1][j], choose[i - 1][j - 1]);
		}
	}
	for (int j = 0; j <= k; ++j) {
		dp[0][0] = 1;
	}
	for (int i = 0; i <= n; ++i) {
		for (int j = 1; j <= k; ++j) {
			for (int c = 0; c <= min(i, a[j]); ++c) {
				dp[i][j] = add(dp[i][j], mul(mul(choose[i][c], choose[n - c][a[j] - c]), dp[i - c][j - 1]));
			}
		}
	}
	cout << dp[n][k] << endl;
}