結果
問題 | No.2211 Frequency Table of GCD |
ユーザー |
![]() |
提出日時 | 2023-02-10 22:15:06 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,078 ms / 2,000 ms |
コード長 | 1,481 bytes |
コンパイル時間 | 4,378 ms |
コンパイル使用メモリ | 255,508 KB |
最終ジャッジ日時 | 2025-02-10 12:54:16 |
ジャッジサーバーID (参考情報) |
judge4 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/all>using namespace std;using namespace atcoder;//using mint = modint1000000007;//const int mod = 1000000007;using mint = modint998244353;const int mod = 998244353;//const int INF = 1e9;//const long long LINF = 1e18;#define rep(i, n) for (int i = 0; i < (n); ++i)#define rep2(i,l,r)for(int i=(l);i<(r);++i)#define rrep(i, n) for (int i = (n-1); i >= 0; --i)#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)#define all(x) (x).begin(),(x).end()#define allR(x) (x).rbegin(),(x).rend()#define endl "\n"#define P pair<int,int>template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }//約数列挙vector<long long> DivisorEnumeration(long long n) {vector<long long> x;for (long long i = 1; i*i <= n; i++) {if (0 == n % i) {x.push_back(i);if (n != i * i) x.push_back(n / i);}}sort(all(x));return x;}int main() {ios::sync_with_stdio(false);cin.tie(nullptr);int n, m; cin >> n >> m;vector<int>c(m + 1);rep(i, n) {int a; cin >> a;auto v = DivisorEnumeration(a);rep(j, v.size())c[v[j]]++;}vector<mint>ans(m + 1);rrep2(i, 1, m + 1) {ans[i] = pow_mod(2, c[i], mod) - 1;for (int j = i * 2; j <= m; j += i) ans[i] -= ans[j];}rep2(i, 1, m + 1)cout << ans[i].val() << endl;return 0;}