#include<iostream>
#include<string>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<list>
#include<cstdint>
#include<bits/stdc++.h>
#include<queue>
#include<set>

#define rep(i, n) for(int i = 0; i < (int)(n); i++)
#define rep_r(i, n) for(int i = n-1; i >= 0; i--)
#define rep1(i, n) for(int i = 1; i <= (int)(n); i++)
#define REP(i,n) for(int i=0, i##_len=(n); i<i##_len; ++i)
#define all(x) (x).begin(),(x).end()
#define SZ(x) ((ll)(x).size())
#define bit(n) (1LL<<(n))
#define UNIQUE(v) v.erase( unique(v.begin(), v.end()), v.end() );

#define INF bit(60)
#define pb push_back
#define mod 1000000007

template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }

using namespace std;
using uif=uint_fast64_t;
using ll=long long int;

int qp(int a,ll b){int ans=1;do{if(b&1)ans=1ll*ans*a%mod;a=1ll*a*a%mod;}while(b>>=1);return ans;}
int qp(int a,ll b,int mo){int ans=1;do{if(b&1)ans=1ll*ans*a%mo;a=1ll*a*a%mo;}while(b>>=1);return ans;}

int gcd(int a,int b){return b?gcd(b,a%b):a;}
int dx[4]={1,0,-1,0};
int dy[4]={0,1,0,-1};


#define FACSIZE 2002

ll invfac[FACSIZE];
ll fac[FACSIZE];

ll mod_pow(ll ,ll );
ll mod_add(ll ,ll );
ll mod_mul(ll ,ll );
ll mod_div(ll ,ll );
void make_fact(ll );
void make_invfact(ll );

ll n,m;

struct Node{
  //unordered_set<ll> children;
  vector<ll> children;
};

//vector< set<ll,greater<ll>> > influence;
vector<Node> graph;
vector<double> dp_origin;
vector<double> dp;

bool findPath(ll start,ll target){
  if(start==target) return true;
  queue<ll> que;
  que.push(start);
  while(!que.empty()){
    auto q=que.front();
    que.pop();
    for( auto x:graph[q].children ){
      if( x==target ) return true;
      else{
        if(x<target){
          que.push(x);
        }
      }
    }

  }
  return false;
}

double e(ll x){
  if(x==n) return dp_origin[x]=0.0;
  if( dp_origin[x]<1e290 ) return dp_origin[x];

  ll num=graph[x].children.size();
  double expe=0.0;
  ll cnt=0;
  for(auto c:graph[x].children){
    double val=dp_origin[c];
    if(val<1e290){
      expe += 1.0 * ( val +1.0);
      cnt++;
    }
    else return 1e300;
  }

  return dp_origin[x]=expe/cnt;

}


double e(ll x,pair<ll,ll> p){
  if(x==n) return dp_origin[x]=0.0;
  if( dp_origin[x]<1e290 ) return dp_origin[x];

  ll num=graph[x].children.size();
  double expe=0.0;
  ll cnt=0;
  for(auto c:graph[x].children){
    double val=dp_origin[c];
    if( p==make_pair(x,c) ){

    }
    else if(val<1e290){
      expe += 1.0 * ( val +1.0);
      cnt++;
    }
    else return 1e300;
  }
  if(cnt==0) return 1e300;

  return dp_origin[x]=expe/cnt;

}

vector<double> dp_cut;

double dfs(ll ,ll);

int main(void){

  cin>>n>>m;
  vector< pair<ll,ll> > P(m);

  graph.resize(n+1);
  dp_origin.resize(n+1,1e295);
  rep(i,m){
    ll s,t;
    cin>>s>>t;
    //graph[s].children.insert(t);
    graph[s].children.push_back(t);
    
  }
  //double ans=e(1);

  for(ll i=n;i>=0;i--){
    e(i);
  }
  double ans=dp_origin[1];

  unordered_map<ll,ll> um;
  for(ll i=1;i<n;i++){
    for(auto c:graph[i].children){
      if(um.find(i)==end(um)) um[i]=c;
      else{
        if( dp_origin[um[i]] < dp_origin[c])  um[i]=c;
      }
    }
  }
  

  for(auto p:um)
  {
    ll x=p.first;
    for(ll i=0;i<n+1;i++) dp_origin[i]=1e295;
    //graph[x].children.erase(p.second);
    for(ll i=n;i>=0;i--){
      e(i,p);
    }
    ans=min(ans,dp_origin[1]);
    //graph[x].children.insert(p.second);
  }

  cout << setprecision(20) << ans << endl;

  return 0;

}



ll mod_pow(ll x,ll r){
  if(r==0) return 1;
  else if(r==1) return x%mod;
  else if(r%2==0){
    ll t=mod_pow(x,r/2)%mod;
    return t*t %mod;
  } 
  else{
    ll t=mod_pow(x,r/2)%mod;
    ll k= t*t %mod;
    return k*x%mod;
  } 
}
ll mod_add(ll a,ll b){
  return (a+b)%mod;
}

ll mod_mul(ll a,ll b){
  return (a*b)%mod;
}

ll mod_div(ll a,ll b){
  return mod_mul(a,mod_pow(b,mod-2));
}

void make_fact(ll n){
  fac[0]=1;
  rep(i,n){
    fac[i+1]=mod_mul(fac[i],i+1);
  }  
}

void make_invfact(ll n){
  invfac[n]= mod_pow( fac[n] ,mod-2);
  for(int i=n-1;i>=0;i--){
    invfac[i]=mod_mul(invfac[i+1],i+1);
  }

}


Measured memory_bytes: 5376