結果
問題 | 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;}