2 10 1 1 1
1
10
1
1 5 1 8
 4 4 1 0 1
1

1
1 2 1 4
 2 5 1 0 1
2

1
1 4 1 5
 10 10 1 8 4
10
2 3 4 5 6 7 8 9
10
1 1 3 1
2 1 7 1
1 1 9 1
7 1 4 1
10 1 7 1
2 1 7 1
3 1 2 1
5 1 2 1
10 1 5 1
6 1 9 1
 2 4 1 1 1
1
4
1
1 2 1 3
 2 2 0 1 1

1
1
1 2 2 2
 1000 1000 1 1 10
1
2
1
1 900 1 1000
 5 5 1 1 1
3
2
1
1 5 1 1
 2 4 1 1 1
1
2
1
2 3 2 4
 4 4 1 0 1
4

5
1 1 2 2
1 3 2 2
3 3 4 3
3 2 2 2
1 2 2 3

3
 2
 1
 3
4
4
3
3
4
3
3
4
3
 1
 3
 100
 4
 1
 6
4
3
5
4

#include <bits/stdc++.h>
const int INF = 1000000000;
using namespace std;
const int maxn = 100005;
int st[maxn], el[maxn];
int main() {
  int n, m, cl, cv, v, q;
  scanf("%d%d%d%d%d", &n, &m, &cl, &cv, &v);
  for (int i = 1; i <= cl; i++) scanf("%d", &st[i]);
  for (int i = 1; i <= cv; i++) scanf("%d", &el[i]);
  scanf("%d", &q);
  while (q--) {
    int x, y, x1, y1;
    long long ans = INF;
    scanf("%d%d%d%d", &x, &y, &x1, &y1);
    if (x == x1) {
      printf("%d\n", abs(y - y1));
      continue;
    }
    int p = lower_bound(st + 1, st + 1 + cl, y) - st;
    if (p <= cl)
      ans = ans > abs(st[p] - y) + abs(st[p] - y1) + abs(x - x1)
                ? abs(st[p] - y) + abs(st[p] - y1) + abs(x - x1)
                : ans;
    if (p > 1)
      ans = ans > abs(st[p - 1] - y) + abs(st[p - 1] - y1) + abs(x - x1)
                ? abs(st[p - 1] - y) + abs(st[p - 1] - y1) + abs(x - x1)
                : ans;
    int t = lower_bound(el + 1, el + 1 + cv, y) - el;
    if (t <= cv)
      ans = ans > abs(el[t] - y) + abs(el[t] - y1) + ceil(abs(x - x1) * 1.0 / v)
                ? abs(el[t] - y) + abs(el[t] - y1) + ceil(abs(x - x1) * 1.0 / v)
                : ans;
    if (t > 1)
      ans = ans > abs(el[t - 1] - y) + abs(el[t - 1] - y1) +
                        ceil(abs(x - x1) * 1.0 / v)
                ? abs(el[t - 1] - y) + abs(el[t - 1] - y1) +
                      ceil(abs(x - x1) * 1.0 / v)
                : ans;
    p = lower_bound(st + 1, st + 1 + cl, y1) - st;
    if (p <= cl)
      ans = ans > abs(st[p] - y) + abs(st[p] - y1) + abs(x - x1)
                ? abs(st[p] - y) + abs(st[p] - y1) + abs(x - x1)
                : ans;
    if (p > 1)
      ans = ans > abs(st[p - 1] - y) + abs(st[p - 1] - y1) + abs(x - x1)
                ? abs(st[p - 1] - y) + abs(st[p - 1] - y1) + abs(x - x1)
                : ans;
    t = lower_bound(el + 1, el + 1 + cv, y1) - el;
    if (t <= cv)
      ans = ans > abs(el[t] - y) + abs(el[t] - y1) + ceil(abs(x - x1) * 1.0 / v)
                ? abs(el[t] - y) + abs(el[t] - y1) + ceil(abs(x - x1) * 1.0 / v)
                : ans;
    if (t > 1)
      ans = ans > abs(el[t - 1] - y) + abs(el[t - 1] - y1) +
                        ceil(abs(x - x1) * 1.0 / v)
                ? abs(el[t - 1] - y) + abs(el[t - 1] - y1) +
                      ceil(abs(x - x1) * 1.0 / v)
                : ans;
    printf("%I64d\n", ans);
  }
  return 0;
}
