editing
approved
isA181604 := proc(p)
if isprime(p) and (isprime(p-2) or isprime(p+2)) then
if modp(p, 10) = 3 then
true;
else
false;
end if ;
end if;
end proc: # R. J. Mathar, Feb 14 2017
Select[Prime@Range@700, Mod[ #, 10] == 3 && (PrimeQ[ # - 2] || PrimeQ[ # + 2]) &] [From _]) &] (* _Robert G. Wilson v_, Nov 06 2010] *)
O. Omar E. Pol, <a href="http://www.polprimos.com">Determinacion geometrica de los numeros primos y perfectos</a>
OEIS Server: https://oeis.org/edit/global/2123
_Omar E. Pol (info(AT)polprimos.com), _, Nov 01 2010
OEIS Server: https://oeis.org/edit/global/157
Select[Prime@Range@700, Mod[ #, 10] == 3 && (PrimeQ[ # - 2] || PrimeQ[ # + 2]) &] [From _Robert G. Wilson v (rgwv(AT)rgwv.com), _, Nov 06 2010]
More terms from _Robert G. Wilson v (rgwv(AT)rgwv.com), _, Nov 06 2010
OEIS Server: https://oeis.org/edit/global/156
Twin primes ending in 3.
3, 13, 43, 73, 103, 193, 283, 313, 433, 463, 523, 643, 823, 883, 1033, 1063, 1093, 1153, 1303, 1453, 1483, 1723, 1873, 1933, 2083, 2113, 2143, 2383, 2593, 2713, 2803, 3253, 3373, 3463, 3583, 3673, 3823, 3853, 4003, 4093, 4243, 4273, 4423, 4483, 4723, 4933
1,1
O. E. Pol, <a href="http://www.polprimos.com">Determinacion geometrica de los numeros primos y perfectos</a>
Select[Prime@Range@700, Mod[ #, 10] == 3 && (PrimeQ[ # - 2] || PrimeQ[ # + 2]) &] [From Robert G. Wilson v (rgwv(AT)rgwv.com), Nov 06 2010]
Cf. A001097, A181603, A181605, A181606.
base,nonn,new
Omar E. Pol (info(AT)polprimos.com), Nov 01 2010
More terms from Robert G. Wilson v (rgwv(AT)rgwv.com), Nov 06 2010