login
Primes p such that sum of digits + 9 is prime.
2

%I #12 Sep 08 2022 08:46:08

%S 2,11,13,17,19,31,37,53,59,71,73,101,103,107,109,127,149,163,167,181,

%T 211,233,239,251,257,271,293,307,347,383,389,419,431,433,479,491,499,

%U 503,509,521,523,541,563,569,587,613,617,631,653,659,677,701,743

%N Primes p such that sum of digits + 9 is prime.

%H Robert Israel and Vincenzo Librandi, <a href="/A243588/b243588.txt">Table of n, a(n) for n = 1..10000</a> (n = 1..1000 from Vincenzo Librandi)

%p filter:= n -> isprime(convert(convert(n,base,10),`+`)+9):

%p select(isprime and filter,[$1..1000]); # _Robert Israel_, Jun 08 2014

%t Select[Prime[Range[200]], PrimeQ[Plus@@IntegerDigits[#] + 9] &]

%o (Magma) [p: p in PrimesUpTo(1500) | IsPrime(q) where q is 9+&+Intseq(p)];

%Y Cf. similar sequences listed in A243586.

%K nonn,base

%O 1,1

%A _Vincenzo Librandi_, Jun 07 2014