結果

問題 No.1340 おーじ君をさがせ
コンテスト
ユーザー Shinya Fujita
提出日時 2025-01-19 21:44:42
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 912 ms / 2,000 ms
コード長 994 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,030 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 80,896 KB
最終ジャッジ日時 2025-01-19 21:45:07
合計ジャッジ時間 15,437 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 59
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

N, M, T = map(int, input().split())


def eye(n):
    res = [[0]*n for _ in range(n)]
    for i in range(n):
        res[i][i] = 1
    return res

def prod_matrix(m1, m2):
    # return m1 * m2
    n = len(m1)
    res = [[0]*n for _ in range(n)]
    for i in range(n):
        for j in range(n):
            for k in range(n):
                res[i][j] |= m1[i][k] * m2[k][j]
    return res

def exp_matrix(matrix, times):
    res = eye(len(matrix))
    while times:
        if times & 1:
            res = prod_matrix(res, matrix)
        matrix = prod_matrix(matrix, matrix)
        times >>= 1
    return res

def prod_matrix_vex(matrix, vec):
    res = [0] * len(vec)
    for i, m in enumerate(matrix):
        for mi, vi in zip(m, vec):
            res[i] |= mi * vi
    return res


mat = [[0]*N for _ in range(N)]
for _ in range(M):
    a, b = map(int, input().split())
    mat[b][a] = 1

mat = exp_matrix(mat, T)

s0 = [0] * N
s0[0] = 1
state = prod_matrix_vex(mat, s0)
print(sum(state))
0