結果
問題 | No.2206 Popcount Sum 2 |
ユーザー |
|
提出日時 | 2023-02-04 13:45:01 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 2,293 ms / 4,000 ms |
コード長 | 1,464 bytes |
コンパイル時間 | 247 ms |
コンパイル使用メモリ | 82,380 KB |
実行使用メモリ | 138,872 KB |
最終ジャッジ日時 | 2024-07-03 11:33:12 |
合計ジャッジ時間 | 27,456 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 18 |
ソースコード
import sysinput = lambda :sys.stdin.readline()[:-1]ni = lambda :int(input())na = lambda :list(map(int,input().split()))yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")#######################################################################mod = 998244353nn = 10**6fact = [1] * nnfor i in range(nn - 1):fact[i + 1] = fact[i] * (i + 1) % modinvfact = [1] * nninvfact[nn - 1] = pow(fact[nn - 1], mod - 2, mod)for i in range(nn - 1)[::-1]:invfact[i] = invfact[i + 1] * (i + 1) % moddef binom(x, y):if x < 0 or y < 0 or x - y < 0:return 0return fact[x] * invfact[y] % mod * invfact[x - y] % modque = []t = ni()B = 500ans = [None] * tfor _ in range(t):n,m = na()que.append((n//B, m, n, _))que = sorted(que)sm, sn = 1, 1sum = 1z = pow(2, mod-2, mod)for q in que:_, m, n, i = q#print(sum)if sm < m:for x in range(sm, m):sum += binom(sn-1, x)sum %= modelse:for x in range(sm, m, -1):sum -= binom(sn-1, x-1)sum %= modif sn < n:for x in range(sn, n):sum = (sum * 2 - binom(x-1, m-1))%modelse:for x in range(sn, n, -1):sum = (sum + binom(x-2, m-1))*z %modsum %= modans[i] = sum * (pow(2, n, mod)-1)%modsn, sm = n, mprint(*ans, sep = "\n")