#include <bits/stdc++.h>

using namespace std;

#ifdef DGC
#include "debug.h"
#else
#define debug(...) 9715
#endif
typedef long long ll;
typedef long double ld;
typedef complex<double> point;
#define F first
#define S second

const int N = 1005, mod = 1e9+7;
int fac[N], ifac[N], dp[25][N];

int Pow(int a, int b)
{
	int x = 1;
	for (; b > 0; b >>= 1)
	{
		if (b&1) x = (ll)x * a % mod;
		a = (ll)a * a % mod;
	}
	return x;
}

int comb(int n, int k)
{
	return (ll)fac[n] * ifac[n-k] % mod * ifac[k] % mod;
}

int main()
{
	#ifdef DGC
		//freopen("a.in", "r", stdin);
		//freopen("b.out", "w", stdout);
	#endif

	ios_base::sync_with_stdio(0), cin.tie(0);

	fac[0] = ifac[0] = 1;
	for (int i = 1; i < N; ++i)
	{
		fac[i] = (ll)fac[i-1] * i % mod;
		ifac[i] = Pow(fac[i], mod-2);
	}

	int n, k;
	cin >> n >> k;
	vector<int> a(k);
	for (auto &i : a) cin >> i;

	dp[0][0] = 1;
	for (int d = 0; d < k; ++d)
		for (int i = 0; i <= n; ++i)
			for (int j = 0; i+j <= n && j <= a[d]; ++j)
				dp[d+1][i+j] = (dp[d+1][i+j] + (ll)dp[d][i] * comb(n-i, j) % mod * comb(n-j, a[d]-j)) % mod;

	cout << dp[k][n] << "\n";

	return 0;
}
