On the characterization of p-harmonic functions on the Heisenberg group by mean value properties