結果
問題 | 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*N for i in range(N): for j in range(N): ij = i * N + j for k in range(N): temp[ij] += A[i*N+k] * B[k*N+j] temp[ij] %= mod return temp def Matpow_Linear(A, M, mod, N): Mat = [0] * N*N for i in range(N): Mat[i*N+i] = 1 while M: if M & 1: Mat = Matprod(Mat, A, mod, N) A = Matprod(A, A, mod, N) M >>= 1 return Mat C, N, M = map(int, input().split()) mod = 10 ** 9 + 7 inv = pow(C, mod-2, mod) a1 = (C - 1) * inv % mod a2 = (C - 2) * inv % mod A = [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])%mod ans = 1 - pow(1 - p, M, mod) print(ans%mod)