#include<iostream>
#include<iomanip>
#include<vector>
using namespace std;
int N,M;
vector<int>G[666];
double dp[666];
double x[666][666];
int main()
{
	cin>>N>>M;
	for(int i=0;i<M;i++)
	{
		int s,t;cin>>s>>t;
		s--,t--;
		G[s].push_back(t);
	}
	dp[N-1]=0;
	x[N-1][N-1]=1;
	for(int i=N-2;i>=0;i--)
	{
		double sum=0;
		for(int v:G[i])
		{
			sum+=dp[v];
			for(int j=0;j<N;j++)x[i][j]+=x[v][j]/G[i].size();
		}
		x[i][i]=1;
		dp[i]=sum/G[i].size()+1;
	}
	double ans=dp[0];
	for(int i=0;i<N-1;i++)
	{
		if(G[i].size()==1)continue;
		double sum=0,m=0;
		for(int v:G[i])
		{
			sum+=dp[v];
			m=max(m,dp[v]);
		}
		double ch=(sum-m)/(G[i].size()-1)+1-dp[i];
		ans=min(ans,dp[0]+ch*x[0][i]);
	}
	cout<<fixed<<setprecision(16)<<ans<<endl;
}
