#include <bits/stdc++.h>
const double PI = acos(-1.0);
using namespace std;
const long long mod = 1000000007;
const long long MAXR = 1e6 + 5;
const long long MAXE = 25;
long long r, n;
long long f[MAXR][MAXE];
bool vis[MAXR];
vector<long long> prime;
void setP() {
  long long m = sqrt(MAXR);
  for (long long i = 2; i <= m; i++)
    if (!vis[i])
      for (long long j = i * i; j < MAXR; j += i) vis[j] = true;
  for (long long i = 2; i < MAXR; i++)
    if (!vis[i]) prime.push_back(i);
}
void setD(long long r) {
  for (long long ri = 0; ri <= r; ri++) f[ri][0] = 1;
  for (long long ei = 1; ei < MAXE; ei++) f[0][ei] = 2;
  for (long long ri = 1; ri <= r; ri++)
    for (long long ei = 1; ei < MAXE; ei++)
      f[ri][ei] = (f[ri][ei - 1] + f[ri - 1][ei]) % mod;
}
vector<long long> pf;
vector<long long> exps;
void setPE(long long x) {
  for (long long i = 0; i < prime.size(); i++) {
    long long pi = prime[i];
    if (x % pi == 0) {
      pf.push_back(pi);
      long long ei = 0;
      while (x % pi == 0) {
        x /= pi;
        ei++;
      }
      exps.push_back(ei);
    }
    if (prime[i] * prime[i] > x) break;
  }
  if (x > 1) pf.push_back(x), exps.push_back(1);
}
vector<long long> vr, vn;
int main() {
  long long q;
  setP();
  scanf("%lld", &q);
  long long mr = 0;
  vr.clear();
  vn.clear();
  for (long long i = 0; i < q; i++) {
    long long tr, tn;
    scanf("%lld%lld", &tr, &tn);
    vr.push_back(tr);
    vn.push_back(tn);
    mr = max(mr, vr[i]);
  }
  setD(mr + 3);
  for (long long i = 0; i < q; i++) {
    r = vr[i];
    n = vn[i];
    pf.clear();
    exps.clear();
    setPE(n);
    long long ans = 1;
    for (long long i = 0; i < exps.size(); i++)
      ans = (ans * f[r][exps[i]]) % mod;
    printf("%lld\n", ans);
  }
  return 0;
}
