OFFSET
1,1
COMMENTS
All primes == 1 mod 18.
MATHEMATICA
Select[Prime[Range[756]], Mod[#, 9]==1&&PrimeQ[#-2]&] (* James C. McMahon, Nov 30 2025 *)
PROG
(PARI) lista(nn) = forprime(p=2, nn, if (isprime(p-2) && ((p % 9) == 1), print1(p, ", "))); \\ Michel Marcus, Nov 12 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Jan 09 2009
STATUS
approved
