#include <cstdio>

#include <cmath>

#include <cstring>

#include <iostream>

#include <algorithm>

#include <queue>

#include <vector>

#include <map>



using namespace std;

#define pb push_back
#define mp make_pair
#define F first
#define S second

const int N=100005;

typedef long long LL;
pair<double, double> p[N];
int idx[N];
int a[N];
inline double sqr(double x) {
    return x * x;
}
double inline dis(int a, int b) {
    return sqrt(sqr(p[a].F - p[b].F) + sqr(p[a].S - p[b].S));
}
bool cmpx(const int &a,const int &b) {
    return p[a].F < p[b].F;
}
bool cmpy(const int &a,const int &b) {
    return p[a].S < p[b].S;
}
struct Node
{
    int x,y;
    double d;
};
Node find(int l, int r) {
    if (l == r - 1)    
    {
        Node ans;
        ans.x=idx[l];
        ans.y=idx[r];
        ans.d=dis(ans.x,ans.y);
        return ans;
    }
    if (l == r - 2)    
    {
        Node dq1,dq2,dq3;
        dq1.x=idx[l];dq1.y=idx[l+1];
        dq2.x=idx[l];dq2.y=idx[l+2];
        dq3.x=idx[l+1];dq3.y=idx[l+2];
        dq1.d=dis(dq1.x,dq1.y);
        dq2.d=dis(dq2.x,dq2.y);
        dq3.d=dis(dq3.x,dq3.y);
        if (dq1.d<dq2.d)
        {
            if (dq1.d<dq3.d) return dq1;
            else return dq3;
        }
        else
        {
            if (dq2.d<dq3.d) return dq2;
            else return dq3;
        }
    }
    int mid = l + r >> 1;
    Node dq1,dq2,t;
    dq1=find(l,mid);
    dq2=find(mid+1,r);
    if (dq1.d<dq2.d) t=dq1;else t=dq2;
    int cnt = 0;
    for (int i = l; i <= r; ++i) {
        if (abs(p[idx[i]].F - p[idx[mid]].F) <= t.d)    a[++cnt] = idx[i];
    }
    sort(a + 1, a + cnt + 1, cmpy);
    for (int i = 1; i <= cnt; ++i)
        for (int j = i + 1; j <= cnt; ++j) {
            if (p[a[j]].S - p[a[i]].S >= t.d)    break;
            dq1.x=a[i];dq1.y=a[j];
            dq1.d=dis(dq1.x,dq1.y);
            if (dq1.d<t.d) t=dq1;
        }
    return t;
}
double x[N],y[N];
int main() {
    int n;
    int T;
    scanf("%d",&T);
    while (T--) 
    {
        scanf("%d",&n);
        for (int i = 1; i <= n; ++i) {
            scanf("%lf%lf", &x[i], &y[i]);
            p[i] = mp(x[i], y[i]);
            idx[i] = i;
        }
        sort(idx + 1, idx + n + 1, cmpx);
        Node ans = find(1, n);
        LL sum,dq;
        sum=0;

        dq=p[ans.x].F-p[ans.y].F;
        dq=dq*dq;
        sum+=dq;
        dq=p[ans.x].S-p[ans.y].S;
        dq=dq*dq;
        sum+=dq;
        sum=sum*(n-2);
        int num=0;
        for (int i=1;i<=n;++i)
        {
            if (i==ans.x) continue;
            ++num;p[num]=mp(x[i],y[i]);
            idx[num]=num;
        }
        sort(idx + 1, idx + num + 1, cmpx);
        
        Node ans1 = find(1, num);
        dq=p[ans1.x].F-p[ans1.y].F;
        dq=dq*dq;
        sum+=dq;
        dq=p[ans1.x].S-p[ans1.y].S;
        dq=dq*dq;
        sum+=dq;
        num=0;
        for (int i=1;i<=n;++i)
        {
            if (i==ans.y) continue;
            ++num;p[num]=mp(x[i],y[i]);
            idx[num]=num;
        }
        sort(idx + 1, idx + num + 1, cmpx);
        Node ans2 = find(1, num);
        dq=p[ans2.x].F-p[ans2.y].F;
        dq=dq*dq;
        sum+=dq;
        dq=p[ans2.x].S-p[ans2.y].S;
        dq=dq*dq;
        sum+=dq;
        printf("%I64d\n",sum);
    }
    return 0;
}