結果
問題 | No.2075 GCD Subsequence |
ユーザー |
|
提出日時 | 2022-09-24 03:03:45 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 285 ms / 4,000 ms |
コード長 | 1,200 bytes |
コンパイル時間 | 1,963 ms |
コンパイル使用メモリ | 197,332 KB |
最終ジャッジ日時 | 2025-02-07 14:10:16 |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
#include <bits/stdc++.h> using namespace std; const int mod=998244353; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int N; cin>>N; vector<int>A(N); for(int &i:A)cin>>i; vector<unsigned>cnt(1000001); unsigned ans=0; for(int i=0;i<N;i++){ if(A[i]==1){ ++ans; continue; } int p=A[i]; vector<int>d; for(int j=2;j*j<=p;j++){ if(p%j==0){ d.push_back(j); while(p%j==0)p/=j; } } if(p!=1)d.push_back(p); int s=1; for(int j=1;j<1<<d.size();j++){ int k=1; for(int l=0;l<d.size();l++){ if(j>>l&1)k*=d[l]; } if(__popcount(j)%2==1)s=min(s+cnt[k],s+cnt[k]-mod); else s=min(s+mod-cnt[k],s-cnt[k]); } for(int j=1;j<1<<d.size();j++){ int k=1; for(int l=0;l<d.size();l++){ if(j>>l&1)k*=d[l]; } cnt[k]=min(s+cnt[k],s+cnt[k]-mod); } } vector<int>hurui(1000001); for(int i=2;i<=1000000;i++){ if(hurui[i]!=0)continue; for(int j=i;j<=1000000;j+=i)++hurui[j]; } for(int i=2;i<=1000000;i++){ if(hurui[i]%2)ans=min(ans+cnt[i],ans+cnt[i]-mod); else ans=min(ans+mod-cnt[i],ans-cnt[i]); } cout<<ans<<'\n'; }