結果
問題 | No.2211 Frequency Table of GCD |
ユーザー |
|
提出日時 | 2023-03-17 17:49:01 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,050 ms / 2,000 ms |
コード長 | 1,469 bytes |
コンパイル時間 | 1,101 ms |
コンパイル使用メモリ | 126,956 KB |
最終ジャッジ日時 | 2025-02-11 12:26:54 |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <iostream>#include <string>#include <vector>#include <algorithm>#include <functional>#include <cmath>#include <iomanip>#include <stack>#include <queue>#include <numeric>#include <map>#include <unordered_map>#include <set>#include <fstream>#include <chrono>#include <random>#include <bitset>//#include <atcoder/all>#define rep(i,n) for(int i=0;i<(n);i++)#define all(x) x.begin(), x.end()#define rall(x) x.rbegin(), x.rend()#define sz(x) ((int)(x).size())#define pb push_backusing ll = long long;using namespace std;template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }ll gcd(ll a, ll b) {return b?gcd(b,a%b):a;}ll lcm(ll a, ll b) {return a/gcd(a,b)*b;}const ll mod = 998244353;ll mpow(ll a, ll x){ll res = 1;while(x > 0){if(x & 1) res = (res * a) % mod;a = (a * a) % mod;x >>= 1;}return res;}int main(){ll N,M; cin >> N >> M;vector<ll> A(N);rep(i,N) cin >> A[i];vector<ll> ans(M+1,0);vector<ll> cnt(202000,0);rep(i,N){for(ll a=1;a*a<=A[i];a++){if(A[i]%a==0){cnt[a]++;if(a*a!=A[i]) cnt[A[i]/a]++;}}}for(int m=M;m>0;m--){ll a = (mpow(2,cnt[m])+mod-1);for(int mm=m*2;mm<=M;mm+=m){a += mod-ans[mm];}ans[m] = a%mod;}for(int i=1;i<=M;i++){cout << ans[i] << '\n';}return 0;}