結果

問題 No.867 避難経路
ユーザー aaaaaaaaaa2230
提出日時 2022-05-21 12:48:57
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,580 bytes
コンパイル時間 513 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 371,288 KB
最終ジャッジ日時 2024-09-20 11:34:40
合計ジャッジ時間 28,570 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 3
other AC * 28 TLE * 3 -- * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline
from heapq import heappush,heappop

# def divmod(x,mod):
    
#     a = x//mod
#     return a,x-a*mod
h,w = map(int,input().split())
gx,gy = [int(x)-1 for x in input().split()]
A = [list(map(int,input().split())) for i in range(h)]

inf = 10**20
dp = [[[inf]*(500) for i in range(w)] for j in range(h)]
dist = [[inf]*w for i in range(h)]
step = [[inf]*w for i in range(h)]
dist[gx][gy] = A[gx][gy]
dp[gx][gy][0] = A[gx][gy]
step[gx][gy] = 0
H = []
dx = [1,0,-1,0]
dy = [0,1,0,-1]
mod2 = 1<<16
heappush(H,(A[gx][gy]<<16)+gx*w+gy)
for d in range(500):
    nH = []

    while H:
        num = heappop(H)
        c,num = divmod(num,mod2)
        x,y = divmod(num,w)
        if dp[x][y][d] != c:
            continue
        # print(d,c,x,y)
        for i in range(4):
            nx = x+dx[i]
            ny = y+dy[i]
            if nx < 0 or nx >= h or ny < 0 or ny >= w:
                continue
            nc = c + A[nx][ny]
            if dist[nx][ny] <= nc:
                continue
            if d < 499 and dp[nx][ny][d+1] > nc:
                dp[nx][ny][d+1] = nc
                dist[nx][ny] = nc
                step[nx][ny] = d+1
                nnum = (nc<<16)+nx*w+ny
                heappush(nH,nnum)
    if nH:
        H = nH
    else:
        break


q = int(input())
for _ in range(q):
    x,y,k = map(int,input().split())
    x -= 1
    y -= 1

    ans = inf
    k = k**2
    nc = k
    ndp = dp[x][y]

    for i in range(step[x][y]+1):
        if ans > ndp[i]+nc:
            ans = ndp[i]+nc
        nc += k

    print(ans)
0