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>
using namespace std;
const int mod = 1e9 + 7;
const long long INF = 2e9 + 10;
int n, m, cl, ce, v;
long long floo[100010], eleva[100010];
int main() {
  cin >> n >> m >> cl >> ce >> v;
  floo[1] = eleva[1] = -INF;
  ++cl;
  ++ce;
  for (int i = 2; i <= cl; ++i) scanf("%lld", &floo[i]);
  floo[++cl] = INF + 1;
  for (int i = 2; i <= ce; ++i) scanf("%lld", &eleva[i]);
  eleva[++ce] = INF + 1;
  int q;
  scanf("%d", &q);
  while (q--) {
    long long x1, x2, y1, y2;
    scanf("%lld%lld%lld%lld", &x1, &y1, &x2, &y2);
    swap(x1, y1);
    swap(x2, y2);
    long long ans = INF;
    int r = lower_bound(floo + 1, floo + 1 + cl, x1) - floo;
    int l = r - 1;
    ans = min(ans, min(abs(floo[r] - x1) + abs(floo[r] - x2),
                       abs(x1 - floo[l]) + abs(x2 - floo[l])) +
                       abs(y2 - y1));
    r = lower_bound(floo + 1, floo + 1 + cl, x2) - floo;
    l = r - 1;
    ans = min(ans, min(abs(floo[r] - x1) + abs(floo[r] - x2),
                       abs(x1 - floo[l]) + abs(x2 - floo[l])) +
                       abs(y2 - y1));
    r = lower_bound(eleva + 1, eleva + 1 + ce, x1) - eleva;
    l = r - 1;
    ans = min(ans, min(abs(x1 - eleva[r]) + abs(x2 - eleva[r]),
                       abs(x1 - eleva[l]) + abs(x2 - eleva[l])) +
                       (abs(y2 - y1) + v - 1) / v);
    r = lower_bound(eleva + 1, eleva + 1 + ce, x2) - eleva;
    l = r - 1;
    ans = min(ans, min(abs(x1 - eleva[r]) + abs(x2 - eleva[r]),
                       abs(x1 - eleva[l]) + abs(x2 - eleva[l])) +
                       (abs(y2 - y1) + v - 1) / v);
    if (y1 == y2) ans = min(ans, abs(x1 - x2));
    cout << ans << endl;
  }
}
