#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

struct vec{
    ll x,y;
    vec operator-(vec a){
        return {x-a.x,y-a.y};
    }
};

ll cross(vec a,vec b){
    return abs(a.x*b.y-a.y*b.x);
}

const int N=510;
int n;
int ptr[N],idx[N];
ll f[N];
ll dp[N][N];
vec a[N];
vector<int> p[N];

int nxt(int x){return (x+1)%n;}
int pre(int x){return (x+n-1)%n;}

bool in(int l,int r,int p){
    if(r>=l)return r>=p&&p>=l;
    else return p>=r||p<=l;
}

bool in(int l,int r,int pl,int pr){
    if(pr>=r&&r>=l&&l>=pl)return 1;
    if(r>=l&&l>=pl&&pl>=pr)return 1;
    if(l>=pl&&pl>=pr&&pr>=r)return 1;
    if(pl>=pr&&pr>=r&&r>=l)return 1;
    return 0;
}

void solve(){
    cin>>n;

    map<int,int> mp;
    int cnt=0;
    for(int i=0;i<n;i++){
        cin>>f[i];
        if(mp.find(f[i])==mp.end())mp[f[i]]=cnt++;
        f[i]=mp[f[i]];

        p[i].clear();
    }

    for(int i=0;i<n;i++){
        idx[i]=p[f[i]].size();
        p[f[i]].push_back(i);
    }

    for(int i=0;i<n;i++)cin>>a[i].x>>a[i].y;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            if(f[i]==f[j])dp[i][j]=0;
            else dp[i][j]=-1;
        }
    }
    for(int l=1;l<n;l++){
        for(int i=0;i<n;i++)ptr[i]=(idx[i]+1)%p[f[i]].size();
        for(int i=0;i<n;i++){
            int j=(i+l-1)%n;
            if(dp[i][j]==-1)continue;

            for(int nxti=nxt(j);nxti!=i;nxti=nxt(nxti)){
                int sz=p[f[nxti]].size();
                if(in(i,j,nxti,p[f[nxti]][(idx[nxti]+1)%sz])){
                    ptr[nxti]=(idx[nxti]+1)%sz;
                }
                while (ptr[nxti]!=idx[nxti]&&!in(i,j,nxti,p[f[nxti]][ptr[nxti]])){
                    ptr[nxti]=(ptr[nxti]+1)%sz;
                }
                int nxtj=p[f[nxti]][ptr[nxti]];
                if(nxtj!=nxti){
                    dp[nxti][nxtj]=max(dp[nxti][nxtj],
                        dp[i][j]+cross(a[j]-a[nxtj],a[i]-a[nxtj])+cross(a[i]-a[nxtj],a[nxti]-a[nxtj]));
                }
            }
            
        }
    }
    ll ans=0;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            ans=max(ans,dp[i][j]);
            // cout<<dp[i][j]<<" ";
        }
        // cout<<"\n";
    }
    cout<<ans<<"\n";
}

int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);cout.tie(0);
    int n=1;
    cin>>n;
    while(n--)solve();
    return 0;
}