Skip to main content
4 of 6
added 50 characters in body
Carlo Beenakker
  • 188.3k
  • 18
  • 448
  • 651

The $p$-th guest receives a share of $$s_{p,N}=p N^{-p-1} (N-p+1) \frac{\Gamma(N+1)}{\Gamma(N-p+2)}$$ For large $N$ I substitute the asymptotic expansion of the Gamma function, $$\Gamma(z)\approx e^{-z}z^z(2\pi/z)^{1/2},$$ to arrive at the approximation $$s_{p,N}\approx p(N-p+1) e^{1-p} \sqrt{\frac{1}{N+1}}\sqrt{\frac{1}{N-p+2}} (N+1)^{N+1} N^{-p-1} (N-p+2)^{-N+p-1}$$ which agrees very well with the exact value (compare the dots with the continuous curve in the plot below, for $N=1000$)

Now the desired number $p$ of fortunate guests should follow by solving $s_{p,N}=1/N$ for $p$. I expect an answer $p\approx \sqrt{N}$, but will need to work a bit more for a precise result.

Carlo Beenakker
  • 188.3k
  • 18
  • 448
  • 651