結果
問題 |
No.147 試験監督(2)
|
ユーザー |
![]() |
提出日時 | 2017-10-18 09:50:56 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,158 bytes |
コンパイル時間 | 271 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 29,772 KB |
最終ジャッジ日時 | 2024-11-18 14:27:26 |
合計ジャッジ時間 | 15,982 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | TLE * 4 |
ソースコード
def get_inner_prod(raw, col): return sum([x*y for (x,y) in zip(raw,col)]) def get_matrix_trans(matrix): # malloc n*m matrix matrix_t = [[0 for j in range(len(matrix))] for i in range(len(matrix[0]))] for i in range(len(matrix)): for j in range(len(matrix[i])): matrix_t[j][i] = matrix[i][j] return matrix_t def matmult(matrix1, matrix2): # malloc A_m*B_n matrix result = [[0 for j in range(len(matrix2[0]))] for i in range(len(matrix1))] matrix2_t = get_matrix_trans(matrix2) for i in range(len(result)): for j in range(len(result[0])): result[i][j] = get_inner_prod(matrix1[i], matrix2_t[j]) return result #二乗法の累乗 def matpow(A,p): n=len(A) A=list(A) R=[[0 for i in range(n)] for j in range(n)] for i in range(n): R[i][i]=1 while p: if p%2: R = matmult(A,R) A=matmult(A,A) p >>= 1 return R def pat(c): A=[[1,1],[1,0]] B=matpow(A,c-1) return (B[0][0]+B[0][1]+B[1][0]+B[1][1])%mo N=int(input()) ret = 1 mo=1000000007 for i in range(N): C,D=map(int,input().split()) if D>mo-1: D %= mo-1 ret = ret * pow(pat(C),D,mo) % mo print(ret)