We consider the following perturbed critical Dirichlet problem involving Hardy-Schr\"odinger operator on a smooth bounded domain $\Omega \subset \mathbb{R}^N$, $N\geq 3$, with $0 \in \Omega$: $$ \left\{ \begin{array}{ll}-\Delta u-\gamma \frac{u}{|x|^2}-\epsilon u=|u|^{\frac{4}{N-2}}u &\hbox{in }\Omega u=0 & \hbox{on }\partial \Omega, \end{array}\right. when $\epsilon>0$ is small and $\gamma< {(...