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

Author

Leroy Quet, Dec 05 2005

Keywords

• nonn

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

Crossrefs

Cf. A008475, A114519, A114520, A114522.

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