#include <cstdio>
#define size 26
using namespace std;
int n, k, i, j, a;
float m[size][size], v[size], c=1;

int main()
{
    freopen("dezastru.in", "r", stdin);
    freopen("dezastru.out", "w", stdout);
    scanf("%d %d", &n, &k);
    m[0][0]=1;
    for(i=1;i<=n;i++)
        scanf("%f", &v[i]), m[i][0]=1;

    for(i=1;i<=n;i++)
        for(j=1;j<=i;j++)
            m[i][j]=m[i-1][j]+m[i-1][j-1]*v[i];

    a=n-k;
    for(i=k+1;i<=n;i++)
    {
        c*=i/a;
        if(a>1)
            a--;
    }

    printf("%.6f\n", (float)m[n][k]/c);
    return 0;
}