#include<cstdio>

#include<algorithm>

#include<cstring>

using namespace std;

const int maxn=;

struct edge{int from,v;}e[maxn*];

int n,m,a,b,ins[maxn],first[maxn],d[maxn],q[],p[maxn][maxn],tot;

double f[maxn][maxn];

void insert(int u,int v)

{

    tot++;e[tot].v=v;e[tot].from=first[u];first[u]=tot;ins[u]++;

    tot++;e[tot].v=u;e[tot].from=first[v];first[v]=tot;ins[v]++;

}

double dp(int a,int b)

{

    if(f[a][b])return f[a][b];

    if(a==b)return f[a][b]=;

    if(p[a][b]==b||p[p[a][b]][b]==b)return f[a][b]=;

    double ans=dp(p[p[a][b]][b],b);

    for(int i=first[b];i;i=e[i].from)

     ans+=dp(p[p[a][b]][b],e[i].v);

    return f[a][b]=ans/(1.0*ins[b]+)+;

}

void bfs(int x)

{

    memset(d,-,sizeof(d));

    int head=,tail=;d[x]=;p[x][x]=;

    for(int i=first[x];i;i=e[i].from)p[x][e[i].v]=e[i].v,d[e[i].v]=,q[tail++]=e[i].v;

    while(head!=tail)

     {

         int y=q[head++];if(head>)head=;int num=p[x][y];

         for(int i=first[y];i;i=e[i].from)

          if(d[e[i].v]==-||((d[y]+==d[e[i].v])&&num<p[x][e[i].v]))

           {

               d[e[i].v]=d[y]+;

               p[x][e[i].v]=num;

               q[tail++]=e[i].v;

               if(tail>)tail=;

           }

     }

}

int main()

{

    scanf("%d%d%d%d",&n,&m,&a,&b);

    for(int i=;i<=m;i++)

     {

         int u,v;

         scanf("%d%d",&u,&v);

         insert(u,v);

     }

    for(int i=;i<=n;i++)bfs(i);

    printf("%.3lf",dp(a,b));

    return ;

}