結果
問題 |
No.833 かっこいい電車
|
ユーザー |
![]() |
提出日時 | 2025-06-05 05:14:27 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 408 ms / 2,000 ms |
コード長 | 5,632 bytes |
コンパイル時間 | 666 ms |
コンパイル使用メモリ | 82,680 KB |
実行使用メモリ | 92,124 KB |
最終ジャッジ日時 | 2025-06-05 05:14:39 |
合計ジャッジ時間 | 10,729 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 30 |
ソースコード
import typing import sys input = sys.stdin.readline def _ceil_pow2(n: int) -> int: x = 0 while (1 << x) < n: x += 1 return x class SegTree: def __init__(self, op: typing.Callable[[typing.Any, typing.Any], typing.Any], e: typing.Any, v: typing.Union[int, typing.List[typing.Any]]) -> None: self._op = op self._e = e if isinstance(v, int): v = [e] * v self._n = len(v) self._log = _ceil_pow2(self._n) self._size = 1 << self._log self._d = [e] * (2 * self._size) for i in range(self._n): self._d[self._size + i] = v[i] for i in range(self._size - 1, 0, -1): self._update(i) def set(self, p: int, x: typing.Any) -> None: assert 0 <= p < self._n p += self._size self._d[p] = x for i in range(1, self._log + 1): self._update(p >> i) def get(self, p: int) -> typing.Any: assert 0 <= p < self._n return self._d[p + self._size] def prod(self, left: int, right: int) -> typing.Any: assert 0 <= left <= right <= self._n sml = self._e smr = self._e left += self._size right += self._size while left < right: if left & 1: sml = self._op(sml, self._d[left]) left += 1 if right & 1: right -= 1 smr = self._op(self._d[right], smr) left >>= 1 right >>= 1 return self._op(sml, smr) def all_prod(self) -> typing.Any: return self._d[1] def max_right(self, left: int, f: typing.Callable[[typing.Any], bool]) -> int: assert 0 <= left <= self._n assert f(self._e) if left == self._n: return self._n left += self._size sm = self._e first = True while first or (left & -left) != left: first = False while left % 2 == 0: left >>= 1 if not f(self._op(sm, self._d[left])): while left < self._size: left *= 2 if f(self._op(sm, self._d[left])): sm = self._op(sm, self._d[left]) left += 1 return left - self._size sm = self._op(sm, self._d[left]) left += 1 return self._n def min_left(self, right: int, f: typing.Callable[[typing.Any], bool]) -> int: assert 0 <= right <= self._n assert f(self._e) if right == 0: return 0 right += self._size sm = self._e first = True while first or (right & -right) != right: first = False right -= 1 while right > 1 and right % 2: right >>= 1 if not f(self._op(self._d[right], sm)): while right < self._size: right = 2 * right + 1 if f(self._op(self._d[right], sm)): sm = self._op(self._d[right], sm) right -= 1 return right + 1 - self._size sm = self._op(self._d[right], sm) return 0 def _update(self, k: int) -> None: self._d[k] = self._op(self._d[2 * k], self._d[2 * k + 1]) class FenwickTree: def __init__(self, n: int): self.data = [0] * (n+10) self.n = (n+10) def add(self, p: int, x: int): assert 0 <= p < self.n p += 1 while p < len(self.data): self.data[p] += x p += p & -p def sum(self, p: int) -> int: """区間 [0, p] の和""" assert 0 <= p < self.n p += 1 s = 0 while p > 0: s += self.data[p] p -= p & -p return s def rangesum(self, l: int, r: int) -> int: """区間 [l, r] の和""" assert 0 <= l <= r < self.n s = self.sum(r) if l > 0: s -= self.sum(l-1) return s N, Q = map(int, input().split()) A = list(map(int, input().split())) ft = FenwickTree(N) for i in range(N): ft.add(i, 1) segt = SegTree(lambda a, b: a+b, 0, A) def get_low(p: int): lo = 0 hi = N-1 res = hi while lo <= hi: m = (lo + hi) // 2 if ft.sum(m) >= p+1: res = min(res, m) hi = m - 1 else: lo = m + 1 return res def get_high(p: int): lo = 0 hi = N-1 res = 0 while lo <= hi: m = (lo + hi) // 2 if ft.sum(m) <= p+1: res = max(res, m) lo = m + 1 else: hi = m - 1 return res def dump(): xs = [ft.rangesum(i, i) for i in range(N)] print(f'{xs=}') ys = [segt.get(i) for i in range(N)] print(f'{ys=}') for _ in range(Q): q, x = map(int, input().split()) #print() #print(f'{q=} {x=}') #dump() match q: case 1: # 連結 x -= 1 if ft.rangesum(x+1, x+1) == 1: ft.add(x+1, -1) case 2: # 切り離す x -= 1 if ft.rangesum(x+1, x+1) == 0: ft.add(x+1, 1) case 3: # inc x -= 1 segt.set(x, segt.get(x) + 1) case 4: # 区間和 x -= 1 p = ft.sum(x)-1 l = get_low(p) r = get_high(p) res = segt.prod(l, r+1) #print(f'{x=} {l=} {r=} {res=}') print(res)