結果
問題 | No.2211 Frequency Table of GCD |
ユーザー |
![]() |
提出日時 | 2023-02-10 23:54:08 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 501 ms / 2,000 ms |
コード長 | 789 bytes |
コンパイル時間 | 186 ms |
コンパイル使用メモリ | 82,176 KB |
実行使用メモリ | 138,044 KB |
最終ジャッジ日時 | 2024-07-07 18:11:44 |
合計ジャッジ時間 | 9,725 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#模範解答def prime_factors(n):#2~nまでの数を素因数分解して、各々の素因数リストを返します。O(nlog(n))d=[[] for _ in range(n+1)]for i in range(1,n+1):for j in range(i,n+1,i):d[j].append(i)return(d)N,M=map(int,input().split())mod=998244353A=list(map(int,input().split()))ans=[0]*(M+1)#ans[i]:最大公約数がians_ex=[0]*(M+1)#ans_ex[i]:最大公約数がiの倍数num=[0]*(M+1)d=prime_factors(M)for i in range(N):for n in d[A[i]]:num[n]+=1for i in range(1,M+1):ans_ex[i]=pow(2,num[i],mod)-1for i in range(M,0,-1):ima=ans_ex[i]%modima_ind=2*iwhile ima_ind<=M:ima-=ans[ima_ind]ima_ind+=iima%=modans[i]=imafor i in range(1,M+1):print(ans[i])