OFFSET
1,1
COMMENTS
If n = p^q, where p is prime and q > 0, then p+q is in the sequence.
If n is not of that form, omit the term.
Might be a good "puzzle" sequence - guess the rule given the first ten or so terms.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
n = 3 = 3^1, so 3+1 = 4 is a term; n = 4 = 2^2, so 2+2 = 4 is again a term; n = 5 = 5^1, so we get 5+1 =6.
But 6 is not a prime power, so we skip it.
MATHEMATICA
fsum[a_] := Total[Flatten[FactorInteger[a]]]; fsum/@Select[Range[242], PrimePowerQ](* James C. McMahon, Jun 08 2024 *)
PROG
(PARI) for(n=1, 300, fac=factor(n); if(matsize(fac)[1]==1, print1(fac[1, 1]+fac[1, 2], ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Dec 04 2005
EXTENSIONS
Edited and extended by Klaus Brockhaus, Jan 21 2006
Further edited by N. J. A. Sloane, Nov 19 2018
STATUS
approved