結果
問題 | No.2441 行列累乗 |
ユーザー |
|
提出日時 | 2023-08-25 21:21:54 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 37 ms / 2,000 ms |
コード長 | 771 bytes |
コンパイル時間 | 99 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 10,752 KB |
最終ジャッジ日時 | 2024-12-24 07:40:58 |
合計ジャッジ時間 | 1,626 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
def mat_mul(a, b): """ a: 行列(2次元配列)I*J b: 行列(2次元配列)J*K """ MOD = 1000000007 I, J, K = len(a), len(b[0]), len(b) c = [[0] * J for _ in range(I)] for i in range(I): for j in range(J): for k in range(K): c[i][j] += a[i][k] * b[k][j] return c def pow_mat(A, n): """ A: 正方行列(2次元配列) n: 累乗指数 """ p = [[0]*len(A) for _ in range(len(A))] # 基本行列にする for i in range(len(A)): p[i][i] = 1 while n > 0: if n & 1: p = mat_mul(A, p) A = mat_mul(A, A) n >>= 1 return p M = [list(map(int, input().split())) for _ in range(2)] M = pow_mat(M, 3) for c in M: print(*c)