#pragma GCC target ("avx2") #pragma GCC optimize ("O3") #pragma GCC optimize("Ofast") #pragma GCC optimize ("unroll-loops") // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native") #include #include using namespace __gnu_pbds; using namespace std; #define io ios_base::sync_with_stdio(false);cin.tie(NULL);cout.tie(NULL) #define endl '\n' typedef long long ll; #define mod1 (ll)998244353 #define int ll #define pll pair typedef long double lb; typedef tree< pair, null_type, less>, rb_tree_tag, tree_order_statistics_node_update> ordered_set; #define eps (lb)(1e-9) struct custom_hash { static uint64_t splitmix64(uint64_t x) { x += 0x9e3779b97f4a7c15; x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9; x = (x ^ (x >> 27)) * 0x94d049bb133111eb; return x ^ (x >> 31); } size_t operator()(uint64_t x) const { static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count(); return splitmix64(x + FIXED_RANDOM); } }; // Operator overloads template // cin >> pair istream& operator>>(istream &istream, pair &p) { return (istream >> p.first >> p.second); } template // cin >> vector istream& operator>>(istream &istream, vector &v) { for (auto &it : v) cin >> it; return istream; } template // cout << pair ostream& operator<<(ostream &ostream, const pair &p) { return (ostream << p.first << " " << p.second); } template // cout << vector ostream& operator<<(ostream &ostream, const vector &c) { for (auto &it : c) cout << it << " "; return ostream; } // Utility functions template void print(T &&t) { cout << t << "\n"; } template void print(T &&t, Args &&... args) { cout << t << " "; print(forward(args)...); } mt19937 rng(chrono::steady_clock::now().time_since_epoch().count()); ll random(ll p){ // gives random number in [0,p] return uniform_int_distribution(1<<19, p)(rng); } ll power(ll a, ll b, ll c=mod1){ if(a==0){ return 0;} else if(b==0||a==1){ return 1;} ll p=power(a,b/2,c); p*=p; p%=c; if(b%2){p*=a;} p%=c; return p;} ll modinv(ll a, ll c=mod1){ return power(a,c-2,c); } ll n,m; vectorfreq,dp,cnt; void solve(); signed main() { io; ll t=1,n=1; // cin>>t; while (t--){ solve(); } return 0; } void solve(){ cin>>n>>m; freq.resize(m+1); freq=vector(m+1,0); dp=cnt=freq; vectorv(n); cin>>v; for(auto i:v){freq[i]++;} for(ll i(1);i<=m;++i){ for(ll j(i);j<=m;j+=i){ dp[i]+=freq[j]; } } // cout<0;--i){ cnt[i]=power(2,dp[i])-1; for(ll j(2*i);j<=m;j+=i){ cnt[i]-=cnt[j]; } cnt[i]%=mod1; cnt[i]+=mod1; cnt[i]%=mod1; } for(ll i(1);i<=m;++i){ cout<