login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A350573
Base-10 version of A038103 interpreted as base 3.
3
0, 1, 2, 3, 6, 7, 11, 12, 21, 22, 23, 24, 31, 41, 43, 62, 145, 146, 149, 152, 2626, 2627, 4370, 8463, 8466, 13562, 20841, 42320, 43171, 46609, 47749, 48150, 52723, 55179, 55180, 55181
OFFSET
1,3
LINKS
MATHEMATICA
B3[n_] := FromDigits[IntegerDigits[n, 3]]
SP[n_] := (b = B3[n]; StringPosition[ToString[B3[b]], ToString[b]])
n=0; t={}; While[n<10^6, If[SP[n]!={}, AppendTo[t, n]]; n++]; t
PROG
(Python)
from sympy.ntheory.digits import digits
from itertools import count, islice, product
def agen(): # generator of terms
yield 0
for d in count(1):
for first in "12":
for rest in product("012", repeat=d-1):
s = first + "".join(rest)
if s in "".join(str(d) for d in digits(int(s), 3)[1:]):
yield int(s, 3)
print(list(islice(agen(), 36))) # Michael S. Branicky, Jan 06 2022
(Python)
from itertools import count, islice
from gmpy2 import digits
def A350573_gen(): return (n for n in count(0) if (s:=digits(n, 3)) in digits(int(s), 3))
A350573_list = list(islice(A350573_gen(), 30)) # Chai Wah Wu, Jan 09 2022
CROSSREFS
Sequence in context: A075995 A102432 A024561 * A274546 A366641 A113545
KEYWORD
nonn,base
AUTHOR
Hans Havermann, Jan 06 2022
STATUS
approved