Efficient simulation of non-Markovian path integrals via imaginary time evolution of an effective Hamiltonian