OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
MATHEMATICA
f[e_] := If[OddQ[e], e+2, e+1]; fun[p_, e_] := 1 + (p^f[e] - p)/(p^2 - 1); s[1] = 1; s[n_] := Times @@ (fun @@@ FactorInteger[n]); Select[Range[2*10^5], s[#] == s[#+1] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 11 2021
STATUS
approved