Operator splitting methods for numerical integration of weakly perturbed N-body systems