Multi-value numerical methods for Hamiltonian systems