結果

問題 No.801 エレベーター
ユーザー stngstng
提出日時 2022-07-03 19:54:49
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 1,387 bytes
コンパイル時間 579 ms
コンパイル使用メモリ 82,276 KB
実行使用メモリ 67,456 KB
最終ジャッジ日時 2024-05-07 03:39:58
合計ジャッジ時間 3,438 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

# 行列の乗算(mod)
from regex import R


def mat_mul(a, b):
    I, K, J = len(a), len(b), len(b[0])
    c = [[0 for j in range(J)] for i in range(I)]
    for i in range(I):
        for k in range(K):
            for j in range(J):
                c[i][j] += a[i][k] * b[k][j]
                c[i][j] %= mod
    return c

# 行列の累乗(mod)
def mat_pow(a, n):
    b = [[0 for j in range(len(a))] for i in range(len(a))]
    for i in range(len(a)):
        b[i][i] = 1
    while n > 0:
        if n & 1:
            b = mat_mul(b, a)
        a = mat_mul(a, a)
        n >>= 1
    return b

class Bit:
    def __init__(self, n):
        self.size = n
        self.tree = [0] * (n + 1)
 
    def sum(self, i):
        s = 0
        while i > 0:
            s += self.tree[i]
            i -= i & -i
        return s
 
    def add(self, i, x):
        while i <= self.size:
            self.tree[i] += x
            i += i & -i

mod = 10**9+7
n,m,k = map(int,input().split())

st = [[0]*n for i in range(n)]

lr = [[int(i) for i in input().split()] for j in range(m)]

for i in range(1,n+1):
    ki = Bit(n+1)
    for j in range(m):
        l,r = lr[j]
        if l <= i <= r:
            #print(i,j,2)
            ki.add(l,1)
            ki.add(r+1,-1)
    for j in range(1,n+1):
        st[i-1][j-1] = ki.tree[j]
#print(st)
#print(ki)
#exit()
A = mat_pow(st,k)

print(A[0][n-1]%mod)
0