#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cstring>
#include<queue>
#include<map>
#define ll long long
using namespace std;
int tt,n,a[2005];
struct node{
	ll x,y;
}p[500005];
ll f[1005][1005][3];
node ms(node a,node b){
	node c;
	c.x=a.x-b.x,c.y=a.y-b.y;
	return c;
}
ll cross(node a,node b){
	return abs(a.x*b.y-a.y*b.x);
}
ll area(int p1,int p2,int p3){
	p1=(p1-1)%n+1;
	p2=(p2-1)%n+1;
	p3=(p3-1)%n+1;
	return cross(ms(p[p1],p[p2]),ms(p[p1],p[p3]));
}
int main(){
	scanf("%d",&tt);
	while(tt--){
		scanf("%d",&n);
		for(int i=1;i<=n;++i) scanf("%d",&a[i]),a[i+n]=a[i];
		for(int i=1;i<=n;++i){
			scanf("%lld%lld",&p[i].x,&p[i].y);
		}
		for(int i=1;i<=n;++i)
			for(int j=i+n-1;j>=i;--j) f[i][j][0]=f[i][j][1]=f[i][j][2]=-1;
		for(int i=1;i<=n;++i){
			for(int j=i+n-1;j>=1;--j){
				if(a[i]==a[j]) f[i][j][0]=max(f[i][j][0],0LL);
				for(int k=i+1;k<j;++k){
					if(f[i][j][0]!=-1) f[k][j][1]=max(f[k][j][1],f[i][j][0]+area(i,j,k));
					if(f[i][j][0]!=-1) f[i][k][2]=max(f[i][k][1],f[i][j][0]+area(i,j,k));
					if(f[i][j][1]!=-1&&a[k]==a[i]) f[i][k][0]=max(f[i][k][0],f[i][j][1]+area(i,j,k));
					if(f[i][j][2]!=-1&&a[k]==a[j]) f[k][j][0]=max(f[k][j][1],f[i][j][2]+area(i,j,k));
				}
			}
		}
		ll ans=0;
		for(int i=1;i<=n;++i)
			for(int j=i+n-1;j>=i;--j){
				ans=max(ans,f[i][j][0]);
				ans=max(ans,f[i][j][1]);
				ans=max(ans,f[i][j][2]);
			}
		printf("%lld\n",ans);
	}
}