結果
| 問題 | No.1049 Zero (Exhaust) | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2022-04-17 17:05:31 | 
| 言語 | PyPy3 (7.3.15) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 41 ms / 2,000 ms | 
| コード長 | 644 bytes | 
| コンパイル時間 | 365 ms | 
| コンパイル使用メモリ | 81,976 KB | 
| 実行使用メモリ | 52,480 KB | 
| 最終ジャッジ日時 | 2024-12-26 10:46:26 | 
| 合計ジャッジ時間 | 2,426 ms | 
| ジャッジサーバーID (参考情報) | judge3 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 22 | 
ソースコード
MOD = 10 ** 9 + 7
def matpow(A, B, w):
    l = len(A)
    while w:
        if w & 1:
            C = [0] * l
            for i in range(l):
                for j in range(l):
                    C[i] += A[i][j] * B[j]
                    C[i] %= MOD
            B = C
        C = [[0] * l for _ in range(l)]
        for i in range(l):
            for j in range(l):
                for k in range(l):
                    C[i][j] += A[i][k] * A[k][j]
                    C[i][j] %= MOD
        A = C
        w >>= 1
    return B
p, k = map(int, input().split())
B = [1, 0]
A = [[p + 1, p - 1], [2, 2 * p - 2]]
B = matpow(A, B, k)
print(B[0])
            
            
            
        