Greedy algorithm for Egyptian fractions: Difference between revisions

m
m (→‎{{header|Phix}}: bigatom -> mpfr)
Line 1,693:
procedure efrac(integer num, denom)
string prefixfraction = sprintf("%d/%d",{num,denom}),
prefix = ""
if num>=denom then
integer whole = floor(num/denom)
Line 1,700 ⟶ 1,701:
end if
string e = join(egyptian(num, denom)," + ")
printf(1,"%d/%ds -> %s%s\n",{num,denomfraction,prefix,e})
end procedure
Line 1,754 ⟶ 1,755:
43/48 -> 1/2 + 1/3 + 1/16
5/121 -> 1/25 + 1/757 + 1/763309 + 1/873960180913 + 1/1527612795642093418846225
82014/59 -> [34] + 1/8 + 1/95 + 1/14947 + 1/670223480
 
for proper fractions with 1 to 2 digits
7,818

edits