#include <bits/stdc++.h>

using namespace std;
#define int long long
const int MOD = 1000000007;
template <typename Int, Int MOD, int N>
struct comb_util {
	std::array<Int, N + 1> fc, ifc;

	comb_util() {
		fc[0] = 1;
		for (int i = 1; i <= N; i++) fc[i] = fc[i - 1] * i % MOD;
		ifc[N] = inv(fc[N]);
		for (int i = N - 1; i >= 0; i--) ifc[i] = ifc[i + 1] * (i + 1) % MOD;
	}

	Int fact(Int n) { return fc[n]; }

	Int inv_fact(Int n) { return ifc[n]; }

	Int inv(Int n) { return pow(n, MOD - 2); }

	Int pow(Int n, Int a) {
		Int res = 1, exp = n;
		for (; a; a /= 2) {
			if (a & 1) res = res * exp % MOD;
			exp = exp * exp % MOD;
		}
		return res;
	}

	Int perm(Int n, Int r) {
		if (r < 0 || n < r)
			return 0;
		else
			return fc[n] * ifc[n - r] % MOD;
	}

	Int binom(Int n, Int r) {
		if (n < 0 || r < 0 || n < r) return 0;
		return fc[n] * ifc[r] % MOD * ifc[n - r] % MOD;
	}

	Int homo(Int n, Int r) {
		if (n < 0 || r < 0) return 0;
		return r == 0 ? 1 : binom(n + r - 1, r);
	}
};

using comb = comb_util<long long, 1000000007, 300000>;
comb cb;

int N, K;
vector<int> A;
vector<vector<int> > dp;
int rec(const int day, const int sum) {
	//cerr << day<<" "<<sum<<endl;
	if (dp[day][sum] == -1) {
		dp[day][sum] = 0;
		int mx = min(A[day], N - sum);
		for (int i = 0; i <= mx; i++) {
			int coef = cb.binom(N - i, A[day] - i);
			//coef = (coef * i) % MOD;
			coef = (coef * cb.binom(N - sum, i)) % MOD;
			dp[day][sum] = (dp[day][sum] + coef * rec(day + 1, sum + i)) % MOD;
		}
		//cerr << day << " " << sum << " " << dp[day][sum] << endl;
	}
	return dp[day][sum];
}
signed main() {
	cin.tie(0);
	ios::sync_with_stdio(false);
	cin >> N >> K;
	A.resize(K);
	dp.resize(K + 1, vector<int>(N + 1, -1));
	for (int i = 0; i < N; i++) {
		dp[K][i] = 0;
	}
	dp[K][N] = 1;
	for (int i = 0; i < K; i++) {
		cin >> A[i];
	}
	cout << rec(0, 0) << endl;
}