#include<bits/stdc++.h>
using namespace std;
const int N=4e5;
int n;
int tot;
struct E{
    int v,id;
    E(int v=0,int id=0):v(v),id(id){}
};
vector<E> e[N+5];
unordered_map<int,int> mpx,mpy;
void adde(int u,int v,int id){
    e[u].push_back(E(v,id));
    e[v].push_back(E(u,id));
}
int dfn[N+5],dfns;
bool use[N+5];
vector<pair<int,int>> ans;
void dfs(int u,int tf){
    int fid=-1;
    dfn[u]=++dfns;
    for(auto it:e[u]){
        if(it.v==tf) fid=it.id;
        if(!dfn[it.v])
            dfs(it.v,u);
    }
    int lst=-1;
    for(auto [v,id]:e[u]){
        if(v!=tf&&!use[id]){
            if(lst==-1) lst=id;
            else{
                ans.push_back({lst,id});
                use[lst]=use[id]=1;
                lst=-1;
            }
        }
    }
    if(lst!=-1&&fid!=-1){
        ans.push_back({lst,fid});
        use[lst]=use[fid]=1;
    }
}
void solve()
{
    tot=0;
    scanf("%d",&n);
    for(int i=1,x,y;i<=n*2;i++){
        scanf("%d %d",&x,&y);
        if(!mpx[x])
            mpx[x]=++tot;
        if(!mpy[y])
            mpy[y]=++tot;
        adde(mpx[x],mpy[y],i);
    }
    for(int i=1;i<=tot;i++)
        if(!dfn[i]){
            dfs(i,0);
        }
    printf("%d\n",ans.size());
    for(auto [a,b]:ans){
        printf("%d %d\n",a,b);
    }
    for(int i=1,j=-1;i<=n*2;i++)
        if(!use[i]){
            if(j!=-1){
                printf("%d %d\n",i,j);
                j=-1;
            }
            else j=i;
        }
    fill(use,use+n*2+1,0);
    ans.clear();
    fill(dfn,dfn+tot+1,0); dfns=0;
    mpx.clear(); mpy.clear();
    for(int i=1;i<=tot;i++) e[i].clear();
    tot=0;
}
int main()
{
    int t;
    scanf("%d",&t);
    while (t--)
        solve();
    return 0;
}