We present a methodology for accelerating the estimation of the free energy from path integral Monte Carlo simulations by considering an intermediate artificial reference system where interactions are inexpensive to evaluate numerically. Using the spherically averaged Ewald interaction as this intermediate reference system for the uniform electron gas, the interaction contribution for the free energy was evaluated up to 18 times faster than the Ewald-only method. Furthermore, an extrapolation technique with respect to the quantum statistics was tested and applied to alleviate the Fermion sign problem. Combining these two techniques enabled the evaluation of the free energy for a system of 1000 electrons, where both finite-size and statistical errors are below chemical accuracy. The general procedure can be applied to systems relevant for planetary and inertial confinement fusion modeling with low to moderate levels of quantum degeneracy.