OFFSET
0,1
COMMENTS
Not to be confused with the Pisot T(2,32) sequence as defined in A008776, which is A013776. - R. J. Mathar, Feb 13 2016
LINKS
FORMULA
Conjecture: a(n) = 16*a(n-1)-8*a(n-3). G.f.: -(x^2-2) / (8*x^3-16*x+1). - Colin Barker, Sep 18 2015
a(n+1) = ceiling(a(n)^2/a(n-1))-1 for all n > 0. a(n+1)/a(n) ~ 15.968627... as n -> oo. - M. F. Hasler, Feb 11 2016
PROG
(PARI) a=[2, 32]; for(n=2, 2000, a=concat(a, ceil(a[n]^2/a[n-1])-1)); A022028(n)=a[n+1] \\ M. F. Hasler, Feb 11 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved