OFFSET
0,1
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..250
Eric Weisstein's World of Mathematics, Pierce Expansion
Wikipedia, Engel Expansion
EXAMPLE
1/phi = Pi/6 + Pi^2/(6*18) + Pi^3/(6*18*95) + ...
MATHEMATICA
EngelExp[A_, n_] := Join[Array[Pi &, Floor[A]], First@Transpose@
NestList[{Ceiling[Pi/Expand[#[[1]] #[[2]] - 1]], Expand[#[[1]] #[[2]] - 1]/Pi} &, {Ceiling[Pi/(A - Floor[A])], (A - Floor[A])/Pi}, n - 1]]; EngelExp[N[2/(1 + Sqrt[5]), 7!], 20]
CROSSREFS
KEYWORD
nonn
AUTHOR
G. C. Greubel, Dec 25 2016
STATUS
approved