#include <bits/stdc++.h>
#define NMAX 26
using namespace std;
ifstream f("dezastru.in");
ofstream g("dezastru.out");
int n,k,v[NMAX];
long long comb[NMAX][NMAX];
double prob[NMAX],sol,dp[NMAX][NMAX];
int main()
{
    f>>n>>k;
    for(int i=1;i<=n;i++)
    {
        f>>prob[i];
        comb[i][1]=i;
        comb[i][i]=1;
    }
    for(int i=0;i<=n;i++)
        dp[i][0]=1;
    for(int i=2;i<=n;i++)
    {
        for(int j=2;j<i;j++)
            comb[i][j]=comb[i-1][j]+comb[i-1][j-1];
    }
    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]*prob[i]);
    }
    //cout<<dp[n][k];
    sol=dp[n][k]/comb[n][k];
    g<<fixed<<setprecision(7)<<sol;
    return 0;
}