結果
問題 |
No.980 Fibonacci Convolution Hard
|
ユーザー |
![]() |
提出日時 | 2020-02-29 13:19:43 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 1,297 ms / 2,000 ms |
コード長 | 364 bytes |
コンパイル時間 | 145 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 103,300 KB |
最終ジャッジ日時 | 2024-10-13 19:57:40 |
合計ジャッジ時間 | 27,303 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 17 |
ソースコード
import sys, numpy as np m = 10**9+7 p, _ = int(input()), input() q = np.genfromtxt(sys.stdin, int) - 2 M = np.asarray([[-p, 0], [-p, 2]], 'i8')*pow(4+p**2, -1, m)%m N = np.matrix([[0, 1], [1, p]], 'i8') T = np.eye(2, 2**21, -1, 'i8') for i in range(21): T[:, 2**i:2**(i+1)] = N * T[:, :2**i] % m N = N**2%m C = M @ T[:, q] % m print(*((q*C[1]+C[0])%m), sep='\n')