#include<bits/stdc++.h>
#define eps 1e-6
#define N 205
#define x1 drasgjdlslsh
#define x2 drasgjdlslhh
#define Ms(a,b) memset(a,b,sizeof a)
#define db(x) cout<<#x<<"="<<x<<endl;
#define db2(x,y) cout<<#x<<"="<<x<<" "<<#y<<"="<<y<<endl;
#define db3(x,y,z) cout<<#x<<"="<<x<<" "<<#y<<"="<<y<<" "<<#z<<"="<<z<<endl;
using namespace std;
struct node {
    int x1,x2,y1,y2;
} a[N];
int T,n,x[N<<1],y[N<<1];
int sx,sy,ex,ey;
int mp[2][N<<1][N<<1];
struct Node {
    int x,y;
    double dis;
    bool operator <(const Node &b)const {
        return dis>b.dis;
    }
} now,nxt;
priority_queue<Node>q;
double dis[N<<1][N<<1];
//int dx[]= {1,1,1,-1,-1,-1,0,0};
//int dy[]= {1,-1,0,1,-1,0,1,-1};
int dx[]= {1,-1,0,0};
int dy[]= {0,0,1,-1};
int dd[]= {0,0,1,1};
int rx[]= {0,-1,0,0};
int ry[]= {0,0,0,-1};
int lx,ly;
bool check(Node &x) {
    return x.x>=1&&x.y>=1&&x.x<=lx&&x.y<=ly;
}
int main() {
    scanf("%d",&T);
    while(T--) {
        lx=ly=0,Ms(mp,0);
        scanf("%d",&n);
        for(int i=1; i<=n; i++) {
            scanf("%d%d%d%d",&a[i].x1,&a[i].y1,&a[i].x2,&a[i].y2);
//            a[i].x2--,a[i].y2--;
            x[++lx]=a[i].x1,x[++lx]=a[i].x2;
            y[++ly]=a[i].y1,y[++ly]=a[i].y2;
        }
        scanf("%d%d%d%d",&sx,&sy,&ex,&ey);
//        ex--,ey--;
        x[++lx]=sx,x[++lx]=ex;
        y[++ly]=sy,y[++ly]=ey;
        sort(x+1,x+1+lx),lx=unique(x+1,x+1+lx)-x-1;
        sort(y+1,y+1+ly),ly=unique(y+1,y+1+ly)-y-1;
//        for(int i=1; i<=lx; i++)printf("%d ",x[i]);
//        puts("");
//        for(int i=1; i<=ly; i++)printf("%d ",y[i]);
//        puts("");

        sx=lower_bound(x+1,x+1+lx,sx)-x;
        sy=lower_bound(y+1,y+1+ly,sy)-y;
        ex=lower_bound(x+1,x+1+lx,ex)-x;
        ey=lower_bound(y+1,y+1+ly,ey)-y;
        for(int i=1; i<=n; i++) {
            a[i].x1=lower_bound(x+1,x+1+lx,a[i].x1)-x;
            a[i].x2=lower_bound(x+1,x+1+lx,a[i].x2)-x;

            a[i].y1=lower_bound(y+1,y+1+ly,a[i].y1)-y;
            a[i].y2=lower_bound(y+1,y+1+ly,a[i].y2)-y;

            for(int j=a[i].y1; j<=a[i].y2; j++)mp[0][a[i].x1][j]++,mp[0][a[i].x2][j]--;
            for(int j=a[i].x1; j<=a[i].x2; j++)mp[1][j][a[i].y1]++,mp[1][j][a[i].y2]--;
        }
//        for(int i=1; i<=lx; i++) {
//            for(int j=1; j<=ly; j++)
//                printf("%d ",mp[0][i][j]);
//            puts("");
//        }
//        puts("");
        for(int i=1; i<=lx; i++) {
            for(int j=1; j<=ly; j++) {
                mp[0][i][j]+=mp[0][i-1][j];
                mp[1][i][j]+=mp[1][i][j-1];
                dis[i][j]=2e18;
            }
        }
//        for(int i=1; i<=lx; i++) {
//            for(int j=1; j<=ly; j++)
//                printf("%d ",mp[0][i][j]);
//            puts("");
//        }
        
        dis[sx][sy]=0,q.push((Node)<%sx,sy,0%>);
        while(!q.empty()) {
            now=q.top(),q.pop();
            if(fabs(dis[now.x][now.y]-now.dis)>eps)continue;
//            puts("------");
//            db3(now.x,now.y,now.dis)
//            puts("------");
            for(int i=0; i<4; i++) {
                nxt=now;
                nxt.x+=dx[i],nxt.y+=dy[i];
                if(!check(nxt))continue;
//                if(now.x==2&&now.y==2&&nxt.x==3&&nxt.y==2){
//                    db2(x[nxt.x],x[now.x])
//                }
                nxt.dis+=1.0*(abs(x[nxt.x]-x[now.x])+abs(y[nxt.y]-y[now.y]))/(mp[dd[i]][now.x+rx[i]][now.y+ry[i]]+1);
                if(nxt.dis<dis[nxt.x][nxt.y]) {
//                    db3(nxt.x,nxt.y,nxt.dis)
                    dis[nxt.x][nxt.y]=nxt.dis;
                    q.push((Node)<%nxt.x,nxt.y,nxt.dis%>);
                }
            }
        }
        printf("%.5lf\n",dis[ex][ey]);
    }
    return 0;
}