[Numpy-discussion] Efficient way to draw multinomial distribution random samples