A Feynman Path Centroid Dynamics Approach for the Computation of Time Correlation Functions Involving Nonlinear Operators