A dynamic relaxation method with operator splitting and random-choice strategy for SPH