We present a state-of-the-art calculation of the next-to-leading-orderelectroweak corrections to ZZ production, including the leptonic decays of theZ bosons into $\mu^+\mu^-\mathrm{e}^+\mathrm{e}^-$ or $\mu^+\mu^-\mu^+\mu^-$final states. We use complete leading-order and next-to-leading-order matrixelements for four-lepton production, including contributions of virtual photonsand all off-shell effects of Z bosons, where the finite Z-boson width is takeninto account using the complex-mass scheme. The matrix elements are implementedinto Monte Carlo programs allowing for the evaluation of arbitrary differentialdistributions. We present integrated and differential cross sections for theLHC at 13 TeV both for an inclusive setup where only lepton identification cutsare applied, and for a setup motivated by Higgs-boson analyses in thefour-lepton decay channel. The electroweak corrections are divided intophotonic and purely weak contributions. The former show the well-knownpronounced tails near kinematical thresholds and resonances; the latter aregenerically at the level of $\sim-5\%$ and reach several $-10\%$ in thehigh-energy tails of distributions. Comparing the results for$\mu^+\mu^-\mathrm{e}^+\mathrm{e}^-$ and $\mu^+\mu^-\mu^+\mu^-$ final states,we find significant differences mainly in distributions that are sensitive tothe $\mu^+\mu^-$ pairing in the $\mu^+\mu^-\mu^+\mu^-$ final state. Differencesbetween $\mu^+\mu^-\mathrm{e}^+\mathrm{e}^-$ and $\mu^+\mu^-\mu^+\mu^-$channels due to interferences of equal-flavour leptons in the final state canreach up to $10\%$ in off-shell-sensitive regions. Contributions induced byincoming photons, i.e. photon-photon and quark-photon channels, are included,but turn out to be phenomenologically unimportant.