結果
| 問題 | No.2571 列辞書順列 |
| コンテスト | |
| ユーザー |
AngrySadEight
|
| 提出日時 | 2023-12-01 08:05:42 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 612 ms / 3,000 ms |
| コード長 | 1,594 bytes |
| 記録 | |
| コンパイル時間 | 419 ms |
| コンパイル使用メモリ | 81,920 KB |
| 実行使用メモリ | 97,664 KB |
| 最終ジャッジ日時 | 2024-09-27 03:40:44 |
| 合計ジャッジ時間 | 17,541 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 29 |
ソースコード
mod = 998244353
class BIT:
def __init__(self, n):
self.n = n
self.lst = [0 for _ in range(n + 1)]
def add(self, p, x):
while True:
if p > self.n:
break
self.lst[p] += x
self.lst[p] %= 998244353
p += (p & -p)
def sum(self, p):
res = 0
while True:
if p <= 0:
break
res += self.lst[p]
res %= 998244353
p -= (p & -p)
return res
def get(self, p):
return (self.sum(p) - self.sum(p - 1)) % 998244353
N, M, K, Q = map(int, input().split())
A = list(map(int, input().split()))
for i in range(N):
A[i] -= 1
inv_m = pow(M - 1, mod - 2, mod)
pow_m = [1]
for i in range(N):
add_res = (pow_m[-1] * M) % mod
pow_m.append(add_res)
bit1 = BIT(N)
bit2 = BIT(N)
for i in range(N):
bit1.add(i + 1, (pow(M, K - i, mod) * A[i]) % mod)
bit2.add(i + 1, A[i])
for i in range(Q):
t, l, r = map(int, input().split())
if t == 1:
# do nothing
l -= 1
r -= 1
ans = (r - l + 1)
prd1 = bit1.sum(r + 1) - bit1.sum(l)
prd1 = (prd1 * pow_m[l]) % mod
prd2 = bit2.sum(r + 1) - bit2.sum(l)
sums = (prd1 - prd2) % mod
sums = (sums * inv_m) % mod
ans = (ans + sums) % mod
print(ans)
else:
l -= 1
r -= 1
now_bit1 = bit1.get(l + 1)
now_bit2 = bit2.get(l + 1)
bit1.add(l + 1, (-now_bit1 + pow(M, K - l, mod) * r) % mod)
bit2.add(l + 1, (-now_bit2 + r) % mod)
AngrySadEight