We present a formalism for the fermionic quasiparticle propagator in a superfluid fermionic system. Starting from a general many-body Hamiltonian confined by the two-body instantaneous interaction, the equation of motion for the fermionic propagator is obtained in the Dyson form. Before making any approximation, the interaction kernel is found to be decomposed into static and dynamical (time-dependent) contributions, where the latter translates to the energy-dependent and the former maps to the energy-independent terms in the energy domain. The three-fermion correlation function, being the heart of the dynamical part of the kernel, is factorized into the two-fermion and one-fermion ones. With the relaxed particle number constraint, the normal propagator is coupled to the anomalous one via both the static and dynamical kernels, which is formalized by introducing the generalized quasiparticle propagator of the Gor'kov type. The dynamical kernel in the factorized form is associated with the quasiparticle-vibration coupling (QVC), with the vibrations unifying both the normal and pairing phonons. The QVC vertices are related to the variations of the Hamiltonian of the Bogoliubov quasiparticles, which can be obtained by the finite amplitude method.