#include<cstdio>
#include<cmath>
using namespace std;
typedef double db;
const int N=305;
int T,n,m,cnt;
struct node{
	int x,y,next;
}edge[N*N];
int h[N];
void add(int i,int x,int y){
	edge[i].x=x;edge[i].y=y;
	edge[i].next=h[x];h[x]=i;
}
db px[N],py[N],dp[N][N][2];
bool vis[N][N][2];
db max1(db a,db b){
	return a>b?a:b;
}
db w(int i,int j){
	return sqrt((px[i]-px[j])*(px[i]-px[j])+(py[i]-py[j])*(py[i]-py[j]));
}
db f(int i,int j,int v){
	i=i%n;j=j%n;
	if(i==0) i=n;
	if(j==0) j=n;
//	printf("%d %d\n",i,j);
	if(i==j) return 0;
	if(vis[i][j][v]) return dp[i][j][v];
	vis[i][j][v]=true;
	if(i<j){
		for(int k=h[i],y;k;k=edge[k].next){
			y=edge[k].y;
			if(v==0){
				if(i<y and y<=j){
					dp[i][j][v]=max1(dp[i][j][v],w(i,y)+max1(f(y,j,v),f(y,i+1,v^1)));
				}
			}else{
				if(y<i or j<=y){
					dp[i][j][v]=max1(dp[i][j][v],w(i,y)+max1(f(y,j,v),f(y,i-1,v^1)));
				}
			}
		}
	}else{
		for(int k=h[i],y;k;k=edge[k].next){
			y=edge[k].y;
			if(v==1){
				if(j<=y and y<i){
					dp[i][j][v]=max1(dp[i][j][v],w(i,y)+max1(f(y,j,v),f(y,i-1,v^1)));
				}
			}else{
				if(y<=j or i<y){
					dp[i][j][v]=max1(dp[i][j][v],w(i,y)+max1(f(y,j,v),f(y,i+1,v^1)));
				}
			}
		}
	}
	return dp[i][j][v];
} 
int main(){
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n);
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				dp[i][j][0]=dp[i][j][1]=0;
				vis[i][j][0]=vis[i][j][1]=false;
			}
			h[i]=0;
		}
		for(int i=1;i<=n;i++){
			scanf("%lf%lf",&px[i],&py[i]);
		}
		scanf("%d",&m);
		cnt=0;
		for(int i=1,x,y;i<=m;i++){
			scanf("%d%d",&x,&y);
			add(++cnt,x,y);
			add(++cnt,y,x);
		}
//		for(int i=1;i<=cnt;i++) printf("%d--%d\n",edge[i].x,edge[i].y);
//		puts(""); 
		db ans=0;
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				ans=max1(ans,max1(f(i,j,0),f(i,j,1)));
//				printf("%d %d: %lf %lf\n",i,j,dp[i][j][0],dp[i][j][1]);
			}
		}
		printf("%.10lf\n",ans);
	}
	return 0;
}