結果
問題 | 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 = 998244353N = int(input())ans = 0ln = [0] * 61dv = [0] * 61for i in range(61):ln[i] = pow((1 << i) + 1, MOD - 2, MOD)dv[i] = ln[i] * ln[i] * ((MOD + 1) // 2) % MODfor k in range(100):if N <= (1 << k) - 1:ans += dv[k]ans %= MODbreakl, r = 0, kp, q = 1, 1now = Nwhile 0 < now:if l <= r:now -= qq *= 2r -= 1else:now -= pp *= 2l -= 1if l <= r:ans += (dv[k] - dv[k + 1]) * pow(p, MOD - 2, MOD)ans %= MODelse:ans += (ln[k] - ln[k + 1] - (dv[k] - dv[k + 1])) * pow(q, MOD - 2, MOD)ans %= MODprint(ans * 2 % MOD)