OFFSET
1,3
COMMENTS
A018847 lists primes in this sequence. - M. F. Hasler, May 05 2012
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)
PROG
(PARI) is_A018846(n, t=Vec("012..59.86"))={ apply(x->t[eval(x)+1], n=Vec(Str(n)))==vecextract(n, "-1..1") } \\ M. F. Hasler, May 05 2012
(Python)
from itertools import count, islice, product
def ud(s): return s[::-1].translate({ord('6'):ord('9'), ord('9'):ord('6')})
def A018846gen(): # generator of terms
yield from [0, 1, 2, 5, 8]
for d in count(2):
for first in "125689":
for rest in product("0125689", repeat=d//2-1):
left = first + "".join(rest)
for mid in [[""], ["0", "1", "2", "5", "8"]][d%2]:
yield int(left + mid + ud(left))
print(list(islice(A018846gen(), 54))) # Michael S. Branicky, Jul 09 2022
CROSSREFS
KEYWORD
nonn,base,easy,nice
AUTHOR
STATUS
approved