#include <bits/stdc++.h>
using namespace std;
long long int mul(long long int a, long long int b) {
  return ((a % int(1e9 + 7)) * (b % int(1e9 + 7))) % int(1e9 + 7);
}
long long int add(long long int a, long long int b) {
  return ((a % int(1e9 + 7)) + (b % int(1e9 + 7))) % int(1e9 + 7);
}
long long int sub(long long int a, long long int b) {
  long long int ans = ((a % int(1e9 + 7)) - (b % int(1e9 + 7))) % int(1e9 + 7);
  ans %= int(1e9 + 7);
  ans = (ans + int(1e9 + 7)) % int(1e9 + 7);
  return ans;
}
long long int mpow(long long int a, long long int b) {
  long long int ans = 1;
  long long int po = a;
  while (b != 0) {
    if (b % 2) {
      ans = ((ans % int(1e9 + 7)) * (po % int(1e9 + 7))) % int(1e9 + 7);
    }
    po = ((po % int(1e9 + 7)) * (po % int(1e9 + 7))) % int(1e9 + 7);
    b /= 2;
  }
  return ans;
}
long long int fact[1010];
void pre() {
  fact[0] = 1;
  for (long long int i = 0; i < 1000; i++) {
    fact[i + 1] = (fact[i] * (i + 1)) % int(1e9 + 7);
  }
}
void solve() {
  long long int n, m;
  cin >> n >> m;
  vector<long long int> a(m);
  for (long long int i = 0; i < m; i++) cin >> a[i];
  sort(a.begin(), a.end());
  long long int ans = fact[n - m];
  for (long long int i = 1; i < m; i++) {
    if (a[i] == a[i - 1] + 1) continue;
    ans = (ans * mpow(2, (a[i] - a[i - 1] - 2))) % int(1e9 + 7);
    ans = (ans * mpow(fact[a[i] - a[i - 1] - 1], int(1e9 + 7) - 2)) %
          int(1e9 + 7);
  }
  ans = (ans * mpow(fact[a[0] - 1], int(1e9 + 7) - 2)) % int(1e9 + 7);
  ans = (ans * mpow(fact[n - a[m - 1]], int(1e9 + 7) - 2)) % int(1e9 + 7);
  cout << ans << '\n';
}
signed main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  cout.tie(0);
  long long int t;
  pre();
  solve();
}
