#include #include "unistd.h" using namespace std; #define WA cout<<"char134217728";exit(0) #define RE throw #define TLE sleep(1000) #define FOR(i,a,b) for(int i=(a),i##formax=(b);i=i##formax;i--) #define pb push_back #define mp make_pair #define fi first #define se second #define pcnt __builtin_popcount #define sz(x) (int)(x).size() #define maxs(x,y) x=max(x,y) #define mins(x,y) x=min(x,y) #define show(x) cout<<#x<<" = "<(l,r)(rnd)) #define randDouble(l,r) (uniform_real_distribution(l,r)(rnd)) #define dout(d) printf("%.12f\n",d) typedef long long ll; typedef __int128_t lll; typedef pair pii; typedef pair pll; templateusing V=vector; templateusing VV=V>; templateostream& operator<<(ostream& o,const pair& p){return o<<"("<void Fill(A (&array)[N],const T&v){fill((T*)array,(T*)(array+N),v);} lll gcd(lll a,lll b,lll &x,lll &y){if(!b){x=1;y=0;return a;}lll d=gcd(b,a%b,y,x);y-=a/b*x;return d;} ll gcd(ll a,ll b){lll x=0,y=0;return gcd(a,b,x,y);} ll modInv(ll a,ll m){lll x,y;gcd(a,m,x,y);return (x%m+m)%m;} ll modPow(lll a,lll n,ll m){if(!a)return a;lll p=1;for(;n>0;n>>=1,a=a*a%m)if(n&1)p=p*a%m;return(ll)p;} bool isPrime(ll n){if(n<2||n%2==0)return n==2;lll t=n-1,d=t/(t&-t);for(lll a:{2,325,9375,28178,450775,9780504,1795265022})if(a%n){for(t=d,a=modPow(a,t,n);t!=n-1&&a!=1&&a!=n-1;a=a*a%n,t=t*2%n);if(a!=n-1&&t%2==0)return 0;}return 1;} const int IINF = 1e9+6; const ll LINF = 1e18; const int MOD = 1e9+7; const double PI = acos(-1); const double EPS = 1e-10; static random_device rd; static mt19937 rnd(rd()); const int N = 1e5; struct Dijk{ typedef long long ll; const ll tinf = 1e18; int n; vector d; void solve(const vector>>&E, int S, int G=-1){ n = E.size(); d = vector(n, tinf); d[S] = 0; priority_queue,vector>,greater>> Q; Q.push(pair(0, S)); while(!Q.empty()){ pair t=Q.top(); Q.pop(); if(t.second == G) return; if(t.first > d[t.second]) continue; for(auto itr = E[t.second].begin(); itr != E[t.second].end(); itr++){ ll k = t.first + itr->second; int m = itr->first; if(k < d[itr->first]){ d[m] = k; Q.push(pair(k, m)); } } } } } dijk; int n, m, l; V t; VV> e; main(){ cin.tie(0); ios::sync_with_stdio(false); cin >> n >> m >> l; l--; FOR(i, 0, n){ int _t; cin >> _t; t.pb(_t); } e.resize(n); FOR(i, 0, m){ int a, b, c; cin >> a >> b >> c; a--; b--; e[a].pb(mp(b, c)); e[b].pb(mp(a, c)); } ll ans = LINF; dijk.solve(e, l, -1); V dl = dijk.d; FOR(i, 0, n){ dijk.solve(e, i, -1); ll s = 0, ma = -LINF; FOR(j, 0, n){ if(!t[j]) continue; s += dijk.d[j]*t[j]; maxs(ma, dijk.d[j]-dl[j]); } mins(ans, s*2-ma); } cout << ans << endl; }