OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..53
Amiram Eldar, Table of n, a(n), A332038(a(n)) for n = 1..53
EXAMPLE
There are 3 solutions to isigma(x) = 24: isigma(14) = isigma(15) = isigma(23) = 24. For all m < 24 there are 2 or fewer solutions to isigma(x) = m, thus 24 is in the sequence.
MATHEMATICA
fun[p_, e_] := Module[{b = IntegerDigits[e, 2], m}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; isigma[1] = 1; isigma[n_] := Times @@ (fun @@@ FactorInteger[n]); m = 10000; v = Table[0, {m}]; Do[i = isigma[k]; If[i <= m, v[[i]]++], {k, 1, m}]; s = {}; vm = -1; Do[If[v[[k]] > vm, vm = v[[k]]; AppendTo[s, k]], {k, 1, m}]; s
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 05 2020
STATUS
approved