結果
| 問題 |
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)