#include<bits/stdc++.h>
using namespace std;
const int MAXN=4e5+5;
struct info { int len,w,wc,wl,wr,wq,pr; };
struct tag { int w,mw,pr; };
info operator +(const info &x,const info &y) {
	info z{0,0,0,0,0,0,0};
	z.len=x.len+y.len,z.w=min(x.w,y.w);
	z.wc=(x.w==z.w?x.wc:0)+(y.w==z.w?y.wc:0);
	z.wl=max((x.w==z.w?x.wl:0),(x.w==y.w&&x.wl==x.len?x.len+y.wl:0));
	z.wr=max((y.w==z.w?y.wr:0),(x.w==y.w&&y.wr==y.len?y.len+x.wr:0));
	z.wq=max({(x.w==z.w?x.wq:0),(y.w==z.w?y.wq:0),(x.w==y.w?x.wr+y.wl:0)});
	z.pr=(x.w==y.w?(x.pr==y.pr?x.pr:-1):(x.w==z.w?x.pr:y.pr));
	return z;
}
info operator *(const info &x,const tag &y) {
	info z=x;
	if(x.w+y.mw==0&&~y.pr) z.pr=y.pr;
	z.w+=y.w;
	return z;
}
tag operator *(const tag &x,const tag &y) {
	tag z=y;
	z.mw+=x.w;
	if(~x.pr&&(x.mw<z.mw||(x.mw==z.mw&&z.pr==-1))) z.pr=x.pr;
	z.mw=min(z.mw,x.mw),z.w+=x.w;
	return z;
}
map <array<int,2>,array<int,2>> ANS;
void gans(int xl,int xr,int cnt,int mxq) {
	auto &z=ANS[{xl,xr}]; z[0]+=cnt,z[1]=max(z[1],mxq);
}
struct colu {
	info x; int sf; bool tg;
	void clr() { if(tg) gans(x.pr,sf,x.wc,x.wq); }
	void add(const info &y) {
		if(!y.w&&y.pr) {
			if(tg) x.pr==y.pr?(x=x+y):(clr(),x=y);
			else tg=true,x=y;
		} else clr(),tg=false;
	}
};
struct scl { int x,yl,yr; bool c; } a[MAXN];
int n,m,vals[MAXN];
struct SegmentTree {
	info tr[MAXN<<2]; tag tg[MAXN<<2];
	void adt(int p,tag k) { tg[p]=tg[p]*k,tr[p]=tr[p]*k; }
	void psd(int p) { adt(p<<1,tg[p]),adt(p<<1|1,tg[p]),tg[p]={0,0,-1}; }
	void psu(int p) { tr[p]=tr[p<<1]+tr[p<<1|1]; }
	void init(int l=1,int r=m,int p=1) {
		tg[p]={0,0,-1};
		if(l==r) {
			int z=vals[l+1]-vals[l];
			return tr[p]={z,0,z,z,z,z,0},void();
		}
		int mid=(l+r)>>1;
		init(l,mid,p<<1),init(mid+1,r,p<<1|1),psu(p);
	}
	void upd(int ul,int ur,tag k,int l=1,int r=m,int p=1) {
		if(ul<=l&&r<=ur) return adt(p,k);
		int mid=(l+r)>>1; psd(p);
		if(ul<=mid) upd(ul,ur,k,l,mid,p<<1);
		if(mid<ur) upd(ul,ur,k,mid+1,r,p<<1|1);
		psu(p);
	}
	void qry(colu &c,int l=1,int r=m,int p=1) {
		if(~tr[p].pr||tr[p].w) return c.add(tr[p]);
		int mid=(l+r)>>1; psd(p);
		qry(c,l,mid,p<<1),qry(c,mid+1,r,p<<1|1);
	}
}	T;
void solve() {
	cin>>n,ANS.clear();
	for(int i=1,xl,xr,yl,yr;i<=n;++i) {
		cin>>xl>>yl>>xr>>yr,++xr,++yr;
		vals[2*i-1]=yl,vals[2*i]=yr;
		a[2*i-1]={xl,yl,yr,0},a[2*i]={xr,yl,yr,1};
	}
	sort(vals+1,vals+2*n+1),m=unique(vals+1,vals+2*n+1)-vals-1;
	for(int i=1;i<=2*n;++i) {
		a[i].yl=lower_bound(vals+1,vals+m+1,a[i].yl)-vals;
		a[i].yr=lower_bound(vals+1,vals+m+1,a[i].yr)-vals-1;
	}
	sort(a+1,a+2*n+1,[&](auto i,auto j){ return i.x==j.x?i.c<j.c:i.x<j.x; });
	--m,T.init();
	vector <int> cur;
	for(int i=1;i<=2*n;++i) {
		if(a[i].c) T.upd(a[i].yl,a[i].yr,{-1,-1,-1}),cur.push_back(i);
		else T.upd(a[i].yl,a[i].yr,{1,0,a[i].x});
		if(i==2*n||a[i+1].x>a[i].x) {
			colu c{{},a[i].x-1,0};
			T.qry(c),c.clr();
			for(int j:cur) T.upd(a[j].yl,a[j].yr,{0,0,0});
			cur.clear();
		}
	}
	cout<<ANS.size()<<"\n";
	for(auto o:ANS) cout<<o.first[0]<<" "<<o.first[1]<<" "<<o.second[0]<<" "<<o.second[1]<<"\n";
}
signed main() {
	ios::sync_with_stdio(false);
	int Q; cin>>Q;
	while(Q--) solve();
	return 0;
}