結果

問題 No.2211 Frequency Table of GCD
ユーザー chankei271828chankei271828
提出日時 2023-02-10 23:50:14
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 770 bytes
コンパイル時間 358 ms
コンパイル使用メモリ 87,108 KB
実行使用メモリ 112,376 KB
最終ジャッジ日時 2023-09-22 00:44:24
合計ジャッジ時間 8,185 ms
ジャッジサーバーID
(参考情報)
judge11 / judge15
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
権限があれば一括ダウンロードができます

ソースコード

diff #

#模範解答
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=998244353
A=list(map(int,input().split()))
ans=[0]*(M+1)#ans[i]:最大公約数がi
ans_ex=[0]*(M+1)#ans_ex[i]:最大公約数がiの倍数
num=[0]*(M+1)
for i in range(N):
    for n in d[A[i]]:
        num[n]+=1
for i in range(1,M+1):
    ans_ex[i]=pow(2,num[i],mod)-1
for i in range(M,0,-1):
    ima=ans_ex[i]%mod
    ima_ind=2*i
    while ima_ind<=M:
        ima-=ans[ima_ind]
        ima_ind+=i
        ima%=mod
    ans[i]=ima
for i in range(1,M+1):
    print(ans[i])
0