Least k such that exp(log(1) - log(2)*n) > (exp((lngamma(k+1) + k) - (log(sqrt(2*Pi)) + log(k)*(k + 0.5))) - 1)