結果
| 問題 |
No.2512 Mountain Sequences
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-09-07 00:05:05 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 1,360 ms / 3,000 ms |
| コード長 | 1,707 bytes |
| コンパイル時間 | 390 ms |
| コンパイル使用メモリ | 82,432 KB |
| 実行使用メモリ | 98,520 KB |
| 最終ジャッジ日時 | 2024-06-24 18:26:01 |
| 合計ジャッジ時間 | 25,574 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 29 |
ソースコード
import io,os,sys
#https://recruit.gmo.jp/engineer/jisedai/blog/python_standard_io_with_bytesio/
buf = io.BytesIO()
max_K = 5*10**5
MOD = 998244353
def main():
fact =[-1] * (max_K+1) #MOD下での階乗
inv =[-1] * (max_K+1) #MOD下での逆元
fact_inv=[-1] * (max_K+1) #MOD下での階乗の逆元
fact[0]=fact[1]=1
inv[1] = 1
fact_inv[0]=fact_inv[1] = 1
for i in range(2,max_K+1):
fact[i] = (fact[i-1]*i)%MOD
inv[i] = (-inv[MOD%i]*(MOD//i))%MOD
fact_inv[i] = (fact_inv[i-1]*inv[i])%MOD
def c(n,r):
if n < 0 or r < 0 or r > n:
return 0
return (fact[n]*fact_inv[r])%MOD*fact_inv[n-r]%MOD
input = sys.stdin.readline
T = int(input())
block = 1000
Qli = [[] for i in range(200001//block+1)]
for i in range(T):
l,r = map(int,input().split())
Qli[l//block].append(r*(200001)*(200001)+l*(200001)+i)
l = 0
r = 0
ans = 0
answers = [0]*T
f = 0
f2 = fact_inv[2]
for li in Qli:
if len(li)==0:
continue
for x in sorted(li,reverse = f):
R = x//((200001)*(200001))
L = x//(200001)-R*(200001)
i = x%(200001)
while l > L:
l -= 1
ans =(c(2*r, l+1) - 2*ans)%MOD
while l < L:
ans = ((c(2*r, l+1) - ans) * fact_inv[2])%MOD
l += 1
while r > R:
ans =( ans - c(2*r-2, l-1))%MOD
r -= 1
while r < R:
r += 1
ans = (ans + c(2*r-2, l-1))%MOD
answers[i] = ans
f^=1
for i in range(T):
print(answers[i])
main()