#include using namespace std; using ll = long long; const ll mod = 998'244'353; //const ll mod = 1'000'000'007; //const ll mod = 67'280'421'310'721; struct mint{ long long x; mint(long long x=0):x((x%mod+mod)%mod){} mint operator-() const{ return mint(-x); } mint& operator+=(const mint& a){ if((x+=a.x)>=mod)x-=mod; return *this; } mint& operator-=(const mint& a){ if((x+=mod-a.x)>=mod)x-=mod; return *this; } mint& operator*=(const mint& a){ (x *= a.x) %= mod; return *this; } mint operator+(const mint& a) const{ mint res(*this); return res+=a; } mint operator-(const mint& a) const{ mint res(*this); return res-=a; } mint operator*(const mint& a) const{ mint res(*this); return res*=a; } mint pow(long long n) const { assert(0 <= n); mint a = *this, r = 1; while (n) { if (n & 1) r *= a; a *= a; n >>= 1; } return r; } mint inv() const{ return pow(mod-2); } mint& operator/=(const mint& a){ return (*this)*=a.inv(); } mint operator/(const mint& a) const { mint res(*this); return res/=a; } friend ostream& operator<<(ostream& os, const mint& m){ os << m.x; return os; } bool operator==(const mint& a) const { return x == a.x; } bool operator<(const mint& a) const{ return x < a.x; } }; ll calc(ll n,ll p){ ll cnt = 0; while(n){ cnt += n / p; n /= p; } return cnt; } const unsigned long long numset[] = {2, 325, 9375, 28178, 450775, 9780504, 1795265022ULL}; unsigned long long mod_mul(unsigned long long a, unsigned long long b, unsigned long long mod) { long long ret = a * b - mod * (unsigned long long)((long double)(a) * (long double)(b) / (long double)(mod)); return ret + mod * (ret < 0) - mod * (ret >= (ll)mod); } unsigned long long mod_pow(unsigned long long x, unsigned long long k, unsigned long long mod){ unsigned long long res = 1; while(k){ if(k & 1) res = mod_mul(res, x, mod); x = mod_mul(x, x, mod); k >>= 1; } return res; } bool check(unsigned long long n){ if(n < 2 || ((n % 6 != 1) && (n % 6 != 5))) return (n == 2) || (n == 3); unsigned long long d = n - 1, s = 0; while(d % 2 == 0){ d /= 2; s++; } for(unsigned long long a : numset){ if(a % n == 0) return true; unsigned long long res = mod_pow(a, d, n); if(res == 1) continue; bool ok = true; for(unsigned int r = 0; r < s; r++){ if(res == n-1){ ok = false; break; } res = mod_mul(res, res, n); } if(ok) return false; } return true; } const int mx = 1e6; vector p; int cnt[mx+1]; int main(){ int n; cin>>n; while(n--){ ll now; cin>>now; if(check(now)) cout<>n>>k; for(int i = 0;i=0); if(cnt) ans *= mint(cnt+1); } for(ll i = n;i>n-k;i--) if(i>=mx){ if(check(i)) { ans *= calc(n,i) + 1; } } cout<