結果

問題 No.2211 Frequency Table of GCD
ユーザー HydruHydru
提出日時 2023-02-10 22:03:15
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,848 ms / 2,000 ms
コード長 881 bytes
コンパイル時間 532 ms
コンパイル使用メモリ 82,296 KB
実行使用メモリ 122,132 KB
最終ジャッジ日時 2024-07-07 16:18:28
合計ジャッジ時間 22,562 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
52,340 KB
testcase_01 AC 35 ms
52,364 KB
testcase_02 AC 33 ms
52,688 KB
testcase_03 AC 124 ms
89,424 KB
testcase_04 AC 700 ms
97,152 KB
testcase_05 AC 1,074 ms
111,032 KB
testcase_06 AC 778 ms
91,764 KB
testcase_07 AC 1,105 ms
107,644 KB
testcase_08 AC 67 ms
82,940 KB
testcase_09 AC 60 ms
76,864 KB
testcase_10 AC 97 ms
104,780 KB
testcase_11 AC 78 ms
94,144 KB
testcase_12 AC 114 ms
111,420 KB
testcase_13 AC 611 ms
92,664 KB
testcase_14 AC 681 ms
90,108 KB
testcase_15 AC 463 ms
86,984 KB
testcase_16 AC 618 ms
89,804 KB
testcase_17 AC 967 ms
107,200 KB
testcase_18 AC 1,585 ms
121,744 KB
testcase_19 AC 1,848 ms
121,928 KB
testcase_20 AC 1,603 ms
121,736 KB
testcase_21 AC 1,573 ms
122,064 KB
testcase_22 AC 1,571 ms
121,792 KB
testcase_23 AC 81 ms
94,304 KB
testcase_24 AC 1,293 ms
122,132 KB
testcase_25 AC 80 ms
112,648 KB
testcase_26 AC 34 ms
52,928 KB
testcase_27 AC 1,791 ms
122,020 KB
testcase_28 AC 1,335 ms
122,112 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def divisor_enumation(n):
    i=1
    result=set()
    while(i*i<n):
        if n%i==0:
            result.add(i)
            result.add(n//i)
        i+=1
    if i*i==n:
        result.add(i)
    return result

mod=998244353

n,m=map(int,input().split())
A=list(map(int,input().split()))

pow_2=[1]
for _ in range(n+1):
    pow_2.append(pow_2[-1]*2%mod)

cnt=[0 for _ in range(m+1)]

for i in range(n):
    a=A[i]
    n=1

    while(n*n<a):
        if a%n==0:
            cnt[n]+=1
            cnt[a//n]+=1
        n+=1
    if n*n==a:
        cnt[n]+=1

dem=[0 for _ in range(m+1)]
ans=[]

for i in range(m,0,-1):
    
    ans.append((pow_2[cnt[i]]-1-dem[i])%mod)

    if cnt[i]==0:
        continue
    
    X=divisor_enumation(i)
    for u in X:
        if u==i:
            continue
        dem[u]+=pow_2[cnt[i]]-1-dem[i]
        dem[u]%=mod

for x in ans[::-1]:
    print(x)
0