Advertisement
Guest User

Untitled

a guest
May 26th, 2018
75
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 0.84 KB | None | 0 0
  1. #include <iostream>
  2. #include <math.h>
  3. #include <bits/stdc++.h>
  4.  
  5. using namespace std;
  6. int BC(int n, int k)
  7. {
  8.     int C[k+1];
  9.     memset(C, 0, sizeof(C));
  10.  
  11.     C[0] = 1;  // nC0 is 1
  12.  
  13.     for (int i = 1; i <= n; i++)
  14.     {
  15.         for (int j = min(i, k); j > 0; j--)
  16.             C[j] = C[j] + C[j-1];
  17.     }
  18.     return C[k];
  19. }
  20.  
  21. double sBern (int n, int k, double p) // n > k
  22. {
  23.     double q = 1 - p;
  24.     double B_C = BC(n,k);
  25.     return B_C * pow(p,k) * pow(q,n-k);
  26. }
  27.  
  28. int main()
  29. {
  30.     ios_base::sync_with_stdio(0);
  31.     double n,s,d;
  32.  
  33.     while(cin >> s >> d >> n)
  34.     {
  35.         //d = double(10)/3;
  36.         s /= d;
  37.         s /= d;
  38.         double p=n/2, P=0;
  39.         for(int i=n; i >= p; i--)
  40.         {
  41.             P += sBern(n,i,s);
  42.         }
  43.         cout << setprecision(16) << P << "\n";
  44.     }
  45. }
  46.  
  47.  
  48. // 0.5 3.333333 10
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement