結果

問題 No.1968 Distance
ユーザー shiomusubi496
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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)
0