#include <bits/stdc++.h>
using namespace std;
const int MOD = 1000000007;
long long POW(long long a, long long b, long long MMM = MOD) {
  long long ret = 1;
  for (; b; b >>= 1, a = (a * a) % MMM)
    if (b & 1) ret = (ret * a) % MMM;
  return ret;
}
long long gcd(long long a, long long b) { return b ? gcd(b, a % b) : a; }
long long lcm(long long a, long long b) {
  if (a == 0 || b == 0) return a + b;
  return a * (b / gcd(a, b));
}
int dx[] = {0, 1, 0, -1, 1, 1, -1, -1}, dy[] = {1, 0, -1, 0, 1, -1, 1, -1};
int p[1000001], d[1000001][21];
int main() {
  for (int(i) = (2); (i) <= (1000000); (i) += (1)) {
    for (int(j) = (i); (j) <= (1000000); (j) += (i))
      if (!p[j]) p[j] = i;
  }
  for (int(r) = (0); (r) <= (1000000); (r) += (1)) d[r][0] = 1;
  for (int(i) = (1); (i) <= (20); (i) += (1)) d[0][i] = 2;
  for (int(i) = (1); (i) <= (1000000); (i) += (1))
    for (int(j) = (1); (j) <= (20); (j) += (1))
      d[i][j] = (d[i][j - 1] + d[i - 1][j]) % MOD;
  int q;
  scanf("%d", &q);
  while (q--) {
    int r, n;
    scanf("%d%d", &r, &n);
    long long ans = 1;
    vector<pair<int, int> > v;
    while (n > 1) {
      int x = p[n];
      if (!v.empty() && v.back().first == x)
        v.back().second++;
      else
        v.push_back(make_pair(x, 1));
      n /= x;
    }
    for (pair<int, int> p : v) ans = (ans * d[r][p.second]) % MOD;
    printf("%lld\n", ans);
  }
}
