OFFSET
1,3
COMMENTS
This sequence is equivalent to numbers k such that sigma(k) > sigma(m) for all m < k where m and k are not multiples of 4 as no terms are multiples of 4 and A000118(n) = 8*sigma(n) for nonmultiples of 4. - David A. Corneth, May 19 2023
LINKS
David A. Corneth, Table of n, a(n) for n = 1..392 (terms <= 10^8)
MATHEMATICA
f[n_] := SquaresR[4, n];
Module[{record = -1, n}, Reap[For[n = 0, n <= 2000, n++, If[f[n] > record, record = f[n]; Sow[n]]]][[2, 1]]] (* Jean-François Alcover, May 18 2023 *)
DeleteDuplicates[Table[{n, SquaresR[4, n]}, {n, 0, 3600}], GreaterEqual[#1[[2]], #2[[2]]]&][[;; , 1]] (* Harvey P. Dale, Jun 20 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Vos Post, Mar 21 2007
EXTENSIONS
Missing term 378 added by Giovanni Resta, Jun 18 2016
More terms from Jean-François Alcover, May 18 2023
More terms from David A. Corneth, May 18 2023
STATUS
approved