#include<bits/stdc++.h>
using namespace std;
int st[35],k,n,nr;
float v[35],sum,sol,dp[35][35];
int main()
{
    freopen("dezastru.in","r",stdin);
    freopen("dezastru.out","w",stdout);
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++)
    {
        scanf("%f",&v[i]);
    }
    for(int i=1;i<=n;i++) dp[i][0]=1;
 //   dp[1][1]=1.0*v[1];
    for(int i=2;i<=n;i++)
    {
        for(int j=1;j<=i;j++)
        {
            dp[i][j]=dp[i-1][j]+dp[i-1][j-1]*v[i];
        }
    }

    long long  comb=1;
    for(int i=n;i>=(n-k+1);i--)
    {
        comb=comb*1LL*i;
    }
    for(int i=1;i<=k;i++)
    {
        comb=comb/i;
    }
    sol=(2*dp[n][k])/(1.0*comb);
        printf("%.6f\n",sol);
    return 0;
}