結果
問題 | 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 + 7res = [[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] %= modreturn resdef mat_pow(A, n):mod = 10 ** 9 + 7res = [[0] * 6 for _ in range(6)]for i in range(6):res[i][i] = 1while n:if n & 1:res = mat_mul(res, A)A = mat_mul(A, A)n >>= 1return resmod = 10 ** 9 + 7N = int(input())p = [1]for i in range(5):tmp = sum(p) * pow(6, mod - 2, mod) % modp.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] = 1p.reverse()An = mat_pow(A, N - 5)ans = 0for i in range(6):ans += p[i] * An[0][i]ans %= modprint(ans)