#include <bits/stdc++.h>
typedef long double ld;
const ld inf=1e18;
using namespace std;
int n,m,x[610],y[610];
ld f[610][610][2];
bool e[610][610];
ld dis(int i,int j){
	return sqrtl(1ll*(x[i]-x[j])*(x[i]-x[j])+1ll*(y[i]-y[j])*(y[i]-y[j]));
}
void add(int u,int v){
	if(u>v) swap(u,v);
	e[u][v]=e[v][u]=1;
	f[u][v][0]=f[u][v][1]=dis(u,v); 
}
void sol(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d%d",&x[i],&y[i]);
	for(int i=n+1;i<=2*n;i++) x[i]=x[i-n],y[i]=y[i-n];
	for(int i=1;i<=2*n;i++) for(int j=1;j<=2*n;j++) e[i][j]=0,f[i][j][0]=f[i][j][1]=-(i!=j)*inf;
	scanf("%d",&m);
	for(int i=1,x,y;i<=m;i++) scanf("%d%d",&x,&y),add(x,y),add(x+n,y),add(x,y+n),add(x+n,y+n);
	ld ans=0;
	for(int len=1;len<=n;len++) for(int l=1;l+len-1<=n;l++){
		int r=l+len-1;
		for(int i=l+1;i<r;i++){
			if(e[l][r]){
				f[l][r][0]=max(f[l][r][0],f[i][r][1]+dis(l,r));
				f[l][r][1]=max(f[l][r][1],f[l][i][0]+dis(l,r));
			}
			if(e[i][r]) f[l][r][1]=max(f[l][r][1],f[l][i][1]+dis(i,r));
			if(e[l][i]) f[l][r][0]=max(f[l][r][0],f[i][r][0]+dis(l,i));
		}
		ans=max({ans,f[l][r][0],f[l][r][1]});
	}
	printf("%.12Lf\n",ans);
}
int main(){
	int T;
	scanf("%d",&T);
	while(T--) sol();
}