login
Sum of digits of n in fractional base 3/2.
18

%I #38 Mar 24 2023 15:54:48

%S 0,1,2,2,3,4,3,4,5,3,4,5,5,6,7,4,5,6,5,6,7,7,8,9,5,6,7,5,6,7,7,8,9,8,

%T 9,10,5,6,7,7,8,9,6,7,8,7,8,9,9,10,11,9,10,11,5,6,7,7,8,9,8,9,10,6,7,

%U 8,8,9,10,8,9,10,9,10,11,11,12,13,10,11,12,5

%N Sum of digits of n in fractional base 3/2.

%C The base 3/2 expansion is unique, and thus the sum of digits function is well-defined.

%C Fixed point starting with 0 of the two-block substitution a,b -> a,a+1,a+2 for a = 0,1,2,... and b = 0,1,2,.... - _Michel Dekking_, Sep 29 2022

%H Reinhard Zumkeller, <a href="/A244040/b244040.txt">Table of n, a(n) for n = 0..10000</a>

%H Michel Dekking, <a href="https://arxiv.org/abs/2301.13563">The Thue-Morse sequence in base 3/2</a>, arXiv:2301.13563 [math.CO], 2023. See also J. Int. Seq., Vol. 26 (2023), <a href="https://cs.uwaterloo.ca/journals/JIS/VOL26/Dekking/dek25.html">Article 23.2.3</a>.

%F a(0)=0, a(3n+r) = a(2n)+r for n >= 0 and r = 0, 1, 2. - _David Radcliffe_, Aug 21 2021

%e In base 3/2 the number 7 is represented by 211 and so a(7) = 2 + 1 + 1 = 4.

%t a[n_]:= a[n]= If[n==0, 0, a[2*Floor[n/3]] + Mod[n,3]]; Table[a[n], {n, 0, 85}] (* _G. C. Greubel_, Aug 20 2019 *)

%o (Sage)

%o def base32sum(n):

%o L, i = [n], 1

%o while L[i-1]>2:

%o x=L[i-1]

%o L[i-1]=x.mod(3)

%o L.append(2*floor(x/3))

%o i+=1

%o return sum(L)

%o [base32sum(n) for n in [0..85]]

%o (Haskell)

%o a244040 0 = 0

%o a244040 n = a244040 (2 * n') + t where (n', t) = divMod n 3

%o -- _Reinhard Zumkeller_, Sep 05 2014

%o (Python) a244040 = lambda n: a244040((n // 3) * 2) + (n % 3) if n else 0 # _David Radcliffe_, Aug 21 2021

%Y Cf. A024629, A007953, A000120, A053735, A244041, A246435.

%K nonn,base

%O 0,3

%A _James Van Alstine_, Jun 17 2014