%I #15 Mar 15 2023 15:06:36
%S 21609,27889,54289,1461681,6610041,8037225,25250625,39778249,50908225,
%T 60637369,171321921,207504025,218123361,307265841,315808441,380367009,
%U 435598641,437353569,528126361,561073969,581340321,708890625,877285161
%N Squares in A037159.
%t 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 *)
%Y Cf. A037159.
%K nonn,nice
%O 1,1
%A _Naohiro Nomoto_
%E Better description from _Jud McCranie_, Dec 24 2000
%E More terms from Sam Handler (sam_5_5_5_0(AT)yahoo.com), Aug 09 2006
|