結果
問題 | No.1065 電柱 / Pole (Easy) |
ユーザー |
![]() |
提出日時 | 2020-05-29 21:54:18 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 357 ms / 2,000 ms |
コード長 | 3,227 bytes |
コンパイル時間 | 3,169 ms |
コンパイル使用メモリ | 209,152 KB |
最終ジャッジ日時 | 2025-01-10 17:04:20 |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 46 |
ソースコード
#include<bits/stdc++.h> using namespace std; typedef long long ll; #define F first #define S second #define pii pair<int, int> #define pli pair<ll, int> #define pil pair<int, ll> #define pll pair<ll, ll> #define eb emplace_back #define all(v) v.begin(), v.end() #define rep(i, n) for (int i = 0; i < (n); ++i) #define rep3(i, l, n) for (int i = l; i < (n); ++i) #define sz(v) (int)v.size() #define endl '\n' const int inf = 1000000007; const ll INF = 1e18; int mod = 998244353; // int mod = 1000000007; #define abs(x) (x >= 0 ? x : -(x)) #define lb(v, x) (int)(lower_bound(all(v), x) - v.begin()) #define ub(v, x) (int)(upper_bound(all(v), x) - v.begin()) template<typename T1, typename T2> inline bool chmin(T1 &a, T2 b) { if (a > b) { a = b; return 1; } return 0; } template<typename T1, typename T2> inline bool chmax(T1 &a, T2 b) { if (a < b) { a = b; return 1; } return 0; } ll gcd(ll a, ll b) { if (b == 0) return a; return gcd(b, a % b); } ll lcm(ll a, ll b) { return a / gcd(a, b) * b; } template<typename T, typename U> T pow_(T a, U b) { return b ? pow_(a * a, b / 2) * (b % 2 ? a : 1) : 1; } ll modpow(ll a, ll b, ll _mod) { return b ? modpow(a * a % _mod, b / 2, _mod) * (b % 2 ? a : 1) % _mod : 1; } template<class T, class U> ostream& operator << (ostream& os, const pair<T, U>& p) { os << p.F << " " << p.S; return os; } template<class T> ostream& operator << (ostream& os, const vector<T>& vec) { rep(i, sz(vec)) { if (i) os << " "; os << vec[i]; } return os; } template<typename T> inline istream& operator >> (istream& is, vector<T>& v) { rep(j, sz(v)) is >> v[j]; return is; } template<class T, class T2> inline void add(T &a, T2 b) { a += b; if (a >= mod) a -= mod; } template<class T> void operator += (vector<T>& v, vector<T> v2) { rep(i, sz(v2)) v.eb(v2[i]); } void solve(); int main() { ios::sync_with_stdio(false); cin.tie(0); cout.tie(0); cout << fixed << setprecision(11); int T; T = 1; while (T--) solve(); } void solve() { int n, m; cin >> n >> m; int x, y; cin >> x >> y; x--; y--; vector<int> p(n), q(n); rep(i, n) cin >> p[i] >> q[i]; using pid = pair<int, double>; vector<vector<pid> > G(n); auto d = [&](int a, int b) { return sqrt((p[a] - p[b]) * (p[a] - p[b]) + (q[a] - q[b]) * (q[a] - q[b])); }; rep(i, m) { int a, b; cin >> a >> b; a--; b--; G[a].eb(b, d(a, b)); G[b].eb(a, d(a, b)); } auto dijkstra = [&](int s) { using pdi = pair<double, int>; vector<pdi> d(n, pdi(INF, 0)); // 最短距離, その本数 vector<int> vis(n); priority_queue<pdi, vector<pdi>, greater<pdi> > q; d[s] = pdi(0, 1); q.push(pdi(0, s)); while (!q.empty()) { pdi p = q.top(); q.pop(); int v = p.S; if (vis[v]) continue; vis[v] = 1; for (pid np : G[v]) { int nv; double cost; tie(nv, cost) = np; pdi tmp = pdi(d[v].F + cost, d[v].S); if (d[nv].F != tmp.F) { if (chmin(d[nv], tmp)) { } } else d[nv] = pdi(d[nv].F, (d[nv].S + d[v].S) % mod); // 最短距離の本数増やす q.push(pli(d[nv].F, nv)); } } return d; }; auto di = dijkstra(x); cout << di[y].F << endl; }