OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := (p^(e + 1) - 1)/(p - 1) - p; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; q[n_] := (s1 = s[n]) > 1 && s1 == s[n + 1]; Select[Range[10^5], q]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 07 2021
STATUS
approved