Problem: There is an infinite pond that can be represented with a number line. There are n rocks in the pond, numbered from 1 to n. The i-th rock is located at an integer coordinate a_i. The coordinates of the rocks are pairwise distinct. The rocks are numbered in the increasing order of the coordinate, so a_1 < a_2 < ... < a_n.

A robot frog sits on the rock number s. The frog is programmable. It has a base jumping distance parameter d. There also is a setting for the jumping distance range. If the jumping distance range is set to some integer k, then the frog can jump from some rock to any rock at a distance from d - k to d + k inclusive in any direction. The distance between two rocks is an absolute difference between their coordinates.

You are assigned a task to implement a feature for the frog. Given two integers i and k determine if the frog can reach a rock number i from a rock number s performing a sequence of jumps with the jumping distance range set to k. The sequence can be arbitrarily long or empty.

You will be given q testcases for that feature, the j-th testcase consists of two integers i and k. Print "Yes" if the i-th rock is reachable and "No" otherwise.

You can output "YES" and "NO" in any case (for example, strings "yEs", "yes", "Yes" and 'YES"' will be recognized as a positive answer).

Input

The first line contains four integers n, q, s and d (1 ≤ n, q ≤ 2 ⋅ 10^5; 1 ≤ s ≤ n; 1 ≤ d ≤ 10^6) — the number of rocks, the number of testcases, the starting rock and the base jumping distance parameter.

The second line contains n integers a_1, a_2, ..., a_n (1 ≤ a_i ≤ 10^6) — the coordinates of the rocks. The coordinates of the rocks are pairwise distinct. The rocks are numbered in the increasing order of distance from the land, so a_1 < a_2 < ... < a_n.

Each of the next q lines contains two integers i and k (1 ≤ i ≤ n; 1 ≤ k ≤ 10^6) — the parameters to the testcase.

Output

For each of the testcases print an answer. If there is a sequence of jumps from a rock number s to a rock number i with the jumping distance range set to k, then print "Yes". Otherwise, print "No".

Examples

Input


7 4 4 5
1 5 10 13 20 22 28
4 1
7 2
7 3
3 2


Output


Yes
No
Yes
Yes


Input


10 8 6 11
1 2 4 7 8 9 11 13 19 20
2 13
5 8
8 1
6 15
1 15
2 7
7 6
8 9


Output


Yes
Yes
No
Yes
Yes
Yes
Yes
Yes


Input


6 9 6 6
1 2 4 9 18 19
2 17
1 18
5 4
2 11
5 17
6 8
4 3
3 3
6 6


Output


Yes
Yes
Yes
Yes
Yes
Yes
No
No
Yes


Input


4 1 1 10
1 8 10 19
2 1


Output


Yes

Note

Explanation of the first example:

In the first testcase the destination rock is the same as the starting rock, thus no jumps are required to reach it.

In the second testcase the frog can jump any distance in the range [5 - 2; 5 + 2]. Thus, it can reach rock number 5 (by jumping 7 to the right) and rock number 3 (by jumping 3 to the left). From rock number 3 it can reach rock number 2 (by jumping 5 to the left). From rock number 2 it can reach rock number 1 (by jumping 4 to the left). However, there is no way to reach rock number 7.

In the third testcase the frog can jump any distance in the range [5 - 3; 5 + 3]. Thus, it can reach rock number 7 by jumping to rock 5 first and to 7 afterwards.

The fourth testcase is shown in the explanation for the second testcase.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
#define rep(i, x, y) for (int i = x; i <= y; i++)
#define per(i, x, y) for (int i = x; i >= y; i--)
#define vi vector<int>
#define pb push_back
#define mkp make_pair
#define fi first
#define se second
#define gc getchar
#define db double
#define bs bitset
#define ll long long
#define lll __int128
#define ui unsigned int
#define ull unsigned long long
#define pii pair<int, int>
#define clr(x) memset(x, 0, sizeof(x))
#define lbit(x) ((x) & (-(x)))
#define heap priority_queue<int>
#define SI set<int>::iterator
#define MSI multiset<int>::iterator
#define Go(i, x, v) for (int i = 0, x = (i == v.size() ? 0 : v[i]); i < v.size(); i++, x = (i == v.size() ? 0 : v[i]))
using namespace std;
template <class T> void rd(T &x) {
  x = 0; char ch = gc(); int f = 1;
  for (; !isdigit(ch); ch = gc()) if (ch == '-') f = -1;
  for (; isdigit(ch); ch = gc()) x = (x << 1) + (x << 3) + ch - '0';
  x *= f;
}
template <class T> void cmin(T &x, T y) { x = min(x, y); }
template <class T> void cmax(T &x, T y) { x = max(x, y); }
const int mod = 998244353;
template <class T> T add(T x, T y) { return x + y >= mod ? x + y - mod : x + y; }
template <class T> T sub(T x, T y) { return x - y < 0 ? x - y + mod : x - y; }
template <class T> T mul(T x, T y) { return 1ll * x * y % mod; }
template <class T> void Add(T &x, T y) { x = (x + y >= mod ? x + y - mod : x + y); }
template <class T> void Sub(T &x, T y) { x = (x - y < 0 ? x - y + mod : x - y); }
template <class T> void Mul(T &x, T y) { x = 1ll * x * y % mod; }
template <class T> T qpow(T a, T b = mod - 2) { T ret = 1; for (; b; b >>= 1, Mul(a, a)) if (b & 1) Mul(ret, a); return ret; }
/* ============ Header Template ============ */

const int N = 2e5 + 10, M = 1e6 + 10, inf = 0x3f3f3f3f;
int n, Q, s, d, a[N], mn[N], bel[N], idx[M];
vi id[N];
vector<pii> g[N];

struct edge {
  int x, y, w;
};
bool operator < (edge a, edge b) {
  if (a.w ^ b.w) return a.w < b.w;
  return min(a.x, a.y) == min(b.x, b.y) ? max(a.x, a.y) < max(b.x, b.y) : min(a.x, a.y) < min(b.x, b.y);
}

bool unite(int x, int y) {
  x = bel[x], y = bel[y];
  if (x == y) return 0;
  if (id[x].size() < id[y].size()) swap(x, y);
  Go(i, w, id[y]) {
    bel[w] = x;
    id[x].pb(w);
  } id[y].clear();
  return 1;
}

void dfs(int x, int fa, int d) {
  mn[x] = d;
  for (int i = 0; i < g[x].size(); i++) {
    int y = g[x][i].fi, w = g[x][i].se;
    if (y != fa) dfs(y, x, max(d, w));
  }
}

void boruvka() {
  set<int> pos;
  rep(i, 1, n) pos.insert(a[i]);
  int cnt = n;
  while (cnt > 1) {
    // printf("cnt = %d\n", cnt);
    vector<edge> e;
    rep(i, 1, n) if (id[i].size()) {
      // printf("%d ..\n", i);
      Go(o, x, id[i]) pos.erase(a[x]);
      edge cur = (edge){-1, -1, inf};
      Go(o, x, id[i]) {
        rep(j, -1, 1) if (j) {
          int dx = d * j;
          SI it = pos.lower_bound(a[x] + dx);
          if (it != pos.end())
            cmin(cur, (edge){x, idx[*it], abs(d - abs(a[x] - (*it)))});
          if (it != pos.begin()) {
            --it;
            cmin(cur, (edge){x, idx[*it], abs(d - abs(a[x] - (*it)))});
          }
        }
      }
      Go(o, x, id[i]) pos.insert(a[x]);
      e.pb(cur);
      // printf("%d : %d %d %d\n", i, cur.x, cur.y, cur.w);
    }
    for (int i = 0; i < e.size(); i++) {
      edge t = e[i];
      if (unite(t.x, t.y)) {
        --cnt;
        g[t.x].pb(mkp(t.y, t.w));
        g[t.y].pb(mkp(t.x, t.w));
      }
    }
  }
}

signed main() {
  rd(n), rd(Q), rd(s), rd(d);
  rep(i, 1, n) rd(a[i]), idx[a[i]] = i, bel[i] = i, id[i].pb(i);
  boruvka();
  dfs(s, -1, 0);
  while (Q--) {
    int i, k; rd(i), rd(k);
    puts(mn[i] <= k ? "Yes" : "No");
  }
  return 0;
}