|
MATHEMATICA
|
seq[nmax_] := Module[{v = Table[0, {nmax}], i, s = {}, vmax = -1}, Do[i = DivisorSigma[1, n] - n; If[0 < i <= nmax, v[[i]]++], {n, 1, nmax^2}]; Do[If[v[[n]] > vmax, vmax = v[[n]]; If[v[[n - 1]] == 0 && v[[n + 1]] == 0, AppendTo[s, n]]], {n, 2, nmax - 1}]; s]; seq[8000]
|