#include<bits/stdc++.h>
using namespace std;
//#define int long long
#define debug(x) cout<<x<<"--------------------------------------\n" ;
const int N = 1e5 + 10;
typedef long long ll;
const int INF = 1e9 ;
const int MOD = 1e9 + 7 ;
typedef pair<int,int> PII ;
int gcd(int a,int b){
    return b ? gcd(b,a%b) : a ;
}
ll qpow(int a,int k){
    ll res = 1 ;
    while (k){
        if(k&1) res = 1ll * res * a % MOD ;
        a = 1ll * a * a % MOD ;
        k >>= 1 ;
    }
    return res ;
}
#define double long long
const double eps=0;
const double PI=acos(-1);
inline int sgn(double x){
    return fabs(x)<=eps?0:(x<0?-1:1);
}
double ans2 = 0;
struct point3{
    double x,y,z;
    point3(double a=0,double b=0,double c=0):x(a),y(b),z(c){}
    point3 operator +(const point3 &A)const{
        return point3(x+A.x,y+A.y,z+A.z);
    }
    point3 operator -(const point3 &A)const{
        return point3(x-A.x,y-A.y,z-A.z);
    }
    point3 operator *(const double v)const{
        return point3(x*v,y*v,z*v);
    }
    point3 operator /(const double v)const{
        return point3(x/v,y/v,z/v);
    }
    bool operator ==(const point3 &A)const{
        return sgn(x-A.x)==0&&sgn(y-A.y)==0&&sgn(z-A.z)==0;
    }
    double length(){
        return sqrt(x*x+y*y+z*z);
    }
    point3 unit(){
        return *this/length();
    }
};
point3 det(point3 a,point3 b){
    return point3(a.y*b.z-a.z*b.y,a.z*b.x-a.x*b.z,a.x*b.y-a.y*b.x);
}
ll dot(point3 a,point3 b){
    return a.x*b.x+a.y*b.y+a.z*b.z;
}
//除以6就是a,b,c三个向量所构成的四面体的体积
double mix(point3 a,point3 b,point3 c){
    return dot(a,det(b,c));
}
struct plane{
    point3 v[3];
    plane() {}
    plane(point3 uu,point3 vv,point3 ww) {v[0]=uu,v[1]=vv,v[2]=ww;}
    point3 normal() {return det((v[1]-v[0]),(v[2]-v[0]));}
};
ll ans;
ll Volume(point3 a,point3 b,point3 c,point3 d){
    return dot(det(b-a,c-a) ,(d-a));
}
namespace CH3{
    vector<vector<point3>> ret;
    set<pair<int,int> > eg;
    set<int> pt;
    int n;
    vector<point3> p,q;

    void wrap(int a,int b){
        if (eg.find({a,b})==eg.end()){
            int c=-1;
            for (int i=0;i<n;i++)if (i!=a && i!=b){
                    if (c==-1 || Volume(q[c],q[a],q[b],q[i])>0)
                        c=i;
                }
            if (c!=-1){
                ans2 += det(p[b]-p[a],p[b]-p[c]).length()/2;
//                debug(1)
//                cout << ans2 << '\n';
//                cout << p[a].x << ' ' << p[a].y << ' ' << p[a].z << '\n';
//                cout << p[b].x << ' ' << p[b].y << ' ' << p[b].z << '\n';
//                cout << p[c].x << ' ' << p[c].y << ' ' << p[c].z << '\n';
                ans+=Volume(point3(0,0,0),p[a],p[b],p[c]);
                ret.push_back({p[a],p[b],p[c]});
                pt.insert(a); pt.insert(b); pt.insert(c);
                eg.insert({a,b});eg.insert({b,c});eg.insert({c,a});
                wrap(c,b);wrap(a,c);
            }
        }
    };

    vector<point3> convexHull3d(vector<point3> _p){
        p = q = _p; n = p.size();
        ret.clear(); eg.clear(); pt.clear();
        //ans=0;
        for (int i=1;i<n;i++)if (q[i].x<q[0].x)swap(p[0],p[i]),swap(q[0],q[i]);
        for (int i=2;i<n;i++)if (
                    (q[i].x-q[0].x)*(q[1].y-q[0].y)>
                    (q[i].y-q[0].y)*(q[1].x-q[0].x)) swap(q[1],q[i]),swap(p[1],p[i]);
        wrap(0,1);
        vector<int> pc(pt.begin(),pt.end());
        reverse(pc.begin(),pc.end());
        for (auto x:pc) {
            q.erase(q.begin()+x);
        }
//        cout << ans2 << '\n';
        return q;
    }
}
void solve(){
    int n;cin >> n;
    ans = 0;
    vector<point3> p;
    for(int i=0;i<n;i++){
        double x , y , z; cin >> x >> y >> z;
        p.push_back({x,y,z});
    }
    while (p.size() >= 4){
        p = CH3::convexHull3d(p);
    }
//    cout << fixed << setprecision(3);
//    cout << ans2 << '\n';
    cout << ans << '\n';
}
signed main(){
    ios::sync_with_stdio(false) ;
    cin.tie(0) ;
    int tt = 1 ;
    cin >> tt ;
    while (tt--) solve() ;
    return 0 ;
}