結果
問題 | No.1810 RGB Biscuits |
ユーザー | lilictaka |
提出日時 | 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 |
ソースコード
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')