結果
問題 | No.2645 Sum of Divisors? |
ユーザー |
|
提出日時 | 2024-02-20 01:16:12 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 126 ms / 2,000 ms |
コード長 | 450 bytes |
コンパイル時間 | 179 ms |
コンパイル使用メモリ | 82,040 KB |
実行使用メモリ | 92,156 KB |
最終ジャッジ日時 | 2024-09-29 03:21:00 |
合計ジャッジ時間 | 3,778 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 32 |
ソースコード
import mathN = int(input())def zeta(x):return 1/xzeta_c = [0] * (10**6 + 1)for i in range(1,10**6 + 1):zeta_c[i] = zeta_c[i-1] + 1/idef zeta_cumsum(x):if x <= 10**6:return zeta_c[x]else:return math.log(x) + 0.5772156649M = int(N ** 0.5)ans = 0for i in range(1, M+1):ans += zeta_cumsum(N//i) * zeta(i)for j in range(1, M+1):ans += (zeta_cumsum(N//j) - zeta_cumsum(M))* zeta(j)print(ans)