結果
問題 |
No.2211 Frequency Table of GCD
|
ユーザー |
![]() |
提出日時 | 2023-02-10 22:28:18 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,284 ms / 2,000 ms |
コード長 | 1,056 bytes |
コンパイル時間 | 436 ms |
コンパイル使用メモリ | 82,312 KB |
実行使用メモリ | 140,572 KB |
最終ジャッジ日時 | 2024-07-07 18:05:59 |
合計ジャッジ時間 | 12,139 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
def div(n): if n <= 0: return [] S = set() i = 1 while i * i <= n: if n % i == 0: S.add(i) S.add(n // i) i += 1 return list(S) N, M = map(int, input().split()) A = list(map(int, input().split())) mod = 998244353 L = [0] * (M + 1) D = [0] * (M + 1) for a in A: D[a] += 1 dic = dict() for i in range(1, M + 1): if D[i] == 0: continue dic[i] = div(i) for d in dic[i]: L[d] += D[i] two = [1] * 202020 for i in range(202000): two[i + 1] = two[i] * 2 two[i + 1] %= mod ans = [0] * (M + 1) for i in range(M + 1): ans[i] = two[L[i]] - 1 for i in range(M, 0, -1): if ans[i] == 0: continue if i in dic: for d in dic[i]: if d == i: continue ans[d] -= ans[i] ans[d] %= mod else: for d in div(i): if d == i: continue ans[d] -= ans[i] ans[d] %= mod for i in range(1, M + 1): print(ans[i])