結果
問題 | No.2075 GCD Subsequence |
ユーザー |
![]() |
提出日時 | 2022-09-04 08:30:34 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 995 ms / 4,000 ms |
コード長 | 1,046 bytes |
コンパイル時間 | 812 ms |
コンパイル使用メモリ | 73,488 KB |
最終ジャッジ日時 | 2025-02-07 02:31:46 |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
#include<cassert> #include<iostream> #include<vector> using namespace std; typedef long long lint; typedef vector<int>vi; #define rep(i,n)for(int i=0;i<(int)(n);++i) const lint mod=998244353; #define N 1000001 vi fac[N]; bool pr[N]; void pr_init(void){ rep(i,N)if(i>=2)pr[i]=1; rep(i,N){ if(i<=1)continue; if(!pr[i])continue; fac[i].push_back(i); for(lint j=lint(i)*2;j<N;j+=i){ pr[j]=0; fac[j].push_back(i); } } } void add(lint&a,lint b){ a+=b; if(a>=mod)a-=mod; } lint dp[N]; int main(){ pr_init(); int n;cin>>n; for(int i=2;i<N;i++)dp[i]=1; lint ans=0; rep(i,n){ int a;cin>>a; if(a==1){ add(ans,1); continue; } int b=fac[a].size(); lint val=0; rep(c,1<<b){ if(c==0)continue; int d=1,e=0; rep(f,b)if(c&1<<f)d*=fac[a][f],e++; if(e%2)add(val,dp[d]); else add(val,mod-dp[d]); } rep(c,1<<b){ int d=1; rep(f,b)if(c&1<<f)d*=fac[a][f]; add(dp[d],val); } add(ans,val); } cout<<ans<<endl; }