OFFSET
1,1
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = A038597(n)^2. - M. F. Hasler, Oct 05 2013
PROG
(PARI) is(n)=for(k=sqrtnint(n, 3)+1, (sqrtint(12*n-3)+3)\6, if(ispower(n-k^3, 3), return(issquare(n)))); 0 \\ Charles R Greathouse IV, Oct 28 2013
(PARI) mm=820188; cb=vector(mm); for(i=1, mm, cb[i]=i^3); mb=1420608; v=vector(mb); mx=mb^2; for(i=1, mm-1, for(j=i+1, mm, d=cb[j]-cb[i]; if(d<=mx, if(issquare(d, &r), v[r]=1), next(2)))); c=0; for(n=1, mb, if(v[n]==1, c++; write("b038596.txt", c " " n^2))) \\ Donovan Johnson, Oct 31 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Jud McCranie
STATUS
approved