OFFSET
0,1
COMMENTS
This coincides with the linearly recurrent sequence defined by the expansion of (6 - 5*x^2)/(1 - 5*x - x^2 + 4*x^3) only up to n <= 69. - Bruno Berselli, Feb 11 2016
LINKS
Colin Barker, Table of n, a(n) for n = 0..1000
FORMULA
a(n+1) = floor(a(n)^2/a(n-1))+1 for all n > 0. - M. F. Hasler, Feb 10 2016
MAPLE
A022023 := proc(n)
option remember;
if n <= 1 then
op(n+1, [6, 30]) ;
else
a := procname(n-1)^2/procname(n-2) ;
if type(a, 'integer') then
a+1 ;
else
ceil(a) ;
fi;
end if;
end proc: # R. J. Mathar, Feb 10 2016
PROG
(PARI) a=[6, 30]; for(n=2, 30, a=concat(a, a[n]^2\a[n-1]+1)); a \\ M. F. Hasler, Feb 10 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Double-checked and extended to 3 lines of data by M. F. Hasler, Feb 10 2016
STATUS
approved