結果
問題 | No.1073 無限すごろく |
ユーザー |
![]() |
提出日時 | 2020-06-05 21:58:35 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 36 ms / 2,000 ms |
コード長 | 895 bytes |
コンパイル時間 | 207 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 11,008 KB |
最終ジャッジ日時 | 2024-12-17 14:45:49 |
合計ジャッジ時間 | 2,104 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
def mat_mul(A, B): mod = 10 ** 9 + 7 res = [[0]*6 for _ in range(6)] for i in range(6): for k in range(6): for j in range(6): res[i][j] += A[i][k] * B[k][j] res[i][j] %= mod return res def mat_pow(A, n): mod = 10 ** 9 + 7 res = [[0] * 6 for _ in range(6)] for i in range(6): res[i][i] = 1 while n: if n & 1: res = mat_mul(res, A) A = mat_mul(A, A) n >>= 1 return res mod = 10 ** 9 + 7 N = int(input()) p = [1] for i in range(5): tmp = sum(p) * pow(6, mod - 2, mod) % mod p.append(tmp) if N < 6: print(p[N]) exit() A = [[0] * 6 for _ in range(6)] for i in range(6): A[0][i] = p[1] for i in range(1, 6): A[i][i-1] = 1 p.reverse() An = mat_pow(A, N - 5) ans = 0 for i in range(6): ans += p[i] * An[0][i] ans %= mod print(ans)