結果
問題 |
No.2571 列辞書順列
|
ユーザー |
![]() |
提出日時 | 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)