Cod sursa(job #84927)

Utilizator DastasIonescu Vlad Dastas Data 18 septembrie 2007 17:33:02
Problema Dezastru Scor 30
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.03 kb
#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;
}