  We give a construction of completely integrable $(2m)$-dimensional
Hamiltonian systems with cubic Hamilton functions. The construction depends on
a constant skew-Hamiltonian matrix $A$, that is, a matrix satisfying $A^{\rm
T}J=JA$, where $J$ is a non-degenerate skew-symmetric matrix defining the
standard symplectic structure on the phase space $\mathbb R^{2m}$. Applying to
any such system the so called Kahan-Hirota-Kimura discretization scheme, we
arrive at a birational $(2m)$-dimensional map. We show that this map is
symplectic with respect to a symplectic structure that is a perturbation of the
standard symplectic structure on $\mathbb R^{2m}$, and possesses $m$
independent integrals of motion, which are perturbations of the original
Hamilton functions and are in involution with respect to the invariant
symplectic structure. Thus, this map is completely integrable in the
Liouville-Arnold sense. Moreover, under a suitable normalization of the
original $m$-tuples of commuting vector fields, their Kahan-Hirota-Kimura
discretizations also commute and share the invariant symplectic structure and
the $m$ integrals of motion.
