結果

問題 No.867 避難経路
ユーザー risujirohrisujiroh
提出日時 2019-08-16 23:08:12
言語 C++14
(gcc 11.2.0 + boost 1.78.0)
結果
AC  
実行時間 5,085 ms / 6,000 ms
コード長 2,690 bytes
コンパイル時間 2,001 ms
使用メモリ 29,972 KB
最終ジャッジ日時 2022-11-22 12:40:56
合計ジャッジ時間 107,199 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 1 ms
3,436 KB
testcase_01 AC 1 ms
3,432 KB
testcase_02 AC 2 ms
3,620 KB
testcase_03 AC 2 ms
3,556 KB
testcase_04 AC 2 ms
3,536 KB
testcase_05 AC 1 ms
3,536 KB
testcase_06 AC 2 ms
3,580 KB
testcase_07 AC 2 ms
3,528 KB
testcase_08 AC 2 ms
3,528 KB
testcase_09 AC 2 ms
3,496 KB
testcase_10 AC 2 ms
3,492 KB
testcase_11 AC 2 ms
3,568 KB
testcase_12 AC 2 ms
3,480 KB
testcase_13 AC 2 ms
3,480 KB
testcase_14 AC 2 ms
3,552 KB
testcase_15 AC 2,232 ms
29,840 KB
testcase_16 AC 2,625 ms
29,832 KB
testcase_17 AC 2,443 ms
29,964 KB
testcase_18 AC 2,370 ms
29,848 KB
testcase_19 AC 2,345 ms
29,780 KB
testcase_20 AC 434 ms
29,840 KB
testcase_21 AC 394 ms
29,848 KB
testcase_22 AC 606 ms
29,972 KB
testcase_23 AC 665 ms
29,788 KB
testcase_24 AC 611 ms
29,848 KB
testcase_25 AC 5,078 ms
29,896 KB
testcase_26 AC 5,081 ms
29,848 KB
testcase_27 AC 5,079 ms
29,832 KB
testcase_28 AC 5,085 ms
29,836 KB
testcase_29 AC 5,081 ms
29,948 KB
testcase_30 AC 5,070 ms
29,872 KB
testcase_31 AC 5,071 ms
29,880 KB
testcase_32 AC 5,069 ms
29,940 KB
testcase_33 AC 5,074 ms
29,952 KB
testcase_34 AC 5,071 ms
29,728 KB
testcase_35 AC 5,071 ms
29,556 KB
testcase_36 AC 5,071 ms
29,564 KB
testcase_37 AC 5,066 ms
28,880 KB
testcase_38 AC 5,061 ms
29,056 KB
testcase_39 AC 5,070 ms
29,024 KB
testcase_40 AC 5,061 ms
29,088 KB
testcase_41 AC 1 ms
3,448 KB
testcase_42 AC 2 ms
3,508 KB
testcase_43 AC 1 ms
3,560 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using lint = long long;
template<class T = int> using V = vector<T>;
template<class T = int> using VV = V< V<T> >;

using namespace chrono;
auto start_t = steady_clock::now();
void start() { start_t = steady_clock::now(); }
int elapsed_time() { return duration_cast<milliseconds>(steady_clock::now() - start_t).count(); }

struct Edge { int to; lint cost; };

template<class T, class Edge> V<T> dijkstra(const VV<Edge>& g, int s = 0) {
  V<T> dist(g.size(), numeric_limits<T>::max());
  using P = pair<T, int>;
  priority_queue< P, V<P>, greater<P> > pque;
  pque.emplace(dist[s] = 0, s);
  while (!pque.empty()) {
    T d; int v;
    tie(d, v) = pque.top(); pque.pop();
    if (d > dist[v]) continue;
    for (const auto& e : g[v]) if (dist[v] + e.cost < dist[e.to]) {
      pque.emplace(dist[e.to] = dist[v] + e.cost, e.to);
    } 
  }
  return dist;
}

int main() {
  cin.tie(nullptr); ios::sync_with_stdio(false);
  int h, w; cin >> h >> w;
  int gi, gj; cin >> gi >> gj, --gi, --gj;
  auto in = [&](int i, int j) -> int { return i * w + j; };
  auto out = [&](int i, int j) -> int { return (h + i) * w + j; };
  VV<Edge> g(2 * h * w);
  VV<> a(h, V<>(w));
  for (int i = 0; i < h; ++i) for (int j = 0; j < w; ++j) {
    cin >> a[i][j];
    g[in(i, j)].emplace_back(Edge{out(i, j), a[i][j]});
    for (int di = -1; di <= 1; ++di) for (int dj = -1; dj <= 1; ++dj) if (abs(di) + abs(dj) == 1) {
      int ni = i + di, nj = j + dj;
      if (ni < 0 or ni >= h or nj < 0 or nj >= w) continue;
      g[out(i, j)].emplace_back(Edge{in(ni, nj), 0});
    }
  }

  int q; cin >> q;
  struct Q { int t, i, j; lint k; };
  V<Q> qs(q);
  for (int t = 0; t < q; ++t) {
    qs[t].t = t;
    cin >> qs[t].i >> qs[t].j >> qs[t].k, --qs[t].i, --qs[t].j;
  }
  sort(begin(qs), end(qs), [](const auto& l, const auto& r) { return l.k < r.k; });

  V<lint> res(q);
  auto itr = begin(qs);
  while (itr != end(qs)) {
    lint k = itr->k;
    if (elapsed_time() >= 5000) {
      for (int i = 0; i < h; ++i) for (int j = 0; j < w; ++j) {
        g[in(i, j)][0].cost = a[i][j] + k * k;
      }
      auto d = dijkstra<lint>(g, in(gi, gj));
      while (itr != end(qs)) {
        int dist = abs(itr->i - gi) + abs(itr->j - gj) + 1;
        res[itr->t] = d[out(itr->i, itr->j)] + dist * (itr->k * itr->k - k * k);
        ++itr;
      }
      break;
    }

    for (int i = 0; i < h; ++i) for (int j = 0; j < w; ++j) {
      g[in(i, j)][0].cost = a[i][j] + k * k;
    }
    auto d = dijkstra<lint>(g, in(gi, gj));
    while (itr != end(qs) and itr->k == k) {
      res[itr->t] = d[out(itr->i, itr->j)];
      ++itr;
    }
  }
  for (lint e : res) cout << e << '\n';
}
0