No material da aula, há um erro na forma como a função factorial é chamada. O correto é utilizar math.factorial(k) após importar o módulo math, e não np.math.factorial(k).
Não sei o que eu estaria fazendo de errado, eu tive que importar a biblioteca math e coloquei o seguinte código para funcionar:
import math
probabilidade = (np.e ** (-media)) * (media ** k) / (math.factorial(k))
print('%0.8f' % probabilidade)
Conte com o apoio da comunidade Alura na sua jornada. Abraços e bons estudos!