#include<bits/stdc++.h>
using namespace std;
#define FOR(i,s,t) for(int i=(s),_t=(t); i<=_t; ++i)
typedef long long ll;
const int N=1e5+50;
const ll INF=1e18;
vector<int> E[N];
ll h[N];
int n;
void init(int s,ll *dis) {
    priority_queue<pair<int,int>> Q;
    FOR(i,1,n) dis[i]=INF;
    dis[s]=0;
    Q.push({0,s});
    while(Q.size()) {
        auto [d,x]=Q.top();
        Q.pop();
        if(-d>dis[x]) continue;
        for(auto y: E[x]) {
            ll v=-d+(h[y]-h[x])*(h[y]-h[x]);
            if(v<dis[y]) {
                dis[y]=v;
                Q.push({-v,y});
            }
        }
    }
}
ll dis1[N],dis2[N];
struct seg {
    ll k,b;
    ll value(int x) const{
        return k*x+b;
    }
};
struct node {
    int tag;
    seg s;
} T[N<<2];
#define ls p<<1
#define rs p<<1|1
#define lson L,mid,ls
#define rson mid+1,R,rs
void Down(int p) {
    if(T[p].tag) {
        T[ls]={1,{0,INF}};
        T[rs]={1,{0,INF}};
        T[p].tag=0;
    }
}
void ins(int L,int R,int p,seg s) {
    int mid=L+R>>1;
    if(s.value(mid)<T[p].s.value(mid)) swap(T[p].s,s);
    if(L==R) return ;
    Down(p);
    if(s.value(L)<T[p].s.value(L)) ins(lson,s);
    if(s.value(R)<T[p].s.value(R)) ins(rson,s);
}
ll qry(int L,int R,int p,int x) {
    if(L==R) return T[p].s.value(x);
    Down(p);
    int mid=L+R>>1;
    if(x<=mid) return min(T[p].s.value(x),qry(lson,x));
    else return min(T[p].s.value(x),qry(rson,x));
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout << fixed << setprecision(9);
    int m,s,t;
    cin >> n >> m >> s >> t;
    FOR(i,1,n) cin >> h[i];
    FOR(i,1,m) {
        int u,v;
        cin >> u >> v;
        E[u].push_back(v);
        E[v].push_back(u);
    }
    init(s,dis1);
    init(t,dis2);
    int mx=1e5;
    double res=dis1[t];
    for(auto j: E[s]) {
        res=min<double>(res,dis2[j]);
    }
    for(auto j: E[t]) {
        res=min<double>(res,dis1[j]);
    }
    FOR(i,1,n) {
        T[1]={1,{0,INF}};
        for(auto j: E[i]) {
            ins(0,mx,1,{-h[j]*2,dis1[j]*2+h[j]*h[j]});
        }
        for(auto j: E[i]) {
            res=min(res,0.5*(dis2[j]*2+h[j]*h[j]+qry(0,mx,1,h[j])));
        }
    }
    cout << res << '\n';
    return 0;
}