A friend designed a drinking game with a lucky wheel of 30 distinct icons. When playing, each one takes turn to spin the wheel, and write down the items until the first one who gets the item that has already appeared then he drinks. And he asked what the average (expectation) of turns is. The formula is not so difficult to get but it is quite unpractical to compute the value. So I was thinking about an asymptotic formula.

Luckily, I found this here: See section 5.6, "average number of people".

My question is how to prove the Ramanujan's formula given there,

$$ Q(M)=\sum_{k=1}^{M} \frac{M!}{(M-k)! M^k}\sim\sqrt{\frac{\pi M}{2}}-\frac{1}{3}+\frac{1}{12}\sqrt{\frac{\pi}{2M}}-\frac{4}{135M}+\cdots$$

Any good reference available?