結果
問題 | No.659 徘徊迷路 |
ユーザー |
![]() |
提出日時 | 2023-10-18 12:47:56 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,573 bytes |
コンパイル時間 | 280 ms |
コンパイル使用メモリ | 81,968 KB |
実行使用メモリ | 79,752 KB |
最終ジャッジ日時 | 2024-09-18 09:08:24 |
合計ジャッジ時間 | 11,077 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 11 TLE * 1 |
ソースコード
import sys input = sys.stdin.readline def Matprod(A, B, 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] return temp def Matpow_Linear(A, M, N): Mat = [0] * N*N for i in range(N): Mat[i*N+i] = 1 pre = Mat[f(gx, gy) * N + f(sx, sy)] while M: if M & 1: Mat = Matprod(Mat, A, N) nex = Mat[f(gx, gy) * N + f(sx, sy)] if abs(pre - nex) < 1e-7 and nex > 1e-7: return Mat pre = nex A = Matprod(A, A, N) M >>= 1 return Mat H, W, T = map(int, input().split()) sx, sy = map(int, input().split()) gx, gy = map(int, input().split()) N = H * W M = [0] * N * N B = [] for i in range(H): B.extend(list(input().rstrip())) dx = [-1, 0, 1, 0] dy = [0, -1, 0, 1] def f(h, w): return h * W + w for i in range(H): for j in range(W): v = f(i,j) if B[v] == "#": continue cnt = 0 for k in range(4): x = i + dx[k] y = j + dy[k] if B[f(x,y)] == "#": continue cnt += 1 if cnt == 0: M[v * N + v] = 1 else: for k in range(4): x = i + dx[k] y = j + dy[k] if B[f(x,y)] == "#": continue M[f(x,y) * N + v] += 1 / cnt M = Matpow_Linear(M, T, N) print(M[f(gx, gy) * N + f(sx, sy)])