A114518 Numbers n such that A008475(n) is prime.
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
Keywords
Examples
24 = 2^3 * 3 and 2^3 + 3 = 11, which is prime. So 24 is included.
Links
- Nathaniel Johnston, Table of n, a(n) for n = 1..10000
Programs
-
Mathematica
f[n_] := Plus @@ Power @@@ FactorInteger[n]; Select[Range[165], PrimeQ[f[ # ]] &] (* Ray Chandler, Dec 07 2005 *)
-
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)
Extensions
Extended by Ray Chandler and Lambert Herrgesell (zero815(AT)googlemail.com), Dec 07 2005