Path integral Monte Carlo approach for weakly bound van der Waals complexes with rotations: Algorithm and benchmark calculations