Letting $$ a_n = \sqrt{2\pi n}\cdot e^{-n}\sum_{k=0}^{n-1}\frac{n^k}{k!}\cdot{\frac1{n-k}}, $$ does anyone know of a simple asymptotic equivalent for $a_n$? Numerical experimentation suggests that $$ a_n \stackrel{?}\sim \tfrac12\log n. $$ If we get rid of the $\frac1{n-k}$, there is a simple answer, because $$ \lim_{n\to\infty} e^{-n}\sum_{k=0}^{n-1}\frac{n^k}{k!}=\frac12. $$ as shown here. Throwing in the factor $\frac1{n-k}$ seems to change the asymptotics entirely.
The context is that $a_n$ is result of applying Stirling's approximation to $\sum_{k=1}^n \frac{n_{(k)}}{kn^k},$ which is the expected number of "cycles" that a random function on a set of $n$ elements to itself has, where a cycle is a set of distinct numbers $\{x_1,x_2,\dots,x_k\}$ for which $f(x_i)=x_{i+1}$. If you replace "random function" with "random bijection," the expected number of cycles is about $\log n$, so if my guess is correct then random functions tend to have half as many cycles as permutations.
Here we follow exercise 5.2 from Asymptopia by Joel Spencer and show OP's assumption is correct.
A convenient approach is to split the index range of $k$ into three parts:
a small range: $\qquad\quad k<\frac{\sqrt{n}}{\ln n}$
a middle range: $\quad\quad \frac{\sqrt{n}}{\ln n} < k < \sqrt{n}\ln n$
and a large range: $\ \quad k>\sqrt{n}\ln n$.