結果
| 問題 | No.1050 Zero (Maximum) | 
| コンテスト | |
| ユーザー |  sasa8uyauya | 
| 提出日時 | 2024-08-31 11:47:27 | 
| 言語 | PyPy3 (7.3.15) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 141 ms / 2,000 ms | 
| コード長 | 461 bytes | 
| コンパイル時間 | 482 ms | 
| コンパイル使用メモリ | 82,292 KB | 
| 実行使用メモリ | 76,360 KB | 
| 最終ジャッジ日時 | 2024-08-31 11:47:30 | 
| 合計ジャッジ時間 | 2,934 ms | 
| ジャッジサーバーID (参考情報) | judge5 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 15 | 
ソースコード
m,k=map(int,input().split())
M=10**9+7
def seki(a,b):
  c=[[0]*m for i in range(m)]
  for i in range(m):
    for j in range(m):
      for k in range(m):
        c[i][j]+=a[i][k]*b[k][j]
        c[i][j]%=M
  return c
A=[[0]*m for i in range(m)]
for i in range(m):
  for j in range(m):
    A[(i+j)%m][i]+=1
    A[(i*j)%m][i]+=1
B=[[int(j==i) for j in range(m)] for i in range(m)]
for i in range(30):
  if (k>>i)&1:
    B=seki(A,B)
  A=seki(A,A)
print(B[0][0])
            
            
            
        