結果
問題 | No.2354 Poor Sight in Winter |
ユーザー |
![]() |
提出日時 | 2023-06-16 22:20:26 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 199 ms / 2,000 ms |
コード長 | 2,653 bytes |
コンパイル時間 | 1,381 ms |
コンパイル使用メモリ | 123,996 KB |
最終ジャッジ日時 | 2025-02-14 06:19:42 |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <iostream> #include <algorithm> #include <iomanip> #include <vector> #include <queue> #include <deque> #include <set> #include <map> #include <tuple> #include <cmath> #include <numeric> #include <functional> #include <cassert> #define debug_value(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << #x << "=" << x << endl; #define debug(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << x << endl; template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } using namespace std; typedef long long ll; template<typename T> vector<vector<T>> vec2d(int n, int m, T v){ return vector<vector<T>>(n, vector<T>(m, v)); } template<typename T> vector<vector<vector<T>>> vec3d(int n, int m, int k, T v){ return vector<vector<vector<T>>>(n, vector<vector<T>>(m, vector<T>(k, v))); } template<typename T> void print_vector(vector<T> v, char delimiter=' '){ if(v.empty()) { cout << endl; return; } for(int i = 0; i+1 < v.size(); i++) cout << v[i] << delimiter; cout << v.back() << endl; } // calc ceil(a/b) template<typename T> T ceil_div(T a, T b){ return (a+b-1)/b; } const ll inf = 1e18; using P = pair<ll, ll>; int main(){ ios::sync_with_stdio(false); cin.tie(0); cout << setprecision(10) << fixed; int n, k; cin >> n >> k; ll sx, sy, gx, gy; cin >> sx >> sy >> gx >> gy; vector<P> ps = {{sx, sy}}; for(int i = 0; i < n; i++){ ll x, y; cin >> x >> y; ps.push_back({x, y}); } ps.push_back({gx, gy}); n = ps.size(); auto d = vec2d(n, n, 0ll); for(int i = 0; i < n; i++){ for(int j = 0; j < n; j++){ auto [xi, yi] = ps[i]; auto [xj, yj] = ps[j]; d[i][j] = abs(xi-xj)+abs(yi-yj); } } auto f = [&](ll x){ vector<ll> dist(n, inf); priority_queue<P, vector<P>, greater<P>> que; auto set_dist = [&](int v, ll d){ if(chmin(dist[v], d)) que.push({d, v}); }; set_dist(0, 0); while(!que.empty()){ auto [dd, v] = que.top(); que.pop(); if(dd != dist[v]) continue; for(int i = 0; i < n; i++){ if(i == v) continue; ll cnt = ceil_div(d[v][i], x)-1; set_dist(i, dd+cnt); } } return dist[n-1] <= k; }; ll l = 0, r = inf; while(r-l > 1){ ll x = (l+r)/2; if(f(x)) { r = x; } else l = x; } cout << r << endl; }