Multiple small solutions for Schrödinger equations involving positive quasilinear term

We consider the multiplicity of solutions of a class of quasilinear Schrödinger equations involving the $p$-Laplacian: \begin{equation*} -\Delta_{p} u+V(x)|u|^{p-2}u+\Delta_{p}(u^{2})u=K(x)f(x,u),\qquad x\in \mathbb{R}^{N}, \end{equation*} where $\Delta_{p} u=\operatorname{div}(|\nabla u|^{p-2}\nabla u)$, $1<p