https://stats-learn.com/en/python/methods-to-virtue-the-multinomial-distribution-in-python/