Das deutsche Python-Forum

Code: Alles auswählen

import scipy.special

n = 6
k = 3
p = 1/6

ergebnis = scipy.special.binom(n, k) * p**k * (1-p)**(n-k)
print(f"Wahrscheinlichkeit 3er z.B. '555':\t  {ergebnis}")
print(f"Wahrscheinlichkeit aller 3er Pasch:\t  {ergebnis * 6}")

Code: Alles auswählen

Wahrscheinlichkeit 3er z.B. '555':	  0.05358367626886146
Wahrscheinlichkeit aller 3er Pasch:	  0.32150205761316875

Code: Alles auswählen

import itertools
products = list(itertools.product("123456", repeat=6))
# products = list(itertools.product("123", repeat=4))

def convertTuple(tup):
    str = ''.join(tup)
    return str

kartesisches_produkt = len(list(products))
# print(f"kartesisches_produkt:\t  {kartesisches_produkt}")
treffer_pasch = 0
# Print 
for tup in list(products): 
  tup = tuple(sorted(tup))  
  # print(tup)
  tup_as_string = convertTuple(tup)
  # print(f"tup_as_string:\t  {tup_as_string}")
  if '555' in tup_as_string:
    treffer_pasch += 1
    
print(f"treffer_pasch:\t  {treffer_pasch}")
print("Treffer geteilt durch Gesamtzahl", (treffer_pasch / kartesisches_produkt))
print("multiplizier mit 6 Würfelseiten: ", (treffer_pasch / kartesisches_produkt) * 6)

Code: Alles auswählen

treffer_pasch:	  2906
Treffer geteilt durch Gesamtzahl 0.062285665294924554
multiplizier mit 6 Würfelseiten:  0.37371399176954734

Code: Alles auswählen

In [1]: from itertools import product

In [2]: trials = list(trial for trial in product(range(1,7), repeat=6))

In [3]: sum(trial.count(5) == 3 for trial in trials) / len(trials) * 6
Out[3]: 0.3215020576131687

In [4]: sum(trial.count(5) >= 3 for trial in trials) / len(trials) * 6
Out[4]: 0.37371399176954734