OFFSET
1,2
COMMENTS
a(8) has 87 decimal digits.
A subsequence of A003215 (centered hexagonal numbers: 3n(n+1)+1, also first differences of A000578). - Klaus Brockhaus, Mar 20 2010
a(11) has 693 digits and is the last term in the b-file. a(12) has 1386 digits and is too large to include in the b-file. - Harvey P. Dale, Jul 31 2019
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..11
FORMULA
a(n) = 1 + 3*a(n-1)*(a(n-1) + 1). - Zak Seidov, Jun 25 2010
EXAMPLE
n = 2: for k = 7, a(1)^3+k = 1^3+7 = 9 = 2^3 is a cube; 7 is the smallest such k, therefore a(2) = 7.
n = 4: for k = 86191, a(3)^3+k = 169^3+86191 = 4913000 = 170^3 is a cube; 86191 is the smallest such k, therefore a(4) = 86191.
MATHEMATICA
NestList[3#^2+3#+1&, 1, 6] (* Harvey P. Dale, Jul 31 2019 *)
PROG
(Magma) /* inefficient, uses definition */ a:=1; S:=[a]; for n in [2..4] do k:=0; flag:= true; while flag do k+:=1; if IsPower(a^3+k, 3) then Append(~S, k); a:=k; flag:=false; end if; end while; end for; S;
/* uses formula from R. J. Mathar, see A172028 */ [ n eq 1 select 1 else 1+3*Self(n-1)*(Self(n-1)+1): n in [1..8] ]; // Klaus Brockhaus, Mar 16 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Vincenzo Librandi, Jan 23 2010
EXTENSIONS
Edited and a(7) added by Klaus Brockhaus, Mar 16 2010
STATUS
approved