#ifdef NACHIA #define _GLIBCXX_DEBUG #else // disable assert #define NDEBUG #endif #include #include #include #include using namespace std; using ll = long long; const ll INF = 1ll << 60; #define REP(i,n) for(ll i=0; i using V = vector; template void chmax(A& l, const B& r){ if(l < r) l = r; } template void chmin(A& l, const B& r){ if(r < l) l = r; } #include #include namespace nachia{ // ax + by = gcd(a,b) // return ( x, - ) std::pair ExtGcd(long long a, long long b){ long long x = 1, y = 0; while(b){ long long u = a / b; std::swap(a-=b*u, b); std::swap(x-=y*u, y); } return std::make_pair(x, a); } } // namespace nachia namespace nachia{ template struct StaticModint{ private: using u64 = unsigned long long; unsigned int x; public: using my_type = StaticModint; template< class Elem > static Elem safe_mod(Elem x){ if(x < 0){ if(0 <= x+MOD) return x + MOD; return MOD - ((-(x+MOD)-1) % MOD + 1); } return x % MOD; } StaticModint() : x(0){} StaticModint(const my_type& a) : x(a.x){} StaticModint& operator=(const my_type&) = default; template< class Elem > StaticModint(Elem v) : x(safe_mod(v)){} unsigned int operator*() const { return x; } my_type& operator+=(const my_type& r) { auto t = x + r.x; if(t >= MOD) t -= MOD; x = t; return *this; } my_type operator+(const my_type& r) const { my_type res = *this; return res += r; } my_type& operator-=(const my_type& r) { auto t = x + MOD - r.x; if(t >= MOD) t -= MOD; x = t; return *this; } my_type operator-(const my_type& r) const { my_type res = *this; return res -= r; } my_type operator-() const noexcept { my_type res = *this; res.x = ((res.x == 0) ? 0 : (MOD - res.x)); return res; } my_type& operator*=(const my_type& r){ x = (u64)x * r.x % MOD; return *this; } my_type operator*(const my_type& r) const { my_type res = *this; return res *= r; } bool operator==(const my_type& r) const { return x == r.x; } my_type pow(unsigned long long i) const { my_type a = *this, res = 1; while(i){ if(i & 1){ res *= a; } a *= a; i >>= 1; } return res; } my_type inv() const { return my_type(ExtGcd(x, MOD).first); } unsigned int val() const { return x; } int hval() const { return int(x > MOD/2 ? x - MOD : x); } static constexpr unsigned int mod() { return MOD; } static my_type raw(unsigned int val) { auto res = my_type(); res.x = val; return res; } my_type& operator/=(const my_type& r){ return operator*=(r.inv()); } my_type operator/(const my_type& r) const { return operator*(r.inv()); } }; } // namespace nachia using Mint = nachia::StaticModint<998244353>; namespace nachia { template struct FenwickTree { public: FenwickTree() : _n(0){} explicit FenwickTree(int n, T ZERO = T(0)) : _n(n), a(n, ZERO){} void add(int p, T x){ assert(0 <= p && p < _n); p++; while(p <= _n){ a[p-1] += T(x); p += p & -p; } } T sum(int r){ return sumr(r); } T sum(int l, int r){ return sumr(r) - sumr(l); } private: int _n; std::vector a; T sumr(int r){ T s = 0; while(r > 0){ s += a[r-1]; r -= r & -r; } return s; } }; } // namespace nachia namespace nachia { template using PointAddRangeSumFwt = FenwickTree; } // namespace nachia void testcase(){ ll N, X; cin >> N >> X; V> st; V A(N); REP(i,N) cin >> A[i]; A.insert(A.begin(), X); N++; ll lim = N; auto getat = [&](ll a) -> ll { ll ok = -1, ng = st.size(); while(ok + 1 < ng){ ll w = ok + (ng - ok) / 2; if(a % st[w].first == 0) ok = w; else ng = w; } if(ng == st.size()) return lim; return st[ng].second; }; nachia::PointAddRangeSumFwt ds(N+1); // for(auto a : A) cout << a << " "; cout << endl; ds.add(N, 1); for(ll i=N-1; i>=0; i--){ ll p = getat(A[i]); ds.add(i, ds.sum(0, p+1)); // cout << "i = " << i << " , p = " << p << endl; if(A[i] != 1){ st.insert(st.begin(), { 1, i }); for(ll j=(ll)st.size()-1; j>=0; j--){ if(INF / st[j].first < A[i]){ lim = st[j].second; st.pop_back(); } else { st[j].first *= A[i]; } } } } cout << ds.sum(0, 1).val() << "\n"; } int main(){ cin.tie(0)->sync_with_stdio(0); testcase(); return 0; }