A method is presented that allows exact calculations of fragment multiplicity distributions for a canonical ensemble of non-interacting clusters. Fragmentation properties are shown to depend on only a few parameters. Fragments are shown to be copiously produced above the transition temperature. At this transition temperature, the calculated multiplicity distributions broaden and become strongly super-Poissonian. This behavior is compared to predictions from a percolation model. A corresponding microcanonical formalism is also presented.