Pagini recente » Cod sursa (job #1108156) | Cod sursa (job #1421621) | Cod sursa (job #3156058) | Cod sursa (job #933806) | Cod sursa (job #1142032)
#include <iostream>
#include <fstream>
#include <iomanip>
#include <ctime>
using namespace std;
ifstream f("dezastru.in");
ofstream g("dezastru.out");
int n, k, C[26][26];
double prod, prob[26], rez[26][26];
void read()
{
f >> n >> k;
for(int i = 1; i <= n; ++i)
f >> prob[i];
f.close();
}
void back(int niv, int x, double aux)
{
if(niv == k)
{
prod += aux;
return;
}
for(int i = x; i+k <= n+niv; ++i)
{
back(niv + 1, i+1, aux * prob[i]);
}
}
void solve()
{
int i, j;
prod = 0;
//back(0, 0, 1.);
C[0][0] = rez[0][0] = 1;
for (i = 1; i <= n; ++i)
for (C[i][0] = 1, rez[i][0] = 1, j = 1; j <= k; ++j)
rez[i][j] = rez[i-1][j] + rez[i-1][j-1] * prob[i], C[i][j] = C[i-1][j]+C[i-1][j-1];
g << fixed << setprecision(6) << (double)rez[n][k] / (double)C[n][k];
/*do {
prod += prob[nrp/k];
for(int i = 0; i < perm.size(); ++i)
cout << perm[i] << " ";
cout << '\n';
nrp++;
} while(next_permutation(perm.begin(), perm.end()));*/
}
int main()
{
//double tclock = clock();
read();
solve();
//cout << (clock() - tclock) / CLOCKS_PER_SEC;
return 0;
}