A037160 Squares in A037159.
21609, 27889, 54289, 1461681, 6610041, 8037225, 25250625, 39778249, 50908225, 60637369, 171321921, 207504025, 218123361, 307265841, 315808441, 380367009, 435598641, 437353569, 528126361, 561073969, 581340321, 708890625, 877285161
Offset: 1
Crossrefs
Cf. A037159.
Programs
-
Mathematica
f[0] = 0; f[n_ /; 0 < n < 10^10] := 3*n - DivisorSigma[1, n]; f[] = -1; Reap[For[k = 1, k < 3*10^4, k++, If[FixedPoint[f, k^2] == 0, Print[k^2]; Sow[k^2]]]][[2, 1]] (* _Jean-François Alcover, Sep 17 2013 *)
Extensions
Better description from Jud McCranie, Dec 24 2000
More terms from Sam Handler (sam_5_5_5_0(AT)yahoo.com), Aug 09 2006