結果

問題 No.867 避難経路
ユーザー kwm_tkwm_t
提出日時 2023-05-01 23:04:03
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2,339 ms / 6,000 ms
コード長 2,641 bytes
コンパイル時間 3,087 ms
コンパイル使用メモリ 224,964 KB
実行使用メモリ 128,768 KB
最終ジャッジ日時 2024-04-30 20:58:17
合計ジャッジ時間 65,047 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
6,820 KB
testcase_01 AC 4 ms
6,944 KB
testcase_02 AC 4 ms
6,944 KB
testcase_03 AC 4 ms
6,944 KB
testcase_04 AC 4 ms
6,940 KB
testcase_05 AC 5 ms
6,944 KB
testcase_06 AC 4 ms
6,944 KB
testcase_07 AC 4 ms
6,944 KB
testcase_08 AC 5 ms
6,944 KB
testcase_09 AC 4 ms
6,940 KB
testcase_10 AC 3 ms
6,944 KB
testcase_11 AC 3 ms
6,940 KB
testcase_12 AC 4 ms
6,940 KB
testcase_13 AC 4 ms
6,944 KB
testcase_14 AC 3 ms
6,940 KB
testcase_15 AC 2,339 ms
128,640 KB
testcase_16 AC 2,260 ms
128,640 KB
testcase_17 AC 2,334 ms
128,768 KB
testcase_18 AC 2,258 ms
128,768 KB
testcase_19 AC 2,210 ms
128,640 KB
testcase_20 AC 2,123 ms
128,768 KB
testcase_21 AC 2,092 ms
128,768 KB
testcase_22 AC 2,131 ms
128,640 KB
testcase_23 AC 2,231 ms
128,640 KB
testcase_24 AC 2,096 ms
128,768 KB
testcase_25 AC 2,117 ms
128,768 KB
testcase_26 AC 2,120 ms
128,768 KB
testcase_27 AC 2,138 ms
128,768 KB
testcase_28 AC 2,172 ms
128,768 KB
testcase_29 AC 2,247 ms
128,640 KB
testcase_30 AC 1,784 ms
128,640 KB
testcase_31 AC 1,826 ms
128,640 KB
testcase_32 AC 1,809 ms
128,640 KB
testcase_33 AC 1,805 ms
128,640 KB
testcase_34 AC 1,860 ms
127,232 KB
testcase_35 AC 2,014 ms
127,360 KB
testcase_36 AC 1,989 ms
125,312 KB
testcase_37 AC 1,814 ms
119,552 KB
testcase_38 AC 1,794 ms
119,424 KB
testcase_39 AC 1,796 ms
119,552 KB
testcase_40 AC 1,846 ms
119,424 KB
testcase_41 AC 2 ms
6,940 KB
testcase_42 AC 3 ms
6,940 KB
testcase_43 AC 3 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
//#include <atcoder/all>
using namespace std;
//using namespace atcoder;
//using mint = modint1000000007;
//const int mod = 1000000007;
//using mint = modint998244353;
//const int mod = 998244353;
//const int INF = 1e9;
const long long LINF = 1e18;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i,l,r)for(int i=(l);i<(r);++i)
#define rrep(i, n) for (int i = (n-1); i >= 0; --i)
#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)
#define all(x) (x).begin(),(x).end()
#define allR(x) (x).rbegin(),(x).rend()
#define endl "\n"
#define P pair<long long,pair<int,int>>
template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }
const int dx[] = { 1,0,-1,0 };
const int dy[] = { 0,-1,0,1 };
int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int h, w; cin >> h >> w;
	int gx, gy; cin >> gx >> gy;
	gx--, gy--;
	vector a(h, vector<int>(w));
	rep(i, h)rep(j, w)cin >> a[i][j];
	vector dp(250, vector(h, vector<long long>(w, LINF)));
	rep(i, 250) {
		long long k = i;
		dp[i][gx][gy] = k * k + a[gx][gy];
		priority_queue<P, vector<P>, greater<P>> q;
		q.push({ dp[i][gx][gy], {gx, gy} });
		while (!q.empty()) {
			auto tmp = q.top();
			q.pop();
			long long cost = tmp.first;
			int x = tmp.second.first;
			int y = tmp.second.second;
			if (cost > dp[i][x][y]) continue;
			rep(j, 4) {
				int nx = x + dx[j];
				int ny = y + dy[j];
				if ((nx < 0) || (nx >= h) || (ny < 0) || (ny >= w))continue;
				long long ncost = cost + k * k + a[nx][ny];
				if (chmin(dp[i][nx][ny], ncost))q.push({ ncost, { nx, ny } });
			}
		}
	}
	vector dp2(h, vector<long long>(w, LINF));
	{
		set<pair<int, int>>st;
		dp2[gx][gy] = a[gx][gy];
		st.insert({ gx,gy });
		while (!st.empty()) {
			set<pair<int, int>>st2;
			for (auto e : st) {
				int x = e.first;
				int y = e.second;
				rep(j, 4) {
					int nx = x + dx[j];
					int ny = y + dy[j];
					if ((nx < 0) || (nx >= h) || (ny < 0) || (ny >= w))continue;
					if (LINF == dp2[nx][ny])st2.insert({ nx, ny });
					if (!st2.count({ nx,ny })) continue;
					long long ncost = dp2[x][y] + a[nx][ny];
					chmin(dp2[nx][ny],ncost);
				}
			}
			swap(st2, st);
		}
	}
	int q; cin >> q;
	while (q--) {
		int x, y; cin >> x >> y;
		x--, y--;
		long long k; cin >> k;
		if (k >= 250) {
			long long ans = dp2[x][y];
			long long d = abs(gx - x) + abs(gy - y) + 1;
			ans += k * k * d;
			cout << ans << endl;
		}
		else {
			cout << dp[k][x][y] << endl;
		}
	}
	return 0;
}
0