|
%I #13 Feb 18 2019 13:44:36
%S 3,5,7,11,21,23,25,29,41,43,47,49,61,65,67,83,85,89,101,111,113,115,
%T 119,131,133,137,139,151,155,157,173,175,179,191,193,197,199,201,203,
%U 205,209,221,223,227,229,241,245,247,263,265,269,281,283,287,289,311
%N Odd numbers whose decimal digits sum to a prime.
%H Nathaniel Johnston, <a href="/A186773/b186773.txt">Table of n, a(n) for n = 1..10000</a>
%p A186773 := proc(n) option remember: local d,k: if(n=1)then return 3: fi: for k from procname(n-1)+2 by 2 do d:=convert(k,base,10): if(isprime(add(d[j],j=1..nops(d))))then return k: fi: od: end: seq(A186773(n),n=1..56); # _Nathaniel Johnston_, May 27 2011
%t Select[Range[3,311,2],PrimeQ[Total[IntegerDigits[#]]]&] (* _Harvey P. Dale_, Oct 20 2011 *)
%Y Cf. A186647, A119450.
%K nonn,base,easy
%O 1,1
%A _Giovanni Teofilatto_, Feb 26 2011
%E Corrected by _Nathaniel Johnston_, May 27 2011
|
