In this paper an efficient formulation of the Path integral (PI) approach is developed for determining the response probability density functions (PDFs) and first-passage statistics of nonlinear oscillators subject to stationary and time-modulated external Gaussian white noise excitations. Specifically, the evolution of the response PDF is obtained in short time steps, by using a discrete version of the Chapman-Kolmogorov equation and assuming a Gaussian form for the conditional response PDF. Next, the technique involves proceeding to treating the problem via an analytical asymptotic expansion procedure, namely the Laplace’s method of integration. In this manner, the repetitive double integrals involved in the standard implementation of the PI approach are evaluated in a closed form, while the response and first-passage PDFs are obtained by mundane step-by-step application of the derived approximate analytical expression. It is shown that the herein proposed formulation can drastically decrease the associated computational cost by several orders of magnitude, as compared to both the standard PI technique and Monte Carlo solution (MCS) approach. A number of nonlinear oscillators are considered in the numerical examples. Notably, for these systems both response PDFs and first-passage probabilities are presented, whereas comparisons with pertinent MCS data demonstrate the efficiency and accuracy of the technique.