結果
問題 | No.1706 Many Bus Stops (hard) |
ユーザー |
![]() |
提出日時 | 2023-11-27 12:40:43 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 49 ms / 2,000 ms |
コード長 | 816 bytes |
コンパイル時間 | 363 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 60,800 KB |
最終ジャッジ日時 | 2024-09-26 12:20:44 |
合計ジャッジ時間 | 3,325 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 41 |
ソースコード
def Matprod(A, B, mod, N):temp = [0] * N*Nfor i in range(N):for j in range(N):ij = i * N + jfor k in range(N):temp[ij] += A[i*N+k] * B[k*N+j]temp[ij] %= modreturn tempdef Matpow_Linear(A, M, mod, N):Mat = [0] * N*Nfor i in range(N):Mat[i*N+i] = 1while M:if M & 1:Mat = Matprod(Mat, A, mod, N)A = Matprod(A, A, mod, N)M >>= 1return MatC, N, M = map(int, input().split())mod = 10 ** 9 + 7inv = pow(C, mod-2, mod)a1 = (C - 1) * inv % moda2 = (C - 2) * inv % modA = [inv, 0, 0, inv,1, 0, 0, 0,0, a1, inv, a2,0, 0, 1, 0]A = Matpow_Linear(A, N - 1, mod, 4)p = (A[0] * inv + A[1])%modans = 1 - pow(1 - p, M, mod)print(ans%mod)