OFFSET
1,1
COMMENTS
No more terms below 10^10. - Amiram Eldar, Aug 23 2019
EXAMPLE
24 + product of prime factors of 24 = 24 + 2 * 3 = 30; 25 + product of prime factors of 25 = 25 + 5 = 30; hence 24 belongs to the sequence.
MATHEMATICA
s[n_] := n + Apply[Times, Transpose[FactorInteger[n]][[1]]]; Select[Range[2, 10^5], s[ # ] == s[ # + 1] &]
Flatten[Position[Partition[Table[n+Times@@Transpose[FactorInteger[ n]] [[1]], {n, 2, 100000}], 2, 1], _?(#[[1]]==#[[2]]&), {1}, Heads->False]]+1 (* Harvey P. Dale, Apr 20 2014 *)
rad[n_] := Times @@ First /@ FactorInteger[n]; r1 = 1; s = {}; Do[r2 = rad[n]; If[r1 - r2 == 1, AppendTo[s, n-1]]; r1 = r2, {n, 2, 10^5}]; s (* Amiram Eldar, Aug 23 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Joseph L. Pe, Oct 12 2002
EXTENSIONS
a(7)-a(9) from Amiram Eldar, Aug 23 2019
STATUS
approved