# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a086593
Showing 1-1 of 1

%I A086593 #46 Jul 14 2018 17:06:03
%S A086593 2,3,4,5,5,7,7,8,6,9,10,11,9,12,11,13,7,11,13,14,13,17,15,18,11,16,17,
%T A086593 19,14,19,18,21,8,13,16,17,17,22,19,23,16,23,24,27,19,26,25,29,13,20,
%U A086593 23,25,22,29,26,31,17,25,27,30,23,31,29,34,9,15,19,20,21,27,23,28,21,30
%N A086593 Bisection of A086592, denominators of the left-hand half of Kepler's tree of fractions.
%C A086593 Also denominator of alternate fractions in Kepler's tree as shown in A294442. - _N. J. A. Sloane_, Nov 20 2017
%H A086593 Antti Karttunen, <a href="/A086593/b086593.txt">Table of n, a(n) for n = 1..2048</a> (computed from b-file of A020650 provided by _T. D. Noe_)
%F A086593 a(n) = A086592(2n-1) = A020650(4n-2).
%F A086593 a(n+1) = A071585(n) + A071766(n), n >= 0. - _Yosu Yurramendi_, Jun 30 2014
%F A086593 From _Yosu Yurramendi_, Jan 04 2016: (Start)
%F A086593 a(2^(m+1)+k+1) - a(2^m+k+1) = A071585(k),  m >= 0, 0 <= k < 2^m.
%F A086593 a(2^(m+2)-k) = a(2^(m+1)-k) + a(2^m-k),   m > 0, 0 <= k < 2^m-1.
%F A086593 (End)
%F A086593 a(2^n) = A000045(n+3). - _Antti Karttunen_, Jan 29 2016, based on above.
%F A086593 a(n) = A020651(4n-1), a(n+1) = A020651(4n+1), n > 0. - _Yosu Yurramendi_, May 08 2018
%F A086593 a(2^m+k) = A071585(2^(m+1)+k),  m >= 0, 0 <= k < 2^m. - _Yosu Yurramendi_, May 16 2018
%t A086593 (* b = A020650 *) b[1] = 1; b[2] = 2; b[3] = 1; b[n_] := b[n] = Switch[ Mod[n, 4], 0, b[n/2 + 1] + b[n/2], 1, b[(n - 1)/2 + 1], 2, b[(n - 2)/2 + 1] + b[(n - 2)/2], 3, b[(n - 3)/2]]; a[1] = 2; a[n_] := b[4 n - 4]; Array[a, 100] (* _Jean-François Alcover_, Jan 22 2016, after _Yosu Yurramendi_'s formula for A020650 *)
%o A086593 (R)
%o A086593   maxlevel <- 15
%o A086593 d <- c(1,2)
%o A086593 for(m in 0:maxlevel)
%o A086593  for(k in 1:2^m) {
%o A086593    d[2^(m+1)    +k] <- d[k] + d[2^m+k]
%o A086593    d[2^(m+1)+2^m+k] <- d[2^(m+1)+k]
%o A086593 }
%o A086593 a <- vector()
%o A086593 for(m in 0:maxlevel) for(k in 0:(2^m-1)) a[2^m+k] <- d[2^(m+1)+k]
%o A086593   a[1:63]
%o A086593 # _Yosu Yurramendi_, May 16 2018
%Y A086593 Cf. A000045, A020650, A020651, A071585, A071766, A086592, A294442.
%K A086593 nonn
%O A086593 1,1
%A A086593 _Antti Karttunen_, Aug 28 2003

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE