結果
問題 | No.840 ほむほむほむら |
ユーザー |
|
提出日時 | 2019-08-20 00:29:07 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 886 ms / 4,000 ms |
コード長 | 864 bytes |
コンパイル時間 | 188 ms |
コンパイル使用メモリ | 82,080 KB |
実行使用メモリ | 77,256 KB |
最終ジャッジ日時 | 2024-10-05 13:26:18 |
合計ジャッジ時間 | 6,916 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 25 |
ソースコード
N, k = list(map(int, input().split())) mod = 998244353 mat_len = k ** 3 def mul(A, B) : C = [[0 for _ in range(mat_len)] for _ in range(mat_len)] for i in range (mat_len): for j in range(mat_len): for k in range(mat_len): C[i][j]+=A[i][k]*B[k][j]%mod C[i][j]%=mod return C B = [[0] * mat_len for _ in range(mat_len)] for i in range(k): for j in range(k): for l in range(k): cur = i*k*k+j*k+l B[i*k*k+j*k+(l+1)%k][cur] += 1 B[i*k*k+(j+l)%k*k+l][cur] += 1 B[(i+j)%k*k*k+j*k+l][cur] += 1 ans = [[0] * mat_len for _ in range(mat_len)] for i in range(mat_len): ans[i][i] = 1 while N: if N & 1: ans = mul(ans, B) B=mul(B, B) N>>=1 ret = 0 for i in range(k**2): ret = (ret + ans[i][0])%mod print(int(ret))