結果
問題 | No.1234 典型RMQ |
ユーザー | convexineq |
提出日時 | 2020-09-19 00:49:25 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 490 ms / 2,000 ms |
コード長 | 7,161 bytes |
コンパイル時間 | 550 ms |
コンパイル使用メモリ | 82,560 KB |
実行使用メモリ | 95,352 KB |
最終ジャッジ日時 | 2024-11-09 02:04:13 |
合計ジャッジ時間 | 12,209 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 45 ms
53,888 KB |
testcase_01 | AC | 46 ms
53,888 KB |
testcase_02 | AC | 44 ms
53,760 KB |
testcase_03 | AC | 47 ms
53,632 KB |
testcase_04 | AC | 43 ms
53,760 KB |
testcase_05 | AC | 44 ms
53,760 KB |
testcase_06 | AC | 453 ms
89,728 KB |
testcase_07 | AC | 402 ms
79,504 KB |
testcase_08 | AC | 490 ms
95,040 KB |
testcase_09 | AC | 441 ms
85,020 KB |
testcase_10 | AC | 447 ms
92,868 KB |
testcase_11 | AC | 467 ms
89,548 KB |
testcase_12 | AC | 425 ms
83,740 KB |
testcase_13 | AC | 420 ms
79,084 KB |
testcase_14 | AC | 436 ms
83,924 KB |
testcase_15 | AC | 432 ms
82,624 KB |
testcase_16 | AC | 469 ms
92,288 KB |
testcase_17 | AC | 434 ms
84,276 KB |
testcase_18 | AC | 386 ms
79,076 KB |
testcase_19 | AC | 483 ms
94,336 KB |
testcase_20 | AC | 289 ms
94,208 KB |
testcase_21 | AC | 444 ms
89,728 KB |
testcase_22 | AC | 436 ms
95,352 KB |
testcase_23 | AC | 429 ms
94,848 KB |
testcase_24 | AC | 416 ms
95,304 KB |
testcase_25 | AC | 413 ms
95,044 KB |
testcase_26 | AC | 415 ms
95,232 KB |
testcase_27 | AC | 42 ms
53,760 KB |
testcase_28 | AC | 46 ms
53,888 KB |
testcase_29 | AC | 44 ms
53,376 KB |
ソースコード
""" 遅延セグメント木(区間演算、区間更新) data[] の要素に モノイド X をもつ lazy[] の要素に Aut(X) をもつ(ただし作用素は「左」から作用とする) アクセスは0-indexed, 内部のツリーは 1-indexed(つまりすべての和は tree[1]) 関数は半開区間 引数: op_X: モノイド演算 (max, min, __add__,ラムダ式,関数定義など) e_X: 単位元 compose: 作用素を合成させる関数[注:普通の関数合成と同じく、左作用](max, min, __add__,ラムダ式,関数定義など) mapping(f,x) = f(x) 関数適用 id_M: 恒等作用素 N: 処理する区間の長さ array: この配列で初期化 """ """ 作用素の例: f \mapsto min(f,-) compose = min funcval = min ID_M = INF = 10**18 f \mapsto max(f,-) compose = max funcval = max ID_M = 0 f \mapsto f (定数関数)(区間代入、最後の操作のみが影響する) compose = lambda f,g: (f if g == ID_M else f) funcval = lambda f,x: (x if f == ID_M else f) ID_M = None #Noneではなく、範囲外の数にすると速くなる """ class LazySegmentTree: def __init__(self, op_X, e_X, mapping, compose, id_M, N, array=None): __slots__ = ["op_X","e_X","mapping","compose","id_M","N","log","N0","data","lazy"] # それぞれ Xの演算、単位元、f(x), f\circ g, Xの恒等変換 self.e_X = e_X; self.op_X = op_X; self.mapping = mapping; self.compose = compose; self.id_M = id_M self.N = N self.log = (N-1).bit_length() self.N0 = 1<<self.log self.data = [e_X]*(2*self.N0) self.lazy = [id_M]*self.N0 if array: assert N == len(array) self.data[self.N0:self.N0+self.N] = array for i in range(self.N0-1,0,-1): self.update(i) # 1点更新 def point_set(self, p, x): p += self.N0 for i in range(self.log, 0,-1): self.push(p>>i) self.data[p] = x for i in range(1, self.log + 1): self.update(p>>i) # 1点取得 def point_get(self, p): p += self.N0 for i in range(self.log, 0, -1): self.push(p>>i) return self.data[p] # 半開区間[L,R)をopでまとめる def prod(self, l, r): if l == r: return self.e_X l += self.N0 r += self.N0 for i in range(self.log, 0, -1): if (l>>i)<<i != l: self.push(l>>i) if (r>>i)<<i != r: self.push(r>>i) sml = smr = self.e_X while l < r: if l & 1: sml = self.op_X(sml, self.data[l]) l += 1 if r & 1: r -= 1 smr = self.op_X(self.data[r], smr) l >>= 1 r >>= 1 return self.op_X(sml, smr) # 全体をopでまとめる def all_prod(s): return self.data[1] # 1点作用 def apply(self, p, f): p += self.N0 for i in range(self.log, 0, -1): self.push(p>>i) self.data[p] = self.mapping(f, self.data[p]) for i in range(1, self.log + 1): self.update(p>>i) # 区間作用 def apply(self, l, r, f): if l == r: return l += self.N0 r += self.N0 for i in range(self.log, 0, -1): if (l>>i)<<i != l: self.push(l>>i) if (r>>i)<<i != r: self.push((r-1)>>i) l2, r2 = l, r while l < r: if l & 1: self.all_apply(l, f) l += 1 if r & 1: r -= 1 self.all_apply(r, f) l >>= 1 r >>= 1 l, r = l2, r2 for i in range(1, self.log + 1): if (l>>i)<<i != l: self.update(l>>i) if (r>>i)<<i != r: self.update((r-1)>>i) """ 始点 l を固定 f(x_l*...*x_{r-1}) が True になる最大の r つまり TTTTFFFF となるとき、F となる最小の添え字 存在しない場合 n が返る f(e_M) = True でないと壊れる """ def max_right(self, l, g): if l == self.N: return self.N l += self.N0 for i in range(self.log, 0, -1): self.push(l>>i) sm = self.e_X while True: while l&1 == 0: l >>= 1 if not g(self.op_X(sm, self.data[l])): while l < self.N0: self.push(l) l *= 2 if g(self.op_X(sm, self.data[l])): sm = self.op_X(sm, self.data[l]) l += 1 return l - self.N0 sm = self.op_X(sm, self.data[l]) l += 1 if l&-l == l: break return self.N """ 終点 r を固定 f(x_l*...*x_{r-1}) が True になる最小の l つまり FFFFTTTT となるとき、T となる最小の添え字 存在しない場合 r が返る f(e_M) = True でないと壊れる """ def min_left(self, r, g): if r == 0: return 0 r += self.N0 for i in range(self.log, 0, -1): self.push((r-1)>>i) sm = self.e_X while True: r -= 1 while r>1 and r&1: r >>= 1 if not g(self.op_X(self.data[r], sm)): while r < self.N0: self.push(r) r = 2*r + 1 if g(self.op_X(self.data[r], sm)): sm = self.op_X(self.data[r], sm) r -= 1 return r + 1 - self.N0 sm = self.op_X(self.data[r], sm) if r&-r == r: break return 0 # 以下内部関数 def update(self, k): self.data[k] = self.op_X(self.data[2*k], self.data[2*k+1]) def all_apply(self, k, f): self.data[k] = self.mapping(f, self.data[k]) if k < self.N0: self.lazy[k] = self.compose(f, self.lazy[k]) def push(self, k): #propagate と同じ if self.lazy[k] == self.id_M: return self.data[2*k ] = self.mapping(self.lazy[k], self.data[2*k]) self.data[2*k+1] = self.mapping(self.lazy[k], self.data[2*k+1]) if 2*k < self.N0: self.lazy[2*k] = self.compose(self.lazy[k], self.lazy[2*k]) self.lazy[2*k+1] = self.compose(self.lazy[k], self.lazy[2*k+1]) self.lazy[k] = self.id_M ################################################################### # ################################################################### import sys readline = sys.stdin.readline n, = map(int, readline().split()) *a, = map(int, readline().split()) q, = map(int, readline().split()) from operator import add e_X = INF = 1<<60 id_M = 0 op_X = min compose = add mapping = add seg = LazySegmentTree(op_X, e_X, mapping, compose, id_M, n, array=a) for _ in range(q): k,l,r,c = map(int, readline().split()) if k==1: seg.apply(l-1,r,c) else: print(seg.prod(l-1,r))