結果

問題 No.1105 Many Triplets
ユーザー 双六双六
提出日時 2020-07-21 21:06:57
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 473 ms / 2,000 ms
コード長 822 bytes
コンパイル時間 74 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 44,500 KB
最終ジャッジ日時 2024-06-09 21:46:05
合計ジャッジ時間 16,058 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 472 ms
43,988 KB
testcase_01 AC 447 ms
44,112 KB
testcase_02 AC 449 ms
44,248 KB
testcase_03 AC 446 ms
43,988 KB
testcase_04 AC 448 ms
43,860 KB
testcase_05 AC 442 ms
44,112 KB
testcase_06 AC 444 ms
44,000 KB
testcase_07 AC 440 ms
44,244 KB
testcase_08 AC 449 ms
44,244 KB
testcase_09 AC 444 ms
44,004 KB
testcase_10 AC 449 ms
43,996 KB
testcase_11 AC 446 ms
44,496 KB
testcase_12 AC 443 ms
44,120 KB
testcase_13 AC 442 ms
43,992 KB
testcase_14 AC 458 ms
44,108 KB
testcase_15 AC 454 ms
44,240 KB
testcase_16 AC 449 ms
44,248 KB
testcase_17 AC 446 ms
44,372 KB
testcase_18 AC 443 ms
44,112 KB
testcase_19 AC 449 ms
44,372 KB
testcase_20 AC 443 ms
44,240 KB
testcase_21 AC 446 ms
43,988 KB
testcase_22 AC 460 ms
44,124 KB
testcase_23 AC 452 ms
43,860 KB
testcase_24 AC 442 ms
44,240 KB
testcase_25 AC 454 ms
44,240 KB
testcase_26 AC 473 ms
44,244 KB
testcase_27 AC 442 ms
44,500 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys; input = sys.stdin.buffer.readline
sys.setrecursionlimit(10**7)
import numpy as np
mod = 10 ** 9 + 7

def getlist():
	return list(map(int, input().split()))

def matrixPow(A, n, N):
	#B:単位行列 演算の種類によって初期化方法を変える必要もある
	B = np.full((N, N), 0)
	for i in range(N):
		B[i, i] = 1

	while n > 0:
		if n & 1 == 1:
			B = A @ B
			B = np.mod(B, mod)
		A = A @ A
		A = np.mod(A, mod)
		n = n >> 1

	return B


#処理内容
def main():
	N = int(input())
	A1, B1, C1 = getlist()
	start = np.array([[A1], [B1], [C1]])
	mul = np.array([[1, -1, 0], [0, 1, -1], [-1, 0, 1]])
	B = matrixPow(mul, N - 1, 3)
	# print(B)
	ans = B @ start
	ans = np.mod(ans, mod)
	# print(ans)
	anslist = [ans[0, 0], ans[1, 0], ans[2, 0]]
	print(*anslist)





if __name__ == '__main__':
	main()
0