fork download
  1. from math import factorial
  2. F = factorial
  3.  
  4. # a: (-1, -1)
  5. # b: (-1, 0)
  6. # c: (-1, 1)
  7. # d: (0, 0)
  8. # e: (0, 1)
  9. # f: (1, 1)
  10.  
  11. N = 1
  12.  
  13. total_probability = 0
  14.  
  15. for a in range(N+1):
  16. for b in range(N-a+1):
  17. for c in range(N-a-b+1):
  18. for d in range(N-a-b-c+1):
  19. for e in range(N-a-b-c-d+1):
  20. f = N-a-b-c-d-e
  21.  
  22. prob = (1/36)**a * (1/18)**b * (1/12)**c * (1/9)**d * (1/6)**e * (1/4)**f
  23. combs = F(N)/(F(a)*F(b)*F(c)*F(d)*F(e)*F(f)) * 2**(b+c+e)
  24.  
  25. if a-c+f == 0:
  26. total_probability += prob * combs
  27.  
  28. print(total_probability)
  29.  
Success #stdin #stdout 0.1s 10088KB
stdin
Standard input is empty
stdout
0.5555555555555556