結果
| 問題 | 
                            No.147 試験監督(2)
                             | 
                    
| コンテスト | |
| ユーザー | 
                             ayaoni
                         | 
                    
| 提出日時 | 2020-12-15 19:39:12 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 1,798 ms / 2,000 ms | 
| コード長 | 1,263 bytes | 
| コンパイル時間 | 134 ms | 
| コンパイル使用メモリ | 82,016 KB | 
| 実行使用メモリ | 82,376 KB | 
| 最終ジャッジ日時 | 2024-09-20 02:02:35 | 
| 合計ジャッジ時間 | 8,002 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge3 / judge4 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 4 | 
ソースコード
import sys
sys.setrecursionlimit(10**7)
def I(): return int(sys.stdin.readline().rstrip())
def MI(): return map(int,sys.stdin.readline().rstrip().split())
def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))
def LI2(): return list(map(int,sys.stdin.readline().rstrip()))
def S(): return sys.stdin.readline().rstrip()
def LS(): return list(sys.stdin.readline().rstrip().split())
def LS2(): return list(sys.stdin.readline().rstrip())
def matrix_multiplication_mod(A,B,p):
    l = len(A)
    m = len(B)
    n = len(B[0])
    C = [[0]*n for _ in range(l)]
    for i in range(l):
        for j in range(n):
            C[i][j] = sum((A[i][k]*B[k][j]) % p for k in range(m)) % p
    return C
def matrix_power_mod(A,n,p):
    l = len(A)
    C = [[0] * l for _ in range(l)]
    for i in range(l):
        C[i][i] = 1
    while n > 0:
        if n % 2 == 1:
            C = matrix_multiplication_mod(C,A,p)
        A = matrix_multiplication_mod(A,A,p)
        n >>= 1
    return C
N = I()
A = [[1,1],[1,0]]
mod = 10**9+7
ans = 1
for _ in range(N):
    C,D = MI()
    if C % mod == 0:
        ans = 0
        break
    X = matrix_power_mod(A,C % (2*(mod+1)),mod)
    a = sum(X[0]) % mod
    ans *= pow(a,D % (mod-1),mod)
    ans %= mod
print(ans)
            
            
            
        
            
ayaoni