login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A245657 Primes p for which none of the concatenations p3, p9, 3p, 9p are primes. 1
3, 107, 113, 179, 317, 443, 487, 599, 641, 653, 751, 773, 937, 977, 991, 1021, 1087, 1103, 1187, 1201, 1213, 1217, 1301, 1409, 1427, 1439, 1483, 1553, 1559, 1579, 1609, 1637, 1693, 1747, 1777, 1787, 1789, 1861, 1949, 1987, 1993, 2081, 2129, 2239, 2281, 2287, 2293, 2351, 2393, 2477 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
LINKS
MATHEMATICA
Select[Prime[Range[400]], NoneTrue[{10#+3, 10#+9, 3*10^IntegerLength[#]+#, 9*10^IntegerLength[ #]+#}, PrimeQ]&] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Sep 06 2020 *)
PROG
(PARI) lista(nn) = {forprime(p=2, nn, if (!isprime(eval(concat(Str(p), Str(3)))) && ! isprime(eval(concat(Str(p), Str(9)))) && ! isprime(eval(concat(Str(3), Str(p)))) && ! isprime(eval(concat(Str(9), Str(p)))), print1(p, ", ")); ); } \\ Michel Marcus, Sep 14 2014
(Python)
import sympy
from sympy import isprime
from sympy import prime
for n in range(1, 10**3):
..p = str(prime(n))
..if not isprime(p+'3') and not isprime(p+'9') and not isprime('3'+p) and not isprime('9'+p):
....print(int(p), end=', ') # Derek Orr, Sep 16 2014
CROSSREFS
Sequence in context: A343214 A234636 A213862 * A139921 A176807 A301369
KEYWORD
nonn,base,easy
AUTHOR
Vladimir Shevelev, Sep 13 2014
EXTENSIONS
More terms from Derek Orr, Sep 16 2014
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 12:15 EDT 2024. Contains 371969 sequences. (Running on oeis4.)