#include using namespace std; using int128 = __int128_t; using int64 = long long; using int32 = int; using uint128 = __uint128_t; using uint64 = unsigned long long; using uint32 = unsigned int; #define ALL(obj) (obj).begin(),(obj).end() template using priority_queue_reverse = priority_queue,greater>; constexpr int64 MOD = 1'000'000'000LL + 7; //' constexpr int64 MOD2 = 998244353; constexpr int64 HIGHINF = 1'000'000'000'000'000'000LL; constexpr int64 LOWINF = 1'000'000'000'000'000LL; //' constexpr long double PI = 3.1415926535897932384626433L; template vector multivector(size_t N,T init){return vector(N,init);} template auto multivector(size_t N,T... t){return vector(N,multivector(t...));} template void corner(bool flg, T hoge) {if (flg) {cout << hoge << endl; exit(0);}} template ostream &operator<<(ostream &o, const map&obj) {o << "{"; for (auto &x : obj) o << " {" << x.first << " : " << x.second << "}" << ","; o << " }"; return o;} template ostream &operator<<(ostream &o, const set&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;} template ostream &operator<<(ostream &o, const multiset&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;} template ostream &operator<<(ostream &o, const vector&obj) {o << "{"; for (int i = 0; i < (int)obj.size(); ++i)o << (i > 0 ? ", " : "") << obj[i]; o << "}"; return o;} template ostream &operator<<(ostream &o, const deque&obj) {o << "{"; for (int i = 0; i < (int)obj.size(); ++i)o << (i > 0 ? ", " : "") << obj[i]; o << "}"; return o;} template ostream &operator<<(ostream &o, const pair&obj) {o << "{" << obj.first << ", " << obj.second << "}"; return o;} void print(void) {cout << endl;} template void print(Head&& head) {cout << head;print();} template void print(Head&& head, Tail&&... tail) {cout << head << " ";print(forward(tail)...);} template void chmax(T& a, const T b){a=max(a,b);} template void chmin(T& a, const T b){a=min(a,b);} vector split(const string &str, const char delemiter) {vector res;stringstream ss(str);string buffer; while( getline(ss, buffer, delemiter) ) res.push_back(buffer); return res;} inline constexpr int msb(int x) {return x?31-__builtin_clz(x):-1;} inline constexpr int64 ceil_div(const int64 a,const int64 b) {return (a+(b-1))/b;}// return ceil(a/b) void YN(bool flg) {cout << (flg ? "YES" : "NO") << endl;} void Yn(bool flg) {cout << (flg ? "Yes" : "No") << endl;} void yn(bool flg) {cout << (flg ? "yes" : "no") << endl;} /* * @title Graph * @docs md/graph/Graph.md */ template class Graph{ private: const size_t N,H,W; public: vector>> edges; Graph(const size_t N):H(-1),W(-1),N(N), edges(N) {} Graph(const size_t H, const size_t W):H(H),W(W),N(H*W), edges(H*W) {} inline void make_edge(size_t from, size_t to, T w) { edges[from].emplace_back(to,w); } //{from_y,from_x} -> {to_y,to_x} inline void make_edge(pair from, pair to, T w) { make_edge(from.first*W+from.second,to.first*W+to.second,w); } inline void make_bidirectional_edge(size_t from, size_t to, T w) { make_edge(from,to,w); make_edge(to,from,w); } inline void make_bidirectional_edge(pair from, pair to, T w) { make_edge(from.first*W+from.second,to.first*W+to.second,w); make_edge(to.first*W+to.second,from.first*W+from.second,w); } inline size_t size(){return N;} inline size_t idx(pair yx){return yx.first*W+yx.second;} }; /** * @url * @est */ int main() { cin.tie(0);ios::sync_with_stdio(false); int N,S,T,K; cin >> N >> S >> T >> K; S--,T--; vector X(N); for(int i=0;i> X[i]; int M; cin >> M; vector>> edge(N); for(int i=0;i> a >> b >> y; --a,--b; edge[a].push_back({b,y}); } auto dp = multivector(N,(K+1),LOWINF); auto dp2 = multivector(N,K+1,make_pair(-1,-1)); dp[S][1] = X[S]; priority_queue_reverse> pq; pq.push({dp[S][1],S,1}); while(pq.size()) { auto tp = pq.top(); pq.pop(); int64 cost = get<0>(tp); int from = get<1>(tp); int from_cnt = get<2>(tp); if(cost > dp[from][from_cnt]) continue; int nxt = min(from_cnt+1,K); for(auto p:edge[from]) { int to = p.first; int64 w = p.second; if(dp[to][nxt] > cost + w + X[to]) { dp[to][nxt] = cost + w + X[to]; dp2[to][nxt] = {from,from_cnt}; pq.push({dp[to][nxt],to,nxt}); } } } int ans_K = min_element(dp[T].begin()+K,dp[T].end())-dp[T].begin(); int64 ans = dp[T][ans_K]; if(ans>=LOWINF) { cout << "Impossible" << endl; } else { cout << "Possible" << endl; cout << ans << endl; vector P; int p = T, q = ans_K; while(p!=-1) { P.push_back(p+1); auto tmp=p,tmq=q; p = dp2[tmp][tmq].first; q = dp2[tmp][tmq].second; } reverse(ALL(P)); cout << P.size() << endl; for(int i=0;i