#include <bits/stdc++.h>
using namespace std;
map<int, int> memo;
vector<long long> w;
int phiF(int n) {
  auto a = memo.find(n);
  if (a != memo.end()) return a->second;
  int nn = n;
  if (n <= 1) {
    memo[n] = n;
    return n;
  }
  int ans = n;
  for (int i = 2; i * i <= n; i++) {
    if (n % i == 0) {
      while (n % i == 0) n /= i;
      ans -= ans / i;
    }
  }
  if (n > 1) ans -= ans / n;
  memo[nn] = ans;
  return ans;
}
long long expMod(long long a, long long exp, long long mod) {
  if (exp == 0) return 1;
  long long pa = expMod(a, exp / 2, mod);
  pa = (pa * pa) % mod;
  if (exp % 2 == 1) pa = (pa * a) % mod;
  return pa;
}
bool isBig(int l, int r, long long lim, long long &res) {
  if (w[l] == 1) {
    res = 1;
    return false;
  }
  if (l == r) {
    res = w[l];
    return w[l] >= lim;
  }
  if (w[l] >= lim) return true;
  long long wp = 1;
  int e = 0;
  while (wp < lim) {
    wp *= w[l];
    e++;
  }
  long long rr;
  if (isBig(l + 1, r, e, rr)) return true;
  res = 1;
  for (int i = 0; i < rr; i++) res *= w[l];
  return false;
}
long long solve(int l, int r, long long mod) {
  if (l == r) return w[l] % mod;
  if (mod == 1) return 0;
  long long res;
  long long lim = 0, mm = mod;
  while (mm) {
    mm >>= 1;
    lim++;
  }
  if (isBig(l + 1, r, lim, res)) {
    long long m2 = phiF(mod);
    long long e = solve(l + 1, r, m2) + m2;
    return expMod(w[l], e, mod);
  } else {
    return expMod(w[l], res, mod);
  }
}
int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int n, m;
  cin >> n >> m;
  w.resize(n);
  for (int i = 0; i < n; i++) cin >> w[i];
  int q;
  cin >> q;
  while (q--) {
    int l, r;
    cin >> l >> r;
    l--;
    r--;
    auto ans = solve(l, r, m);
    cout << ans << endl;
  }
  return 0;
}
