結果

問題 No.659 徘徊迷路
ユーザー H3PO4H3PO4
提出日時 2020-11-06 08:52:12
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
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
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 494 ms
44,384 KB
testcase_01 AC 489 ms
44,512 KB
testcase_02 AC 484 ms
44,764 KB
testcase_03 AC 486 ms
44,508 KB
testcase_04 AC 492 ms
44,600 KB
testcase_05 AC 492 ms
44,832 KB
testcase_06 AC 492 ms
44,772 KB
testcase_07 AC 489 ms
44,252 KB
testcase_08 AC 499 ms
45,020 KB
testcase_09 AC 497 ms
44,124 KB
testcase_10 AC 502 ms
44,508 KB
testcase_11 AC 500 ms
44,516 KB
testcase_12 AC 501 ms
46,088 KB
testcase_13 AC 501 ms
44,508 KB
testcase_14 AC 495 ms
44,324 KB
testcase_15 AC 498 ms
44,636 KB
testcase_16 AC 501 ms
44,640 KB
権限があれば一括ダウンロードができます

ソースコード

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