/* -*- coding: utf-8 -*- * * 3776.cc: No.3776 Magic Square of Product - yukicoder */ #include #include #include #include using namespace std; /* constant */ const int MAX_K = 1000000; const int MAX_S = 20; const int MOD = 998244353; /* typedef */ using vb = vector; using vi = vector; using pii = pair; using vpii = vector; template struct MI { int v; MI(): v() {} MI(int _v): v(_v % MOD) { if (v < 0) v += MOD; } MI(long long _v): v(_v % MOD) { if (v < 0) v += MOD; } explicit operator int() const { return v; } MI operator+(const MI m) const { return MI(v + m.v); } MI operator-(const MI m) const { return MI(v + MOD - m.v); } MI operator-() const { return MI(MOD - v); } MI operator*(const MI m) const { return MI((long long)v * m.v); } MI &operator+=(const MI m) { return (*this = *this + m); } MI &operator-=(const MI m) { return (*this = *this - m); } MI &operator*=(const MI m) { return (*this = *this * m); } bool operator==(const MI m) const { return v == m.v; } bool operator!=(const MI m) const { return v != m.v; } MI pow(int n) const { // a^n % MOD MI pm = 1, a = *this; while (n > 0) { if (n & 1) pm *= a; a *= a; n >>= 1; } return pm; } MI inv() const { return pow(MOD - 2); } MI operator/(const MI m) const { return *this * m.inv(); } MI &operator/=(const MI m) { return (*this = *this / m); } }; using mi = MI; /* global variables */ vb primes; vi pnums; int cs[MAX_S + 1]; /* subroutines */ int gen_primes(int maxp) { primes.assign(maxp + 1, true); primes[0] = primes[1] = false; pnums.clear(); int p; for (p = 2; p * p <= maxp; p++) if (primes[p]) { pnums.push_back(p); for (int q = p * p; q <= maxp; q += p) primes[q] = false; } for (; p <= maxp; p++) if (primes[p]) pnums.push_back(p); return (int)pnums.size(); } vpii prime_decomp(int n) { vpii pds; for (auto pi: pnums) { if (pi * pi > n) { if (n > 1) pds.push_back(pii(n, 1)); break; } if (n % pi == 0) { int fi = 0; while (n % pi == 0) n /= pi, fi++; pds.push_back(pii(pi, fi)); } } return pds; } int countmagic(int s) { int cnt = 0; for (int i0 = 0; i0 <= s; i0++) for (int i1 = 0; i0 + i1 <= s; i1++) { int i2 = s - (i0 + i1); for (int i3 = 0; i0 + i3 <= s; i3++) { int i6 = s - (i0 + i3); int max13 = max(i1, i3); for (int i4 = 0; max13 + i4 <= s; i4++) { int i5 = s - (i3 + i4); int i7 = s - (i1 + i4); if (i2 + i5 <= s && i2 + i5 == i6 + i7) cnt++; } } } //printf(" countmgic(%d) = %d\n", s, cnt); return cnt; } /* main */ int main() { for (int s = 0; s <= MAX_S; s++) cs[s] = countmagic(s); gen_primes(MAX_K); int tn; scanf("%d", &tn); while (tn--) { int k; scanf("%d", &k); auto pds = prime_decomp(k); mi c = 1; for (auto [p, d]: pds) c *= cs[d]; printf("%d\n", (int)c); } return 0; }