A114518 - OEIS

%I #12 Oct 11 2019 16:28:02

%S 2,3,5,6,7,10,11,12,13,17,18,19,22,23,24,28,29,31,34,36,37,40,41,43,

%T 47,48,52,53,54,58,59,61,67,71,72,73,76,79,82,83,88,89,97,100,101,103,

%U 107,108,109,112,113,118,127,131,137,139,142,148,149,151,157,160,162,163

%N Numbers n such that A008475(n) is prime.

%H Nathaniel Johnston, <a href="/A114518/b114518.txt">Table of n, a(n) for n = 1..10000</a>

%e 24 = 2^3 * 3 and 2^3 + 3 = 11, which is prime. So 24 is included.

%t f[n_] := Plus @@ Power @@@ FactorInteger[n]; Select[Range[165], PrimeQ[f[ # ]] &] (* _Ray Chandler_, Dec 07 2005 *)

%o (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)

%Y Cf. A008475, A114519, A114520, A114522.

%K nonn

%O 1,1

%A _Leroy Quet_, Dec 05 2005

%E Extended by _Ray Chandler_ and Lambert Herrgesell (zero815(AT)googlemail.com), Dec 07 2005