Problem: There are n students living in the campus. Every morning all students wake up at the same time and go to wash. There are m rooms with wash basins. The i-th of these rooms contains ai wash basins. Every student independently select one the rooms with equal probability and goes to it. After all students selected their rooms, students in each room divide into queues by the number of wash basins so that the size of the largest queue is the least possible. Calculate the expected value of the size of the largest queue among all rooms.

Input

The first line contains two positive integers n and m (1 ≤ n, m ≤ 50) — the amount of students and the amount of rooms. The second line contains m integers a1, a2, ... , am (1 ≤ ai ≤ 50). ai means the amount of wash basins in the i-th room.

Output

Output single number: the expected value of the size of the largest queue. Your answer must have an absolute or relative error less than 10 - 9.

Examples

Input

1 1
2


Output

1.00000000000000000000


Input

2 2
1 1


Output

1.50000000000000000000


Input

2 3
1 1 1


Output

1.33333333333333350000


Input

7 5
1 1 2 3 1


Output

2.50216960000000070000

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
using namespace std::chrono;
template <class e, class t, int N>
std::basic_istream<e, t>& operator>>(std::basic_istream<e, t>& in,
                                     const e (&literal)[N]) {
  std::array<e, N - 1> buffer;
  in >> buffer[0];
  if (N > 2) in.read(&buffer[1], N - 2);
  if (strncmp(&buffer[0], literal, N - 1))
    in.setstate(in.rdstate() | std::ios::badbit);
  return in;
}
template <class e, class t>
std::basic_istream<e, t>& operator>>(std::basic_istream<e, t>& in,
                                     const e& literal) {
  e buffer;
  in >> buffer;
  if (buffer != literal) in.setstate(in.rdstate() | std::ios::badbit);
  return in;
}
template <class e, class t, int N>
std::basic_istream<e, t>& operator>>(std::basic_istream<e, t>& in,
                                     e (&carray)[N]) {
  return std::operator>>(in, carray);
}
template <class t1, class t2>
istream& operator>>(istream& in, pair<t1, t2>& p) {
  in >> p.first >> p.second;
  return in;
}
template <class t>
istream& operator>>(istream& in, vector<t>& v) {
  for (int i = 0; i < v.size(); i++) {
    in >> v[i];
  }
  return in;
}
template <class t>
ostream& operator<<(ostream& out, vector<t>& a) {
  for (auto x : a) {
    out << x << " ";
  }
  return out;
}
template <class t1, class t2>
ostream& operator<<(ostream& out, map<t1, t2>& a) {
  for (auto x : a) {
    out << x << " ";
  }
  return out;
}
template <class t>
ostream& operator<<(ostream& out, set<t>& a) {
  for (auto x : a) {
    out << x << " ";
  }
  return out;
}
template <class t1, class t2>
ostream& operator<<(ostream& out, const pair<t1, t2>& a) {
  out << a.first << " " << a.second;
  return out;
}
template <class t1, class t2>
pair<t1, t2> operator+(const pair<t2, t1>& a, const pair<t1, t2>& b) {
  return {a.first + b.first, a.second + b.second};
}
template <class t1, class t2>
pair<t1, t2> operator-(const pair<t2, t1>& a, const pair<t1, t2>& b) {
  return {a.first - b.first, a.second - b.second};
}
long long gcd(long long a, long long b) { return b == 0 ? a : gcd(b, a % b); }
long long lcm(long long a, long long b) { return a * (b / gcd(a, b)); }
long long fpow(long long b, long long e, long long m = (long long)(1e9 + 7)) {
  long long r = 1;
  while (e) {
    if (e & 1) r = (r * b) % m;
    e >>= 1;
    b = (b * b) % m;
  }
  return r;
}
long long n, m;
vector<long long> a;
long long c[100][100];
long double dp[100][100][100];
void init() {
  for (long long i = 0; i < (long long)60; i++)
    for (long long j = 0; j < (long long)60; j++)
      for (long long k = 0; k < (long long)70; k++) dp[i][j][k] = -2;
  for (long long i = 0; i < (long long)51; i++) {
    for (long long j = 0; j < (long long)51; j++) {
      if (j == 0)
        c[i][j] = 1;
      else if (i == 0)
        c[i][j] = 0;
      else
        c[i][j] = c[i - 1][j - 1] + c[i - 1][j];
    }
  }
}
long double solve(long long i, long long k, long long mx) {
  if (i == m) return (!k ? mx : 0);
  long double& r = dp[i][k][mx];
  if (r > -1) return r;
  long double p = 1;
  r = 0;
  for (long long j = 0; j < (long long)k + 1; j++) {
    long long len = (j + a[i] - 1) / a[i];
    r += (c[k][j] * p * solve(i + 1, k - j, max(mx, len)));
    p /= m;
  }
  return r;
}
int32_t main(int32_t argc, char const* argv[]) {
  ios::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  ;
  cin >> n >> m;
  a.resize(m);
  cin >> a;
  init();
  cout << setprecision(15) << solve(0, n, 0) << "\n";
  return 0;
}