This note reports on a search for an extended scalar sector of the standard model, where a new CP-even (odd) scalar decays to a Z boson and a lighter CP-odd (even) scalar, which further decays to . The Z boson is reconstructed via its decays to leptons. The analysed data were recorded in proton-proton collisions at \sqrt{s} = 13~\mathrm{TeV}, collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9~\mathrm{fb}^{-1}. Data and predictions from the standard model are in agreement within uncertainties. Upper limits at confidence level are set on the production cross section times branching ratio, with masses of the resonances ranging up to 1000~\mathrm{GeV}. The results are interpreted in the context of the Two-Higgs Doublet Model.