結果
問題 |
No.2211 Frequency Table of GCD
|
ユーザー |
![]() |
提出日時 | 2024-07-07 12:05:09 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 634 ms / 2,000 ms |
コード長 | 1,011 bytes |
コンパイル時間 | 4,141 ms |
コンパイル使用メモリ | 252,272 KB |
最終ジャッジ日時 | 2025-02-22 02:30:57 |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; #define rep(i, n) for(int i=0;i<(n);++i) #define rep1(i, n) for(int i=1;i<=(n);i++) #define ll long long using mint = modint998244353; using P = pair<ll,ll>; using lb = long double; using T = tuple<ll, ll, ll>; #ifdef LOCAL # include <debug_print.hpp> # define dbg(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__) #else # define dbg(...) (static_cast<void>(0)) #endif int main() { int n, m; cin >> n >> m; vector<int> a(n); rep(i,n) cin >> a[i]; vector<mint> ans(m+1); vector<int> cnt(m+1); rep(i,n) { for(int j=1;j*j<=a[i];j++){ if(a[i]%j==0) { cnt[j]++; if(a[i]/j!=j) cnt[a[i]/j]++; } } } for(int k=m;k>=1;k--){ ans[k] = mint(2).pow(cnt[k])-1; for(int j=k+k;j<=m;j+=k) ans[k] -= ans[j]; } for(int i=1;i<=m;i++){ cout << ans[i].val() << endl; } return 0; }