結果

問題 No.1810 RGB Biscuits
ユーザー lilictakalilictaka
提出日時 2022-05-09 17:01:39
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 257 ms / 2,000 ms
コード長 845 bytes
コンパイル時間 289 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 78,448 KB
最終ジャッジ日時 2024-07-16 15:01:07
合計ジャッジ時間 4,645 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #

def mulmod(a,b,mod):
    return a*b%mod
def mat_mul_Pypy(A,B):
    return [[sum([mulmod(A[i][k],B[k][j],MOD) for k in range(len(B))])%MOD for j in range(len(B[0]))] for i in range(len(A))]

def mat_power_Pypy(A, N):#行列累乗
    P = [[1 if i==j else 0 for i in range(len(A[0]))] for j in range(len(A))]
    while N:
        if N & 1:
            P = mat_mul_Pypy(P, A)
        A = mat_mul_Pypy(A, A)
        N >>= 1
    return P
MOD = 10 ** 9 + 7
A,B = map(int,input().split())
N = int(input())
ANS = []
def solve(t):
    M = [[A,B],[1,0]]
    tmp = mat_power_Pypy(M,t//2)
    board = mat_mul_Pypy(tmp,[[1],[1]])
    if t % 2 == 0:
        return (board[0][0] + board[1][0]) % MOD
    else:
        return (board[0][0] * (A+1) + board[1][0] * (B+1)) % MOD
for _ in range(N):
    t = int(input())
    ANS.append(solve(t))
print(*ANS,sep='\n')
0