#include<bits/stdc++.h>
#define fi first
#define se second
using namespace std;

void solve() {
    int n,lena=0,lenb=0; cin>>n;
    vector<int>a((n<<2)+10,0),b((n<<2)+10,0);
    vector<int>la((n<<2)+10,0),lb((n<<2)+10,0);
    vector<bool>vis((n<<2)+10,false),siv((n<<2)+10,false);
    for(int i=1;i<=(n<<1);++i)
        cin>>a[i]>>b[i],la[++lena]=a[i],lb[++lenb]=b[i];
    sort(la.begin()+1,la.begin()+1+lena);
    sort(lb.begin()+1,lb.begin()+1+lenb);
    lena=unique(la.begin()+1,la.begin()+1+lena)-la.begin()-1;
    lenb=unique(lb.begin()+1,lb.begin()+1+lenb)-lb.begin()-1;
    vector<set<pair<int,int>>>g((n<<2)+10);
    for(int i=1;i<=(n<<1);++i) {
        a[i]=lower_bound(la.begin()+1,la.begin()+1+lena,a[i])-la.begin();
        b[i]=lower_bound(lb.begin()+1,lb.begin()+1+lenb,b[i])-lb.begin()+lena;
        g[a[i]].emplace(b[i],i),g[b[i]].emplace(a[i],i);
    }
    vector<pair<int,int>>res; int ret=0,edge=0;
    function<void(int,int)>dfs=[&](int x,int f,int fid=0) {
        vis[x]=siv[x]=true,edge+=g[x].size();
        vector<int>son; vector<pair<int,int>>bac; auto tmp=g[x];
        for(auto [to,id]: tmp) {
            if(to==f) { fid=id; continue; }
            if(vis[to]) { bac.emplace_back(to,id); continue; } dfs(to,x);
            if(g[x].find({to,id})!=g[x].end()) son.push_back(id);
        }
        vis[x]=false;
        if(!(son.size()&1)) {
            for(int i=1;i<son.size();i+=2)
                res.emplace_back(son[i-1],son[i]);
            for(int i=1;i<bac.size();i+=2) {
                res.emplace_back(bac[i-1].se,bac[i].se);
                int fir=bac[i-1].fi,sec=bac[i].fi;
                g[fir].erase({x,bac[i-1].se});
                g[sec].erase({x,bac[i].se});
            }
            if((bac.size()&1)&&fid) {
                int fir=bac.back().fi,sec=bac.back().se;
                res.emplace_back(sec,fid);
                g[fir].erase({x,sec}),g[f].erase({x,fid});
            }
        } else {
            for(int i=1;i<son.size();i+=2)
                res.emplace_back(son[i-1],son[i]);
            if(!bac.empty()) {
                res.emplace_back(son.back(),bac[0].se);
                g[bac[0].fi].erase(make_pair(x,bac[0].se));
                for(int i=2;i<bac.size();i+=2) {
                    res.emplace_back(bac[i-1].se,bac[i].se);
                    int fir=bac[i-1].fi,sec=bac[i].fi;
                    g[fir].erase({x,bac[i-1].se}),g[sec].erase({x,bac[i].se});
                }
                if(!(bac.size()&1)&&fid) {
                    int fir=bac.back().fi,sec=bac.back().se;
                    res.emplace_back(sec,fid);
                    g[fir].erase({x,sec}),g[f].erase({x,fid});
                }
            } else if(fid) {
                res.emplace_back(son.back(),fid);
                g[f].erase({x,fid});
            }
        }
    };
    for(int i=1;i<=lena+lenb;++i)
        if(!siv[i]) edge=0,dfs(i,0),ret+=edge>>2;
    // if(ret!=res.size()) assert(ret>res.size());
    cout<<ret<<"\n"; vector<int>vec; vis.assign((n<<2)+10,false);
    for(auto [x,y]: res) cout<<x<<' '<<y<<"\n",vis[x]=vis[y]=true;
    for(int i=1;i<=(n<<1);++i) if(!vis[i]) vec.push_back(i);
    for(int i=1;i<vec.size();i+=2) cout<<vec[i-1]<<' '<<vec[i]<<"\n";
}

auto main() -> signed {
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    int T; cin>>T;
    while(T--) solve();
    return 0;
}