結果
問題 |
No.1275 綺麗な式
|
ユーザー |
|
提出日時 | 2021-01-04 14:52:33 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
OLE
|
実行時間 | - |
コード長 | 758 bytes |
コンパイル時間 | 328 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 16,128 KB |
最終ジャッジ日時 | 2024-10-15 00:31:50 |
合計ジャッジ時間 | 6,867 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | -- * 2 |
other | OLE * 1 -- * 59 |
ソースコード
mod = 10 ** 9 + 7 def matrix_multiplication(a, b): n = len(a) res = [[0] * n for i in range(n)] for i in range(n): for k in range(n): for j in range(n): print(b[k][j]) res[i][j] += a[i][k] * b[k][j] res[i][j] %= mod return res def matrix_pow(a, k): if k == 1: return a n = len(a) res = [[0] * n for i in range(n)] for i in range(n): res[i][i] = 1 while k: if k & 1: res = matrix_multiplication(res, a) a = matrix_multiplication(a, a) n >>= 1 return res A, B = map(int, input().split()) n = int(input()) M = [[A * 2, B - A ** 2], [1, 0]] M = matrix_pow(M, n - 1) print((M[0][0] * 2 * A + M[0][1] * 2) % mod)