結果

問題 No.3492 区間冪乗加算一点取得
コンテスト
ユーザー tassei903
提出日時 2026-04-03 23:00:41
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
WA  
実行時間 -
コード長 1,990 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 182 ms
コンパイル使用メモリ 85,084 KB
実行使用メモリ 260,556 KB
最終ジャッジ日時 2026-04-03 23:01:53
合計ジャッジ時間 29,646 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4 WA * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import sys
ni = lambda :int(input())
na = lambda :list(map(int,input().split()))
yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")
no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")
#######################################################################

class DualSegmentTree:
    def __init__(self, size, f, default):
        self.n = (size-1).bit_length()
        self.size = 1<<self.n
        self.default = default
        self.lazy = [default]*(self.size*2)
        self.f = f

    def propagate(self, k):
        lazy = self.lazy
        if lazy[k] != self.default:
            lazy[2*k] = self.f(lazy[2*k], lazy[k])
            lazy[2*k+1] = self.f(lazy[2*k+1], lazy[k])
            lazy[k] = self.default

    def thrust(self, k):
        for i in range(self.n,0,-1):
            self.propagate(k>>i)

    def update(self, a, b, x):
        a += self.size
        b += self.size-1
        self.thrust(a)
        self.thrust(b)
        l = a
        r = b + 1
        lazy = self.lazy
        while l < r:
            if l & 1:
                lazy[l] = self.f(lazy[l], x)
                l += 1
            if r & 1:
                r -= 1
                lazy[r] = self.f(lazy[r], x)
            l >>= 1
            r >>= 1
    
    def get(self, k):
        k += self.size
        self.thrust(k)
        return self.lazy[k]

n, b, q = na()
k = 101
def op(x, y):
    return [(x[i] + y[i]) % b for i in range(k)]

a = DualSegmentTree(n, op, [0] * k)

for _ in range(q):
    l, m, r, c, d = na()
    l -= 1
    m -= 1
    c += 1
    f = [0] * k
    f[0] = 1
    for i in range(d):
        nf = [0] * k
        for j in range(i + 1):
            nf[j + 1] = (nf[j+1] + f[j]) % b
            nf[j] = (nf[j] + f[j] * c) % b
        f = nf
        # print(f)
    # print("!", f)
    a.update(l, r, f)
    res = a.get(m)

    z = 1
    r = 0
    for i in range(k):
        r = (r + z * f[i]) % b
        z = m * z % b
    print(r)
0