結果
| 問題 |
No.2060 AND Sequence
|
| コンテスト | |
| ユーザー |
tamato
|
| 提出日時 | 2022-08-26 21:53:29 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 40 ms / 2,000 ms |
| コード長 | 1,054 bytes |
| コンパイル時間 | 645 ms |
| コンパイル使用メモリ | 82,688 KB |
| 実行使用メモリ | 52,736 KB |
| 最終ジャッジ日時 | 2024-10-13 22:27:12 |
| 合計ジャッジ時間 | 3,232 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 43 |
ソースコード
mod = 998244353
def main():
import sys
input = sys.stdin.readline
# comb init
nmax = 30 + 10 # change here
fac = [0] * nmax
finv = [0] * nmax
inv = [0] * nmax
fac[0] = 1
fac[1] = 1
finv[0] = 1
finv[1] = 1
inv[1] = 1
for i in range(2, nmax):
fac[i] = fac[i - 1] * i % mod
inv[i] = mod - inv[mod % i] * (mod // i) % mod
finv[i] = finv[i - 1] * inv[i] % mod
def comb(n, r):
if n < r:
return 0
else:
return (fac[n] * ((finv[r] * finv[n - r]) % mod)) % mod
N, M = map(int, input().split())
popcnt = [0] * 31
cnt1 = 0
for i in range(M.bit_length() - 1, -1, -1):
if M >> i & 1:
K = i
for k in range(K + 1):
popcnt[k + cnt1] = (popcnt[k + cnt1] + comb(K, k)) % mod
cnt1 += 1
popcnt[cnt1] += 1
ans = 0
for i in range(31):
tmp = pow(N, i, mod)
ans = (ans + tmp * popcnt[i])%mod
print(ans)
if __name__ == '__main__':
main()
tamato