結果
問題 |
No.2211 Frequency Table of GCD
|
ユーザー |
![]() |
提出日時 | 2023-02-10 21:41:23 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,814 ms / 2,000 ms |
コード長 | 883 bytes |
コンパイル時間 | 1,864 ms |
コンパイル使用メモリ | 196,308 KB |
最終ジャッジ日時 | 2025-02-10 12:24:53 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define int long long int int_pow(int x, int n, int MOD) { int ret = 1; while (n > 0) { if (n & 1) ret = ret * x % MOD; x = x * x % MOD; n >>= 1; } return ret; } signed main(){ int N,M; cin>>N>>M; vector<int> A(M+1); for(int i=0;i<N;i++){ int x; cin>>x; for(int j=1;j*j<=x;j++){ if(x%j == 0){ A[j]++; if(j*j != x) A[x/j]++; } } } vector<int> S(M+1); int mod = 998244353; for(int i=M;i>0;i--){ S[i] += int_pow(2,A[i],mod)-1; S[i] += 100*mod; S[i] %= mod; for(int j=1;j*j<=i;j++){ if(i%j == 0){ if(j == i) continue; S[j] -= S[i]; if(j*j != i && j != 1) S[i/j] -= S[i]; } } } for(int i=1;i<=M;i++) cout<<(S[i]+1000000*mod)%mod<<endl; }