結果

問題 No.572 妖精の演奏
ユーザー tktk_snsn
提出日時 2021-04-23 12:14:19
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 267 ms / 2,000 ms
コード長 704 bytes
コンパイル時間 134 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 11,776 KB
最終ジャッジ日時 2024-07-04 06:59:37
合計ジャッジ時間 2,481 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline


def matmul(A, B):
    sz = len(A)
    res = [[0]*sz for _ in range(sz)]
    for i in range(sz):
        for j in range(sz):
            for k in range(sz):
                res[i][j] = max(res[i][j], A[i][k]+B[k][j])
    return res


N = int(input())
M = int(input())
A = list(list(map(int,input().split())) for _ in range(M))
D = N.bit_length()
mat = [-1]*D
mat[0] = A
for i in range(1, D):
    mat[i] = matmul(mat[i-1],mat[i-1])

dp = [0]*M
N -= 1
for d in reversed(range(D)):
    if (N>>d)&1:
        ndp = [0]*M
        for i in range(M):
            for j in range(M):
                ndp[j] = max(ndp[j], dp[i] + mat[d][i][j])
        dp = ndp
print(max(dp))
0