Multipactor discharge on a dielectric surface: Statistical theory and simulation results