結果
| 問題 |
No.147 試験監督(2)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-11-25 17:17:17 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,562 bytes |
| コンパイル時間 | 376 ms |
| コンパイル使用メモリ | 82,120 KB |
| 実行使用メモリ | 92,492 KB |
| 最終ジャッジ日時 | 2024-09-26 10:59:53 |
| 合計ジャッジ時間 | 12,581 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | WA * 1 TLE * 3 |
ソースコード
from collections import *
from itertools import *
from functools import *
import math,sys
input = sys.stdin.readline
mod = 10**9 + 7
def matrix_mul(A,B,mod = None):
nA = len(A)
mA = len(A[0])
mB = len(B[0])
tmp = [[0]*mB for _ in range(nA)]
tB = list(zip(*B))
if mod is None:
for i in range(nA):
Ai = A[i]
for j in range(mB):
Bj = tB[j]
tmp[i][j] = sum(Ai[k]*Bj[k] for k in range(mA))
return tmp
for i in range(nA):
Ai = A[i]
for j in range(mB):
Bj = tB[j]
for k in range(mA):
tmp[i][j] += Ai[k]*Bj[k]
tmp[i][j] %= mod
return tmp
def matrix_pow(A,n,mod = None):
nbit = list(str(bin(n))[2:])
nbit = [int(i) for i in nbit]
N = len(A)
C = [[0]*N for _ in range(N)]
B = A
for i in range(N):
C[i][i] = 1
if mod is None:
for i in range(len(nbit)):
if nbit[-1-i] == 1:
C = matrix_mul(C,B)
B = matrix_mul(B,B)
return C
for i in range(len(nbit)):
if nbit[-1-i] == 1:
C = matrix_mul(C,B,mod)
B = matrix_mul(B,B,mod)
return C
N = int(input())
CD = [tuple(map(int,input().split())) for _ in range(N)]
X = [[0,1],[1,1]]
I = [[1],[1]]
ans = 1
for c,d in CD:
Y = matrix_mul(matrix_pow(X,c-1,mod),I,mod)
a,b = Y[0][0],Y[1][0]
ans *= pow(a+b,d%(mod-1),mod)
ans %= mod
print(ans)