Sampling for a multinomial