2026-07-072026-07-07http://salesiana.dossiersoluciones.com/handle/123456789/167138Generalized probability distributions for Maxwell-Boltzmann, Bose-Einstein and Fermi-Dirac statistics, with unequal source probabilities $q_i$ for each level $i$, are obtained by combinatorial reasoning. For equiprobable degenerate sublevels, these reduce to those given by Brillouin in 1930, more commonly given as a statistical weight for each statistic. These distributions and corresponding cross-entropy (divergence) functions are shown to be special cases of the Pólya urn model, involving neither independent nor identically distributed ("ninid") sampling. The most probable Pólya distribution contains the Acharya-Swamy intermediate statistic.8 pages; 4 figuresStatistical Mechanicstext