結果
問題 | No.1073 無限すごろく |
ユーザー |
![]() |
提出日時 | 2020-05-08 00:49:42 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 35 ms / 2,000 ms |
コード長 | 607 bytes |
コンパイル時間 | 72 ms |
コンパイル使用メモリ | 12,928 KB |
実行使用メモリ | 11,008 KB |
最終ジャッジ日時 | 2024-07-03 09:43:33 |
合計ジャッジ時間 | 1,977 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
def mmult(A, B): global mod n, m, l = len(A), len(B), len(B[0]) ret = [[0]*l for _ in range(n)] for i in range(n): for j in range(m): for k in range(l): ret[i][k] = (ret[i][k]+A[i][j]*B[j][k])%mod return ret def mpow(A, n): if n == 0: return [[1 if i==j else 0 for j in range(len(A))] for i in range(len(A))] return mmult(mpow(A, n-1), A) if n % 2 else mpow(mmult(A, A), n//2) mod = 10**9+7 t = (mod + 1) // 6 X = [[0,0,0,0,0,t],[1,0,0,0,0,t],[0,1,0,0,0,t],[0,0,1,0,0,t],[0,0,0,1,0,t],[0,0,0,0,1,t]] N = int(input()) print(mpow(X, N)[-1][-1])