log(n**n * exp(-n) * sqrt(2*PI*n))
  = n*log(n) - n + (log(2 + log(PI) + log(n))/2