OFFSET
2,1
COMMENTS
For almost all n, a(n) = n^2. Exceptions: a(16) = 211, a(36) = 466, a(39) = 1075, a(53) = 1173, a(58) = 919...
LINKS
Jinyuan Wang, Table of n, a(n) for n = 2..1000
PROG
(PARI) a(n) = {k=n+1; while((k^3-1)%(n^3-1)>0, k++); k; } \\ Jinyuan Wang, Jul 29 2019
(Magma) sol:=[]; m:=1; for n in [2..50] do for k in [n+1..3000] do if IsIntegral((k^3-1)/(n^3-1)) then sol[m]:=k; m:=m+1; break; end if; end for; end for; sol; // Marius A. Burtea, Jul 30 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 28 2002
EXTENSIONS
Name clarified by Jon E. Schoenfield, Jul 28 2019
STATUS
approved