結果

問題 No.1599 Hikyaku
ユーザー 👑 NachiaNachia
提出日時 2021-07-10 15:20:59
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,874 bytes
コンパイル時間 950 ms
コンパイル使用メモリ 90,584 KB
実行使用メモリ 65,412 KB
最終ジャッジ日時 2023-09-14 19:32:59
合計ジャッジ時間 41,509 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 209 ms
62,664 KB
testcase_01 AC 253 ms
63,000 KB
testcase_02 AC 258 ms
62,888 KB
testcase_03 AC 261 ms
62,820 KB
testcase_04 AC 1,147 ms
64,664 KB
testcase_05 AC 524 ms
63,044 KB
testcase_06 AC 1,157 ms
64,652 KB
testcase_07 AC 612 ms
64,692 KB
testcase_08 AC 1,239 ms
64,884 KB
testcase_09 AC 519 ms
62,816 KB
testcase_10 AC 549 ms
62,748 KB
testcase_11 AC 431 ms
62,692 KB
testcase_12 AC 519 ms
64,320 KB
testcase_13 AC 372 ms
64,380 KB
testcase_14 AC 210 ms
62,680 KB
testcase_15 AC 230 ms
62,680 KB
testcase_16 AC 513 ms
62,660 KB
testcase_17 WA -
testcase_18 AC 1,125 ms
62,988 KB
testcase_19 AC 1,265 ms
64,936 KB
testcase_20 AC 1,258 ms
64,672 KB
testcase_21 AC 1,273 ms
64,664 KB
testcase_22 AC 1,287 ms
64,640 KB
testcase_23 AC 1,252 ms
64,616 KB
testcase_24 AC 1,248 ms
64,560 KB
testcase_25 AC 1,277 ms
65,412 KB
testcase_26 AC 1,115 ms
64,264 KB
testcase_27 AC 791 ms
64,200 KB
testcase_28 AC 762 ms
63,416 KB
testcase_29 AC 576 ms
63,544 KB
testcase_30 AC 353 ms
62,784 KB
testcase_31 AC 213 ms
62,680 KB
testcase_32 AC 211 ms
62,732 KB
testcase_33 AC 223 ms
62,700 KB
testcase_34 AC 708 ms
62,848 KB
testcase_35 AC 973 ms
63,528 KB
testcase_36 AC 1,236 ms
63,456 KB
testcase_37 AC 1,211 ms
63,364 KB
testcase_38 AC 1,210 ms
63,448 KB
testcase_39 AC 1,223 ms
63,436 KB
testcase_40 AC 1,232 ms
63,712 KB
testcase_41 AC 166 ms
62,644 KB
testcase_42 AC 213 ms
62,724 KB
testcase_43 AC 165 ms
62,644 KB
testcase_44 WA -
testcase_45 AC 536 ms
62,756 KB
testcase_46 WA -
testcase_47 WA -
testcase_48 WA -
testcase_49 AC 211 ms
62,752 KB
testcase_50 AC 370 ms
62,716 KB
testcase_51 AC 529 ms
62,720 KB
testcase_52 AC 566 ms
63,596 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <queue>
using namespace std;
using ll = long long;
using ull = unsigned long long;
#define rep(i,n) for(int i=0; i<(n); i++)

template<class T>
using nega_queue = priority_queue<T,vector<T>,greater<T>>;

const int INF = 1001001001;

const int W = 600;
const int H = 600;
const int NV = W * H; // number of vertices
const int maxV = 7;

struct Start{ int id,v; };

int pos2id(int y,int x){ return y * W + x; }

int X,Y;
int N;
vector<Start> hikyaku;
double dt[maxV];
double V[maxV];
double T[maxV][NV];
vector<int> E[NV];

double D[maxV][NV];

int main(){
  cin >> X >> Y; X--; Y--;
  cin >> N;

  rep(y,H) rep(x,W-1) E[pos2id(y,x)].push_back(pos2id(y,x+1));
  rep(y,H) rep(x,W-1) E[pos2id(y,x+1)].push_back(pos2id(y,x));
  rep(y,H-1) rep(x,W) E[pos2id(y,x)].push_back(pos2id(y+1,x));
  rep(y,H-1) rep(x,W) E[pos2id(y+1,x)].push_back(pos2id(y,x));

  hikyaku.resize(N);
  rep(i,N){
    int x,y,v; cin >> x >> y >> v;
    x--; y--; v--;
    hikyaku[i] = { pos2id(y,x),v };
  }

  rep(v,maxV) V[v] = 1+v;
  rep(v,maxV) dt[v] = 1000.0 / (1+v);
  
  rep(v,maxV) rep(i,NV) T[v][i] = INF;

  {
    vector<vector<int>> I(maxV);
    for(auto a : hikyaku){
      I[a.v].push_back(a.id);
      T[a.v][a.id] = 0;
    }
    rep(v,maxV){
      rep(i,I[v].size()){
        int p = I[v][i];
        for(int e : E[p]) if(T[v][e] >= INF - 1.0){
          T[v][e] = T[v][p] + dt[v];
          I[v].push_back(e);
        }
      }
    }
  }
  
  rep(v,maxV) rep(i,NV) D[v][i] = INF;

  nega_queue<pair<double,pair<int,int>>> Q;
  D[hikyaku[0].v][hikyaku[0].id] = 0.0;
  Q.push(make_pair(0.0,make_pair(hikyaku[0].v,hikyaku[0].id)));

  while(Q.size()){
    double d = Q.top().first;
    int v = Q.top().second.first;
    int p = Q.top().second.second;
    Q.pop();
    if(D[v][p] != d) continue;
    for(int nxv=v+1; nxv<maxV; nxv++) if(D[nxv][p] > max(T[nxv][p],d)){
      D[nxv][p] = max(T[nxv][p],d);
      Q.push(make_pair(D[nxv][p],make_pair(nxv,p)));
    }
    for(int e : E[p]){
      if(D[v][e] > D[v][p] + dt[v]){
        D[v][e] = D[v][p] + dt[v];
        Q.push(make_pair(D[v][e],make_pair(v,e)));
      }
      for(int nxv=v+1; nxv<maxV; nxv++){
        double nxd = d + dt[v];
        if(d < T[nxv][p] && T[nxv][e] <= d + dt[v]){
          double t0 = T[nxv][e] - D[v][p];
          double t = (1000.0+V[nxv]*t0) / (V[v]+V[nxv]);
          double cand = d + t * 2 - t0;
          nxd = min(nxd,max(d+dt[nxv],cand));
          if(D[nxv][e] <= nxd) continue;
          D[nxv][e] = nxd;
          Q.push(make_pair(nxd,make_pair(nxv,e)));
        }
      }
    }
  }

  double ans = INF;
  rep(v,maxV) ans = min(ans,D[v][pos2id(Y,X)]);
  cout.precision(10); cout << fixed;
  cout << ans << endl;
  return 0;
}

struct ios_do_not_sync{
  ios_do_not_sync(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
  }
} ios_do_not_sync_inst;
0