Problem: Bearland has n cities, numbered 1 through n. There are m bidirectional roads. The i-th road connects two distinct cities ai and bi. No two roads connect the same pair of cities. It's possible to get from any city to any other city (using one or more roads).

The distance between cities a and b is defined as the minimum number of roads used to travel between a and b.

Limak is a grizzly bear. He is a criminal and your task is to catch him, or at least to try to catch him. You have only two days (today and tomorrow) and after that Limak is going to hide forever.

Your main weapon is BCD (Bear Criminal Detector). Where you are in some city, you can use BCD and it tells you the distance between you and a city where Limak currently is. Unfortunately, BCD can be used only once a day.

You don't know much about Limak's current location. You assume that he is in one of n cities, chosen uniformly at random (each city with probability <image>). You decided for the following plan:

  1. Choose one city and use BCD there. 
    * After using BCD you can try to catch Limak (but maybe it isn't a good idea). In this case you choose one city and check it. You win if Limak is there. Otherwise, Limak becomes more careful and you will never catch him (you loose). 
  2. Wait 24 hours to use BCD again. You know that Limak will change his location during that time. In detail, he will choose uniformly at random one of roads from his initial city, and he will use the chosen road, going to some other city. 
  3. Tomorrow, you will again choose one city and use BCD there. 
  4. Finally, you will try to catch Limak. You will choose one city and check it. You will win if Limak is there, and loose otherwise. 



Each time when you choose one of cities, you can choose any of n cities. Let's say it isn't a problem for you to quickly get somewhere.

What is the probability of finding Limak, if you behave optimally?

Input

The first line of the input contains two integers n and m (2 ≤ n ≤ 400, <image>) — the number of cities and the number of roads, respectively.

Then, m lines follow. The i-th of them contains two integers ai and bi (1 ≤ ai, bi ≤ n, ai ≠ bi) — cities connected by the i-th road.

No two roads connect the same pair of cities. It's possible to get from any city to any other city.

Output

Print one real number — the probability of finding Limak, if you behave optimally. Your answer will be considered correct if its absolute error does not exceed 10 - 6.

Namely: let's assume that your answer is a, and the answer of the jury is b. The checker program will consider your answer correct if |a - b| ≤ 10 - 6.

Examples

Input

3 3
1 2
1 3
2 3


Output

0.833333333333


Input

5 4
1 2
3 1
5 1
1 4


Output

1.000000000000


Input

4 4
1 2
1 3
2 3
1 4


Output

0.916666666667


Input

5 5
1 2
2 3
3 4
4 5
1 5


Output

0.900000000000

Note

In the first sample test, there are three cities and there is a road between every pair of cities. Let's analyze one of optimal scenarios.

  1. Use BCD in city 1. 
    * With probability <image> Limak is in this city and BCD tells you that the distance is 0. You should try to catch him now and you win for sure. 
    * With probability <image> the distance is 1 because Limak is in city 2 or city 3. In this case you should wait for the second day. 
  2. You wait and Limak moves to some other city. 
    * There is probability <image> that Limak was in city 2 and then went to city 3. 
    * <image> that he went from 2 to 1. 
    * <image> that he went from 3 to 2. 
    * <image> that he went from 3 to 1. 
  3. Use BCD again in city 1 (though it's allowed to use it in some other city). 
    * If the distance is 0 then you're sure Limak is in this city (you win). 
    * If the distance is 1 then Limak is in city 2 or city 3. Then you should guess that he is in city 2 (guessing city 3 would be fine too). 



You loose only if Limak was in city 2 first and then he moved to city 3. The probability of loosing is <image>. The answer is <image>.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <typename T>
bool readInteger(T& x) {
  char c, r = 0, n = 0;
  x = 0;
  while (true) {
    c = getchar();
    if (c < 0 && !r)
      return 0;
    else if (c == '-' && !r)
      n = 1;
    else if (c >= '0' && c <= '9')
      x = x * 10 + c - '0', r = 1;
    else if (r)
      break;
  }
  if (n) x = -x;
  return 1;
}
template <typename T>
vector<T> readIntegral(int n) {
  vector<T> ret(n);
  for (int i = 0; i < n; i++) readInteger(ret[i]);
  return ret;
}
auto readInts = readIntegral<int>;
auto readLongs = readIntegral<long long>;
template <typename T, T MOD>
struct ModInt {
  T value;
  ModInt() : value(0) {}
  ModInt(T x) {
    x %= MOD;
    if (x < 0) x += MOD;
    value = x;
  }

 public:
  ModInt& operator+=(ModInt x) {
    value += x.value;
    if (value >= MOD) value -= MOD;
    return *this;
  }
  ModInt& operator-=(ModInt x) {
    value -= x.value;
    if (value < 0) value += MOD;
    return *this;
  }
  ModInt& operator*=(ModInt x) {
    value *= x.value;
    value %= MOD;
    return *this;
  }
  ModInt& operator/=(ModInt x) {
    x.invert();
    return *this *= x;
  }
  ModInt operator+(ModInt x) const {
    ModInt o = *this;
    o += x;
    return o;
  }
  ModInt operator-(ModInt x) const {
    ModInt o = *this;
    o -= x;
    return o;
  }
  ModInt operator*(ModInt x) const {
    ModInt o = *this;
    o *= x;
    return o;
  }
  ModInt operator/(ModInt x) const {
    ModInt o = *this;
    o /= x;
    return o;
  }
  bool operator==(ModInt x) const { return value == x.value; }
  bool operator!=(ModInt x) const { return !(*this == x); }
  ModInt operator-() const { return ModInt(0) - *this; }
  ModInt operator^(long long x) const {
    ModInt ret = 1 % MOD;
    ModInt mul = *this;
    while (x) {
      if (x & 1) ret *= mul;
      mul *= mul;
      x >>= 1;
    }
    return ret;
  }
  ModInt& operator^=(long long x) { return *this = *this ^ x; }

 private:
  void invert() { *this ^= MOD - 2; }

 public:
  void answer() { std::cout << value << std::endl; }
};
const int MX = 405;
vector<int> E[MX];
int D[MX][MX];
int n, m;
template <typename T>
void uniquify(vector<T>& v) {
  sort(v.begin(), v.end());
  v.resize(unique(v.begin(), v.end()) - v.begin());
}
vector<pair<int, double>> process(vector<pair<int, double>>& v) {
  vector<pair<int, double>> ret;
  sort(v.begin(), v.end());
  for (auto& pr : (v)) {
    if (ret.empty() || pr.first != ret.back().first) {
      ret.push_back(pr);
    } else {
      ret.back().second += pr.second;
    }
  }
  return ret;
}
double dpr[MX];
int stamp[MX];
double f(int x) {
  double base = 1.0 / n;
  vector<double> pdist(n);
  vector<int> countAt(n);
  for (int i = (0); i < int(n); ++i) pdist[D[x][i]] += base;
  for (int i = (0); i < int(n); ++i) ++countAt[D[x][i]];
  vector<vector<pair<int, double>>> possLoc(n);
  for (int i = (0); i < int(n); ++i)
    for (auto& o : (E[i])) possLoc[D[x][i]].push_back({o, base / E[i].size()});
  double tot = 0;
  for (int i = (0); i < int(n); ++i) {
    if (countAt[i] == 0) continue;
    double catchChance = 1.0 / countAt[i] * pdist[i];
    possLoc[i] = process(possLoc[i]);
    for (int j = (0); j < int(n); ++j) stamp[j] = -1;
    for (int j = (0); j < int(n); ++j) {
      double ccatch = 0;
      for (auto& pr : (possLoc[i])) {
        auto& q = dpr[D[j][pr.first]];
        if (stamp[D[j][pr.first]] != j) {
          stamp[D[j][pr.first]] = j;
          q = 0;
        }
        auto r = pr.second;
        if (r > q) {
          ccatch += r - q;
          q = r;
        }
      }
      catchChance = max(catchChance, ccatch);
    }
    tot += catchChance;
  }
  return tot;
}
int main() {
  for (int i = (0); i < int(MX); ++i)
    for (int j = (0); j < int(MX); ++j) D[i][j] = ((int)1e9 + 10);
  readInteger(n);
  readInteger(m);
  for (int i = (0); i < int(n); ++i) D[i][i] = 0;
  for (int k = (0); k < int(m); ++k) {
    int a, b;
    readInteger(a);
    readInteger(b);
    --a;
    --b;
    D[a][b] = min(D[a][b], 1);
    D[b][a] = min(D[b][a], 1);
    E[a].push_back(b);
    E[b].push_back(a);
  }
  for (int i = (0); i < int(n); ++i)
    for (int j = (0); j < int(n); ++j)
      for (int k = (0); k < int(n); ++k)
        D[j][k] = min(D[j][k], D[j][i] + D[i][k]);
  double ans = 0;
  for (int i = (0); i < int(n); ++i) ans = max(ans, f(i));
  printf("%.10f\n", ans);
}