結果
問題 |
No.160 最短経路のうち辞書順最小
|
ユーザー |
|
提出日時 | 2019-03-22 05:30:20 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 28 ms / 5,000 ms |
コード長 | 11,799 bytes |
コンパイル時間 | 2,061 ms |
コンパイル使用メモリ | 207,516 KB |
最終ジャッジ日時 | 2025-01-06 23:39:34 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 26 |
ソースコード
/** * */ // {{{ header #include <bits/stdc++.h> using namespace std; using i8 = int8_t; using u8 = uint8_t; using i16 = int16_t; using u16 = uint16_t; using i32 = int32_t; using u32 = uint32_t; using i64 = int64_t; using u64 = uint64_t; using i128 = __int128_t; using u128 = __uint128_t; using f32 = float; using f64 = double; using f80 = __float80; using f128 = __float128; // }}} constexpr i64 INF = 1'010'000'000'000'000'000LL; constexpr i64 MOD = 1'000'000'007LL; constexpr f64 EPS = 1e-12; constexpr f64 PI = 3.14159265358979323846; // {{{ util #define FOR(i, start, end) for(i64 i = (start), i##_end=(end); i < i##_end; ++i) #define REP(i, n) FOR(i, 0, n) #define ALL(f,c,...) (([&](decltype((c)) cccc) { return (f)(std::begin(cccc), std::end(cccc), ## __VA_ARGS__); })(c)) #define SLICE(f,c,l,r,...) (([&](decltype((c)) cccc, decltype((l)) llll, decltype((r)) rrrr) {\ auto iiii = llll <= rrrr ? std::begin(cccc)+llll : std::end(cccc);\ auto jjjj = llll <= rrrr ? std::begin(cccc)+rrrr : std::end(cccc);\ return (f)(iiii, jjjj, ## __VA_ARGS__);\ })(c,l,r)) #define GENERIC(f) ([](auto&&... args) -> decltype(auto) { return (f)(std::forward<decltype(args)>(args)...); }) template<typename F> class FixPoint { public: explicit constexpr FixPoint(F&& f) : f_(forward<F>(f)) {} template<typename... Args> constexpr decltype(auto) operator()(Args&&... args) const { return f_(*this, forward<Args>(args)...); } private: const F f_; }; template<typename F> decltype(auto) FIX(F&& f) { return FixPoint<F>(forward<F>(f)); } template<typename C> i64 SIZE(const C& c) { return static_cast<i64>(c.size()); } template<typename T, size_t N> i64 SIZE(const T (&)[N]) { return static_cast<i64>(N); } bool is_odd (i64 x) { return x % 2 != 0; } bool is_even(i64 x) { return x % 2 == 0; } template<typename T> i64 cmp(T x, T y) { return (y<x) - (x<y); } template<typename T> i64 sgn(T x) { return cmp(x, T(0)); } // Haskell の divMod と同じ pair<i64,i64> divmod(i64 a, i64 b) { i64 q = a / b; i64 r = a % b; if((b>0 && r<0) || (b<0 && r>0)) { --q; r += b; } return {q,r}; } i64 div_ceil(i64 a, i64 b) { i64 q = a / b; i64 r = a % b; if((b>0 && r>0) || (b<0 && r<0)) ++q; return q; } i64 div_floor(i64 a, i64 b) { return divmod(a,b).first; } i64 modulo(i64 a, i64 b) { return divmod(a,b).second; } bool feq(f64 x, f64 y, f64 eps=EPS) { return fabs(x-y) < eps; } template<typename T, typename U> bool chmax(T& xmax, const U& x) { if(xmax < x) { xmax = x; return true; } return false; } template<typename T, typename U> bool chmin(T& xmin, const U& x) { if(x < xmin) { xmin = x; return true; } return false; } template<typename InputIt> auto SUM(InputIt first, InputIt last) { using T = typename iterator_traits<InputIt>::value_type; return accumulate(first, last, T()); } template<typename ForwardIt, typename UnaryOperation> ForwardIt transform_self(ForwardIt first, ForwardIt last, UnaryOperation op) { return transform(first, last, first, op); } template<typename C> void UNIQ(C& c) { c.erase(ALL(unique,c), end(c)); } template<typename BinaryFunc, typename UnaryFunc> auto ON(BinaryFunc bf, UnaryFunc uf) { return [bf,uf](const auto& x, const auto& y) { return bf(uf(x), uf(y)); }; } template<typename F> auto LT_ON(F f) { return ON(less<>(), f); } template<typename F> auto GT_ON(F f) { return ON(greater<>(), f); } char digit_chr(i64 n) { return static_cast<char>('0' + n); } i64 digit_ord(char c) { return c - '0'; } char lower_chr(i64 n) { return static_cast<char>('a' + n); } i64 lower_ord(char c) { return c - 'a'; } char upper_chr(i64 n) { return static_cast<char>('A' + n); } i64 upper_ord(char c) { return c - 'A'; } // 出力は operator<< を直接使わず、このテンプレート経由で行う // 提出用出力とデバッグ用出力を分けるため template<typename T> struct Formatter { static ostream& write_str(ostream& out, const T& x) { return out << x; } static ostream& write_repr(ostream& out, const T& x) { return out << x; } }; template<typename T> ostream& WRITE_STR(ostream& out, const T& x) { return Formatter<T>::write_str(out, x); } template<typename T> ostream& WRITE_REPR(ostream& out, const T& x) { return Formatter<T>::write_repr(out, x); } template<typename InputIt> ostream& WRITE_JOIN_STR(ostream& out, InputIt first, InputIt last, const string& sep) { while(first != last) { WRITE_STR(out, *first++); if(first != last) out << sep; } return out; } template<typename InputIt> ostream& WRITE_JOIN_REPR(ostream& out, InputIt first, InputIt last, const string& sep) { while(first != last) { WRITE_REPR(out, *first++); if(first != last) out << sep; } return out; } template<typename InputIt> ostream& WRITE_RANGE_STR(ostream& out, InputIt first, InputIt last) { return WRITE_JOIN_STR(out, first, last, " "); } template<typename InputIt> ostream& WRITE_RANGE_REPR(ostream& out, InputIt first, InputIt last) { out << "["; WRITE_JOIN_REPR(out, first, last, ", "); out << "]"; return out; } template<typename T> void FROM_STR(const string& s, T& x) { istringstream in(s); in >> x; } template<typename T> string TO_STR(const T& x) { ostringstream out; WRITE_STR(out, x); return out.str(); } template<typename T> string TO_REPR(const T& x) { ostringstream out; WRITE_REPR(out, x); return out.str(); } template<typename InputIt> string RANGE_TO_STR(InputIt first, InputIt last) { ostringstream out; WRITE_RANGE_STR(out, first, last); return out.str(); } template<typename InputIt> string RANGE_TO_REPR(InputIt first, InputIt last) { ostringstream out; WRITE_RANGE_REPR(out, first, last); return out.str(); } template<typename InputIt> string JOIN(InputIt first, InputIt last, const string& sep) { ostringstream out; WRITE_JOIN_STR(out, first, last, sep); return out.str(); } template<typename T> struct Formatter<vector<T>> { static ostream& write_str(ostream& out, const vector<T>& v) { return WRITE_RANGE_STR(out, begin(v), end(v)); } static ostream& write_repr(ostream& out, const vector<T>& v) { out << "vector"; return WRITE_RANGE_REPR(out, begin(v), end(v)); } }; template<typename T1, typename T2> struct Formatter<pair<T1,T2>> { static ostream& write_str(ostream& out, const pair<T1,T2>& p) { WRITE_STR(out, p.first); out << ' '; WRITE_STR(out, p.second); return out; } static ostream& write_repr(ostream& out, const pair<T1,T2>& p) { out << "("; WRITE_REPR(out, p.first); out << ","; WRITE_REPR(out, p.second); out << ")"; return out; } }; template<typename... TS> struct Formatter<tuple<TS...>> { template<size_t I=0, enable_if_t<I == sizeof...(TS), nullptr_t> = nullptr> static ostream& write_str_impl(ostream& out, const tuple<TS...>&) { return out; } template<size_t I=0, enable_if_t<I < sizeof...(TS), nullptr_t> = nullptr> static ostream& write_str_impl(ostream& out, const tuple<TS...>& t) { if(I != 0) out << ' '; WRITE_STR(out, get<I>(t)); return write_str_impl<I+1>(out, t); } template<size_t I=0, enable_if_t<I == sizeof...(TS), nullptr_t> = nullptr> static ostream& write_repr_impl(ostream& out, const tuple<TS...>&) { if(sizeof...(TS) == 0) out << "("; return out << ")"; } template<size_t I=0, enable_if_t<I < sizeof...(TS), nullptr_t> = nullptr> static ostream& write_repr_impl(ostream& out, const tuple<TS...>& t) { if(I == 0) out << "("; else out << ","; WRITE_REPR(out, get<I>(t)); return write_repr_impl<I+1>(out, t); } static ostream& write_str(ostream& out, const tuple<TS...>& t) { return write_str_impl(out, t); } static ostream& write_repr(ostream& out, const tuple<TS...>& t) { return write_repr_impl(out, t); } }; template<typename T> void RD(T& x) { cin >> x; #ifdef PROCON_LOCAL assert(cin); #endif } template<typename T> void RD(vector<T>& v, i64 n) { v.reserve(n); REP(_, n) { T e; RD(e); v.emplace_back(e); } } void PRINT() {} template<typename T, typename... TS> void PRINT(const T& x, const TS& ...args) { WRITE_STR(cout, x); if(sizeof...(args)) { cout << ' '; PRINT(args...); } } template<typename... TS> void PRINTLN(const TS& ...args) { PRINT(args...); cout << '\n'; } [[noreturn]] void EXIT() { #ifdef PROCON_LOCAL cerr.flush(); #endif cout.flush(); _Exit(0); } template<typename T> void DBG_IMPL(i64 line, const char* expr, const T& value) { #ifdef PROCON_LOCAL cerr << "[L " << line << "]: "; cerr << expr << " = "; WRITE_REPR(cerr, value); cerr << "\n"; #endif } template<typename T, size_t N> void DBG_ARRAY_IMPL(i64 line, const char* expr, const T (&ary)[N]) { #ifdef PROCON_LOCAL cerr << "[L " << line << "]: "; cerr << expr << " = "; WRITE_RANGE_REPR(cerr, begin(ary), end(ary)); cerr << "\n"; #endif } template<typename InputIt> void DBG_RANGE_IMPL(i64 line, const char* expr1, const char* expr2, InputIt first, InputIt last) { #ifdef PROCON_LOCAL cerr << "[L " << line << "]: "; cerr << expr1 << "," << expr2 << " = "; WRITE_RANGE_REPR(cerr, first, last); cerr << "\n"; #endif } #define DBG(expr) DBG_IMPL(__LINE__, #expr, (expr)) #define DBG_ARRAY(expr) DBG_ARRAY_IMPL(__LINE__, #expr, (expr)) #define DBG_RANGE(first,last) DBG_RANGE_IMPL(__LINE__, #first, #last, (first), (last)) // }}} // {{{ init struct ProconInit { static constexpr int IOS_PREC = 15; static constexpr bool AUTOFLUSH = false; ProconInit() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(IOS_PREC); #ifdef PROCON_LOCAL cerr << fixed << setprecision(IOS_PREC); #endif if(AUTOFLUSH) cout << unitbuf; } } PROCON_INIT; // }}} //-------------------------------------------------------------------- pair<bool,vector<vector<i64>>> graph_floyd(vector<vector<i64>>& g) { i64 n = SIZE(g); vector<vector<i64>> nex(n, vector<i64>(n,-1)); REP(i, n) REP(j, n) { if(g[i][j] != INF) nex[i][j] = j; } REP(k, n) { REP(i, n) { if(g[i][k] == INF) continue; REP(j, n) { if(g[k][j] == INF) continue; if(chmin(g[i][j], g[i][k] + g[k][j])) { nex[i][j] = nex[i][k]; } else if(i != k && g[i][j] == g[i][k] + g[k][j]) { chmin(nex[i][j], nex[i][k]); } if(i == j && g[i][j] < 0) return { false, nex }; } } } return { true, nex }; } void solve() { i64 N; RD(N); i64 M; RD(M); i64 start; RD(start); i64 goal; RD(goal); vector<vector<i64>> G(N, vector<i64>(N, INF)); REP(i, N) { G[i][i] = 0; } REP(_, M) { i64 from,to; RD(from); RD(to); i64 cost; RD(cost); G[from][to] = cost; G[to][from] = cost; } bool ok; vector<vector<i64>> nex; tie(ok,nex) = graph_floyd(G); assert(ok); vector<i64> route; for(i64 v = start; v != goal; v = nex[v][goal]) { route.emplace_back(v); } route.emplace_back(goal); PRINTLN(route); } signed main() { solve(); EXIT(); }