On the calculation of nonlinear optical macroscopic susceptibilities in solution within the polarizable continuum model