OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..3557 (terms below 10^8)
EXAMPLE
MATHEMATICA
s[n_] := Module[{f = FactorInteger[n], p, e}, p = f[[;; , 1]]; e = f[[;; , 2]]; Times @@ ((p^(e + 1) - 1)/(p - 1)) - Times @@ (p + 1)]; s[1] = 0; seq[max_] := Module[{v = Select[Union[Array[s, max]], 0 < # <= max &], i}, i = Position[Differences[v], 1] // Flatten; v[[i]]]; seq[10^4]
PROG
(PARI) s(n) = {my(f = factor(n), p, e); prod(i = 1, #f~, p = f[i, 1]; e = f[i, 2]; ((p^(e + 1) - 1)/(p - 1))) - prod(i = 1, #f~, f[i, 1] + 1); }
lista(kmax) = {my(v = select(x -> (x < kmax), Set(vector(kmax, k, s(k))))); for(k=1, #v-1, if(v[k+1] - v[k] == 1, print1(v[k], ", "))); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Apr 18 2023
STATUS
approved