#include <bits/stdc++.h>
using namespace std;

template<typename T>
struct Dijkstra{
    const T inf=numeric_limits<T>::max();
    using P=pair<T,int>;
    int n;
    vector<vector<pair<int,T>>> G;
    vector<T> d;
    Dijkstra(int n):n(n),G(n),d(n){}
    void add_edge(int u,int v,T w){
        G[u].emplace_back(v,w);
    }
    vector<T> build(int s){
        fill(d.begin(),d.end(),inf);
        d[s]=0;
        priority_queue<P,vector<P>,greater<P>> pq;
        pq.emplace(d[s],s);
        while(!pq.empty()){
            P p=pq.top(); pq.pop();
            int v=p.second;
            if (d[v]<p.first) continue;
            for (auto &e:G[v]){
                int u=e.first; T c=e.second;
                if (d[v]+c<d[u]){
                    d[u]=d[v]+c;
                    pq.emplace(d[u],u);
                }
            }
        }
        return d;
    }
};

const int MAX_N=2e5+10;

int N,M,X,Y;
double p[MAX_N],q[MAX_N];

double dist(int i,int j){
    return sqrt((p[i]-p[j])*(p[i]-p[j])+(q[i]-q[j])*(q[i]-q[j]));
}

int main(){
    cin.tie(0);
    ios::sync_with_stdio(false);
    cin >> N >> M >> X >> Y;
    for (int i=0;i<N;++i) cin >> p[i] >> q[i];
    Dijkstra<double> D(N);
    for (int i=0;i<M;++i){
        int P,Q; cin >> P >> Q;
        double d=dist(--P,--Q);
        D.add_edge(P,Q,d);
        D.add_edge(Q,P,d);
    }
    cout << fixed << setprecision(10);
    cout << D.build(--X)[--Y] << '\n';
}