結果
| 問題 |
No.147 試験監督(2)
|
| コンテスト | |
| ユーザー |
kmjp
|
| 提出日時 | 2015-02-09 00:44:11 |
| 言語 | PyPy2 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 1,051 ms / 2,000 ms |
| コード長 | 753 bytes |
| コンパイル時間 | 876 ms |
| コンパイル使用メモリ | 77,196 KB |
| 実行使用メモリ | 81,280 KB |
| 最終ジャッジ日時 | 2024-06-23 10:07:35 |
| 合計ジャッジ時間 | 5,877 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 4 |
ソースコード
import sys
import math
def matmult(A,B):
n=len(A)
C=[[0 for i in range(n)] for j in range(n)]
C[0][0]=(A[0][0]*B[0][0]+A[0][1]*B[1][0])%mo
C[0][1]=(A[0][0]*B[0][1]+A[0][1]*B[1][1])%mo
C[1][0]=(A[1][0]*B[0][0]+A[1][1]*B[1][0])%mo
C[1][1]=(A[1][0]*B[0][1]+A[1][1]*B[1][1])%mo
return list(C)
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=input()
ret = 1
mo=1000000007
for i in range(N):
C,D=map(int,raw_input().strip().split(" "))
if D>mo-1:
D %= mo-1
ret = ret * pow(pat(C),D,mo) % mo
print ret
kmjp