#include<bits/stdc++.h>
using namespace std;
#define FOR(i,s,t) for(int i=(s),_t=(t); i<=_t; ++i)
typedef long long ll;
const int N=1005;
const int M=1e5+50;
int f[N][N][4];
ll C[M][4];
void add(int x1,int x2,int y1,int y2,int p) {
    f[x1][y1][p]++;
    f[x2+1][y1][p]--;
    f[x1][y2+1][p]--;
    f[x2+1][y2+1][p]++;
}
void solve() {
    memset(f,0,sizeof(f));
    int n;
    cin >> n;
    FOR(i,0,n) {
        C[i][0]=1;
        FOR(j,1,min(i,3)) C[i][j]=C[i-1][j-1]+C[i-1][j];
    }
    while(n--) {
        int x1,x2,y1,y2;
        cin >> x1 >> y1 >> x2 >> y2;
        add(x1,x1,y1,y1,0);
        add(x1,x1,y1+1,y2,1);
        add(x1+1,x2,y1,y1,2);
        add(x1+1,x2,y1+1,y2,3);
    }
    ll res=0;
    FOR(i,1,1000) FOR(j,1,1000) {
        FOR(k,0,3) f[i][j][k]+=f[i-1][j][k]+f[i][j-1][k]-f[i-1][j-1][k];
        // if(i==2 && j==2) cerr << f[2][2][0] << endl;
        FOR(v0,0,3) FOR(v1,0,3-v0) FOR(v2,0,3-v0-v1) {
            if(v0+v1 && v0+v2) res+=C[f[i][j][0]][v0]*C[f[i][j][1]][v1]*C[f[i][j][2]][v2]*C[f[i][j][3]][3-v0-v1-v2];
        }
    }
    cout << res << '\n';
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int T;
    cin >> T;
    while(T--) solve();
    return 0;
}