// #pragma GCC optimize("O3")
// #pragma GCC optimize("unroll-loops")
// #pragma GCC target("avx2,bmi,bmi2,popcnt,lzcnt")
#include<bits/stdc++.h>
 
using namespace std;
 
#define i64 long long
#define mp make_pair
#define pb push_back
#define all(x) (x).begin(), (x).end()
struct Line {
  mutable i64 k, m, p;
  bool operator<(const Line& o) const { return k < o.k; }
  bool operator<(i64 x) const { return p < x; }
};
struct LineContainer : multiset<Line, less<>> {
  // (for doubles, use inf = 1/.0, div(a,b) = a/b)
  static const i64 inf = 1e18;
  i64 div(i64 a, i64 b) { // floored division
    return a / b - ((a ^ b) < 0 && a % b); }
  bool isect(iterator x, iterator y) {
    if (y == end()) return x->p = inf, 0;
    if (x->k == y->k) x->p = x->m > y->m ? inf : -inf;
    else x->p = div(y->m - x->m, x->k - y->k);
    return x->p >= y->p;
  }
  void add(i64 k, i64 m) {
    auto z = insert({k, m, 0}), y = z++, x = y;
    while (isect(y, z)) z = erase(z);
    if (x != begin() && isect(--x, y)) isect(x, y = erase(y));
    while ((y = x) != begin() && (--x)->p >= y->p)
      isect(x, erase(y));
  }
  i64 query(i64 x) {
    // assert(!empty());
    if (empty()) return -inf;
    auto l = *lower_bound(x);
    return l.k * x + l.m;
  }
};

const int MAXN = 1e5 + 10;
const i64 INF = 1e18;
int N, M, S, T;
i64 H[MAXN], f[3][MAXN];
vector<int> G[MAXN];
void dijkstra(int st, int t) {
  for (int i = 1; i <= N; i ++) f[t][i] = INF;
  f[t][st] = 0;
  priority_queue<pair<i64, int>> pq; pq.push(mp(0, st));
  while (!pq.empty()) {
    auto [d, u] = pq.top(); pq.pop();
    if (f[t][u] != -d) continue;
    for (auto v : G[u]) {
      i64 add = (H[u] - H[v]) * (H[u] - H[v]);
      if (f[t][v] > f[t][u] + add) {
        f[t][v] = f[t][u] + add;
        pq.push(mp(-f[t][v], v));
      }
    }
  }
  // for (int i = 1; i <= N; i ++) {
  //   cout << t << " " << f[t][i] << "\n";
  // }
}
void Solve(void) {
  cin >> N >> M >> S >> T;
  for (int i = 1; i <= N; i ++) cin >> H[i];
  for (int i = 1; i <= M; i ++) {
    int u, v; cin >> u >> v;
    G[u].pb(v);
    G[v].pb(u);
  } dijkstra(S, 0); dijkstra(T, 1);
  long double ans = INF;
  for (int u = 1; u <= N; u ++) {
    ans = min(ans, (long double)(f[0][u] + f[1][u]));
    LineContainer convex;
    for (auto v : G[u]) {
      convex.add(2 * H[v], -(2 * f[0][v] + H[v] * H[v]));
    }
    for (auto v : G[u]) {
      long double cur = convex.query(H[v]);
      // cout << u << " " << v << " " << cur << " " << 2 * f[1][v] << "\n";
      cur = (-cur + 2 * f[1][v] + H[v] * H[v])/2;
      ans = min(ans, cur);
      // cout << u << " " << v << " " << cur << " " << f[1][v] << "\n";
    }
  }
  cout << ans << "\n";
}
signed main() {
  ios_base::sync_with_stdio(false); cin.tie(0);
  cout << fixed << setprecision(10);
  int Tests = 1; // cin >> Tests; 
  while (Tests --) {
    Solve();    
  }
}