#include<bits/stdc++.h>
using namespace std;
#ifdef DEBUG
#include"debug.h"
#else
#define debug(...) void()
#endif
#define all(x) (x).begin(),(x).end()
template<class T>
auto ary(T *a,int l,int r){
	return vector<T>{a+l,a+1+r};
}
using ll=long long;
using ull=unsigned long long;
using ld=long double;
using vec=complex<int>;
int dot(const vec &a,const vec &b){
	return real(a)*real(b)+imag(a)*imag(b);
}
int cross(const vec &a,const vec &b){
	return dot(a,b*vec(0,-1));
}
ll dis2(const vec &a){
	return dot(a,a);
}
ld dis(const vec &a){
	return sqrtl(dis2(a));
}
const int N=5e3+10;
int T,n,m;
vec o,a[N],b[N];
void read(vec &a){
	int x,y;
	scanf("%d%d",&x,&y);
	a=vec(x,y);
}
template<class T>
int sign(T x){
	return x>0?1:(x<0?-1:0);
}
ld calc(vec p,vec a,vec b){
	if(dot(p-a,b-a)<0)return dis(a-p);
	if(dot(p-b,a-b)<0)return INFINITY;
	return abs(cross(p-a,p-b))/dis(a-b);
}
bool atline(vec p,vec a,vec b){
	return !cross(p-a,p-b)&&dot(p-a,p-b)<=0;
}
int inTriangle(vec p,vec a,vec b,vec c){
	if(atline(p,a,b)||atline(p,b,c)||atline(p,c,a))return 0;
	return
		sign(cross(p-a,b-a))*sign(cross(p-a,c-a))<0&&
		sign(cross(p-b,a-b))*sign(cross(p-b,c-b))<0&&
		sign(cross(p-c,a-c))*sign(cross(p-c,b-c))<0?1:-1;
}
ld ans[N];
void work(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)read(a[i]);
	read(o);
	scanf("%d",&m);
	for(int i=1;i<=m;i++)read(b[i]);
	int st=[&](){
		for(int i=1;i<=n;i++){
			if(a[i]==o)return i;
		}
		int p=1;
		ld val=calc(o,a[1],a[2]);
		for(int i=2;i<=n;i++){
			ld dis=calc(o,a[i],a[i%n+1]);
			if(dis<val)val=dis,p=i;
		}
		return p;
	}();
	static vec stk[N];
	static ld val[N];
	int top=1;
	stk[top]=o,val[top]=0;
	vec las=a[st];
	for(int i=1,u=st%n+1;i<=n;){
		debug(u,las,ary(stk,1,top));
		if(a[u]==stk[top]){
			debug("op1");
			if(top==1)break;
			top--;
			for(int j=1;j<=m;j++){
				if(atline(b[j],stk[top],a[u])){
					ans[j]=val[top]+dis(b[j]-stk[top]);
				}
			}
			continue;
		}
		vec cur=a[u];
		int s=(u+n-2)%n+1;
		for(int j=1;j<=n;j++){
			if(cross(las-stk[top],a[j]-stk[top])<=0)continue;
			if(cross(a[j]-stk[top],a[u]-stk[top])<=0)continue;
			if(cross(a[u]-a[s],a[j]-a[s])<0)continue;
			if(cross(a[j]-stk[top],cur-stk[top])>0)cur=a[j];
		}
		bool flag=0;
		if(top>1){
			vec p=stk[top]*2-stk[top-1];
			if(cross(las-stk[top],p-stk[top])>0&&cross(p-stk[top],a[u]-stk[top])>=0){
				if(cross(p-stk[top],cur-stk[top])>=0)cur=p,flag=1;
			}
		}
		debug(stk[top],las,cur);
		for(int j=1;j<=m;j++){
			if(cross(las-stk[top],b[j]-stk[top])<0)continue;
			if(cross(b[j]-stk[top],a[u]-stk[top])<0)continue;
			if(cross(a[u]-a[s],b[j]-a[s])<0)continue;
			debug("ans",ary(stk,1,top),val[top]);
			ans[j]=val[top]+dis(b[j]-stk[top]);
		}
		if(flag){
			las=cur;
			debug("op2");
			top--;
		}else if(cur!=a[u]){
			las=cur*2-stk[top];
			debug("op3");
			stk[++top]=cur;
			val[top]=val[top-1]+dis(stk[top]-stk[top-1]);
		}else{
			debug("op4");
			i++;
			int v=u%n+1;
			if(cross(a[u]-stk[top],a[v]-stk[top])<0){
				las=a[v]*2-a[u];
				stk[++top]=a[u];
				val[top]=val[top-1]+dis(stk[top]-stk[top-1]);
			}else las=a[u]*2-stk[top];
			u=v;
		}
	}
	for(int i=1;i<=m;i++)printf("%.15Lf\n",ans[i]);
	fill(ans+1,ans+1+m,ld(0));
}
int main(){
	for(scanf("%d",&T);T--;)work();
	return 0;
}
#ifdef DEBUG
#include"debug.hpp"
#endif