The partition function of all classical spin models, including all discrete standard statistical models and all Abelian discrete lattice gauge theories (LGTs), is expressed as a special instance of the partition function of the 4D Z(2) LGT. This unifies all classical spin models with apparently very different features in a single complete model. This result is applied to establish a new method to compute the mean-field theory of Abelian discrete LGTs with d >= 4, and to show that computing the partition function of the 4D Z(2) LGT is computationally hard (#P hard). The 4D Z(2) LGT is also proved to be approximately complete for Abelian continuous models. The proof uses techniques from quantum information.