login
Rebase n from 3 to 2. Replace 3^k with 2^k in ternary expansion of n.
19

%I #31 Jan 03 2021 14:36:26

%S 0,1,2,2,3,4,4,5,6,4,5,6,6,7,8,8,9,10,8,9,10,10,11,12,12,13,14,8,9,10,

%T 10,11,12,12,13,14,12,13,14,14,15,16,16,17,18,16,17,18,18,19,20,20,21,

%U 22,16,17,18,18,19,20,20,21,22,20,21,22,22,23,24,24,25,26,24,25,26,26,27

%N Rebase n from 3 to 2. Replace 3^k with 2^k in ternary expansion of n.

%C Notation: (3)[n](2).

%C Fixed point of the morphism 0->0,1,2; 1->2,3,4; 2->4,5,6; ...; n->2n,2n+1,2n+2. - _Philippe Deléham_, Oct 22 2011

%H Harry J. Smith, <a href="/A065361/b065361.txt">Table of n, a(n) for n = 0..1000</a>

%F a(0)=0, a(3n)=2*a(n), a(3n+1)=2*a(n)+1, a(3n+2)=2*a(n)+2. - _Benoit Cloitre_, Dec 21 2002

%F a(n) = 2*a(floor(n/3))+n-3*floor(n/3). - _Benoit Cloitre_, Apr 27 2003

%F a(n) = Sum_{k>=0} A030341(n,k)*2^k. - _Philippe Deléham_, Oct 22 2011

%e 15 = 120 -> 1(4)+2(2)+0(1) = 8 = a(15).

%t t = Table[FromDigits[RealDigits[n, 3], 2], {n, 0, 100}]

%t (* _Clark Kimberling_, Aug 02 2012 *)

%o (PARI) a(n)=if(n<1,0,if(n%3,a(n-1)+1,2*a(n/3)))

%o (PARI) a(n)=if(n<1,0,2*a(floor(n/3))+n-3*floor(n/3))

%o (PARI) Rebase(x, b, c)= { local(d, e=0, f=1); while (x>0, d=x-b*(x\b); x\=b; e+=d*f; f*=c); return(e) } { for (n=0, 1000, write("b065361.txt", n, " ", Rebase(n, 3, 2)) ) } \\ _Harry J. Smith_, Oct 17 2009

%o (Julia)

%o function a(n)

%o m, r, b = n, 0, 1

%o while m > 0

%o m, q = divrem(m, 3)

%o r += b * q

%o b *= 2

%o end

%o r end; [a(n) for n in 0:76] |> println # _Peter Luschny_, Jan 03 2021

%Y Cf. A065362, A030341.

%K base,easy,nonn

%O 0,3

%A _Marc LeBrun_, Oct 31 2001