Problem: Polycarpus has t safes. The password for each safe is a square matrix consisting of decimal digits '0' ... '9' (the sizes of passwords to the safes may vary). Alas, Polycarpus has forgotten all passwords, so now he has to restore them.

Polycarpus enjoys prime numbers, so when he chose the matrix passwords, he wrote a prime number in each row of each matrix. To his surprise, he found that all the matrices turned out to be symmetrical (that is, they remain the same after transposition). Now, years later, Polycarp was irritated to find out that he remembers only the prime numbers pi, written in the first lines of the password matrices.

For each safe find the number of matrices which can be passwords to it.

The number of digits in pi determines the number of rows and columns of the i-th matrix. One prime number can occur in several rows of the password matrix or in several matrices. The prime numbers that are written not in the first row of the matrix may have leading zeros.

Input

The first line of the input contains an integer t (1 ≤ t ≤ 30) — the number of safes. Next t lines contain integers pi (10 ≤ pi ≤ 99999), pi is a prime number written in the first row of the password matrix for the i-th safe. All pi's are written without leading zeros.

Output

Print t numbers, the i-th of them should be the number of matrices that can be a password to the i-th safe. Print the numbers on separate lines.

Examples

Input

4
11
239
401
9001


Output

4
28
61
2834

Note

Here is a possible password matrix for the second safe: 
    
    
      
    239  
    307  
    977  
    

Here is a possible password matrix for the fourth safe: 
    
    
      
    9001  
    0002  
    0002  
    1223   
    

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
double pi = asin(1.0) * 2;
struct point {
  double x, y;
  point() {}
  point(double _x, double _y) {
    x = _x;
    y = _y;
  }
  inline void read() { scanf("%lf%lf", &x, &y); }
  inline void show() { printf("%lf %lf\n", x, y); }
  inline double dis() { return x * x + y * y; }
  inline double dist() { return sqrt(x * x + y * y); }
};
inline point operator+(const point &a, const point &b) {
  return point(a.x + b.x, a.y + b.y);
}
inline point operator-(const point &a, const point &b) {
  return point(a.x - b.x, a.y - b.y);
}
inline point operator*(const point &a, const double &b) {
  return point(a.x * b, a.y * b);
}
inline point operator/(const point &a, const double &b) {
  return point(a.x / b, a.y / b);
}
inline double det(const point &a, const point &b) {
  return a.x * b.y - a.y * b.x;
}
inline double dot(const point &a, const point &b) {
  return a.x * b.x + a.y * b.y;
}
int p[100000], np;
bool f[110000];
int num[110000];
int getans(int n) {
  int ans = 0;
  int a = n;
  if (n < 100) {
    ans += num[n % 10 * 10];
    return ans;
  }
  if (n < 1000) {
    for (int z = lower_bound(p, p + np, n % 100 / 10 * 100) - p;
         z < np && p[z] < n % 100 / 10 * 100 + 100; z++) {
      ans += num[n % 10 * 100 + p[z] % 10 * 10];
    }
    return ans;
  }
  if (n < 10000) {
    int lz = n % 1000 / 100 * 1000, rz = lz + 1000;
    for (int z = lower_bound(p, p + np, lz) - p; z < np && p[z] < rz; z++) {
      int b = p[z];
      int lq = n % 100 / 10 * 1000 + b % 100 / 10 * 100, rq = lq + 100;
      for (int q = lower_bound(p, p + np, lq) - p; q < np && p[q] < rq; q++) {
        int c = p[q];
        ans += num[n % 10 * 1000 + b % 10 * 100 + c % 10 * 10];
      }
    }
    return ans;
  }
  int cnt = 0;
  if (n < 100000) {
    int lz = n % 10000 / 1000 * 10000, rz = lz + 10000;
    for (int z = lower_bound(p, p + np, lz) - p; z < np && p[z] < rz; z++) {
      int b = p[z];
      int lq = n % 1000 / 100 * 10000 + b % 1000 / 100 * 1000, rq = lq + 1000;
      for (int q = lower_bound(p, p + np, lq) - p; q < np && p[q] < rq; q++) {
        int c = p[q];
        int lw =
                n % 100 / 10 * 10000 + b % 100 / 10 * 1000 + c % 100 / 10 * 100,
            rw = lw + 100;
        for (int w = lower_bound(p, p + np, lw) - p; w < np && p[w] < rw; w++) {
          ans += num[n % 10 * 10000 + b % 10 * 1000 + c % 10 * 100 +
                     p[w] % 10 * 10];
        }
      }
    }
    return ans;
  }
}
int main() {
  np = 0;
  f[0] = f[1] = true;
  memset(f, false, sizeof(f));
  for (int i = 2; i < 110000; i++)
    if (!f[i]) {
      p[np++] = i;
      for (int j = i + i; j < 110000; j += i) f[j] = true;
    }
  for (int i = 0; i < 100000; i++)
    if (i % 10 == 0)
      for (int j = 0; j < 10; j++) num[i] += (!f[i + j]);
  int t;
  for (scanf("%d", &t); t--;) {
    int n;
    scanf("%d", &n);
    cout << getans(n) << endl;
  }
  return 0;
}