#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5,M=3e5+5;
int n,m,head[N],cnt,tot,du[N],hh[N];
bool v[M];
struct edge{
	int to,next;
}e[M],e1[M];
void add(int u,int v){
	e[++cnt].to=v;
	e[cnt].next=head[u];
	head[u]=cnt;
}
void add1(int u,int v){
	e1[++tot].to=v;
	e1[tot].next=hh[u];
	hh[u]=tot;
	du[v]++;
}
stack<int> stk;
int dfn[N],low[N],dfs_clock,sccno[N],scc_cnt,siz[N];
void tarjan(int u){
	dfn[u]=low[u]=++dfs_clock;
	stk.push(u);
	for(int i=head[u];i;i=e[i].next){
		int v=e[i].to;
		if(!dfn[v]){
			tarjan(v);
			low[u]=min(low[u],low[v]);
		}else if(!sccno[v])
			low[u]=min(low[u],dfn[v]);
	}
	if(low[u]==dfn[u]){
		scc_cnt++;
		while(1){
			int v=stk.top();
			stk.pop();
			sccno[v]=scc_cnt;
			siz[scc_cnt]++;
			if(u==v)
				break;
		}
	}
}
int main()
{
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
	cin>>n>>m;
	if(n==1){
		cout<<"1.000000";
		return 0;
	}
	int x,y;
	for(int i=1;i<=m;i++){
		cin>>x>>y;
		add(x,y);
	}
	for(int i=1;i<=n;i++)
		if(!dfn[i])
			tarjan(i);
	for(int i=1;i<=n;i++){
		for(int j=head[i];j;j=e[j].next){
			if(sccno[i]==sccno[e[j].to]||v[sccno[e[j].to]])
				continue;
			add1(sccno[i],sccno[e[j].to]);
			v[sccno[e[j].to]]=1;
		}
		for(int j=head[i];j;j=e[j].next)
			v[sccno[e[j].to]]=0;
	}
	cnt=0;
	bool flag=1;
	for(int i=1;i<=scc_cnt;i++){
		if(!du[i])
			cnt++;
		if(du[i]||siz[i]>1)
			continue;
		if(flag){
			bool f=1;
			for(int j=hh[i];j;j=e1[j].next){
				if(du[e1[i].to]<=1){
					f=0;
					break;
				}
				if(f)
					cnt--,flag=0;
			}
		}
	}
	int a=n-cnt;
	printf("%.6lf",(double)a/n);
	return 0;
}