%I #10 Feb 05 2014 17:08:26
%S 4,10,12,14,15,18,20,21,24,26,28,33,35,36,38,39,44,45,48,50,51,52,56,
%T 62,63,65,69,72,75,77,80,86,91,92,93,95,98,104,108,111,112,117,122,
%U 123,129,133,135,144,146,147,148,152,159,161,162,171,172,175,176,188,189
%N n - (sum of prime factors of n) is prime.
%C Prime factors taken without multiplicity. - _Harvey P. Dale_, Feb 04 2014
%H Vincenzo Librandi, <a href="/A071179/b071179.txt">Table of n, a(n) for n = 1..2000</a>
%e 20 = 2^2.5 and 20-(2+5) = 13 is prime hence 20 is in the sequence
%t Select[Range[200],PrimeQ[#-Total[Transpose[FactorInteger[#]][[1]]]]&] (* _Harvey P. Dale_, Feb 04 2014 *)
%o (PARI) for(n=1,100,if(isprime(n-sumdiv(n,d,isprime(d)*d))==1,print1(n,",")))
%K easy,nonn
%O 1,1
%A _Benoit Cloitre_, Jun 10 2002
|