結果

問題 No.673 カブトムシ
ユーザー roaris
提出日時 2020-06-24 04:25:20
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 39 ms / 2,000 ms
コード長 656 bytes
コンパイル時間 186 ms
コンパイル使用メモリ 82,104 KB
実行使用メモリ 53,880 KB
最終ジャッジ日時 2024-07-03 19:52:47
合計ジャッジ時間 1,608 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

def mat_mul(A, B):
    C = [[0]*len(B[0]) for _ in range(len(A))]
    
    for i in range(len(A)):
        for j in range(len(B[0])):
            for k in range(len(A[0])):
                C[i][j] += A[i][k]*B[k][j]
                C[i][j] %= MOD
    
    return C

def mat_pow(A, n):
    res = [[0]*len(A) for _ in range(len(A))]
    
    for i in range(len(A)):
        res[i][i] = 1
    
    while n>0:
        if n&1:
            res = mat_mul(res, A)
        
        n >>= 1
        A = mat_mul(A, A)
    
    return res

B, C, D = map(int, input().split())
MOD = 10**9+7
res = mat_mul(mat_pow([[C, B*C], [0, 1]], D-1), [[B*C], [1]])
print(res[0][0])
0