結果

問題 No.659 徘徊迷路
ユーザー H3PO4
提出日時 2020-11-06 08:52:12
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 502 ms / 2,000 ms
コード長 877 bytes
コンパイル時間 100 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 46,088 KB
最終ジャッジ日時 2024-07-22 11:40:31
合計ジャッジ時間 10,530 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 5
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff #

import numpy as np

R, C, T = map(int, input().split())
Sy, Sx = map(int, input().split())
Gy, Gx = map(int, input().split())
B = [input() for _ in range(R)]

grid2int = lambda y, x: y * C + x

matr = [[0] * (R * C) for _ in range(R * C)]

for y in range(1, R - 1):
    for x in range(1, C - 1):
        if B[y][x] == '#':
            continue

        coor = grid2int(y, x)
        movable = []
        for dy, dx in ((1, 0), (0, 1), (-1, 0), (0, -1)):
            ydy, xdx = y + dy, x + dx
            if B[ydy][xdx] == '.':
                movable.append(grid2int(ydy, xdx))

        if not movable:
            matr[coor][coor] = 1
            continue

        l = len(movable)
        for p in movable:
            matr[coor][p] = 1 / l
matr = np.array(matr)
X = np.zeros(R * C)
X[grid2int(Sy, Sx)] = 1

Y = X @ np.linalg.matrix_power(matr, T)
print(Y[grid2int(Gy, Gx)])
0