#include<bits/stdc++.h>
using namespace std;
const int N=100005;
int x[N],y[N];
int xx[N],yy[N];
int vis[N],ans;
int h[N],ne[N<<1],adj[N<<1],id[N<<1],tp;
int mch[N];
void add(int u,int v,int eid){
	ne[++tp]=h[u];
	h[u]=tp;
	adj[tp]=v;
	id[tp]=eid;
}
int dfs(int u,int fa){
	int l=0;
	vis[u]=1;
	for(int i=h[u];i;i=ne[i]){
		int v=adj[i];
		if(v==fa)continue;
		if(vis[v]){
			if(mch[id[i]])continue;
			if(l){
				mch[l]=id[i];
				mch[id[i]]=l;
				l=0;
				ans++;
			}
			else l=id[i];
			continue;
		}
		int k=dfs(v,u);
		if(k){
			mch[k]=id[i];
			mch[id[i]]=k;
			ans++;
		}
		else{
			if(l){
				mch[l]=id[i];
				mch[id[i]]=l;
				l=0;
				ans++;
			}
			else l=id[i];
		}
	}
	return l;
}
void work(){
	int n;
	cin>>n;n<<=1;
	for(int i=1;i<=n;i++){
		cin>>x[i]>>y[i];
		xx[i]=x[i];
		yy[i]=y[i];
		mch[i]=0;
	}
	sort(xx+1,xx+1+n);
	sort(yy+1,yy+1+n);
	int px=unique(xx+1,xx+1+n)-(xx+1);
	int py=unique(yy+1,yy+1+n)-(yy+1);
	tp=0;
	for(int i=1;i<=px+py;i++)h[i]=0;
	
	for(int i=1;i<=n;i++){
		int nx=lower_bound(xx+1,xx+1+px,x[i])-xx;
		int ny=lower_bound(yy+1,yy+1+py,y[i])-yy;
		add(nx,ny+px,i);
		add(ny+px,nx,i);
	}
	
	for(int i=1;i<=px+py;i++)vis[i]=0;
	ans=0;
	for(int i=1;i<=px+py;i++){
		if(vis[i])continue;
		dfs(i,0);
	}
	cout<<ans<<"\n";
	int l=0;
	for(int i=1;i<=n;i++){
		if(mch[i]){
			if(i<mch[i])cout<<i<<" "<<mch[i]<<"\n";
		}
		else{
			if(l){
				cout<<l<<" "<<i<<"\n";
				l=0;
			}
			else l=i;
		}
	}
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	int t;
	cin>>t;
	while(t--)work();
	return 0;
}