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

typedef long long ll;
const int md = 1e9 + 7;
const int mod = 1e9 + 7;

inline int add(int x, int y) {
	if (x + y >= md) return x + y - md;
	return x + y;
}

inline int sub(int x, int y) {
	if (x < y) return x - y + md;
	return x - y;
}

inline int mul(int x, int y) {
	return 1ull * x * y % md;
}
inline int mul(int x,int y,int z) {return mul(mul(x,y),z);}

inline int qpow(int x, int y) {
	int ans = 1;
	while (y) {
		if (y & 1) ans = mul(ans, x);
		y >>= 1; x = mul(x, x);
	}
	return ans;
}
int Inv(int x) {return qpow(x,mod-2);}
struct Node
{
	int siz;
	vector<pair<int,int>> eg;
	Node(int s=0)
	{
		assert(s<=1);
		siz=s;
	}
};
const int L=200;
int n,dv[1005][1005];
unordered_map<int,vector<int>> f[L+1][L+1],g[L+1];
void init()
{
	for(int i=0;i<=n;i++) for(int j=0;j<=i;j++) f[i][j].clear();
	for(int i=0;i<=n;i++) g[i].clear();
	f[0][0][1]=vector<int>(0);
	for(int i=1;i<=n/2+2;i++)
	{
		for(int j=1;j<=i;j++)
		{
			int adw=dv[n-2][j];
			for(int k=0;k<=j;k++)
			{
				for(auto [w,v]:f[i-j][k])
				{
					int gw=mul(adw,w);
					if(f[i][j].find(gw)==f[i][j].end())
					{
						f[i][j][gw]=f[i-j][k][w];
						f[i][j][gw].push_back(j);
					}
				}
			}
		}
		for(int j=1;j<=i;j++)
		{
			for(auto [w,v]:f[i][j])
			{
				if(g[i].find(w)==g[i].end())
				{
					g[i][w]=v;
				}
			}
		}
	}
}
int Q,q[15],ok[15],sok;
vector<int> ans[15];
int ansn[15];
int main() {
 #ifdef xay5421
	freopen("a.in","r",stdin);
 #endif
	for(int i=1;i<=1000;i++)
	{
		dv[i][0]=1;
		for(int j=1;j<=1000;j++) dv[i][j]=mul(dv[i][j-1],i-j+1);
	}
	cin>>Q;
	for(int i=1;i<=Q;i++)
	{
		cin>>q[i];
		if(q[i]==1) ok[i]=1,sok++;
		if(q[i]==2)
		{
			ok[i]=1,sok++;
		}
	}
	for(n=3;n<=L&&sok<Q;n++)
	{
		fprintf(stderr,"n = %d\n",n);
		init();
		int adw=mul(n,n-1);
		for(int qid=1;qid<=Q;qid++)
		{
			// q[qid] sum=n-2
			for(int big=1;big<=n-2&&!ok[qid];big++)
			{
				int cw=dv[n-2][big];
				for(int L=0;L<=big&&!ok[qid];L++)
				{
					int R=n-2-big-L;
					// printf("%d %d %d\n",big,L,R);
					if(!(R>=0)) continue;
					for(auto [w,v]:g[L])
					{
						int way=mul(adw,w,cw);
						int ned=mul(Inv(way),q[qid]);
						if(g[R].find(ned)!=g[R].end())
						{
							ans[qid]=g[R][ned];
							for(int i:v) ans[qid].push_back(i);
							ans[qid].push_back(big);
							ok[qid]=1;
							ansn[qid]=n;
							sok++;
							break;
						}
					}
				}
			}
		}
	}
	for(int i=1;i<=Q;i++)
	{
		if(q[i]==1) printf("1\n");
		else if(q[i]==2) printf("2\n1 2\n");
		else
		{
			printf("%d\n",ansn[i]);
			while(ans[i].size()<ansn[i]) ans[i].push_back(0);
	//		for(int j=0;j<(int)ans[i].size();j++) ans[i][j]++;
			for(int j=0;j<ansn[i]-1;j++)
			{
				for(int c=0;c<(int)ans[i].size();c++) for(int d=0;d<(int)ans[i].size();d++)
				{
					if(ans[i][c]==0&&ans[i][d]>0)
					{
						ans[i][c]--,ans[i][d]--;
						printf("%d %d\n",c+1,d+1);
					}
				}
			}
			for(int c=0;c<(int)ans[i].size();c++) for(int d=c+1;d<(int)ans[i].size();d++)
			{
				if(ans[i][c]==0&&ans[i][d]==0)
				{
					ans[i][c]--,ans[i][d]--;
					printf("%d %d\n",c+1,d+1);
				}
			}
		}
	}
	return 0;
}