結果
| 問題 | No.1050 Zero (Maximum) | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2022-03-23 23:23:33 | 
| 言語 | PyPy3 (7.3.15) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 135 ms / 2,000 ms | 
| コード長 | 571 bytes | 
| コンパイル時間 | 273 ms | 
| コンパイル使用メモリ | 82,780 KB | 
| 実行使用メモリ | 75,904 KB | 
| 最終ジャッジ日時 | 2024-10-12 02:50:35 | 
| 合計ジャッジ時間 | 2,516 ms | 
| ジャッジサーバーID (参考情報) | judge2 / judge4 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 15 | 
ソースコード
m,k = map(int,input().split())
mod = 10**9+7
def calc(A,B):
    ans = [[0]*m for i in range(m)]
    for i in range(m):
        for j in range(m):
            for t in range(m):
                ans[i][j] += A[i][t]*B[t][j]
                ans[i][j] %= mod
    return ans
A = [[0]*m for i in range(m)]
for i in range(m):
    for j in range(m):
        A[(i+j)%m][j] += 1
        A[(i*j)%m][j] += 1
base = [[0]*m for i in range(m)]
for i in range(m):
    base[i][i] = 1
while k:
    if k & 1:
        base = calc(A,base)
    A = calc(A,A)
    k >>= 1
print(base[0][0])
            
            
            
        