We have \(d\) mutually exclusive outcomes and \(n\) independent trials. This function enumerates all possible vectors of length \(d\) of counts of each outcome for \(n\) trials, i.e., the sample space. The result is output as a matrix with \(d\) columns where each row represents a possible observation.