#include <bits/stdc++.h>
using namespace std;
const long long inf = 1000000000;
int n, m, w[100005], tot, a[105], nx[100005];
bool flag;
int ksm(int x, int y, int mo) {
  int ans = 1;
  bool f1 = 0;
  while (y) {
    if (y % 2) {
      flag += ((long long)ans * x > inf) | f1;
      if (flag)
        ans = (long long)ans * x % mo;
      else
        ans *= x;
    }
    f1 += ((long long)x * x > inf);
    if (f1 || flag)
      x = (long long)x * x % mo;
    else
      x *= x;
    y = y / 2;
  }
  return ans;
}
int phi(int n) {
  int result = n;
  for (int i = 2; i * i <= n; i++) {
    if (n % i == 0) {
      while (n % i == 0) n /= i;
      result -= result / i;
    }
  }
  if (n > 1) result -= result / n;
  return result;
}
int solve(int d, int l, int r) {
  int mo = a[d];
  if (l == r) {
    if (w[l] > inf) {
      flag = 1;
      return w[l] % mo;
    } else
      return w[l];
  }
  if (d == tot && r - l + 1 > 10) {
    flag = 1;
    return 0;
  }
  int phi = a[d + 1];
  int x = solve(min(tot, d + 1), l + 1, r);
  if (flag)
    x += phi;
  else if (x >= phi)
    x = x % phi + phi;
  return ksm(w[l], x, mo);
}
int main() {
  scanf("%d%d", &n, &m);
  a[0] = m;
  while (a[tot] > 1) tot++, a[tot] = phi(a[tot - 1]);
  a[tot + 1] = 1;
  for (int i = 1; i <= n; i++) scanf("%d", &w[i]);
  int q;
  scanf("%d", &q);
  nx[n + 1] = n + 1;
  for (int i = n; i >= 1; i--) nx[i] = (w[i] == 1) ? i : nx[i + 1];
  while (q--) {
    int l, r;
    scanf("%d%d", &l, &r);
    int tmp = nx[l];
    if (tmp <= r) r = tmp - 1;
    if (l > r)
      puts(m > 1 ? "1" : "0");
    else
      flag = 0, printf("%d\n", solve(0, l, r) % m);
  }
  return 0;
}
