#include <fstream.h>
#include <iostream.h>

int x[30],n,p;  float y[15000],prob[55]; int c=0;
long long int prd()
{       long long int prod=1; int i;
	for(i=p+1;i<=n;i++)
		prod*=i;
	return prod;
}

int valid(int t)
{       int i;
	for(i=1;i<t;i++)
		if(x[t]<=x[i])
			return 0;
	return 1;
}
void rez()
{       double sol=1;  int i;
	for(i=1;i<=p;i++)
		sol*=prob[x[i]];
	y[++c]=sol;
}


void back(int k)
{       int i;
	for(i=1;i<=n;i++)
	{	x[k]=i;

		if(valid(k))
			if(k==p)
			    rez();
			else
			   back(k+1);
	}
}

int main()
{
	fstream f("dezastru.in",ios::in), g("dezastru.out",ios::out);
	f>>n>>p;  int i;
	for(i=1;i<=n;i++)
		f>>prob[i];

	back(1);
	double sol=0;
	for(i=1;i<=c;i++)
		sol=sol+y[i]/prd();
	g<<sol;
	g.close();
	f.close();
	return 0;


}