Is the best way to compute the division of large factorials to take the exponential of the log gamma function like so:
$a!/b! = e^{ln(\Gamma(a)) - ln(\Gamma(b))}[\math]<br /> <br /> I'm doing this specifically to avoid overflow and to avoid using any large precision packages. Of course I am using the built in loggamma functions, not separately computing the logarithm of the gamma function.$