#include <bits/stdc++.h>
using namespace std;
long long p = 1e9 + 7;
long long fpow(long long a, int exp, long long mod) {
  if (exp == 0) {
    return 1;
  }
  if (exp % 2 == 0) {
    long long x = fpow(a, exp / 2, mod);
    return (x * x) % mod;
  }
  return (a * fpow(a, exp - 1, mod)) % mod;
}
long long inv(long long x) { return (fpow(x, p - 2, p) + p) % p; }
long long n;
long long k;
long long a[5050];
long long sigma[5050];
int main() {
  ios_base::sync_with_stdio(0);
  cin >> n >> k;
  for (int i = 0; i < n; i++) {
    cin >> a[i];
  }
  sigma[0] = 1;
  for (int i = 0; i <= n; i++) {
    for (int j = n; j > 0; j--) {
      sigma[j] = (sigma[j] + a[i] * sigma[j - 1]) % p;
    }
  }
  long long res = 0;
  long long var = 1;
  for (int i = n - 1; i >= 0; i--) {
    var = (var * (k - (n - i - 1))) % p;
    if (k - n + i >= 0)
      res +=
          (var * fpow(-1, n + 1 - i, p) * fpow(n, k - n + i, p)) % p * sigma[i];
    res %= p;
  }
  res *= fpow(inv(n), k, p);
  res %= p;
  cout << (res + p) % p;
}
