#include "iostream"
#include "climits"
#include "list"
#include "queue"
#include "stack"
#include "set"
#include "functional"
#include "algorithm"
#include "string"
#include "map"
#include "unordered_map"
#include "unordered_set"
#include "iomanip"
#include "cmath"
#include "random"
#include "bitset"
#include "cstdio"
#include "numeric"
#include "cassert"
#include "ctime"

using namespace std;

constexpr long long int MOD = 1000000007;
//constexpr int MOD = 1000000007;
//constexpr int MOD = 998244353;
//constexpr long long int MOD = 998244353;
constexpr double EPS = 1e-9;

//int N, M, K, H, W, L, R;
long long int N, M, K, H, W, L, R;

long long int power(long long int x, long long int n, long long int M) {
	long long int ret = 1;
	long long int by = x;
	while (n) {
		if (n & 1) {
			ret *= by;
			ret %= M;
		}
		by *= by;
		by %= M;
		n >>= 1;
	}
	return ret;
}

int main() {
	ios::sync_with_stdio(false);
	cin.tie(0);

	cin >> N >> K;
	vector<vector<long long int>>dp(K + 1, vector < long long int >(N + 1));
	dp[0][N] = 1;
	vector<long long int>by(N + 1, 1);
	vector<long long int>inv(N + 1, 1);
	for (int i = 2; i <= N; i++) {
		by[i] = by[i - 1] * i%MOD;
		inv[i] = power(by[i], MOD - 2, MOD);
	}
	for (int i = 0; i < K; i++) {
		int num;
		cin >> num;
		for (int j = 0; j <= N; j++) {
			for (int k = 0; k <= min(num, j); k++) {
				long long int add = by[j] * inv[k] % MOD*inv[j - k] % MOD*by[N - k] % MOD*inv[num - k] % MOD*inv[N - num] % MOD;
				dp[i + 1][j - k] += dp[i][j] * add%MOD;
				dp[i + 1][j - k] %= MOD;
			}
		}
	}
	cout << dp.back().front() << endl;
	return 0;
}