The value of factorial() grows so rapidly that you can precompute all the useful values for the range of 32 or 64 bit integers, taking into account numeric overflow, and put the precomputed values into an array.