Neural Contraction Metrics with Formal Guarantees for Discrete-Time Nonlinear Dynamical Systems