OFFSET
1,1
COMMENTS
EXAMPLE
a(2) = 810 because the prime divisors of 810 are {2, 3, 5}, the prime divisors of 810 + 2 = 812 are {2, 7, 29} and 2^2 + 3^2 + 5^2 = 2 + 7 + 29 = 38, hence 810 is in the sequence.
MATHEMATICA
sk[n_]:=Module[{k=1}, While[Plus@@(First@#&/@FactorInteger[k]^2)!=Plus@@(First@#&/@FactorInteger[n+k]), k++]; k]; Array[sk, 65, 1]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Aug 15 2013
STATUS
approved