A062343 - OEIS

%I #23 Mar 09 2022 22:49:12

%S 17,53,71,107,233,251,431,503,521,701,1061,1151,1223,1511,1601,2141,

%T 2213,2411,3023,3041,3203,3221,4013,4211,5003,5021,6011,6101,7001,

%U 10007,10061,10133,10151,10223,10313,10331,10601,11213,11321,11411

%N Primes whose sum of digits is 8.

%H Alois P. Heinz, <a href="/A062343/b062343.txt">Table of n, a(n) for n = 1..42745</a> (first 500 terms from Vincenzo Librandi)

%F Intersection of A000040 (primes) and A052222 (digit sum 8). - _M. F. Hasler_, Mar 09 2022

%e 1151 belongs to the sequence since it is a prime with sum of digits = 8.

%t Select[Prime[Range[500000]], Total[IntegerDigits[#]]==8 &] (* _Vincenzo Librandi_, Jul 08 2014 *)

%o (Magma) [p: p in PrimesUpTo(20000) | &+Intseq(p) eq 8]; // _Vincenzo Librandi_, Jul 08 2014

%o From _M. F. Hasler_, Mar 09 2022: (Start)

%o (PARI) select( {is_A062343(p, s=8)=sumdigits(p)==s&&isprime(p)}, primes([1, 12345])) \\ 2nd optional parameter for similar sequences with other digit sums.

%o (PARI) {A062343_upto_length(L, s=8, a=List(), u=[10^(L-k)|k<-[1..L]])=forvec(d=[[1, L]|i<-[1..s]], isprime(p=vecsum(vecextract(u, d))) && listput(a, p), 1); Vecrev(a)} \\ (End)

%Y Cf. A000040 (primes), A007953 (sum of digits), A052222 (digit sum = 8).

%Y Cf. A062339 (same for digit sum s = 4), A062341 (s = 5), A062337 (s = 7), A107579 (s = 10), and others listed in A244918 (s = 68).

%Y Subsequence of A062342 (primes with digit sum divisible by 8).

%K nonn,base,easy

%O 1,1

%A _Amarnath Murthy_, Jun 21 2001

%E More terms from Larry Reeves (larryr(AT)acm.org), Jul 06 2001