#include <bits/stdc++.h>
#define NMAX 30
using namespace std;
ifstream fin("dezastru.in");
ofstream fout("dezastru.out");
int n, k;
double p[NMAX], dp[NMAX][NMAX];
void C();
double factorial(int n);
int main()
{
    fin>>n>>k;
    for(int i=1; i<=n; i++) fin>>p[i];
    for(int i=0; i<=n; i++) dp[i][0]=1;
    C();
    double rezultat=dp[n][k]*factorial(k)*factorial(n-k)/factorial(n);
    fout<<fixed<<setprecision(6)<<rezultat<<'\n';
    return 0;
}
double factorial(int n)
{
    if(!n)  return 1;
    return n*factorial(n-1);
}
void C()
{
    for(int i=1; i<=n; i++)
        for(int j=1; j<=i; j++)
            dp[i][j]=dp[i-1][j]+dp[i-1][j-1]*p[i];
}