#include #pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") #pragma GCC target("sse,sse2,sse3,ssse3,sse4,fma,abm,mmx,avx,avx2") #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define rrep(i, n) for (int i = (int)(n - 1); i >= 0; i--) #define all(x) (x).begin(), (x).end() #define sz(x) int(x.size()) #define yn(joken) cout<<((joken) ? "Yes" : "No")<<"\n" #define YN(joken) cout<<((joken) ? "YES" : "NO")<<"\n" using namespace std; using ll = long long; using vi = vector; using vl = vector; using vs = vector; using vc = vector; using vd = vector; using vld = vector; using vvi = vector>; using vvl = vector>; using vvs = vector>; using vvc = vector>; using vvd = vector>; using vvld = vector>; using vvvi = vector>>; using vvvl = vector>>; using vvvvi = vector>>>; using vvvvl = vector>>>; using pii = pair; using pll = pair; const int INF = 1e9; const ll LINF = 2e18; template bool chmax(T& a, const T& b) { if (a < b) { a = b; return 1; } return 0; } template bool chmin(T& a, const T& b) { if (b < a) { a = b; return 1; } return 0; } bool ispow2(int i) { return i && (i & -i) == i; } bool ispow2(ll i) { return i && (i & -i) == i; } template vector make_vec(size_t a) { return vector(a); } template auto make_vec(size_t a, Ts... ts) { return vector(ts...))>(a, make_vec(ts...)); } template istream& operator>>(istream& is, vector& v) { for (int i = 0; i < int(v.size()); i++) { is >> v[i]; } return is; } template ostream& operator<<(ostream& os, const vector& v) { for (int i = 0; i < int(v.size()); i++) { os << v[i]; if (i < int(v.size()) - 1) os << ' '; } return os; } static uint32_t RandXor(){ static uint32_t x=123456789; static uint32_t y=362436069; static uint32_t z=521288629; static uint32_t w=88675123; uint32_t t; t=x^(x<<11); x=y; y=z; z=w; return w=(w^(w>>19))^(t^(t>>8)); } static double Rand01(){ return (RandXor()+0.5)*(1.0/UINT_MAX); } template struct Edge{ int from, to; T cost; int idx; Edge() = default; Edge(int from, int to, T cost = 1, int idx = -1) : from(from), to(to), cost(cost), idx(idx) {} operator int() const { return to; } }; template struct Graph{ vector>> g; int es; Graph() = default; explicit Graph(int n) : g(n), es(0) {} size_t size() const{ return g.size(); } void add_directed_edge(int from, int to, T cost = 1){ g[from].emplace_back(from, to, cost, es++); } void add_edge(int from, int to, T cost = 1){ g[from].emplace_back(from, to, cost, es); g[to].emplace_back(to, from, cost, es++); } void read(int M, int padding = -1, bool weighted = false, bool directed = false){ for (int i = 0; i < M; i++){ int a, b; cin >> a >> b; a += padding; b += padding; T c = T(1); if (weighted) cin >> c; if (directed) add_directed_edge(a, b, c); else add_edge(a, b, c); } } inline vector> &operator[](const int &k){ return g[k]; } inline const vector> &operator[](const int &k) const{ return g[k]; } }; template using Edges = vector>; // dijkstra(g,start) とする. 返り値は以下の3つ. // startからの最短距離の配列 dist // 最短経路でその頂点の前に通る頂点の配列 from (startおよび到達不能頂点では-1) // 最短経路でその頂点の前に通る辺の辺番号の配列 idx (startおよび到達不能頂点では-1) template struct ShortestPath{ vector dist; vector from, id; }; template ShortestPath dijkstra(const Graph &g, int s){ const auto INF = numeric_limits::max(); vector dist(g.size(), INF); vector from(g.size(), -1), id(g.size(), -1); using Pi = pair; priority_queue, greater<>> que; dist[s] = 0; que.emplace(dist[s], s); while (!que.empty()){ T cost; int idx; tie(cost, idx) = que.top(); que.pop(); if (dist[idx] < cost) continue; for (auto &e : g[idx]){ auto next_cost = cost + e.cost; if (dist[e.to] <= next_cost) continue; dist[e.to] = next_cost; from[e.to] = idx; id[e.to] = e.idx; que.emplace(dist[e.to], e.to); } } return {dist, from, id}; } void solve(){ using T=tuple; int N,V,sy,sx,gy,gx; cin>>N>>V>>sy>>sx>>gy>>gx; sx--; sy--; gx--; gy--; vvi F(N,vi(N)); rep(i,N) cin>>F[i]; vvi MD(N,vi(N,INF)),MV(N,vi(N,-1)); vi dx={0,-1,0,1},dy={-1,0,1,0}; auto inside=[&](int x,int y){ return 0<=x && x G(N*N); rep(i,N){ rep(j,N){ rep(d,4){ int nx=i+dx[d],ny=j+dy[d]; if(!inside(nx,ny)) continue; G.add_directed_edge(i*N+j,nx*N+ny,F[nx][ny]); } } } auto ret=dijkstra(G,sx*N+sy).dist; if(ret[gx*N+gy]>=V){ cout<<-1<,greater> pq; // d,v,x,y pq.emplace(0,-V,sx,sy); MD[sx][sy]=0; MV[sx][sy]=V; while(sz(pq)){ auto [d,v,x,y]=pq.top(); pq.pop(); v=-v; if(x==gx && y==gy){ cout<MD[x][y] && vMD[nx][ny]){ if(v-F[nx][ny]<=MV[nx][ny]) continue; MD[nx][ny]=d+1; MV[nx][ny]=v-F[nx][ny]; pq.emplace(d+1,-v+F[nx][ny],nx,ny); } else if(d+1==MD[nx][ny]){ if(v-F[nx][ny]>MV[nx][ny]){ MV[nx][ny]=v-F[nx][ny]; pq.emplace(d+1,-v+F[nx][ny],nx,ny); } } if(d+1