#{{{ header import algorithm, sequtils, tables, macros, math, sets, strutils when defined(MYDEBUG): import header proc scanf(formatstr: cstring){.header: "", varargs.} proc getchar(): char {.header: "", varargs.} proc nextInt(): int = scanf("%lld",addr result) proc nextFloat(): float = scanf("%lf",addr result) proc nextString(): string = var get = false result = "" while true: var c = getchar() if int(c) > int(' '): get = true result.add(c) else: if get: break get = false template `max=`*(x,y:typed):void = x = max(x,y) template `min=`*(x,y:typed):void = x = min(x,y) template infty(T): untyped = ((T(1) shl T(sizeof(T)*8-2)) - 1) #}}} #{{{ Graph type Edge = object src,dst:int weight:int rev:int proc newEdge(src,dst:int,weight:int = 1,rev:int = -1):Edge = var e:Edge e.src = src e.dst = dst e.weight = weight e.rev = rev return e proc `<`(l,r:Edge):bool = return l.weight < r.weight type Graph = seq[seq[Edge]] proc newGraph(n:int):Graph= var g:Graph g = newSeqWith(n,newSeq[Edge]()) return g proc addBiEdge(g:var Graph,src,dst:int,weight:int=1):void = g[src].add(newEdge(src,dst,weight,g[dst].len)) g[dst].add(newEdge(dst,src,weight,g[src].len-1)) proc addEdge(g:var Graph,src,dst:int,weight:int=1):void = g[src].add(newEdge(src,dst,weight,-1)) #}}} proc cmp2(l,r:Edge):int = if l.dst < r.dst : -1 else: 1 import heapqueue proc shortest_path(g:Graph, s:int): (seq[int],seq[int]) = var n = g.len dist = newSeqWith(n,int.infty) prev = newSeqWith(n,-1) Q = initHeapQueue[Edge]() dist[s] = 0 Q.push(newEdge(-2,s,0)) while Q.len > 0: var e = Q.pop() if prev[e.dst] != -1: continue prev[e.dst] = e.src; for f in g[e.dst]: var w = e.weight + f.weight; if dist[f.dst] > w: dist[f.dst] = w; Q.push(newEdge(f.src, f.dst, w)); discard return (dist,prev) #//{{{ shortest_path(graph g, int s) : Dijkstra #template # const graph &g; # const int &s; # vector dist; # vector prev; #// dijkstra_t(const graph &g,int s):g(g),s(s),dist(g.size(),g.inf),prev(g.size(),-1){} # dijkstra_t(const graph &g,const int &s):g(g),s(s),dist(g.size(),g.inf),prev(g.size(),-1){ # /* # assert(dist.size()==g.size()); # REP(i,dist.size())assert(dist[i]>=0 and dist[i]+dist[i]>=0); # */ # } # vector path(int t) { # vector path; # for (int u = t; u >= 0; u = prev[u])path.push_back(u); # reverse(ALL(path)); # return path; # } # void dijkstra(){ # dist[s] = 0; # priority_queue > Q; // "e < f" <=> "e.weight > f.weight" # for (Q.push(edge(-2, s, 0)); !Q.empty(); ) { # edge e = Q.top(); Q.pop(); # if (prev[e.dst] != -1) continue; # prev[e.dst] = e.src; # for(auto f:g[e.dst]) { #// assert(e.weight<=g.inf); #// assert(0<=f.weight and f.weight<=g.inf); # Weight w = e.weight + f.weight; # if (dist[f.dst] > w) { # dist[f.dst] = w; # Q.push(edge(f.src, f.dst, w)); # } # } # } # } # /* # void count_dp(){ # vector > d; # REP(u,g.size())d.push_back({dist[u],u}); # sort(ALL(d)); # dp.assign(g.size(),mod_int()); # dp[s] = 1; # REP(i,d.size()){ # int u = d[i].second; # for(auto &&e:g[u]){ # if(dist[e.dst]+e.weight==dist[u]){ # dp[u]+=dp[e.dst]; # } # } # } # } # */ #}; #template #dijkstra_t shortest_path(const graph &g, int s) { # dijkstra_t ret(g,s); # ret.dijkstra(); #// ret.count_dp(); # return ret; #} #//}}} proc main():void = var N = nextInt() M = nextInt() S = nextInt() G = nextInt() a = newSeq[int](M) b = newSeq[int](M) c = newSeq[int](M) g = newGraph(N) for i in 0..M-1: a[i] = nextInt();b[i] = nextInt();c[i] = nextInt() g.addBiEdge(a[i],b[i],c[i]) for u in 0..N-1: g[u].sort(cmp2) var (dist,_) = shortest_path(g,G) var cur = S ans = newSeq[int]() while cur != G: ans.add(cur) for e in g[cur]: if dist[cur] == dist[e.dst] + e.weight: cur = e.dst break ans.add(cur) echo ans.join(" ") discard main()