結果

問題 No.1234 典型RMQ
ユーザー convexineqconvexineq
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

  
"""
遅延セグメント木(区間演算、区間更新)
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))
        
        
0