Try using Mill's ratio. This tells that if $x\geq 0$$$\int^{\infty}_x e^{-u^2/2}du\sim \frac{e^{-x^2/2}}{x}$$in fact, if you put $n=xe^{x^2/2}$ and $f(n)$ such that $f(n)=x$ the limit is constant. You can express also $f$ in terms of Lambert W function.
In the book "The Probabilistic Method", the integral $\int_\lambda^\infty e^{-t^2/2}dt$ is said to be "approximately equal" to $\frac{e^{-\lambda^2/2}}{\lambda}$ for large $\lambda$. I assume what is meant is $\int_\lambda^\infty e^{-t^2/2}dt \sim \frac{e^{-\lambda^2/2}}{\lambda}$, which can be verified using l'hopital, since the ratio of the derivatives is $\frac{\lambda^2}{\lambda^2+1}$, which tends to 1.