OFFSET
1,2
COMMENTS
EXAMPLE
MATHEMATICA
f[n_] := DivisorSum[n, 1 &, IntegerQ[(1 + Sqrt[1 + 24*#])/6] &]; fm = -1; s = {}; Do[If[(fn = f[n]) > fm, fm = fn; AppendTo[s, n]], {n, 1, 10^5}]; s (* Amiram Eldar, Jul 28 2022 *)
PROG
(PARI) lista(nn) = my(m=0); for (n=1, nn, my(new = sumdiv(n, d, ispolygonal(d, 5))); if (new > m, m = new; print1(n, ", ")); ); \\ Michel Marcus, Jul 28 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Bernard Schott, Jul 28 2022
EXTENSIONS
a(23)-a(24) from David A. Corneth, Jul 28 2022
STATUS
approved