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,1

Links

Nathaniel Johnston, Table of n, a(n) for n = 1..10000

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

Cf. A008475, A114519, A114520, A114522.

Sequence in context: A110086 A107746 A024899 * A066940 A305318 A100118

Adjacent sequences: A114515 A114516 A114517 * A114519 A114520 A114521

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