|
|
A032352
|
|
Numbers k such that there is no prime between 10*k and 10*k+9.
|
|
35
|
|
|
20, 32, 51, 53, 62, 84, 89, 107, 113, 114, 126, 133, 134, 135, 141, 146, 150, 164, 167, 168, 171, 176, 179, 185, 189, 192, 196, 204, 207, 210, 218, 219, 232, 236, 240, 248, 249, 251, 256, 258, 282, 294, 298, 305, 309, 314, 315, 317, 323, 324, 326, 328, 342
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Numbers k with property that appending any single decimal digit to k does not produce a prime.
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
m=32: 321=3*107, 323=17*19, 325=5*5*13, 327=3*109, 329=7*47, therefore 32 is a term.
|
|
MAPLE
|
a:=proc(n) if map(isprime, {seq(10*n+j, j=1..9)})={false} then n else fi end: seq(a(n), n=1..350); # Emeric Deutsch, Aug 01 2005
|
|
MATHEMATICA
|
f[n_] := PrimePi[10n + 10] - PrimePi[10n]; Select[ Range[342], f[ # ] == 0 &] (* Robert G. Wilson v, Sep 24 2004 *)
|
|
PROG
|
(Haskell)
a032352 n = a032352_list !! (n-1)
a032352_list = filter
(\x -> all (== 0) $ map (a010051 . (10*x +)) [1..9]) [1..]
(Magma) [n: n in [1..350] | IsZero(#PrimesInInterval(10*n, 10*n+9))]; // Bruno Berselli, Sep 04 2012
(PARI) is(n)=!isprime(10*n+1) && !isprime(10*n+3) && !isprime(10*n+7) && !isprime(10*n+9) \\ Charles R Greathouse IV, Mar 29 2013
(Python)
from sympy import isprime
def aupto(limit):
alst = []
for d in range(2, limit+1):
td = [10*d + j for j in [1, 3, 7, 9]]
if not any(isprime(t) for t in td): alst.append(d)
return alst
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|