OFFSET
1,1
COMMENTS
Square terms in the sequence are 16, 144, 576, 1600, 3600, 7056, ..., that is, A060300 ((2n(n+1))^2) except 0. And their indices, ind(n), are: 6, 34, 100, 220, 410, 686, ..., that is, ind(n) = 2*A132124(.). - Michel Marcus, suggested by Zak Seidov, Nov 26 2013
LINKS
Zak Seidov, Table of n, a(n) for n = 1..1515 (all terms up to 20000)
E. J. Barbeau, Numbers differing from consecutive squares by squares, Canad. Math. Bull. 28(1985), pp. 337-342.
EXAMPLE
8 = 3^2 - 1^2.
MATHEMATICA
Select[Range[300], IntegerQ[Sqrt[(1+Floor[Sqrt[#]])^2-#]]&] (* Zak Seidov, Nov 26 2013 *)
PROG
(PARI) sq1(n) = {for (i=1, n, a = sqrtint(i) + 1; if (issquare(a^2-i), print1(i, ", ")); ); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Oct 30 2012
STATUS
approved