#include <bits/stdc++.h>
#define x first
#define y second
using namespace std;
using ll=long long;
using pii=pair<int,int>;
using pll=pair<ll,ll>;
using vi=vector<int>;
using vl=vector<ll>;
using ld=long double;
#define all(a) begin(a),end(a)
#define pb push_back

const int N=310,MOD=1e9+7;
const char en='\n';
const ll LLINF=1ll<<60;

int mul(int a,int b)
{
	return a*1ll*b%MOD;
}

int pot(int n,int k)
{
	if (k==0) return 1;
	int r=pot(n,k/2);
	r=mul(r,r);
	if (k%2) r=mul(r,n);
	return r;
}

int inv(int x)
{
	return pot(x,MOD-2);
}

int t,rrem,rem,n,ms,efn;
vi mog;

vi cur;
int cuu;
vector<pair<int,vi>> svp[N][N];
bool bil;
vi ans;

void bact(int s,int la)
{
	if (mul(cuu,mog[efn-s])==rem)
	{
		ans=cur;
		ans.pb(efn-s);
		//cout<<s<<en;
		//for (auto x: ans) cout<<x<<' ';
		//cout<<en;
		return;
	}
	if (s>=efn/3)
	{
		if (!bil) svp[efn][s].pb({cuu,cur});
		else return;
	}
	if (s<=(efn-1)/2) for (int ne=la;ne>=1;--ne) if (s+ne<=ms)
	{
		cur.pb(ne);
		int scuu=cuu;
		cuu=mul(cuu,mog[ne]);
		bact(s+ne,ne);
		cuu=scuu;
		cur.pop_back();
	}
}

void out()
{
	cout<<n<<endl;
	cout<<1<<' '<<2<<endl;
	int cu=2;
	for (auto x: ans)
	{
		for (int i=1;i<=x;++i) cout<<cu<<' '<<cu+i<<endl;
		cu+=x;
	}
}

bool pok()
{
	rem=mul(rrem,mul(inv(n),inv(n-1)));
	mog.clear();
	mog.pb(1);
	mog.pb(n-2);
	while (mog.back()>1) mog.pb(mog.back()-1);
	for (int i=1;i<(int)mog.size();++i) mog[i]=mul(mog[i],mog[i-1]);
	//cout<<n<<' '<<rem<<endl;
	//for (auto x: mog) cout<<x<<' ';
	//cout<<endl<<endl;
	ans.clear();
	efn=n-2;
	ms=(efn*2+2)/3;
	cur.clear();
	cuu=1;
	//cout<<n<<endl;
	bil=svp[efn][efn/3].size()>0;
	bact(0,efn);
	int s1=0,s2=0;
	if (!bil) for (int i=efn/3;i<=ms;++i)
	{
		if (i<=efn/2) s1+=svp[efn][i].size();
		else s2+=svp[efn][i].size();
		sort(all(svp[efn][i]));
	}
	//cout<<n<<' '<<s1<<' '<<s2<<endl;
	if (ans.size())
	{
		out();
		return 1;
	}
	int jed=0;
	for (int sz=efn/2;sz<=ms;++sz) for (int i=0;i<(int)svp[efn][sz].size();++i)
	{
		if (i && svp[efn][sz][i].x==svp[efn][sz][i-1].x)
		{
			++jed;
			continue;
		}
		int potr=mul(rem,inv(svp[efn][sz][i].x));
		auto it=lower_bound(all(svp[efn][efn-sz]),make_pair(potr,vi(0)));
		if (it==svp[efn][efn-sz].end()) continue;
		if (it->x!=potr) continue;
		//cout<<sz<<en;
		//for (auto x: svp[sz][i].y) cout<<x<<' ';
		//cout<<en;
		//cout<<efn-sz<<en;
		//for (auto x: it->y) cout<<x<<' ';
		//cout<<en;
		ans=svp[efn][sz][i].y;
		for (auto x: it->y) ans.pb(x);
		out();
		return 1;
	}
	//cout<<jed<<endl;
	return 0;
}

int main()
{
	ios_base::sync_with_stdio(0);
	cin.tie(0);
	cin>>t;
	while (t--)
	{
		cin>>rrem;
		if (rrem==1)
		{
			cout<<1<<en;
			continue;
		}
		if (rrem==2)
		{
			cout<<2<<en<<"1 2"<<en;
			continue;
		}
		for (n=3;;++n)
		{
			if (pok()) break;
		}
	}
}