Density-matrix functional theory of the Hubbard model: An exact numerical study