Python - probabilidad de dados con python

cara=[1,2,3,4,5,6]

total_casos = 6 ** 3

sum = 10

casos_a_favor=0

caso=0

for i in range(1,7):

  for j in range(1,7):

    for k in range(1,7):

      if i + j + k == sum:

        casos_a_favor += 1

print(casos_a_favor)

probabiidad = 1.0*casos_a_favor/total_casos

print(probabiidad)

from math import factorial

m, n = 6, 3

#Casos favorables:

cf = 13

# Casos posibles:

cp = factorial(m+n-1)/(factorial(n)*factorial(m-1))

p = cf / cp

print(p, round(p,3))

cf, cp = [], []

for n1 in range(1,7):

    for n2 in range(1,7):

        for n3 in range(1,7):

            # Casos favorables, cf (suman 10)

            if n1+n2+n3 == 10 and sorted((n1,n2,n3)) not in cf:

                cf.append(sorted((n1,n2,n3)))

            # Casos posibles, cp

            if sorted((n1,n2,n3)) not in cp:

                cp.append(sorted((n1,n2,n3)))

print(f'p(10) = {len(cf)/len(cp)}')

print(f'casos posibles = {cp} \ncasos favorables = {cf}')