fork download
copy 
  1. #include <stdio.h>
  2. #include <math.h>
  3.  
  4. double factorial(unsigned n, unsigned e) {
  5. 	return n > e ? n * factorial(n - 1, e) : 1;
  6. }
  7.  
  8. double combination(unsigned n, unsigned k) {
  9. 	unsigned e = n - k;
  10. 	return k < e ? factorial(n, e) / factorial(k, 0) : factorial(n, k) / factorial(e, 0);
  11. }
  12.  
  13. int main(void) {
  14. 	const double p = 0.01;
  15. 	const unsigned n = 100;
  16.  
  17. 	double exp = 0;
  18.  
  19. 	for (unsigned i = 0; i <= n; i++)
  20. 	{
  21. 		double pi = combination(n, i) * pow(1 - p, n - i) * pow(p, i);
  22. 		exp += i * pi;
  23. 		printf("%u %f %f\n", i, pi, exp);
  24. 	}
  25.  
  26. 	return 0;
  27. }
  28.  
Success #stdin #stdout 0s 9424KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
0 0.366032 0.000000
1 0.369730 0.369730
2 0.184865 0.739459
3 0.060999 0.922457
4 0.014942 0.982224
5 0.002898 0.996713
6 0.000463 0.999493
7 0.000063 0.999933
8 0.000007 0.999992
9 0.000001 0.999999
10 0.000000 1.000000
11 0.000000 1.000000
12 0.000000 1.000000
13 0.000000 1.000000
14 0.000000 1.000000
15 0.000000 1.000000
16 0.000000 1.000000
17 0.000000 1.000000
18 0.000000 1.000000
19 0.000000 1.000000
20 0.000000 1.000000
21 0.000000 1.000000
22 0.000000 1.000000
23 0.000000 1.000000
24 0.000000 1.000000
25 0.000000 1.000000
26 0.000000 1.000000
27 0.000000 1.000000
28 0.000000 1.000000
29 0.000000 1.000000
30 0.000000 1.000000
31 0.000000 1.000000
32 0.000000 1.000000
33 0.000000 1.000000
34 0.000000 1.000000
35 0.000000 1.000000
36 0.000000 1.000000
37 0.000000 1.000000
38 0.000000 1.000000
39 0.000000 1.000000
40 0.000000 1.000000
41 0.000000 1.000000
42 0.000000 1.000000
43 0.000000 1.000000
44 0.000000 1.000000
45 0.000000 1.000000
46 0.000000 1.000000
47 0.000000 1.000000
48 0.000000 1.000000
49 0.000000 1.000000
50 0.000000 1.000000
51 0.000000 1.000000
52 0.000000 1.000000
53 0.000000 1.000000
54 0.000000 1.000000
55 0.000000 1.000000
56 0.000000 1.000000
57 0.000000 1.000000
58 0.000000 1.000000
59 0.000000 1.000000
60 0.000000 1.000000
61 0.000000 1.000000
62 0.000000 1.000000
63 0.000000 1.000000
64 0.000000 1.000000
65 0.000000 1.000000
66 0.000000 1.000000
67 0.000000 1.000000
68 0.000000 1.000000
69 0.000000 1.000000
70 0.000000 1.000000
71 0.000000 1.000000
72 0.000000 1.000000
73 0.000000 1.000000
74 0.000000 1.000000
75 0.000000 1.000000
76 0.000000 1.000000
77 0.000000 1.000000
78 0.000000 1.000000
79 0.000000 1.000000
80 0.000000 1.000000
81 0.000000 1.000000
82 0.000000 1.000000
83 0.000000 1.000000
84 0.000000 1.000000
85 0.000000 1.000000
86 0.000000 1.000000
87 0.000000 1.000000
88 0.000000 1.000000
89 0.000000 1.000000
90 0.000000 1.000000
91 0.000000 1.000000
92 0.000000 1.000000
93 0.000000 1.000000
94 0.000000 1.000000
95 0.000000 1.000000
96 0.000000 1.000000
97 0.000000 1.000000
98 0.000000 1.000000
99 0.000000 1.000000
100 0.000000 1.000000
https://ideone.com/RRsQBt
language:
C (gcc 8.3)
created:
4 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!