結果
問題 |
No.2211 Frequency Table of GCD
|
ユーザー |
![]() |
提出日時 | 2023-02-10 22:46:07 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 85 ms / 2,000 ms |
コード長 | 935 bytes |
コンパイル時間 | 523 ms |
コンパイル使用メモリ | 73,680 KB |
実行使用メモリ | 6,528 KB |
最終ジャッジ日時 | 2024-07-07 18:09:31 |
合計ジャッジ時間 | 3,356 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <iostream> #include <vector> #include <algorithm> using namespace std; using ll = long long; constexpr ll MOD = 998244353; ll modpow(ll a, ll x){ ll res = 1; for(; x > 0; x >>= 1){ if(x & 1){ res *= a; res %= MOD; } a *= a; a %= MOD; } return res; } int main(){ int n, m; cin >> n >> m; vector<int> a(n); for(auto &it: a) cin >> it; const int lim = 200010; vector<int> cnt(lim, 0); for(auto &it: a) cnt[it]++; vector<ll> f(lim, 0); for(int i = lim-10; i > 0; i--){ int tot = 0; for(int j = i; j < lim; j += i) tot += cnt[j]; if(tot == 0) continue; f[i] = (modpow(2, tot)-1+MOD)%MOD; for(int j = i+i; j < lim; j += i){ f[i] -= f[j]; f[i] += MOD; f[i] %= MOD; } } for(int i = 1; i <= m; i++) cout << f[i] << '\n'; return 0; }