login
A114518
Numbers n such that A008475(n) is prime.
6
2, 3, 5, 6, 7, 10, 11, 12, 13, 17, 18, 19, 22, 23, 24, 28, 29, 31, 34, 36, 37, 40, 41, 43, 47, 48, 52, 53, 54, 58, 59, 61, 67, 71, 72, 73, 76, 79, 82, 83, 88, 89, 97, 100, 101, 103, 107, 108, 109, 112, 113, 118, 127, 131, 137, 139, 142, 148, 149, 151, 157, 160, 162, 163
OFFSET
1,1
LINKS
EXAMPLE
24 = 2^3 * 3 and 2^3 + 3 = 11, which is prime. So 24 is included.
MATHEMATICA
f[n_] := Plus @@ Power @@@ FactorInteger[n]; Select[Range[165], PrimeQ[f[ # ]] &] (* Ray Chandler, Dec 07 2005 *)
PROG
(PARI) A008475(n)=local(t); if(n<1, 0, t=factor(n); sum(k=1, matsize(t)[1], t[k, 1]^t[k, 2])); for(i=1, 500, if(isprime(A008475(i)), print1(i, ", "))) (Herrgesell)
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Dec 05 2005
EXTENSIONS
Extended by Ray Chandler and Lambert Herrgesell (zero815(AT)googlemail.com), Dec 07 2005
STATUS
approved