#include <bits/stdc++.h>
using namespace std;
const int N = 3e6 + 10, MOD = 1e9 + 7;
int add(long long a, int b) { return (a + b + 2LL * MOD) % MOD; }
int mul(long long a, int b) { return (a * b) % MOD; }
void sadd(int &a, int b) { a = add(a, b); }
void smul(int &a, int b) { a = mul(a, b); }
int bp(int a, int b) {
  if (b == 0) return 1;
  int x = bp(a, b / 2);
  smul(x, x);
  return b & 1 ? mul(x, a) : x;
}
int inv(int a) { return bp(a, MOD - 2); }
int fact[N], inv_fact[N];
void prep(int n) {
  inv_fact[0] = fact[0] = 1;
  for (int i = 1; i < n; i++) fact[i] = mul(fact[i - 1], i);
  inv_fact[n - 1] = inv(fact[n - 1]);
  for (int i = n - 2; i >= 1; i--) inv_fact[i] = mul(inv_fact[i + 1], i + 1);
}
int C(int m, int n) {
  if (m > n) return 0;
  return mul(fact[n], mul(inv_fact[n - m], inv_fact[m]));
}
int dp[N][3];
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  prep(N);
  int three = inv(3);
  int n, q;
  cin >> n >> q;
  for (int i = 0; i < 3; i++) dp[0][i] = n;
  for (int i = 1; i <= 3 * n; i++) {
    int z = C(i + 1, 3 * n);
    int x = dp[i - 1][1];
    int y = dp[i - 1][0];
    dp[i][0] = mul(add(z, -add(x, mul(2, y))), three);
    dp[i][1] = dp[i][0] + y;
    dp[i][2] = dp[i][1] + x;
  }
  while (q--) {
    int x;
    cin >> x;
    cout << add(dp[x][0], C(x, 3 * n)) << "\n";
  }
}
