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;
namespace IO {
template <class T>
inline void readin(T& x) {
  char c;
  bool f = 0;
  while ((c = getchar()) < '0' || '9' < c) f |= (c == '-');
  for (x = (c ^ 48); '0' <= (c = getchar()) && c <= '9';
       x = (x << 1) + (x << 3) + (c ^ 48))
    ;
  if (f) x = -x;
}
template <class T>
inline T readin(T x) {
  x = 0;
  char c;
  bool f = 0;
  while ((c = getchar()) < '0' || '9' < c) f |= (c == '-');
  for (x = (c ^ 48); '0' <= (c = getchar()) && c <= '9';
       x = (x << 1) + (x << 3) + (c ^ 48))
    ;
  return f ? -x : x;
}
template <class T>
void fwrit(const T x) {
  if (x < 0) return (void)(putchar('-'), fwrit(-x));
  if (x > 9) fwrit(x / 10);
  putchar(x % 10 ^ 48);
}
template <class T>
inline T Max(const T x, const T y) {
  return x < y ? y : x;
}
template <class T>
inline T Min(const T x, const T y) {
  return x < y ? x : y;
}
template <class T>
inline T fab(const T x) {
  return x > 0 ? x : -x;
}
inline int gcd(const int a, const int b) { return b ? gcd(b, a % b) : a; }
inline void getInv(int inv[], const int lim, const int MOD) {
  inv[0] = inv[1] = 1;
  for (int i = 2; i <= lim; ++i)
    inv[i] = 1ll * inv[MOD % i] * (MOD - MOD / i) % MOD;
}
inline long long mulMod(const long long a, const long long b,
                        const long long mod) {
  return ((a * b - (long long)((long double)a / mod * b + 1e-8) * mod) % mod +
          mod) %
         mod;
}
}  // namespace IO
using namespace IO;
const int mod = 998244353;
const int maxm = 50;
const int maxn = 50;
double f[maxm + 5][maxn + 5][maxn + 5], C[maxn + 5][maxn + 5];
int n, m, a[maxm + 5];
inline void Init() {
  scanf("%d%d", &n, &m);
  for (signed i = (1), i_end_ = (m); i <= i_end_; ++i) scanf("%d", &a[i]);
  for (signed i = (0), i_end_ = (n); i <= i_end_; ++i) {
    C[i][0] = C[i][i] = 1;
    for (signed j = (1), j_end_ = (i - 1); j <= j_end_; ++j)
      C[i][j] = C[i - 1][j - 1] + C[i - 1][j];
  }
}
inline void Get_f() {
  f[0][0][0] = 1;
  for (signed i = (1), i_end_ = (m); i <= i_end_; ++i)
    for (signed j = (0), j_end_ = (n); j <= j_end_; ++j)
      for (signed k = (0), k_end_ = (n); k <= k_end_; ++k)
        for (signed l = (0), l_end_ = (n - j); l <= l_end_; ++l)
          f[i][j + l][Max((l + a[i] - 1) / a[i], k)] +=
              f[i - 1][j][k] * C[n - j][l];
}
inline void Get_ans() {
  double all = 0;
  for (signed i = (0), i_end_ = (n); i <= i_end_; ++i) all += f[m][n][i];
  double ans = 0;
  for (signed i = (1), i_end_ = (n); i <= i_end_; ++i)
    ans += i * f[m][n][i] * 1.0 / all;
  printf("%.20f\n", ans);
}
signed main() {
  Init();
  Get_f();
  Get_ans();
  return 0;
}