結果
問題 | No.2211 Frequency Table of GCD |
ユーザー | tktk_snsn |
提出日時 | 2023-02-10 21:51:07 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 202 ms / 2,000 ms |
コード長 | 926 bytes |
コンパイル時間 | 275 ms |
コンパイル使用メモリ | 87,124 KB |
実行使用メモリ | 111,104 KB |
最終ジャッジ日時 | 2023-09-22 01:05:15 |
合計ジャッジ時間 | 6,794 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 85 ms
77,388 KB |
testcase_01 | AC | 78 ms
77,484 KB |
testcase_02 | AC | 79 ms
77,384 KB |
testcase_03 | AC | 155 ms
86,532 KB |
testcase_04 | AC | 160 ms
93,096 KB |
testcase_05 | AC | 178 ms
103,208 KB |
testcase_06 | AC | 164 ms
89,004 KB |
testcase_07 | AC | 186 ms
100,544 KB |
testcase_08 | AC | 95 ms
87,908 KB |
testcase_09 | AC | 91 ms
83,328 KB |
testcase_10 | AC | 109 ms
100,276 KB |
testcase_11 | AC | 102 ms
92,904 KB |
testcase_12 | AC | 112 ms
103,120 KB |
testcase_13 | AC | 149 ms
90,164 KB |
testcase_14 | AC | 166 ms
88,616 KB |
testcase_15 | AC | 156 ms
87,356 KB |
testcase_16 | AC | 160 ms
88,820 KB |
testcase_17 | AC | 171 ms
100,376 KB |
testcase_18 | AC | 198 ms
110,652 KB |
testcase_19 | AC | 201 ms
110,672 KB |
testcase_20 | AC | 198 ms
110,596 KB |
testcase_21 | AC | 202 ms
110,652 KB |
testcase_22 | AC | 201 ms
110,580 KB |
testcase_23 | AC | 156 ms
87,000 KB |
testcase_24 | AC | 192 ms
110,844 KB |
testcase_25 | AC | 110 ms
106,184 KB |
testcase_26 | AC | 78 ms
77,604 KB |
testcase_27 | AC | 191 ms
110,712 KB |
testcase_28 | AC | 202 ms
111,104 KB |
ソースコード
mod = 998244353 U = 2 * 10 ** 5 two = [1] * (U + 1) for i in range(1, U+1): two[i] = two[i-1] * 2 % mod def zeta_div(A, primes): n = len(A) - 1 for p in primes: for i in reversed(range(1, n // p + 1)): A[i] += A[i * p] def moebius_div(A, primes): n = len(A) for p in primes: for i in range(1, n): if i * p >= n: break A[i] -= A[i * p] def main(): N, M = map(int, input().split()) A = list(map(int, input().split())) F = [0] * (M + 1) for a in A: F[a] += 1 for i in range(M + 1): F[i] = two[F[i]] - 1 for i in range(1, M + 1): p = 1 for j in range(i, M+1, i): p = p * (F[j] + 1) % mod F[i] = p - 1 for i in range(1, M + 1)[::-1]: for j in range(i+i, M+1, i): F[i] -= F[j] F[i] %= mod print(*F[1:], sep="\n") main()