結果
問題 |
No.2211 Frequency Table of GCD
|
ユーザー |
![]() |
提出日時 | 2023-02-11 02:29:16 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 302 ms / 2,000 ms |
コード長 | 1,095 bytes |
コンパイル時間 | 1,118 ms |
コンパイル使用メモリ | 107,508 KB |
最終ジャッジ日時 | 2025-02-10 13:58:33 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <iostream> #include <vector> #include <cmath> #include <map> #include <set> #include <iomanip> #include <queue> #include <algorithm> #include <numeric> #include <deque> #include <complex> #include <cassert> using namespace std; const long long modc = 998244353; long long mod_exp(long long b, long long e){ const long long m=modc; if (e > 0 && b == 0) return 0; long long ans = 1; b %= m; while (e > 0){ if ((e & 1LL)) ans = (ans * b) % m; e = e >> 1LL; b = (b*b) % m; } return ans; } int main(){ long long N, M; cin >> N >> M; vector<long long> A(N), c(M+1), cnt(M+1), ans(M+1); for (int i=0; i<N; i++){ cin >> A[i]; c[A[i]]++; } for (int i=1; i<=M; i++){ for (int j=i; j<=M; j+=i) cnt[i] += c[j]; } for (int i=M; i>=1; i--){ ans[i] = (modc+mod_exp(2, cnt[i])-1) % modc; for (int j=2*i; j<=M; j+=i){ ans[i] += modc - ans[j]; ans[i] %= modc; } } for (int i=1; i<=M; i++) cout << ans[i] << endl; return 0; }