#include<bits/stdc++.h>
#define double long double
using namespace std;
inline int read()
{
	int n=0,f=1,ch=getchar();
	while(ch<'0'||ch>'9')
	{
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		n=n*10+ch-'0';
		ch=getchar();
	}
	return n*f;
}
int p[200005],v[200005];
int np[200005],nv[200005],ncnt;
int main()
{
	int n,q,V,x1,y1,x2,y2;
	n=read();
	q=read();
	V=read();
	for(int i=0;i<=n;i++)p[i]=read();
	for(int i=1;i<=n;i++)v[i]=read();
	if(n==0)
	{
		for(int i=1;i<=q;i++)
		{
			x1=read();
			y1=read();
			x2=read();
			y2=read();
			printf("%.9lf\n",(fabs(x2-x1)+fabs(y2-y1))/V);
		}
		return 0;
	}
	double wy=0,nans=0;
	for(int i=1;i<=n;i++)nans+=1.0*(p[i]-p[i-1])/V,wy+=1.0*(p[i]-p[i-1])/V*v[i];
	int maxv=0,minv=0;
	for(int i=1;i<=n;i++)maxv=max(maxv,v[i]);
	for(int i=1;i<=n;i++)minv=min(minv,v[i]);
	for(int i=1;i<=q;i++)
	{
		x1=read();
		y1=read();
		x2=read();
		y2=read();
		if(x1==p[0]&&x2==p[n])
		{
			double ans=nans;
			if(y1+wy<y2)
			{
				ans+=(y2-y1-wy)/(maxv+V);
			}
			else
			{
				ans+=(y1+wy-y2)/(V-minv);
			}
			printf("%.9Lf\n",ans);
			continue;
		}
		if(x1>x2)swap(x1,x2);
		double ans=1000000000;
		int pos1=0,pos2=0;
		ncnt=0;
		if(x1<p[0])
		{
			np[++ncnt]=x1;
			nv[ncnt]=0;
			pos1=ncnt;
		}
		if(x2<p[0])
		{
			np[++ncnt]=x2;
			nv[ncnt]=0;
			pos2=ncnt;
		}
		for(int i=0;i<=n-1;i++)
		{
			np[++ncnt]=p[i];
			nv[ncnt]=v[i];
			if(x1>=p[i]&&x1<p[i+1])
			{
				np[++ncnt]=x1;
				nv[ncnt]=v[i+1];
				pos1=ncnt;
			}
			if(x2>=p[i]&&x2<p[i+1])
			{
				np[++ncnt]=x2;
				nv[ncnt]=v[i+1];
				pos2=ncnt;
			}
		}
		np[++ncnt]=p[n];
		nv[ncnt]=v[n];
		if(x1>=p[n])
		{
			np[++ncnt]=x1;
			nv[ncnt]=0;
			pos1=ncnt;
		}
		if(x2>=p[n])
		{
			np[++ncnt]=x2;
			nv[ncnt]=0;
			pos2=ncnt;
		}
		np[++ncnt]=np[ncnt-1];
		nv[ncnt]=0; 
		//for(int i=1;i<=ncnt;i++)printf("%d %d\n",np[i],nv[i]);
		//printf("%d %d\n",pos1,pos2);
		int maxv=-1000000000,minv=1000000000;
		for(int i=pos1+1;i<=pos2;i++)maxv=max(maxv,nv[i]);
		for(int i=pos1+1;i<=pos2;i++)minv=min(minv,nv[i]);
		double wy=0,nans=0;
		for(int i=pos1+1;i<=pos2;i++)nans+=1.0*(np[i]-np[i-1])/V,wy+=1.0*(np[i]-np[i-1])/V*nv[i];
		if(y1+wy<y2)
		{
			ans=min(ans,nans+(y2-y1-wy)/(maxv+V));
		}
		else
		{
			ans=min(ans,nans+(y1+wy-y2)/(V-minv));
		}
		printf("%.9Lf\n",ans);
		continue;
		int nmaxv=maxv,nminv=minv;
		double nwy=wy,nnans=nans;
		for(int i=pos1;i>=1;i--)
		{
			nmaxv=max(nmaxv,nv[i]);
			nminv=min(nminv,nv[i]);
			if(i!=pos1)
			{
				nwy=nwy+2.0*(np[i+1]-np[i])/V*nv[i+1];
				nnans=nnans+2.0*(np[i+1]-np[i])/V;
			}
			if(y1+nwy<y2)
			{
				ans=min(ans,nnans+(y2-y1-nwy)/(nmaxv+V));
			}
			else
			{
				ans=min(ans,nnans+(y1+nwy-y2)/(V-nminv));
			}
		}
		nmaxv=maxv;
		nminv=minv;
		nwy=wy;
		nnans=nans;
		for(int i=pos2;i<=ncnt-1;i++)
		{
			nmaxv=max(nmaxv,nv[i+1]);
			nminv=min(nminv,nv[i+1]);
			if(i!=pos2)
			{
				nwy=nwy+2.0*(np[i]-np[i-1])/V*nv[i];
				nnans=nnans+2.0*(np[i]-np[i-1])/V;
			}
			if(y1+nwy<y2)
			{
				ans=min(ans,nnans+(y2-y1-nwy)/(nmaxv+V));
			}
			else
			{
				ans=min(ans,nnans+(y1+nwy-y2)/(V-nminv));
			}
		}
		printf("%.9Lf\n",ans);
	}
	
}