結果
| 問題 |
No.1105 Many Triplets
|
| コンテスト | |
| ユーザー |
rlangevin
|
| 提出日時 | 2023-07-03 08:51:03 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 45 ms / 2,000 ms |
| コード長 | 814 bytes |
| コンパイル時間 | 151 ms |
| コンパイル使用メモリ | 82,292 KB |
| 実行使用メモリ | 61,616 KB |
| 最終ジャッジ日時 | 2024-07-17 06:00:21 |
| 合計ジャッジ時間 | 2,384 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 25 |
ソースコード
import sys
input = sys.stdin.readline
def Matprod(A, B, mod, N):
temp = [0] * N*N
for i in range(N):
for j in range(N):
ij = i * N + j
for k in range(N):
temp[ij] += A[i*N+k] * B[k*N+j]
temp[ij] %= mod
return temp
def Matpow_Linear(A, M, mod, N):
Mat = [0] * N*N
for i in range(N):
Mat[i*N+i] = 1
while M:
if M & 1:
Mat = Matprod(Mat, A, mod, N)
A = Matprod(A, A, mod, N)
M >>= 1
return Mat
N = int(input())
X = list(map(int, input().split()))
A = [1, -1, 0, 0, 1, -1, -1, 0, 1]
mod = 10 ** 9 + 7
A = Matpow_Linear(A, N - 1, mod, 3)
ans = [0] * 3
for i in range(3):
for j in range(3):
ans[i] += A[3 * i + j] * X[j]
ans[i] %= mod
print(*ans)
rlangevin