結果

問題 No.1340 おーじ君をさがせ
ユーザー Shinya Fujita
提出日時 2025-01-19 21:38:40
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 993 bytes
コンパイル時間 241 ms
コンパイル使用メモリ 82,708 KB
実行使用メモリ 81,024 KB
最終ジャッジ日時 2025-01-19 21:39:03
合計ジャッジ時間 15,167 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 39 WA * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

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[a][b] = 1

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