#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int mod=1000000007;
int ksm(ll a,int b,int c=1){
    for(;b;b/=2,a=a*a%mod)
        if(b&1)c=c*a%mod;
    return c;
}
const int N=210;

vector<pair<int,int> >res;

int a[N];

int g[N],ig[N];
int last_ans;

const int B=5;

vector<int>resL_list,resR_list;

int ans_flag;
bool dfsL(int s,int x,int f,int n,int l){
    if(s==0){
        if(ans_flag){
            if(f!=ans_flag)
                return 0;
            for(int i=1;i<x;++i)
                res.emplace_back(a[i]+1,res.size()+1);
            return 1;
        }else{
            resL_list.push_back(f);
            return 0;
        }
    }
    for(int i=min(s,l);i>=1;--i){
        a[x]=i;
        if(dfsL(s-i,x+1,(ll)f*ig[i]%mod,n,i))
            return 1;
    }
    return 0;
}

bool dfsR(int s,int x,int f,int n,int l){
    if(s==0){
        if(ans_flag){
            if(f!=ans_flag)
                return 0;
            for(int i=1;i<x;++i)
                res.emplace_back(a[i]+1,res.size()+1);
            return 1;
        }else{
            resR_list.push_back(f);
            return 0;
        }
    }
    for(int i=l;i<=s;++i){
        a[x]=i;
        if(dfsR(s-i,x+1,(ll)f*g[i]%mod,n,i))
            return 1;
    }
    return 0;
}


vector<pair<int,int> >res_list[N];
vector<pair<int,int> >x_list;
void output(vector<pair<int,int> >res){
    int n=res.size();
   // for(auto i:res)
     //   cout<<i.first<<' '<<i.second<<"!\n";
	cout<<n<<"\n";
    for(;;){
        if(res.size()<=1)
            return;
        sort(res.begin(),res.end());
        reverse(res.begin(),res.end());
        cout<<res[0].second<<' '<<res[res.size()-1].second<<'\n';
        res[0].first--;
        res.pop_back();
    }
}

void solve(){
    if(x_list.empty())
        return;
    for(int i=2;;++i){
        assert(i<N);
		int iv=ksm(i*(i-1),mod-2);
        g[0]=1;
        for(int j=1;j<=i;++j)
            g[j]=(ll)g[j-1]*(i-j-1)%mod;
        for(int j=0;j<=i;++j)
            ig[j]=ksm(g[j],mod-2);
        for(int j=0;j<=i-2;++j){
			ans_flag=0;
            resL_list.clear();
            resR_list.clear();
            dfsL(j,1,1,i,B);
            dfsR(i-2-j,1,1,i,B+1);
            sort(resR_list.begin(),resR_list.end());
            for(int k:resL_list){
                for(int x=0;x<(int)x_list.size();++x){
                    int rx=(ll)iv*x_list[x].first%mod;
                    int u=(ll)k*rx%mod;
                    auto it=lower_bound(resR_list.begin(),resR_list.end(),u);
                    if(it!=resR_list.end()&&*it==u){
                        res.clear();
                        ans_flag=k;
                        dfsL(j,1,1,i,B);
                        ans_flag=u;
                        dfsR(i-2-j,1,1,i,B+1);
                        while(res.size()<i)
                            res.emplace_back(0,res.size()+1);
                        res_list[x_list[x].second]=res;
                        x_list.erase(x_list.begin()+x);
                        --x;
                    }
                }
                if(x_list.empty())
                    return;
            }
        }
    }
}
int T;
int main(){
    ios::sync_with_stdio(0);cin.tie(0);
    cin>>T;

    for(int i=1;i<=T;++i){
        int x;cin>>x;
        if(x!=1)
            x_list.emplace_back(x,i);
        else
            res_list[i].emplace_back(0,1);
    }
    solve();
    for(int i=1;i<=T;++i)
        output(res_list[i]);
    return 0;
}