OFFSET
1,1
COMMENTS
All terms are congruent to (17,23) mod 30.
MAPLE
isA247941 := proc(p)
local m ;
for m from 0 to 10 by 2 do
if not isprime(p-m^2) then
return false;
end if;
end do:
return true;
end proc:
for n from 1 to 100000 do
p := ithprime(n) ;
if isA247941(p) then
printf("%d, ", p);
end if;
end do: # R. J. Mathar, Sep 28 2014
MATHEMATICA
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 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Sep 27 2014
STATUS
approved