In reply to @icecolbeveridge
@icecolbeveridge you're right, I only need an approximation!
In reply to @christianp
@icecolbeveridge but Stirling's approximation doesn't seem to be close enough, unless I've misunderstood. And Python doesn't like making massive floats, so I have to be clever twice