#include<stdio.h>
#include<algorithm>
#include<queue>
using namespace std;
int xx[512],yy[512];
int x1[256],x2[256],y1[256],y2[256];
int numx,numy;
int bsx(int x)
{
    int l=0,r=numx-1;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(xx[mid]==x)
        {
            return mid;
        }
        if(x<xx[mid])
        {
            r=mid-1;
        }
        else
        {
            l=mid+1;
        }
    }
    return -1;
}
int bsy(int y)
{
    int l=0,r=numy-1;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(yy[mid]==y)
        {
            return mid;
        }
        if(y<yy[mid])
        {
            r=mid-1;
        }
        else
        {
            l=mid+1;
        }
    }
    return -1;
}
int head[262144],last[1048576],to[1048576],cnt=0;
double w[1048576];
void add(int u,int v,double ww)
{
    cnt++;
    last[cnt]=head[u];
    head[u]=cnt;
    to[cnt]=v;
    w[cnt]=ww;
}
int a[410][410];
void inc(int x1,int y1,int x2,int y2)
{
    a[x1][y1]++;
    a[x1][y2+1]--;
    a[x2+1][y1]--;
    a[x2+1][y2+1]++;
}
double d[262144];
struct node{
    int u;
    double d;
    bool operator<(node b) const
    {
        return d>b.d;
    }
};
priority_queue<node> q;
void dijkstra(int s)
{
    d[s]=0;
    node ss;
    ss.u=s;
    ss.d=0;
    q.push(ss);
    while(!q.empty())
    {
        node u=q.top();
        q.pop();
        if(u.d>d[u.u])
        {
            continue;
        }
        for(int i=head[u.u];i;i=last[i])
        {
            int v=to[i];
            if(u.d+w[i]<d[v])
            {
                d[v]=u.d+w[i];
                node vv;
                vv.u=v;
                vv.d=d[v];
                q.push(vv);
            }
        }
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&n);
        for(int i=0;i<=n;i++)
        {
            scanf("%d%d%d%d",&x1[i],&y1[i],&x2[i],&y2[i]);
            xx[i<<1]=x1[i];
            yy[i<<1]=y1[i];
            xx[(i<<1)+1]=x2[i];
            yy[(i<<1)+1]=y2[i];
        }
        sort(xx,xx+((n+1)<<1));
        sort(yy,yy+((n+1)<<1));
        numx=0;
        numy=0;
        for(int i=1;i<((n+1)<<1);i++)
        {
            if(xx[i]>xx[numx])
            {
                numx++;
                xx[numx]=xx[i];
            }
            if(yy[i]>yy[numy])
            {
                numy++;
                yy[numy]=yy[i];
            }
        }
        numx++;
        numy++;
        for(int i=0;i<=n;i++)
        {
            x1[i]=bsx(x1[i])+1;
            x2[i]=bsx(x2[i])+1;
            y1[i]=bsy(y1[i])+1;
            y2[i]=bsy(y2[i])+1;
        }
        for(int i=1;i<=numx*numy;i++)
        {
            head[i]=0;
        }
        cnt=0;
        for(int i=0;i<=numx+1;i++)
        {
            for(int j=0;j<=numy+1;j++)
            {
                a[i][j]=0;
            }
        }
        for(int i=0;i<n;i++)
        {
            inc(x1[i]+1,y1[i],x2[i],y2[i]);
        }
        for(int i=1;i<=numx;i++)
        {
            for(int j=1;j<=numy;j++)
            {
                a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];
                if(i>1)
                {
                    add((i-2)*numy+j,(i-1)*numy+j,(double)(xx[i-1]-xx[i-2])/(a[i][j]+1));
                    add((i-1)*numy+j,(i-2)*numy+j,(double)(xx[i-1]-xx[i-2])/(a[i][j]+1));
                }
            }
        }
        for(int i=0;i<=numx+1;i++)
        {
            for(int j=0;j<=numy+1;j++)
            {
                a[i][j]=0;
            }
        }
        for(int i=0;i<n;i++)
        {
            inc(x1[i],y1[i]+1,x2[i],y2[i]);
        }
        for(int i=1;i<=numx;i++)
        {
            for(int j=1;j<=numy;j++)
            {
                a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];
                if(j>1)
                {
                    add((i-1)*numy+j-1,(i-1)*numy+j,(double)(yy[j-1]-yy[j-2])/(a[i][j]+1));
                    add((i-1)*numy+j,(i-1)*numy+j-1,(double)(yy[j-1]-yy[j-2])/(a[i][j]+1));
                }
            }
        }
        for(int i=1;i<=numx*numy;i++)
        {
            d[i]=2000000005;
        }
        dijkstra((x1[n]-1)*numy+y1[n]);
        printf("%.5f\n",d[(x2[n]-1)*numy+y2[n]]);
    }
    return 0;
}