An open-chain imaginary-time path-integral sampling approach to the calculation of approximate symmetrized quantum time correlation functions