#include <bits/stdc++.h>
using namespace std;
struct CD {
  double r, i;
  CD(double r = 0, double i = 0) : r(r), i(i) {}
  double real() const { return r; }
  void operator/=(const int c) { r /= c, i /= c; }
};
CD operator*(const CD& a, const CD& b) {
  return CD(a.r * b.r - a.i * b.i, a.r * b.i + a.i * b.r);
}
CD operator+(const CD& a, const CD& b) { return CD(a.r + b.r, a.i + b.i); }
CD operator-(const CD& a, const CD& b) { return CD(a.r - b.r, a.i - b.i); }
const int MX = (1 << 19);
const double PI = acos(-1.0);
CD cp1[MX + 9], cp2[MX + 9];
int R[MX + 9];
void fft(CD* a, int n, bool inverso) {
  for (int i = 0; i < n; i++) {
    if (R[i] < i) {
      swap(a[R[i]], a[i]);
    }
  }
  for (int m = 2; m <= n; m <<= 1ll) {
    double z = 2 * PI / m * (inverso ? -1 : 1);
    CD wi = CD(cos(z), sin(z));
    for (int j = 0; j < n; j += m) {
      CD w(1);
      for (int k = j, k2 = j + m / 2; k2 < j + m; k++, k2++) {
        CD u = a[k];
        CD v = a[k2] * w;
        a[k] = u + v;
        a[k2] = u - v;
        w = w * wi;
      }
    }
  }
  if (inverso) {
    for (int i = 0; i < n; i++) {
      a[i] /= n;
    }
  }
  return;
}
vector<int> multiply(vector<int>& p1, vector<int>& p2) {
  int n = p1.size() + p2.size() + 1;
  int m = 1, cnt = 0;
  while (m <= n) {
    m += m;
    cnt++;
  }
  for (int i = 0; i < m; i++) {
    R[i] = 0;
    for (int j = 0; j < cnt; j++) {
      R[i] = (R[i] << 1ll) | ((i >> j) & 1);
    }
  }
  for (int i = 0; i < m; i++) {
    cp1[i] = cp2[i] = 0;
  }
  for (int i = 0; i < p1.size(); i++) {
    cp1[i] = p1[i];
  }
  for (int i = 0; i < p2.size(); i++) {
    cp2[i] = p2[i];
  }
  fft(cp1, m, false);
  fft(cp2, m, false);
  for (int i = 0; i < m; i++) {
    cp1[i] = cp1[i] * cp2[i];
  }
  fft(cp1, m, true);
  vector<int> result;
  n -= 2;
  for (int i = 0; i < n; i++) {
    result.push_back(int(cp1[i].real() + (cp1[i].real() > 0 ? 0.5 : -0.5)));
  }
  return result;
}
bool vis[MX];
bool used[MX];
void build(int l, int r) {
  if (l >= r) return;
  int med = (l + r) >> 1;
  vector<int> A, B;
  A.resize(med - l + 1, 0);
  B.resize(r - med + 1, 0);
  for (int i = med; i >= l; i--) {
    if (vis[i]) A[med - i] = 1;
  }
  for (int i = med; i <= r; i++) {
    if (vis[i]) B[i - med] = 1;
  }
  vector<int> AB = multiply(A, B);
  for (int i = 0; i < AB.size(); i++) {
    if (AB[i] != 0) used[i] = true;
  }
  build(l, med);
  build(med + 1, r);
  return;
}
int fact[MX];
void extSieve() {
  for (int i = 2; i * i < MX; ++i) {
    if (fact[i] == 0) {
      for (int j = i * i; j < MX; j += i)
        if (fact[j] == 0) fact[j] = i;
    }
  }
  for (int i = 2; i < MX; ++i) {
    if (fact[i] == 0) {
      fact[i] = i;
    }
  }
}
void solve(int n, int h) {
  vector<int> v;
  v.push_back(1);
  while (n > 1) {
    int f = fact[n];
    int exp = 0;
    while (n % f == 0) {
      n /= f;
      exp++;
    }
    int sz = v.size();
    int aux = f;
    for (int i = 0; i < exp; i++) {
      for (int j = 0; j < sz; j++) {
        v.push_back(aux * v[j]);
      }
      if (i != exp - 1) aux *= f;
    }
  }
  int ans = -1;
  for (int xx : v) {
    if (xx > h && used[xx - h]) {
      ans = max(xx * 2, ans);
    }
  }
  cout << ans << " ";
  return;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  cout.tie(NULL);
  cout.precision(10);
  cout << fixed;
  extSieve();
  int n, D, h, x;
  cin >> n >> D >> h;
  for (int i = 0; i <= n; i++) {
    cin >> x;
    vis[x] = true;
  }
  build(0, D);
  int q;
  cin >> q;
  for (int i = 0; i < q; i++) {
    cin >> x;
    x /= 2;
    solve(x, h);
  }
  return 0;
}
