Evaluate binomial coefficients: Difference between revisions
m
→{{header|Phix}}: replaced with gmp version
Alextretyak (talk | contribs) (Added 11l) |
m (→{{header|Phix}}: replaced with gmp version) |
||
Line 1,754:
</pre>
However errors will creep in should any result or interim value exceed 9,007,199,254,740,992 (on 32-bit), so (and using a different algorithm just for kicks):
{{libheader|
<lang Phix>include builtins\
function
for i=1 to k do
if mpz_fdiv_q_ui(r, r, i)!=0 then ?9/0 end if
-- r = ba_divide(ba_multiply(r,n-i+1),i)
end for
return mpz_get_str(r)
end function
?mpz_binom(5,3)
?
?
?
{{out}}
<pre>
|