For a quasi-periodically forced differential equation, response solutions are quasi-periodic ones whose frequency vector coincides with that of the forcing function and they are known to play a fundamental role in the harmonic and synchronizing behaviors of quasi-periodically forced oscillators. These solutions are well-understood in quasi-periodically perturbed nonlinear oscillators either in the presence of large damping or in the non-degenerate cases with small or free damping. In this paper, we consider the existence of response solutions in quasi-periodically perturbed, second order differential equations, including nonlinear oscillators, of the form x ¨ + λ x l = ϵ f (ω t , x , x ˙) , x ∈ R , where λ is a constant, 0 < ϵ ≪ 1 is a small parameter, l > 1 is an integer, ω ∈ R d is a frequency vector, and f : T d × R 2 → R 1 is real analytic and non-degenerate in x up to a given order p ≥ 0 , i.e., [ f (· , 0 , 0) ] = [ ∂ f (· , 0 , 0) ∂ x ] = [ ∂ 2 f (· , 0 , 0) ∂ x 2 ] = ⋯ = [ ∂ p - 1 f (· , 0 , 0) ∂ x p - 1 ] = 0 and [ ∂ p f (· , 0 , 0) ∂ x p ] ≠ 0 , where [ ] denotes the average value of a continuous function on T d . In the case that λ = 0 and f is independent of x ˙ , the existence of response solutions was first shown by Gentile (Ergod Theory Dyn Syst 27:427–457, 2007) when p = 1 . This result was later generalized by Corsi and Gentile (Commun Math Phys 316:489–529, 2012; Ergod Theory Dyn Syst 35:1079–1140, 2015; Nonlinear Differ Equ Appl 24(1):article 3, 2017) to the case that p > 1 is odd. In the case λ ≠ 0 , the existence of response solutions is studied by the authors Si and Yi (Nonlinearity 33(11):6072–6099, 2020) when p = 0 . The present paper is devoted to the study of response solutions of the above quasi-periodically perturbed differential equations for the case λ ≠ 0 by allowing p > 0 . Under the conditions that 0 ≤ p < l / 2 and λ [ ∂ p f (· , 0 , 0) ∂ x p ] > 0 when l - p is even, we obtain a general result which particularly implies the following: (1) If either l is odd and λ < 0 or l is even and [ ∂ p f (· , 0 , 0) ∂ x p ] > 0 , then as ϵ sufficiently small response solutions exist for each ω satisfying a Brjuno-like non-resonant condition; (2) If either l is odd and λ > 0 or l is even and [ ∂ p f (· , 0 , 0) ∂ x p ] < 0 , then there exists an ϵ ∗ > 0 sufficiently small and a Cantor set E ∈ (0 , ϵ ∗) with almost full Lebesgue measure such that response solutions exist for each ϵ ∈ E and ω satisfying a Diophantine condition. Similar results are also obtained in the case λ = ± ϵ which particularly concern the existence of large amplitude response solutions. [ABSTRACT FROM AUTHOR]