#include <bits/stdc++.h>
#include <fstream>
using namespace std;
typedef long long int64;
typedef vector<int> vec;
typedef vector<int64> vec64;
  string __fname = "dezastru";
  ifstream in (__fname + ".in");
  ofstream out (__fname + ".out");
  #define cin in
  #define cout out
#define pii pair <int, int> 
#define pii64 pair <int64, int64>
#define ss cout << " ";
#define nn cout << "\n";
#define ct(x) cout << x;
#define cts(x) cout << x << " ";
#define ctn(x) cout << x << "\n";
#define db(x) cout << "> " << #x << ": " << x << "\n";
#define qr queries();
void solve(int);
void queries(){int n;cin >> n;for (int i = 1; i <= n; i++) solve(i);}
int64 ceildiv(int64 a, int64 b) {return a / b + !!(a % b);}
template<class T>T gcd (T a, T b){return (b ? gcd(b, a % b): a);}
template<class T>T lcm (T a, T b){return a * b / gcd(a, b);}
// // // // // // // // // // // // // // // // // // // // // // 
/*                  TEMPLATE - VANILLA                         */
// // // // // // // // // // // // // // // // // // // // // //
const int maxn = 200200;
const int64 mod = 1000000007;
const long double pi = 3.14159265359;
const int ddx[] = {-1, -1, 0, 1, 1, 1, 0, -1};
const int ddy[] = {0, 1, 1, 1, 0, -1, -1, -1};
const int dx[] = {-1, 0, 1, 0};
const int dy[] = {0, 1, 0, -1};
long double rs = 0, d = 0;
int n,k;
vector <long double> a;

void solve(int id){
    
    return;
}

void back (int step, int sf, long double calc) {
    if (sf == k) {
        rs+=calc;
        d+=1;
        return;
    }
    if (step == n) return;
    for (int i = step + 1; i < n; i++){
        back(i, sf + 1, calc*a[i]);
    }
}


int main(){
    ios_base::sync_with_stdio(0);cin.tie(0);
    cin >> n >> k;
    a.resize(n);
    for (int i = 0; i < n; i++){
        cin >> a[i];
    }
    back(-1, 0, 1);
    // db(rs);
    // db(d);
    cout << fixed << setprecision(12) << (rs / d);
    return 0;
}