OFFSET
1,1
COMMENTS
Or, primes in A033075. - Zak Seidov, Feb 01 2011
REFERENCES
J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 67, p. 23, Ellipses, Paris 2008.
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000 (terms n = 1..1500 from Zak Seidov)
MATHEMATICA
Select[Prime[Range[10000]], # < 10 || Union[Abs[Differences[IntegerDigits[#]]]] == {1} &]
PROG
(Haskell)
a048398 n = a048398_list !! (n-1)
a048398_list = filter ((== 1) . a010051') a033075_list
-- Reinhard Zumkeller, Feb 21 2012, Nov 04 2010
(Python 3.2 or higher)
from itertools import product, accumulate
from sympy import isprime
A048398_list = [2, 3, 5, 7]
for l in range(1, 17):
for d in [1, 3, 7, 9]:
dlist = [d]*l
for elist in product([-1, 1], repeat=l):
flist = [str(d+e) for d, e in zip(dlist, accumulate(elist)) if 0 <= d+e < 10]
if len(flist) == l and flist[-1] != '0':
n = 10*int(''.join(flist[::-1]))+d
if isprime(n):
A048398_list.append(n)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Patrick De Geest, Apr 15 1999
STATUS
approved