Use Inverse Gamma Function in `factorial_inv`

[Rudxain]

This would extend the input and output domain to the Reals instead of just Ints