結果
問題 |
No.1275 綺麗な式
|
ユーザー |
|
提出日時 | 2021-01-04 14:56:54 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
TLE
|
実行時間 | - |
コード長 | 727 bytes |
コンパイル時間 | 381 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 17,700 KB |
最終ジャッジ日時 | 2024-10-15 00:35:26 |
合計ジャッジ時間 | 5,851 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 38 TLE * 1 -- * 21 |
ソースコード
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): 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) k >>= 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)