#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1001
using namespace std;
int n,opl,opr,w;
double ans,x1,x2,y3,y4,has[N<<1],sum1[N<<4],sum2[N<<4];
int f[N<<4];
struct node
{
    double l,r,h;
    int f;
    bool operator < (node p)const
    {
        return h<p.h;
    } 
}e[N<<1];
void up(int k,int l,int r)
{
    if(f[k]>=2) sum2[k]=has[r+1]-has[l];
    else if(f[k]==1) 
    {
        sum1[k]=has[r+1]-has[l];
        sum2[k]=sum1[k<<1]+sum1[k<<1|1];
    }
    else
    {
        sum1[k]=sum1[k<<1]+sum1[k<<1|1];
        sum2[k]=sum2[k<<1]+sum2[k<<1|1];
    }
}
void change(int k,int l,int r)
{
    if(opl<=l && r<=opr)
    {
        f[k]+=w;
        up(k,l,r);
        return;
    }
    int mid=l+r>>1;
    if(opl<=mid) change(k<<1,l,mid);
    if(opr>mid) change(k<<1|1,mid+1,r);
    up(k,l,r);
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y3,&x2,&y4);
            has[i*2-1]=x1; has[i*2]=x2;
            e[i*2-1].l=x1; e[i*2-1].r=x2; e[i*2-1].f=-1; e[i*2-1].h=y4;
            e[i*2].l=x1;   e[i*2].r=x2;   e[i*2].f=1;    e[i*2].h=y3;
        }
        sort(e+1,e+n*2+1);
        sort(has+1,has+n*2+1);
        int cnt=unique(has+1,has+n*2+1)-(has+1);
        for(int i=1;i<=n*2;i++)
        {
            opl=lower_bound(has,has+cnt+1,e[i].l)-has;
            opr=lower_bound(has,has+cnt+1,e[i].r)-has-1;
            w=e[i].f;
            change(1,1,cnt);
            ans+=sum2[1]*(e[i+1].h-e[i].h);
        }
        printf("%.2lf\n",ans);
        memset(f,0,sizeof(f));
        memset(sum1,0,sizeof(sum1));
        memset(sum2,0,sizeof(sum2));
        ans=0;
    }
}