OFFSET
1,2
FORMULA
a(n+1) = floor((sqrt(8*10^floor(log(a(n))/log(10))*a(n)+1)+1)/2)*(floor((sqrt(8*10^floor(log(a(n))/log(10))*a(n)+1)+1)/2)+1)/2. - Benoit Cloitre, Feb 25 2002
PROG
(PARI) /* with 1000 significant digits : */ a(1)=1, a(2)=10 and for n>2 s=10; for(i=1, 10, s=floor((sqrt(8*10^floor(log(s)/log(10))*s+1)+1)/2)*(floor((sqrt(8*10^floor(log(s)/log(10))*s+1)+1)/2)+1)/2; print1(s, ", "); );
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Amarnath Murthy, Feb 19 2002
EXTENSIONS
More terms from Vladeta Jovovic, Feb 24 2002
More terms from Benoit Cloitre, Mar 07 2002
STATUS
approved