結果
問題 | No.2902 ZERO!! |
ユーザー |
|
提出日時 | 2024-09-27 20:58:26 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 477 ms / 2,000 ms |
コード長 | 717 bytes |
コンパイル時間 | 233 ms |
コンパイル使用メモリ | 82,388 KB |
実行使用メモリ | 100,348 KB |
最終ジャッジ日時 | 2024-09-27 20:58:36 |
合計ジャッジ時間 | 9,379 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 41 |
ソースコード
from collections import defaultdictN = int(input())lpf = [-1] * (N + 1)prime = []for i in range(2, N + 1):if lpf[i] == -1:lpf[i] = iprime.append(i)for p in prime:if p * i > N or p > lpf[i]:breaklpf[p * i] = pcnt = []for p in prime:k = pc = 0while k <= N:c += N // kk *= pcnt.append(c)mod = 998244353m = cnt[0]ac = [1] * (m + 2)# 後ろから累積積for i in cnt:for j in range(1, i + 1):if i // j != i // (j + 1):ac[j] *= (i // j + 1) * pow(i // (j + 1) + 1, mod - 2, mod) % modac[j] %= modfor i in range(m + 1, 0, -1):ac[i - 1] *= ac[i]ac[i - 1] %= modans = 0for i in range(1, m + 1):ans += i * (ac[i] - ac[i + 1])ans %= modprint(ans)