OFFSET
0,2
PROG
(Haskell)
a261788 n = a261788_list !! (n-1)
a261788_list = f 1 1 a261786_list' where
f z k (x:xs) | x >= z = k : f (3 * z) (k + 1) xs
| otherwise = f z (k + 1) xs
(PARI) ts(n) = Str(fromdigits(digits(n, 3))); \\ A007089
f(n) = my(s=ts(n), k=1); while (#strsplit(s, ts(k)) != 1, k++); k; \\ A261787
lista(nn) = my(last=1, k=0, p=3^k); for (n=1, nn, if (last >= p, print1(n, ", "); k++; p = 3^k); last += f(last); ); \\ Michel Marcus, Feb 06 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Reinhard Zumkeller, Sep 01 2015
EXTENSIONS
a(11)-a(20) from Michel Marcus, Feb 06 2022
STATUS
approved