#include #include #include #include #include #include #include #include #include #include #include #define fir first #define sec second #define sz(s) (s).size() #define pb push_back #define get(n) scanf("%d",&n); #define gets(s) string s;cin >> (s); #define prfi(n) printf("%d", &n); #define prfd(n) printf("%lf", &n); #define All(s) (s).begin(), (s).end() #define rep(i,j) for(int (i)=0;(i)<(j);(i)++) #define For(i,j,k) for(int (i)=(j);(i)<(k);(i)++) #define drep(i,j) for(int (i)=(j);(i)>=0;(i)--) #define Ford(i,j,k) for(int (i)=(j);i>=(k);i--) #define fore(c,v) for(auto (c): (v)) #define lp for(int __=0;__; using pll = std::pair; using vi = std::vector ; using vvi = std::vector ; using vll = std::vector; using vvll = std::vector; using vd = std::vector ; using vvd = std::vector ; using qi = std::queue ; using vpii = std::vector >; using vpll = std::vector; using namespace std; const int Mod = (1e9) + 7; const int INF = 1e9 + 19; const ll INFL = 1e18 + 19; const int dx[] = {-1, 0, 0, 1}; const int dy[] = {0, -1, 1, 0}; const int dx2[] = {-1, -1, -1, 0, 0, 0, 1, 1, 1}; const int dy2[] = {1, 0, -1, 1, 0, -1, 1, 0, -1}; const double EPS = 1e-10; //_____________________________________Templates_________________________________________// template inline void chmin(T1 &a, T2 b){if(a > b) a = b;} template inline void chmax(T1 &a, T2 b){if(a < b) a = b;} template inline void pri(T a){cout << a << endl;} template using vec = vector; template using min_priority_queue = priority_queue, greater>; //mainly use for dynamic prog template void update(T1 &a, T2 b){ a += b; if(a > Mod) a %= Mod; } inline void IN(void){ return; } template void IN(First& first, Rest&... rest){ cin >> first; IN(rest...); return; } inline void OUT(void){ cout << "\n"; return; } template void OUT(First first, Rest... rest){ cout << first << " "; OUT(rest...); return; } bool pairsort(pll pl, pll pr){ if(pl.first == pr.first)return pl.second > pr.second; return pl.first < pr.first; } int cntbit(ll a,int n,int j){int res = 0;For(i,j,n){if(a>>i & 1){res++;}}return res;} vector make_bit(int a){vector res; for(int i=31;i>=0;i--)if(a&(1< a)return GCD(b,a);if(a%b == 0)return b;else return GCD(b, a%b);} int LCM(int a, int b){return a*b/GCD(a,b);} int roundup(int a, int b){if(a % b == 0)return a/b;else return (a+b)/b;} int rounddown(int a, int b){if(a%b == 0)return a/b;else {return (a-b)/b;}} ll pow(ll a, ll n){ll res = 1;while(n > 0){if(n&1)res *= a; a *= a; n = n >> 1;}return res;} ll GetDiviserCount(ll N)//約数の個数 { ll res = 1; For(i,2,sqrt(N)+1) { ll cnt = 0; while(N%i == 0) { cnt++; N /= i; } res *= (cnt + 1); if(N == 1)break; } if(N != 1)res *= 2; return res; } vll GetDivisor(ll N)//約数列挙 { vll res; for(ll i = 1;i*i <= N;i++) { if(N%i == 0) { res.pb(i); if(i*i != N)res.pb(N/i); } } sort(All(res)); return res; } struct mint { ll x; // typedef long long ll; mint(ll x=0):x((x%Mod+Mod)%Mod){} mint operator-() const { return mint(-x);} mint& operator+=(const mint a) { if ((x += a.x) >= Mod) x -= Mod; return *this; } mint& operator-=(const mint a) { if ((x += Mod-a.x) >= Mod) x -= Mod; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= Mod; return *this; } mint operator+(const mint a) const { mint res(*this); return res+=a; } mint operator-(const mint a) const { mint res(*this); return res-=a; } mint operator*(const mint a) const { mint res(*this); return res*=a; } mint pow(ll t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime mod mint inv() const { return pow(Mod-2); } mint& operator/=(const mint a) { return (*this) *= a.inv(); } mint operator/(const mint a) const { mint res(*this); return res/=a; } friend ostream& operator<<(ostream& os, const mint& a); }; ostream& operator<<(ostream& os, const mint& a) { os << a.x; return os; } using vm = vector; using vvm = vector; struct combination { vector fact, ifact; combination(int n):fact(n+1),ifact(n+1) { assert(n < Mod); fact[0] = 1; for (int i = 1; i <= n; ++i) fact[i] = fact[i-1]*i; ifact[n] = fact[n].inv(); for (int i = n; i >= 1; --i) ifact[i-1] = ifact[i]*i; } mint operator()(int n, int k) { if (k < 0 || k > n) return 0; return fact[n]*ifact[k]*ifact[n-k]; } }; template struct MyMatrix { using mat = MyMatrix; vector> m_dat; int m_h; int m_w; MyMatrix(int h, int w) : m_h(h), m_w(w) ,m_dat(h,vector(w)) {} vector &operator[] (int idx) { return m_dat[idx]; } mat Multiple(mat &a) { mat C(m_h, a.m_w); for(int i=0;i 0) { if(x&1)B = B.Multiple(A); A = A.Multiple(A); x = x >> 1; } return B; } //friend ostream& operator<<(ostream &os, const mat &A); }; /* template ostream& operator<<(ostream &os, Matrix& A) { for(int i=0;i using mat = MyMatrix; //_____________________ following sorce code_________________________// const int max_n = 3 * (1e5) + 1; const int max_m = 83 * (1e5) + 1; int n,m,k; ll N; int h,w; string S; int a,b,c; vi v; int ans; double x,y,z; struct Edge { int to; ll cost; Edge(int to,ll cost) : to(to), cost(cost) {} bool operator>(const Edge &a) const { return cost > a.cost; } }; vector G[1010101]; vll score; vvll dp; void Dijkstra(int root) { min_priority_queue qu; qu.push(Edge(root,0)); // dp[root][root] = 0; while(qu.size()) { auto e = qu.top();qu.pop(); int curr = e.to; ll cost = e.cost; if(dp[root][curr] < e.cost)continue; dp[root][curr] = e.cost; fore(to,G[curr]) { if(dp[root][to.to] < e.cost + to.cost)continue; dp[root][to.to] = e.cost + to.cost; qu.push(Edge(to.to,e.cost + to.cost)); } } } signed main (int argc, char* argv[]) { cin.tie(0); ios::sync_with_stdio(false); IN(n,m,k);k--; score = vll(n); int whi = 0; rep(i,n) { IN(score[i]); if(score[i] == 0)whi++; } if(whi == (n-1)) { pri(-1); return 0; } rep(i,m) { ll c; IN(a,b,c);a--;b--; G[a].push_back(Edge(b,c)); G[b].push_back(Edge(a,c)); } dp = vvll(n+1,vll(n+1,INFL)); ll ans = INFL; rep(i,n)Dijkstra(i); rep(i,n) { ll sum = 0; ll del = INFL; rep(j,n) { if(score[j]) { sum += score[j]*2*dp[i][j]; chmin(del,dp[k][j]-dp[j][i]); } } chmin(ans,sum+del); } pri(ans); //pri(ans); //for(auto c : ans){cout << c << endl;} //cout << fixed << setprecision(15) << ans << endl; return 0; }