%I #13 Sep 16 2019 17:23:58

%S 167,227,677,1217,5843,13163,15683,15923,24107,122267,148403,148727,

%T 157307,186023,198413,227597,229253,313997,338267,344273,360293,

%U 372833,451937,464483,469367,517613,525257,541547,603917,732233,742073,991073,1006253,1196873,1219847,1328927

%N Primes p such that all of p - m^2, m = 2, 4, 6, 8, 10, are (positive) primes.

%C All terms are congruent to (17,23) mod 30.

%p isA247941 := proc(p)

%p local m ;

%p for m from 0 to 10 by 2 do

%p if not isprime(p-m^2) then

%p return false;

%p end if;

%p end do:

%p return true;

%p end proc:

%p for n from 1 to 100000 do

%p p := ithprime(n) ;

%p if isA247941(p) then

%p printf("%d,",p);

%p end if;

%p end do: # _R. J. Mathar_, Sep 28 2014

%t Select[Prime[Range[25,103000]],AllTrue[#-(2Range[5])^2,PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Sep 16 2019 *)

%Y Cf. A046132, A246873, A246874 (supersequence).

%K nonn

%O 1,1

%A _Zak Seidov_, Sep 27 2014