Every factorial is a power

The previous post mentioned that 24! 10²⁴ and 25! 10²⁵.

For every n, there is some base b such that n! = bⁿ. For example, 30! 12³⁰.

It's easy to find b [1]:

What's interesting is that b is very nearly a linear function of n.

In hindsight it's clear that this should be the case-it follows easily from Stirling's approximation-but I didn't anticipate this before I plotted it.

Now fix n and find b such that n! = bⁿ. Since the relationship between n and b(n) is nearly linear, this suggests

which is true. It follows from the multiplication identity for the gamma function:

Let z = n + 1/2 so that the left side is (2n)!. On the right side, (z + 1/2) = n! and (z) is not too different from n!. The rest of the right side is 2²ⁿ/.

So our observation that b(n) is nearly linear gave us a hint of Gauss's multiplication formula.

[1] Numerically you would probably evaluate this function by calling a routine that computes log (n + 1) directly without computing(n + 1) first. This avoids overflow for large n.

This is why mathematical libraries will have not only gamma functions but also loggamma functions. The latter seems redundant, but it's not.

The post Every factorial is a power first appeared on John D. Cook.