Problem: Nauuo is a girl who loves coding.

One day she was solving a problem which requires to calculate a sum of some numbers modulo p.

She wrote the following code and got the verdict "Wrong answer".

<image>

She soon discovered the bug — the ModAdd function only worked for numbers in the range [0,p), but the numbers in the problem may be out of the range. She was curious about the wrong function, so she wanted to know the result of it.

However, the original code worked too slow, so she asked you to help her.

You are given an array a_1,a_2,…,a_n and a number p. Nauuo will make m queries, in each query, you are given l and r, and you have to calculate the results of Sum(a,l,r,p). You can see the definition of the Sum function in the pseudocode above.

Note that the integers won't overflow in the code above.

Input

The first line contains three integers n, m, p (1 ≤ n ≤ 10^6, 1 ≤ m ≤ 2 ⋅ 10^5, 1 ≤ p ≤ 10^9) — the length of the given array, the number of queries and the modulus. Note that the modulus is used only in the ModAdd function.

The second line contains n integers a_1,a_2,…,a_n (-10^9≤ a_i≤10^9) — the given array.

In the following m lines, each line contains two integers l, r (1≤ l≤ r≤ n) — you have to calculate the result of Sum(a,l,r,p).

Output

The output contains m integers to answer the queries in the given order.

Example

Input


4 5 6
7 2 -3 17
2 3
1 3
1 2
2 4
4 4


Output


-1
0
3
10
11

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
char ibuf[1 << 25], *ih = ibuf, obuf[1 << 23], *oh = obuf;
inline void read(int& x) {
  int f = 1;
  for (; !isdigit(*ih); ++ih)
    if (*ih == '-') f = -1;
  for (x = 0; isdigit(*ih); x = x * 10 + *ih++ - 48)
    ;
  x *= f;
}
inline void write(long long x) {
  if (x < 0) *oh++ = '-', x = -x;
  if (x == 0) return void(*oh++ = '0');
  static int buf[30], xb;
  for (xb = 0; x; x /= 10) buf[++xb] = x % 10;
  for (; xb;) *oh++ = buf[xb--] | 48;
}
inline unsigned int R() {
  static unsigned int seed = 19260817;
  return seed ^= seed >> 5, seed ^= seed << 17, seed ^= seed >> 13;
}
int C, C1, C2;
const int N = 1e6 + 5;
struct node {
  int l, r, fa, sz;
  long long v, t, mx, mn;
  inline void add(const long long& z) {
    v += z;
    mx += z;
    mn += z;
    t += z;
  }
} t[N];
inline void pushdown(int x) {
  if (t[x].t) t[t[x].l].add(t[x].t), t[t[x].r].add(t[x].t), t[x].t = 0;
}
inline void maintain(int x) {
  t[x].sz = t[t[x].l].sz + t[t[x].r].sz + 1;
  t[x].mx = t[x].r ? t[t[x].r].mx + t[x].t : t[x].v;
  t[x].mn = t[x].l ? t[t[x].l].mn + t[x].t : t[x].v;
}
void split(int x, int& p, int& q, const long long& v) {
  ++C1;
  if (!x) {
    p = q = 0;
    return;
  }
  if (t[x].mx < v) {
    p = x;
    q = 0;
    return;
  }
  if (v <= t[x].mn) {
    p = 0;
    q = x;
    return;
  }
  pushdown(x);
  if (t[x].v < v)
    p = x, split(t[x].r, t[p].r, q, v), t[t[p].r].fa = p, maintain(p);
  else
    q = x, split(t[x].l, p, t[q].l, v), t[t[q].l].fa = q, maintain(q);
}
int ins(int x, int y) {
  ++C;
  if (!(R() % (t[x].sz + 1))) {
    split(x, t[y].l, t[y].r, t[y].v);
    t[t[y].l].fa = t[t[y].r].fa = y;
    maintain(y);
    return y;
  }
  pushdown(x);
  int& z = t[x].v <= t[y].v ? t[x].r : t[x].l;
  z = ins(z, y);
  t[z].fa = x;
  maintain(x);
  return x;
}
int getsz(int x) { return x > 0 ? getsz(t[x].l) + getsz(t[x].r) + 1 : 0; }
int merge(int x, int y) {
  ++C;
  if (!x || !y) return x | y;
  pushdown(x);
  pushdown(y);
  if (R() % (t[x].sz + t[y].sz) < t[x].sz)
    return t[x].r = merge(t[x].r, y), t[t[x].r].fa = x, maintain(x), x;
  else
    return t[y].l = merge(x, t[y].l), t[t[y].l].fa = y, maintain(y), y;
}
int join(int x, int y) {
  ++C2;
  if (!x || !y) return x | y;
  if (t[x].mx <= t[y].mn) return merge(x, y);
  if (t[x].mn >= t[y].mx) return merge(y, x);
  if (R() % (t[x].sz + t[y].sz) < t[x].sz) swap(x, y);
  pushdown(x);
  pushdown(y);
  int p, q;
  split(x, p, q, t[y].v);
  t[y].l = join(p, t[y].l);
  t[y].r = join(q, t[y].r);
  t[t[y].l].fa = t[t[y].r].fa = y;
  maintain(y);
  return y;
}
inline long long ask(int x) {
  long long ans = t[x].v;
  for (; x = t[x].fa;) ans += t[x].t;
  return ans;
}
int n, m, p, a[N], l, r, rt, i;
vector<pair<int, int>> qu[N];
long long ans[N];
int main() {
  scanf("%d%d%d", &n, &m, &p);
  fread(ibuf, 1, 1 << 25, stdin);
  for (i = 1; i <= n; ++i) read(a[i]);
  for (i = 1; i <= m; ++i) read(l), read(r), qu[r].push_back(make_pair(l, i));
  for (i = 1; i <= n; ++i) {
    rt = ins(rt, i);
    t[rt].add(a[i]);
    split(rt, l, r, p);
    t[l].fa = t[r].fa = 0;
    t[r].add(-p);
    rt = join(l, r);
    t[rt].fa = 0;
    for (auto u : qu[i]) ans[u.second] = ask(u.first);
  }
  for (i = 1; i <= m; ++i) write(ans[i]), *oh++ = '\n';
  return fwrite(obuf, 1, oh - obuf, stdout), 0;
}