Measurement of the electroweak production of Z$γ$ and two jets in proton-proton collisions at $\sqrt{s} =$ 13 TeV and constraints on anomalous quartic gauge couplings