#include <bits/stdc++.h>
const int P = 1000000007, MAXN = 1000031, MAXL = 21;
bool ip[MAXN] = {1, 1};
int p[MAXN], l[MAXN] = {0, 1}, a[MAXN] = {0, 0};
void sieve() {
  for (int i = 2, c = 0; i < MAXN; i++) {
    if (!ip[i]) {
      p[c++] = i;
      l[i] = 1;
      a[i] = 1;
    }
    for (int j = 0; j < c && i * p[j] < MAXN; j++) {
      ip[i * p[j]] = 1;
      if (i % p[j]) {
        l[i * p[j]] = i;
        a[i * p[j]] = 1;
      } else {
        l[i * p[j]] = l[i];
        a[i * p[j]] = a[i] + 1;
        break;
      }
    }
  }
}
int f[MAXN][MAXL];
void pre() {
  std::fill(f[0], f[0] + MAXL, 2);
  f[0][0] = 1;
  for (int i = 1; i < MAXN; i++) {
    f[i][0] = 1;
    for (int j = 1; j < MAXL; j++) {
      f[i][j] = (f[i - 1][j] + f[i][j - 1]) % P;
    }
  }
}
int solve(int r, int n) {
  long long ans = 1;
  while (n > 1) {
    ans = ans * f[r][a[n]] % P;
    n = l[n];
  }
  return ans;
}
int main() {
  sieve();
  pre();
  int q = 0;
  scanf("%d", &q);
  while (q--) {
    int r = 0, n = 0;
    scanf("%d%d", &r, &n);
    printf("%d\n", solve(r, n));
  }
  return 0;
}
