結果

問題 No.567 コンプリート
ユーザー kurimupy
提出日時 2020-06-18 02:05:47
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 473 ms / 2,000 ms
コード長 494 bytes
コンパイル時間 87 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 44,748 KB
最終ジャッジ日時 2024-07-03 12:37:33
合計ジャッジ時間 8,996 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

import numpy as np

N=int(input())
def mat(A, N):
    if N == 0:
        return np.eye(6)
    elif N == 1:
        return A
    else:
        if N % 2 == 0:
            X = mat(A, N//2)
            return np.dot(X, X)
        else:
            X = mat(A, N//2)
            return np.dot(A, np.dot(X, X))


A = np.zeros((6, 6),dtype="float")
for i in range(5):
    A[i][i] += i+1
    A[i][i+1] += 5-i
A[5][5]=6
A /= 6
X0 = np.array([1, 0, 0, 0, 0, 0])
X0 = np.dot(X0, mat(A, N-1))

print(X0[-1])
0