A Feynman path centroid dynamics approach for the computation of time correlation functions involving nonlinear operators