Randomized exponential integrators for modulated nonlinear Schrödinger equations