結果

問題 No.510 二次漸化式
ユーザー maspymaspy
提出日時 2020-04-03 01:16:06
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 759 ms / 3,000 ms
コード長 2,150 bytes
コンパイル時間 303 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 152,320 KB
最終ジャッジ日時 2024-06-28 15:40:20
合計ジャッジ時間 18,866 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
52,352 KB
testcase_01 AC 43 ms
52,608 KB
testcase_02 AC 418 ms
89,472 KB
testcase_03 AC 412 ms
89,088 KB
testcase_04 AC 405 ms
89,088 KB
testcase_05 AC 412 ms
89,472 KB
testcase_06 AC 566 ms
104,832 KB
testcase_07 AC 571 ms
105,088 KB
testcase_08 AC 588 ms
104,448 KB
testcase_09 AC 561 ms
104,448 KB
testcase_10 AC 253 ms
78,804 KB
testcase_11 AC 264 ms
79,640 KB
testcase_12 AC 257 ms
79,228 KB
testcase_13 AC 268 ms
78,976 KB
testcase_14 AC 259 ms
79,360 KB
testcase_15 AC 259 ms
79,248 KB
testcase_16 AC 416 ms
100,736 KB
testcase_17 AC 425 ms
100,864 KB
testcase_18 AC 434 ms
101,760 KB
testcase_19 AC 423 ms
100,608 KB
testcase_20 AC 422 ms
101,248 KB
testcase_21 AC 425 ms
101,376 KB
testcase_22 AC 410 ms
101,120 KB
testcase_23 AC 683 ms
122,624 KB
testcase_24 AC 673 ms
121,856 KB
testcase_25 AC 669 ms
121,856 KB
testcase_26 AC 670 ms
121,856 KB
testcase_27 AC 670 ms
121,856 KB
testcase_28 AC 661 ms
122,880 KB
testcase_29 AC 700 ms
122,256 KB
testcase_30 AC 690 ms
122,240 KB
testcase_31 AC 745 ms
152,320 KB
testcase_32 AC 745 ms
152,192 KB
testcase_33 AC 759 ms
152,192 KB
testcase_34 AC 460 ms
83,712 KB
testcase_35 AC 405 ms
83,456 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#!/usr/bin/ python3.8
import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines
MOD = 10 ** 9 + 7


def mat(x, y):
    return (1, x, 0, 0, 0, y * y, 2 * y, 1, 0, 0, y, 1, 0, 0, 0, 1)


class SegTree:
    """ segment tree with point modification and range product access. """
    data_unit = mat(0, 0)
    @classmethod
    def data_f(cls, A, B):
        C = [0] * 16
        for i in range(4):
            for j in range(i, 4):
                for k in range(j, 4):
                    C[4 * i + k] += B[4 * i + j] * A[4 * j + k]
        return tuple(x % MOD for x in C)

    def __init__(self, N):
        self.N = N
        self.data = [self.data_unit] * (N + N)

    def build(self, raw_data):
        data = self.data
        f = self.data_f
        N = self.N
        data[N:] = raw_data[:]
        for i in range(N - 1, 0, -1):
            data[i] = f(data[i << 1], data[i << 1 | 1])

    def set_val(self, i, x):
        data = self.data
        f = self.data_f
        i += self.N
        data[i] = x
        while i > 1:
            i >>= 1
            data[i] = f(data[i << 1], data[i << 1 | 1])

    def fold(self, L, R):
        """ compute for [L, R) """
        vL = vR = self.data_unit
        data = self.data
        f = self.data_f
        L += self.N
        R += self.N
        while L < R:
            if L & 1:
                vL = f(vL, data[L])
                L += 1
            if R & 1:
                R -= 1
                vR = f(data[R], vR)
            L >>= 1
            R >>= 1
        return f(vL, vR)


N = int(readline())
X = [0] * N
Y = [0] * N
Q = int(readline())

seg = SegTree(N)


for query in readlines()[:Q]:
    if query.startswith(b'a'):
        i = int(query[2:])
        if i == 0:
            print(1)
        else:
            A = seg.fold(0, i)
            print(sum(A[:4]) % MOD)
    elif query.startswith(b'x'):
        i, v = map(int, query[2:].split())
        X[i] = v
        seg.set_val(i, mat(v, Y[i]))
    elif query.startswith(b'y'):
        i, v = map(int, query[2:].split())
        Y[i] = v
        seg.set_val(i, mat(X[i], v))
0