#include <stdio.h>

using namespace std;

int n,k;
float a[27][27],b[27];
double r;

int main()
{
    FILE *f=fopen("dezastru.in","r");
    fscanf(f,"%d%d",&n,&k);
    for(int i = 1;i<=n;i++)
        fscanf(f,"%f",&b[i]);

    fclose(f);
    a[0][0]=1;
    for(int i = 1;i<=n;i++)
    {
        a[i][0]=1;
        for(int j = 1;j<=i;j++)
        {
            a[i][j] = a[i-1][j-1]*b[i]+a[i-1][j];
        }
    }
    r=1;
    if(k<n/2)
    {
        for(int i =n-k+1;i<=n;i++)
            r*=i;
        for(int i = 1;i<=k;i++)
            r/=i;
    }
    else
    {
        for(int i = k+1;i<=n;i++)
            r*=i;
        for(int i = 1;i<k;i++)
            r/=i;
    }


    FILE*g =fopen("dezastru.out","w");
    fprintf(g,"%g",a[n][k]/r);
    fclose(g);



    return 0;
}