#include <fstream>
#include <iomanip>
using namespace std;

float dp[30][30];

int comb(int n,int k) {
    if(n==k)
        return 1;
    if(k==1)
        return n;
    return comb(n,n-k);
}

int main()
{
    ifstream fin("dezastru.in");
    ofstream fout("dezastru.out");
    int n,k;
    fin>>n>>k;
    float p;
    dp[0][0]=1;
    for(int i=1;i<=n;i++) {
        fin>>p;
        dp[i][0]=1;
        for(int j=1;j<=k;j++)
            dp[i][j]=dp[i-1][j]+dp[i-1][j-1]*p;
    }
    fout<<setprecision(6)<<fixed<<dp[n][k]/comb(n,k);
    return 0;
}