OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..333 (terms below 10^10)
EXAMPLE
The prime divisors of 9936 are {2, 3, 23} and the prime divisors of 9937 are {19, 523} => 2^2 + 3^2 + 23^2 = 19 + 523 = 542, hence 9936 is in the sequence.
MATHEMATICA
fQ[n_] := Plus @@ (First@# & /@ FactorInteger[n]^2) == Plus @@ (First@# & /@ FactorInteger[n + 1]); Select[ Range@ 100000, fQ]
PROG
(Magma) [k:k in [2..2500000]| &+PrimeDivisors(k+1) eq &+[PrimeDivisors(k)[i]^2: i in [1..#PrimeDivisors(k)]]]; // Marius A. Burtea, Feb 18 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Aug 15 2013
STATUS
approved