#include<bits/stdc++.h>
using namespace std;
const int N=200005;
int x[N],y[N];
int xx[N],yy[N];
int idx,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){
	vector<int> a;
	vis[u]=++idx;
	int fid=0;
	for(int i=h[u];i;i=ne[i]){
		int v=adj[i];
		if(v==fa){
			fid=id[i];
			continue;
		}
		if(vis[v]){
			if(vis[u]>vis[v])a.push_back(id[i]);
			continue;
		}
		if(dfs(v,u))a.push_back(id[i]);
	}
	int l=0;
	for(auto i:a){
		if(l){
			mch[i]=l;
			mch[l]=i;
			l=0;
			ans++;
		}
		else l=i;
	}
	if(l&&fid){
		mch[fid]=l;
		mch[l]=fid;
		ans++;
		return 0;
	}
	return 1;
}
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];
	}
	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);
	
	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);
	}
	
	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;
		}
	}
	
	idx=0;
	for(int i=1;i<=px+py;i++)vis[i]=0;
	for(int i=1;i<=n;i++)mch[i]=0;
	for(int i=1;i<=px+py;i++)h[i]=0;
	for(int i=1;i<=tp;i++)ne[i]=0;
	tp=0;
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	int t;
	cin>>t;
	while(t--)work();
	return 0;
}