//header{{{ #include using namespace std; #define rep(i,n) for(int i=0;i<(n);++i) #define reps(i,n) for(int i=1;i<=(n);++i) #define all(x) (x).begin(),(x).end() #define Fixed fixed << setprecision(12) #define int int64_t using pii = pair; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr int mod = 1e9+7; constexpr int MOD = 998244353; template inline bool chmax(A &a, const B &b) { return b > a && (a = b, true); } template inline bool chmin(A &a, const B &b) { return b < a && (a = b, true); } template using min_heap = priority_queue,greater >; template using max_heap = priority_queue; template using umap = unordered_map; int gcd(int a,int b){ return b ? gcd(b,a % b) : a;} int lcm(int a,int b){ return a / gcd(a,b) * b;} constexpr int dx[] = {1,0,-1,0,1,1,-1,-1}; constexpr int dy[] = {0,-1,0,1,1,-1,-1,1}; //}}} signed main(){ cin.tie(0); ios::sync_with_stdio(false); int n,m,S,G; cin >> n >> m >> S >> G; vector > g(n); rep(i,m){ int a,b,c; cin >> a >> b >> c; g[a].emplace_back(b,c); g[b].emplace_back(a,c); } vector min_cost(n,LINF),res(n,99999); min_heap pq; res[S] = S; min_cost[S] = 0; pq.emplace(0,S); while(!pq.empty()){ int now,cost; tie(cost,now) = pq.top(); pq.pop(); if(min_cost[now] < cost) continue; rep(i,g[now].size()){ int next,ncost; tie(next,ncost) = g[now][i]; if(chmin(min_cost[next],min_cost[now] + ncost)){ pq.emplace(min_cost[next],next); res[next] = now; }else if(min_cost[next] == min_cost[now] + ncost && chmin(res[next],now)){ pq.emplace(min_cost[next],next); } } } int now = G; vector ans; ans.emplace_back(G); for(int i = 0;res[now] != now;i++){ ans.emplace_back(res[now]); now = res[now]; } reverse(all(ans)); rep(i,ans.size()){ cout << (i == 0 ? "" : " ") << ans[i]; } cout << '\n'; return 0; }