Pagini recente » Cod sursa (job #1544984) | Cod sursa (job #1769856) | Cod sursa (job #2498979) | Cod sursa (job #2281640) | Cod sursa (job #84927)
Cod sursa(job #84927)
#include <cstdio>
#define maxn 26
FILE *in = fopen("dezastru.in","r"), *out = fopen("dezastru.out","w");
int n, k;
double nr[maxn];
double a[maxn][maxn]; // a[i][j] == suma tuturor produselor de j factori alesi din primele i numere
void read_init()
{
fscanf(in, "%d %d", &n, &k);
for ( int i = 1; i <= n; ++i )
fscanf(in, "%lf", &nr[i]);
for ( int i = 0; i <= n; ++i )
a[i][0] = 1;
}
long long fact(int n)
{
long long rez = 1;
for ( int i = 2; i <= n; ++i )
rez *= i;
return rez;
}
long long cnk(int n, int k)
{
return (fact(n)/fact(n-k))/fact(k);
}
int main()
{
read_init();
for ( int i = 1; i <= n; ++i )
for ( int j = 1; j <= k; ++j )
a[i][j] = a[i-1][j] + a[i-1][j-1]*nr[i];
// for ( int i = 0; i < n; ++i )
// {
// for ( int j = 0; j < k; ++j )
// printf("%lf ", a[i][j]);
// printf("\n");
// }
fprintf(out, "%.6lf\n", a[n][k]/(double)cnk(n,k));
return 0;
}