OFFSET
1,1
FORMULA
See A006784 for definition of Engel expansion.
PROG
(PARI) { my(x=vector(31), a=vector(#x-1)); x[1]=0.62923321313860758778445822033246963679442280615921;
for (i=1, #a, a[i]=ceil(1/x[i]); x[i+1]=x[i]*a[i]-1); a }
CROSSREFS
KEYWORD
nonn
AUTHOR
Jon Perry, Feb 08 2003
STATUS
approved