#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <limits>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

using namespace std;

using Int = long long;

template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> ostream &operator<<(ostream &os, const vector<T> &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; }
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }
#define COLOR(s) ("\x1b[" s "m")


using Double = double;
constexpr Double INF = 1e+100;

int N, Q;
Int V;
vector<Int> P;
vector<Int> W;
vector<Int> X0, Y0, X1, Y1;

int indexOf(Int x) {
  return (upper_bound(P.begin(), P.end(), x) - P.begin()) - 1;
}

namespace brute {
vector<Int> fs;
Int calc(Int x) {
  const int i = indexOf(x);
  if (i < 0) return 0;
  if (i >= N) return fs[N];
  return fs[i] + W[i] * (x - P[i]);
}
vector<Double> run() {
cerr<<"[brute::run]"<<endl;
  fs.assign(N + 1, 0);
  for (int i = 0; i < N; ++i) {
    fs[i + 1] = fs[i] + W[i] * (P[i + 1] - P[i]);
  }
  vector<Double> ans(Q, INF);
  for (int q = 0; q < Q; ++q) {
    const Int f0 = calc(X0[q]);
    const Int f1 = calc(X1[q]);
    auto check = [&](int i, Int x, Int f) -> void {
      const Int w = (0 <= i && i < N) ? W[i] : 0;
      Int dx = 0, dy = 0;
      dx += abs(x - X0[q]);
      dx += abs(x - X1[q]);
      dy += (X0[q] <= x) ? (f - f0) : (f0 - f);
      dy += (X1[q] <= x) ? (f - f1) : (f1 - f);
      /*
        yoko: (dx/V) sec
        tate: katteni (dy/V)
      */
      Double cost = 0.0;
      cost += (Double)dx / V;
      if (V * (Y1[q] - Y0[q]) >= dy) {
        cost += ((Y1[q] - Y0[q]) - (Double)dy / V) / (w + V);
      } else {
        cost += ((Y1[q] - Y0[q]) - (Double)dy / V) / (w - V);
      }
// cerr<<"check i="<<i<<" x="<<x<<" f="<<f<<": w="<<w<<", dx="<<dx<<", dy="<<dy<<", cost="<<cost<<endl;
      chmin(ans[q], cost);
    };
    const int i0 = indexOf(X0[q]);
    check(i0, X0[q], f0);
    for (int i = 0; i <= N; ++i) {
      check(i - 1, P[i], fs[i]);
      check(i    , P[i], fs[i]);
    }
  }
  return ans;
}
}  // brute

int main() {
  for (; ~scanf("%d%d%lld", &N, &Q, &V); ) {
    P.resize(N + 1);
    for (int i = 0; i <= N; ++i) {
      scanf("%lld", &P[i]);
    }
    W.resize(N);
    for (int i = 0; i < N; ++i) {
      scanf("%lld", &W[i]);
    }
    X0.resize(Q);
    Y0.resize(Q);
    X1.resize(Q);
    Y1.resize(Q);
    for (int q = 0; q < Q; ++q) {
      scanf("%lld%lld%lld%lld", &X0[q], &Y0[q], &X1[q], &Y1[q]);
    }
    
    const auto ans = brute::run();
    for (int q = 0; q < Q; ++q) {
      printf("%.8f\n", ans[q]);
    }
  }
  return 0;
}