The effect of the electron–phonon interaction on an electron bound to a hydrogenic impurity in a three-dimensional (3D) anisotropic quantum dot (QD) is studied theoretically. We use the Landau–Pekar variational approach to calculate the binding energy of ground state (GS) and first-excited state (ES) with considering electron–phonon interaction. The expressions of the GS and ES energies under investigation depict a rich variety of dependent relationship with the variational parameters in three different limiting cases. Numerical calculations were performed for ZnSe QDs with different confinement lengths in the xy-plane and the z-direction, respectively. It is illustrated that binding energies of impurity polarons corresponding to each level are larger in small QDs. Furthermore, the contribution to binding energy from phonon is about 15% of the total binding energy.