login
A213912
a(1) = 1; for n>1, a(n) = floor(sqrt(a(n-1))) if that number is not already in the sequence, otherwise a(n) = 3*a(n-1).
3
1, 3, 9, 27, 5, 2, 6, 18, 4, 12, 36, 108, 10, 30, 90, 270, 16, 48, 144, 432, 20, 60, 7, 21, 63, 189, 13, 39, 117, 351, 1053, 32, 96, 288, 864, 29, 87, 261, 783, 2349, 7047, 83, 249, 15, 45, 135, 11, 33, 99, 297, 17, 51, 153, 459, 1377, 37, 111, 333, 999, 31
OFFSET
1,2
COMMENTS
Permutation of natural numbers with inverse A213913.
PROG
(Haskell)
a213912 n = a213912_list !! (n-1)
a213912_list = 1 : f [1] where
f xs@(x:_) = y : f (y : xs) where
y = if z `notElem` xs then z else 3 * x where z = a000196 x
(PARI) lista(nn) = {my(k, v=vector(nn)); v[1]=1; for(n=2, nn, if(vecsearch(vecsort(v), k=sqrtint(v[n-1])), v[n]=3*v[n-1], v[n]=k)); v; } \\ Jinyuan Wang, Jun 26 2020
CROSSREFS
Cf. A114183.
Sequence in context: A114181 A036134 A317502 * A070360 A248131 A070346
KEYWORD
nonn,look
AUTHOR
Reinhard Zumkeller, Mar 05 2013
STATUS
approved