結果
問題 |
No.1968 Distance
|
ユーザー |
|
提出日時 | 2022-06-03 16:47:42 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 34 ms / 2,000 ms |
コード長 | 725 bytes |
コンパイル時間 | 93 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 10,752 KB |
最終ジャッジ日時 | 2024-09-21 02:20:16 |
合計ジャッジ時間 | 1,905 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 25 |
ソースコード
MOD = 998244353 N = int(input()) ans = 0 ln = [0] * 61 dv = [0] * 61 for i in range(61): ln[i] = pow((1 << i) + 1, MOD - 2, MOD) dv[i] = ln[i] * ln[i] * ((MOD + 1) // 2) % MOD for k in range(100): if N <= (1 << k) - 1: ans += dv[k] ans %= MOD break l, r = 0, k p, q = 1, 1 now = N while 0 < now: if l <= r: now -= q q *= 2 r -= 1 else: now -= p p *= 2 l -= 1 if l <= r: ans += (dv[k] - dv[k + 1]) * pow(p, MOD - 2, MOD) ans %= MOD else: ans += (ln[k] - ln[k + 1] - (dv[k] - dv[k + 1])) * pow(q, MOD - 2, MOD) ans %= MOD print(ans * 2 % MOD)