#include <algorithm>
#include <cstdint>
#include <iostream>
#include <numeric>
#include <vector>

using namespace std;

constexpr int mod = 1e9 + 7;

auto& mul(auto&& a, auto b) { return a = a * (uint64_t) b % mod; }
int inv(int a) { for (int r = 1, b = mod - 2; ; b /= 2, mul(a, a)) if (!b) return r; else if (b % 2) mul(r, a); }

using namespace std;

template <typename K, typename V, int mlf = 8>
struct HT {
  int capacity;
  vector<pair<K, int>> knx;
  vector<V> values;
  vector<int> head;

  HT(): capacity{4 * mlf}, head(capacity, -1) {
  }

  void emplace(K k, V v) {
//    if (~find(k)) return;
    knx.push_back({k, head[k & capacity - 1]});
    values.push_back(v);
    head[k & capacity - 1] = knx.size() - 1;
    maybe_rehash();
  }

  int find(K k) {
    for (int i = head[k & capacity - 1]; ~i; i = knx[i].second) if (knx[i].first == k) return i;
    return -1;
  }

  void clear() {
    knx.clear();
    values.clear();
    head.assign(capacity, -1);
  }

  void maybe_rehash() {
    if (knx.size() * mlf <= capacity) return;
    head.assign(capacity *= 2, -1);
    for (int i = 0; i < knx.size(); ++i) {
      knx[i].second = head[knx[i].first & capacity - 1];
      head[knx[i].first & capacity - 1] = i;
    }
  }
};

vector<HT<int, vector<int>>> even(1);
vector<int> inv_(1);

void gen_even(int n0, int n, int c, int p, auto&& cur) {
  even[n0 - n].emplace(p, cur);
  for (; c > 0; c -= 2) if (c <= n) {
    int np = p;
    for (int i = 0; i < c; ++i) mul(np, inv_[n0 - i]);
    cur.push_back(c);
    gen_even(n0, n - c, c, np, cur);
    cur.pop_back();
  }
}

vector<int> qs;
vector<vector<int>> ans;

void solve(int n0, int n, int c, int p, auto&& cur) {
  if (n % 2 == 0)
  for (int i = 0; i < qs.size(); ++i) if (ans[i].empty() && ~even[n].find(mul(+p, qs[i]))) {
    ans[i] = cur;
    auto& t = even[n].values[even[n].find(mul(+p, qs[i]))];
    ans[i].insert(ans[i].end(), t.begin(), t.end());
  }
  for (; c > 0; c -= 2) if (c <= n) {
    int np = p;
    for (int i = 0; i < c; ++i) mul(np, n0 - i);
    cur.push_back(c);
    solve(n0, n - c, c, np, cur);
    cur.pop_back();
  }
}

int main() {
  int tc; cin >> tc;
  qs.resize(tc);
  ans.resize(tc);
  for (int i = 0; i < tc; ++i) {
    cin >> qs[i];
    if (qs[i] == 1) ans[i] = {-1};
    else if (qs[i] == 2) ans[i] = {-1};
    qs[i] = inv(qs[i]);
  }
  for (int n = 1, done = 0; !done++; ++n) {
    inv_.push_back(inv(n));
    even.push_back({});
    for (auto& x: even) x.clear();
    gen_even(n, n, n + n % 2, 1, vector<int>{});
    solve(n, n, n - n % 2 + 1, mul(n + 1, n + 2), vector<int>{});
    for (int i = 0; i < tc; ++i) done &= !ans[i].empty();
  }
  for (int i = 0; i < tc; ++i) {
    if (qs[i] == 1) cout << "1\n";
    else {
      if (ans[i][0] == -1) ans[i] = {};
      cout << accumulate(ans[i].begin(), ans[i].end(), 0) + 2 << '\n' << "1 2\n";
      for (int j = 0, k = 2; j < ans[i].size(); ++j)
      while (ans[i][j]--) cout << j + 1 << ' ' << ++k << '\n';
    }
  }
}