#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1e6 + 10, MOD = 1e9 + 7;
long long f[MAXN * 2], fr[MAXN + 100], ans;
vector<pair<int, int> > pi;
int q, r, n, lst[MAXN];
long long rev(long long a) {
  int t = MOD - 2;
  long long res = 1;
  for (t; t; t /= 2) {
    if (t % 2) res *= a, res %= MOD;
    a *= a, a %= MOD;
  }
  return res;
}
long long ch(int n, int k) { return f[n] * fr[n - k] % MOD * fr[k] % MOD; }
void getd(int num) {
  pi.clear();
  while (num != 1) {
    int cl = lst[num], c = 0;
    while (num % cl == 0) num /= cl, c++;
    pi.push_back({cl, c});
  }
}
void solve() {
  cin >> r >> n;
  ans = 1;
  getd(n);
  if (r == 0) ans = (1 << pi.size());
  for (int i = 0; i < pi.size() && r > 0; i++) {
    ans *= (ch(pi[i].second - 1 + r, r) * 2 + ch(pi[i].second + r - 1, r - 1)) %
           MOD;
    ans %= MOD;
  }
  cout << ans << '\n';
}
int main() {
  ios::sync_with_stdio(false);
  cin.tie(0), cout.tie(0);
  f[0] = 1;
  for (int i = 1; i < 100 + MAXN; i++) f[i] = f[i - 1] * i, f[i] %= MOD;
  for (int i = 0; i < 100 + MAXN; i++) fr[i] = rev(f[i]);
  lst[1] = 1;
  for (int j = 2; j < MAXN; j += 2) lst[j] = 2;
  for (int i = 3; i < MAXN; i += 2)
    for (int j = i; j < MAXN; j += 2 * i)
      if (lst[j] == 0) lst[j] = i;
  cin >> q;
  while (q--) solve();
}
