#include<algorithm>
using namespace std;
short n,k;
double a[26],b[26][26];
long calc(){
    short i;
    long p;
    for(i=k+1,p=1; i<=n; p*=i++);
    for(i=1; i<=n-k; p/=i++);
    return p;}
void solve(){
    short i,j;
    int cnk;
    scanf("%hd%hd",&n,&k);
    for(i=1; i<=n; ++i)
        scanf("%lf",&a[i]);
    for(i=0; i<=n; ++i)
        b[i][0]=1;
    for(i=1; i<=n; ++i)
        for(j=1; j<=k; ++j)
            b[i][j]=b[i-1][j]+b[i-1][j-1]*a[i];
    cnk=calc();
    printf("%lf",b[n][k]/cnk);}
int main(){
    freopen("dezastru.in","r",stdin);
    freopen("dezastru.out","w",stdout);
    solve();
    return 0;}