結果
| 問題 |
No.2211 Frequency Table of GCD
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-02-10 22:38:53 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 302 ms / 2,000 ms |
| コード長 | 768 bytes |
| コンパイル時間 | 3,864 ms |
| コンパイル使用メモリ | 250,320 KB |
| 最終ジャッジ日時 | 2025-02-10 13:06:30 |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 26 |
ソースコード
//借りました
//https://kmjp.hatenablog.jp/entry/2017/04/30/0930
#include<bits/stdc++.h>
using namespace std;
#include<atcoder/all>
using namespace atcoder;
using ll = long long;
int N,A;
ll mo=998244353;
int num[201010];
ll p2[201010];
ll pat[201010];
void solve() {
int i,j,k,l,r,x,y; string s;
p2[0]=1;
for(int i = 0;i<201000;i++) p2[i+1]=p2[i]*2%mo;
cin>>N >> k;
for(int i = 0;i<N;i++){
cin>>x;
num[x]++;
}
for(i=1;i<=200000;i++) {
for(j=i*2;j<=200000;j+=i) num[i]+=num[j];
pat[i]=(p2[num[i]]+mo-1)%mo;
}
for(i=200000;i>=1;i--) {
for(j=i*2;j<=200000;j+=i) pat[i]-=pat[j];
pat[i]=(pat[i]%mo+mo)%mo;
}
for(int i = 1;i<=k;i++){
cout<<pat[i]<<endl;
}
}
signed main(){
cin.tie(nullptr);
ios::sync_with_stdio(false);
solve();
}