We consider the problem of constructing cosmological solutions of the Einstein-Maxwell equations that contain multiple charged black holes. By considering the field equations as a set of constraint and evolution equations, we construct exact initial data for N charged black holes on a hypersphere. This corresponds to the maximum of expansion of a cosmological solution, and provides sufficient information for a unique evolution. We then consider the specific example of a universe that contains eight charged black holes, and show that the existence of non-zero electric charge reduces the scale of the cosmological region of the space. These solutions generalize the Majumdar-Papapetrou solutions away from the extremal limit of charged black holes, and provide what we believe to be some of the first relativistic calculations of the effects of electric charge on cosmological backreaction.