We present the first analytical results for the O(αα_{s}) corrections to the total cross section for the inclusive production of an on-shell Z boson at hadron colliders. We include the complete set of contributions, with photonic and massive weak gauge boson effects, which have been computed in analytical form and expressed in terms of polylogarithmic and elliptic integrals. We present numerical results, relevant for the precision studies at the LHC. These corrections increase the accuracy of the predictions and contribute to the reduction of the QCD component of the theoretical uncertainty.