#include using namespace std; namespace { typedef double real; typedef long long ll; template ostream& operator<<(ostream& os, const vector& vs) { if (vs.empty()) return os << "[]"; os << vs[0]; for (int i = 1; i < vs.size(); i++) os << " " << vs[i]; return os; } template istream& operator>>(istream& is, vector& vs) { for (auto it = vs.begin(); it != vs.end(); it++) is >> *it; return is; } struct Edge { int from, to, cost; Edge(int from, int to, int cost) : from(from), to(to), cost(cost) {} }; int N, M, S, G; vector> Graph; void input() { cin >> N >> M >> S >> G; Graph.clear(); Graph.resize(N); for (int i = 0; i < M; i++) { int a, b, c; cin >> a >> b >> c; Graph[a].push_back(Edge(a, b, c)); Graph[b].push_back(Edge(b, a, c)); } } typedef pair> C; const int INF = 1<<28; const C INF_C = make_pair(INF, vector(1, 1000)); struct State { int v; C cost; State() {} State(int v, C cost) : v(v), cost(cost) {} }; bool operator<(const State& a, const State& b) { if (a.cost.first == b.cost.first) { a.cost.second > b.cost.second; } else { a.cost.first > b.cost.first; } } struct State0 { int v, cost; State0(int v, int cost) : v(v), cost(cost) {} }; bool operator<(const State0& a, const State0& b) { return a.cost > b.cost; } int calc(vector& D, vector>& P) { priority_queue PQ; D[S] = 0; PQ.push(State0(S, D[S])); while (not PQ.empty()) { State0 s = PQ.top(); PQ.pop(); int v = s.v; int cost = s.cost; for (int i = 0; i < Graph[v].size(); i++) { const Edge& e = Graph[v][i]; if (D[e.to] > cost + e.cost) { D[e.to] = cost + e.cost; PQ.push(State0(e.to, cost + e.cost)); P[e.to].clear(); P[e.to].push_back(e); } if (D[e.to] == cost + e.cost) { P[e.to].push_back(e); } } } } vector D; vector> P; ostream& operator<<(ostream& os, const Edge& e) { return os << "(" << e.from << "->" << e.to << ")"; } vector> memo; vector dfs(int v) { if (v == S) return vector(1, S); if (memo[v][0] >= 0) return memo[v]; vector ret; ret.push_back(INF); for (int i = 0; i < P[v].size(); i++) { // e.to縺後>縺セ縺ソ縺ヲ繧久 // e.from縺後″縺溘∪縺。 const Edge& e = P[v][i]; auto p = dfs(e.from); p.push_back(v); ret = min(ret, p); } return memo[v] = ret; } void solve() { D.clear(); D.resize(N, INF); P.clear(); P.resize(N); calc(D, P); /* for (int i = 0; i < N; i++) { cout << i << " -> " << P[i] << endl; } cout << endl; */ memo.clear(); memo.resize(N, vector(1, -1)); cout << dfs(G) << endl; } } int main() { vector a = {4, 1}; vector b = {4, 3}; input(); solve(); return 0; }