結果

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

ソースコード

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 = []

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 dx,dy in ((1,0),(0,1),(0,-1),(-1,0)):
            nx = x+dx
            ny = y+dy
            if 0 <= nx < h and 0 <= ny < w:
                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