結果
問題 | No.2075 GCD Subsequence |
ユーザー |
![]() |
提出日時 | 2022-09-04 08:30:34 |
言語 | C++17 (gcc 12.2.0 + boost 1.81.0) |
結果 |
AC
|
実行時間 | 712 ms / 4,000 ms |
コード長 | 1,046 bytes |
コンパイル時間 | 1,002 ms |
コンパイル使用メモリ | 75,188 KB |
実行使用メモリ | 67,828 KB |
最終ジャッジ日時 | 2023-08-11 16:40:14 |
合計ジャッジ時間 | 22,236 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 437 ms
67,472 KB |
testcase_01 | AC | 440 ms
67,828 KB |
testcase_02 | AC | 439 ms
67,560 KB |
testcase_03 | AC | 434 ms
67,528 KB |
testcase_04 | AC | 433 ms
67,668 KB |
testcase_05 | AC | 437 ms
67,564 KB |
testcase_06 | AC | 438 ms
67,600 KB |
testcase_07 | AC | 435 ms
67,472 KB |
testcase_08 | AC | 517 ms
67,504 KB |
testcase_09 | AC | 570 ms
67,512 KB |
testcase_10 | AC | 522 ms
67,668 KB |
testcase_11 | AC | 549 ms
67,564 KB |
testcase_12 | AC | 526 ms
67,588 KB |
testcase_13 | AC | 507 ms
67,548 KB |
testcase_14 | AC | 543 ms
67,572 KB |
testcase_15 | AC | 505 ms
67,536 KB |
testcase_16 | AC | 511 ms
67,540 KB |
testcase_17 | AC | 563 ms
67,788 KB |
testcase_18 | AC | 700 ms
67,568 KB |
testcase_19 | AC | 703 ms
67,652 KB |
testcase_20 | AC | 709 ms
67,564 KB |
testcase_21 | AC | 704 ms
67,516 KB |
testcase_22 | AC | 700 ms
67,496 KB |
testcase_23 | AC | 709 ms
67,488 KB |
testcase_24 | AC | 712 ms
67,564 KB |
testcase_25 | AC | 706 ms
67,548 KB |
testcase_26 | AC | 705 ms
67,472 KB |
testcase_27 | AC | 705 ms
67,496 KB |
testcase_28 | AC | 695 ms
67,504 KB |
testcase_29 | AC | 432 ms
67,556 KB |
testcase_30 | AC | 429 ms
67,704 KB |
ソースコード
#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; }