cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 13 results. Next

A294444 Number of distinct numbers appearing as denominators in row n of Kepler's triangle A294442.

Original entry on oeis.org

1, 1, 1, 2, 3, 6, 10, 16, 29, 51, 83, 148, 246, 402, 650, 1084, 1740, 2803, 4458
Offset: 0

Views

Author

N. J. A. Sloane, Nov 20 2017

Keywords

Comments

It would be nice to have a formula or recurrence.

Examples

			Row 4 of A294442 contains eight fractions, 1/5, 4/5, 3/7, 4/7, 2/7, 2/7, 3/8, 5/8.
There are three distinct denominators, so a(4) = 3.
		

Crossrefs

Cf. A294442.
See A293160 for the number of distinct numerators, or for numerators or denominators in the Stern-Brocot triangle A002487.

Programs

  • Maple
    # S[n] is the list of fractions, written as pairs [i, j], in row n of Kepler's triangle; nc is the number of distinct numerators, and dc the number of distinct denominators
    S[0]:=[[1,1]]; S[1]:=[[1,2]];
    nc:=[1,1]; dc:=[1,1];
    for n from 2 to 18 do
    S[n]:=[];
    for k from 1 to nops(S[n-1]) do
    t1:=S[n-1][k];
    a:=[t1[1],t1[1]+t1[2]];
    b:=[t1[2],t1[1]+t1[2]];
    S[n]:=[op(S[n]),a,b];
    od:
    listn:={};
    for k from 1 to nops(S[n]) do listn:={op(listn), S[n][k][1]}; od:
    c:=nops(listn);  nc:=[op(nc),c];
    listd:={};
    for k from 1 to nops(S[n]) do listd:={op(listd), S[n][k][2]}; od:
    c:=nops(listd);  dc:=[op(dc),c];
    od:
    nc; # A293160
    dc; # this sequence

A020651 Denominators in recursive bijection from positive integers to positive rationals (the bijection is f(1) = 1, f(2n) = f(n)+1, f(2n+1) = 1/(f(n)+1)).

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 3, 1, 4, 3, 4, 2, 5, 3, 5, 1, 5, 4, 5, 3, 7, 4, 7, 2, 7, 5, 7, 3, 8, 5, 8, 1, 6, 5, 6, 4, 9, 5, 9, 3, 10, 7, 10, 4, 11, 7, 11, 2, 9, 7, 9, 5, 12, 7, 12, 3, 11, 8, 11, 5, 13, 8, 13, 1, 7, 6, 7, 5, 11, 6, 11, 4, 13, 9, 13, 5, 14, 9, 14, 3, 13, 10, 13, 7, 17, 10, 17, 4, 15, 11, 15, 7, 18, 11
Offset: 1

Views

Author

Keywords

Comments

If we insert an initial 1, this is the sequence of numerators in left-hand half of Kepler's tree of fractions. Form a tree of fractions by beginning with 1/1 and then giving every node i/j two descendants labeled i/(i+j) and j/(i+j). See A086592 for denominators. See A294442 for the Kepler tree itself.
Level n of the tree consists of 2^n nodes: 1/2; 1/3, 2/3; 1/4, 3/4, 2/5, 3/5; 1 /5, 4/5, 3/7, 4/7, 2/7, 5/7, 3/8, 5/8; ... Fibonacci numbers occur at the right edge this tree, i.e., a(A000225(n)) = A000045(n+1). The fractions are given in their reduced form, thus gcd(A020650(n), A020651(n)) = 1 and gcd(A020651(n), A086592(n)) = 1 for all n. - Antti Karttunen, May 26 2004
A generalization which includes the "rabbit tree" (A226080) and "all rationals tree" (A226130) follows. Suppose that a,b,c,d,e,f,g,h are complex numbers. Let S be the set of numbers defined by these rules: (1) 1 is in S; (2) if x is in S and cx+d is not 0, then U(x) = (ax+b)/(cx+d) is in S; (3) if x is in S and gx+h is not 0, then D(x) = (ex+f)/(gx+h) is in S. If an infinite path in the resulting tree has convergent nodes, then there is some node after which the path is "updown zigzag" ((UoD)o(UoD)o ...) or "downup zigzag" (DoU)o(DoU)o ...). If ag+ch is not 0, then the updown zigzag limit is invariant of x and equals [ae + cf - bg - dh + sqrt(X)]/(2(ag + ch)), where X = (ae + cf - bg - dh)^2 + 4(be + df + ag + ch). If ce + dg is not 0, then the downup zigzag limit is invariant of x and equals [ae + bg - cf - dh + sqrt(Y)]/(2(ce + dg)), where Y = (ae + bg - cf - dh)^2 + 4(af + bh)(ce + dg)) = X. Thus, for the tree A020651, the updown zigzag limit is -1 + sqrt(2) and the downup zigzag limit, sqrt(2). - Clark Kimberling, Nov 10 2013
From Yosu Yurramendi, Jul 13 2014: (Start)
If the terms (n > 0) are written as an array (left-aligned fashion) with rows of length 2^m, m = 0,1,2,3,...
1,
1,2,
1,3,2,3,
1,4,3,4,2,5,3,5,
1,5,4,5,3,7,4,7,2, 7,5, 7,3, 8,5, 8,
1,6,5,6,4,9,5,9,3,10,7,10,4,11,7,11,2,9,7,9,5,12,7,12,3,11,8,11,5,13,8,13,
then the sum of the m-th row is 3^m (m = 0,1,2,), and each column is an arithmetic sequence. The differences of the arithmetic sequences, except the first on the left, give the sequence A093873 (Numerators in Kepler's tree of harmonic fractions) (a(2^(m+1)-1-k) - a(2^m-1-k) = A093873(k), m = 0,1,2,..., k = 0,1,2,...,2^m-1).
If the rows are written in a right-aligned fashion:
1,
1, 2,
1, 3,2, 3,
1, 4,3, 4,2, 5,3, 5,
1,5,4,5,3, 7,4, 7,2, 7,5, 7,3, 8,5, 8,
1,6,5,6,4,9,5,9,3,10,7,10,4,11,7,11,2,9,7,9,5,12,7,12,3,11,8,11,5,13,8,13,
then each column k is a Fibonacci sequence. (End)
For m >= 0, a(2^m) = 1 and a(3*2^m) = 2. For n >= 0, a(A070875(n)) = 3 (for m >= 0, a(5*2^m) = 3 and a(7*2^m) = 3). - Yosu Yurramendi, Jun 02 2016

Examples

			1, 2, 1/2, 3, 1/3, 3/2, 2/3, 4, 1/4, 4/3, ...
		

Crossrefs

See A294442 and A093873/A093875 for two different versions of the Kepler tree.

Programs

  • Haskell
    import Data.List (transpose); import Data.Ratio (denominator)
    a020651_list = map denominator ks where
       ks = 1 : concat (transpose [map (+ 1) ks, map (recip . (+ 1)) ks])
    -- Reinhard Zumkeller, Feb 22 2014
    
  • Maple
    A020651 := n -> `if`((n < 2),n,`if`(type(n,even), A020651(n/2), A020650(n-1)));
  • Mathematica
    f[1] = 1; f[n_?EvenQ] := f[n] = f[n/2]+1; f[n_?OddQ] := f[n] = 1/(f[(n-1)/2]+1); a[n_] := Denominator[f[n]]; Table[a[n], {n, 1, 94}] (* Jean-François Alcover, Nov 22 2011 *)
  • R
    N <- 25 # arbitrary
    a <- c(1,1,2)
    for(n in 1:N){
      a[4*n]   <- a[2*n]
      a[4*n+1] <- a[2*n] + a[2*n+1]
      a[4*n+2] <-          a[2*n+1]
      a[4*n+3] <- a[2*n] + a[2*n+1]
    }
    a
    # Yosu Yurramendi, Jul 13 2014

Formula

a(1) = 1, a(2n) = a(n), a(2n+1) = A020650(2n). - Antti Karttunen, May 26 2004
a(2n) = A020650(2n+1). - Yosu Yurramendi, Jul 17 2014
a(2^m + k) = A093873(2^(m+1) + k) = A093873(2^(m+1) + 2^m + k), m >= 0, 0 <= k < 2^m. - Yosu Yurramendi, May 18 2016
a(2^m + 2^r + k) = A093873(2^r + k)*(m-(r-1)) + A093873(k), m >= 0, r <= m-1, 0 <= k < 2^r. For k=0 A093873(0) = 0 is needed. - Yosu Yurramendi, Jul 30 2016
a((2n+1)*2^m) = A086592(n), m >= 0, n > 0. For n = 0 A086592(0) = 1 is needed. - Yosu Yurramendi, Feb 14 2017
a(4n+2) = a(4n+1) - a(4n) = a(2n+1) = a(4n+1) - a(n), n > 0. - Yosu Yurramendi, May 08 2018
a(1) = 1, a(n+1) = 2*floor(1/a(n))+1-1/a(n). - Jan Malý, Jul 30 2019
a(n) = A002487(A231551(n)), n > 0. - Yosu Yurramendi, Jul 15 2021

Extensions

Entry revised by N. J. A. Sloane, May 24 2004

A093873 Numerators in Kepler's tree of harmonic fractions.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 3, 2, 3, 1, 3, 2, 3, 1, 4, 3, 4, 2, 5, 3, 5, 1, 4, 3, 4, 2, 5, 3, 5, 1, 5, 4, 5, 3, 7, 4, 7, 2, 7, 5, 7, 3, 8, 5, 8, 1, 5, 4, 5, 3, 7, 4, 7, 2, 7, 5, 7, 3, 8, 5, 8, 1, 6, 5, 6, 4, 9, 5, 9, 3, 10, 7, 10, 4, 11, 7, 11, 2, 9, 7, 9, 5, 12, 7, 12, 3, 11, 8, 11, 5, 13, 8, 13, 1, 6
Offset: 1

Views

Author

Keywords

Comments

Form a tree of fractions by beginning with 1/1 and then giving every node i/j two descendants labeled i/(i+j) and j/(i+j).

Examples

			The first few fractions are:
1 1 1 1 2 1 2 1 3 2 3 1 3 2 3 1 4 3 4 2 5 3 5 1 4 3 4 2 5 3 5
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ...
1 2 2 3 3 3 3 4 4 5 5 4 4 5 5 5 5 7 7 7 7 8 8 5 5 7 7 7 7 8 8
		

Crossrefs

The denominators are in A093875. Usually one only considers the left-hand half of the tree, which gives the fractions A020651/A086592. See A086592 for more information, references to Kepler, etc.
See A294442 for another version of Kepler's tree of fractions.

Programs

  • Haskell
    {-# LANGUAGE ViewPatterns #-}
    import Data.Ratio((%), numerator, denominator)
    rat :: Rational -> (Integer,Integer)
    rat r = (numerator r, denominator r)
    data Harmony = Harmony Harmony Rational Harmony
    rows :: Harmony -> [[Rational]]
    rows (Harmony hL r hR) = [r] : zipWith (++) (rows hL) (rows hR)
    kepler :: Rational -> Harmony
    kepler r = Harmony (kepler (i%(i+j))) r (kepler (j%(i+j)))
               where (rat -> (i,j)) = r
    -- Full tree of Kepler's harmonic fractions:
    k = rows $ kepler 1 :: [[Rational]] -- as list of lists
    h = concat k :: [Rational] -- flattened
    a093873 n = numerator $ h !! (n - 1)
    a093875 n = denominator $ h !! (n - 1)
    a011782 n = numerator $ (map sum k) !! n -- denominator == 1
    -- length (k !! n) == 2^n
    -- numerator $ (map last k) !! n == fibonacci (n + 1)
    -- denominator $ (map last k) !! n == fibonacci (n + 2)
    -- numerator $ (map maximum k) !! n == n
    -- denominator $ (map maximum k) !! n == n + 1
    -- eop.
    -- Reinhard Zumkeller, Oct 17 2010
  • Maple
    M:= 8: # to get a(1) .. a(2^M-1)
    gen[1]:= [1];
    for n from 2 to M do
      gen[n]:= map(t -> (numer(t)/(numer(t)+denom(t)),
          denom(t)/(numer(t)+denom(t))), gen[n-1]);
    od:
    seq(op(map(numer,gen[i])),i=1..M): # Robert Israel, Jan 11 2016
  • Mathematica
    num[1] = num[2] = 1; den[1] = 1; den[2] = 2; num[n_?EvenQ] := num[n] = num[n/2]; den[n_?EvenQ] := den[n] = num[n/2] + den[n/2]; num[n_?OddQ] := num[n] = den[(n-1)/2]; den[n_?OddQ] := den[n] = num[(n-1)/2] + den[(n-1)/2]; A093873 = Table[num[n], {n, 1, 97}] (* Jean-François Alcover, Dec 16 2011 *)

Formula

a(n) = a([n/2])*(1 - n mod 2) + A093875([n/2])*(n mod 2).
a(A029744(n-1)) = 1; a(A070875(n-1)) = 2; a(A123760(n-1)) = 3. - Reinhard Zumkeller, Oct 13 2006
A011782(k) = SUM(a(n)/A093875(n): 2^k<=n<2^(k+1)), k>=0. [Reinhard Zumkeller, Oct 17 2010]
a(1) = 1. For all n>0 a(2n) = a(n), a(2n+1) = A093875(n). - Yosu Yurramendi, Jan 09 2016
a(4n+3) = a(4n+1), a(4n+2) = a(4n+1) - a(4n), a(4n+1) = A071585(n). - Yosu Yurramendi, Jan 11 2016
G.f. G(x) satisfies G(x) = x + (1+x) G(x^2) + Sum_{k>=2} x (1+x^(2^(k-1))) G(x^(2^k)). - Robert Israel, Jan 11 2016
a(2^(m+1)+k) = a(2^(m+1)+2^m+k) = A020651(2^m+k), m>=0, 0<=k<2^m. - Yosu Yurramendi, May 18 2016
a(k) = A020651(2^(m+1)+k) - A020651(2^m+k), m>0, 0Yosu Yurramendi, Jun 01 2016
a(2^(m+1)+k) - a(2^m+k) = a(k) , m >=0, 0 <= k < 2^m. For k=0 a(0)=0 is needed. - Yosu Yurramendi, Jul 22 2016
a(2^(m+2)-1-k) = a(2^(m+1)-1-k) + a(2^m-1-k), m >= 1, 0 <= k < 2^m. - Yosu Yurramendi, Jul 22 2016
a(2^m-1-(2^r -1)) = A000045(m-r), m >= 1, 0 <= r <= m-1. - Yosu Yurramendi, Jul 22 2016
a(2^m+2^r) = m-r, , m >= 1, 0 <= r <= m-1 ; a(2^m+2^r+2^(r-1)) = m-(r-1), m >= 2, 0 <= r <= m-1. - Yosu Yurramendi, Jul 22 2016
A093875(2n) - a(2n) = A093875(n), n > 0; A093875(2n+1) - a(2n+1) = a(n), n > 0. - Yosu Yurramendi, Jul 23 2016

A293160 Number of distinct terms in row n of Stern's diatomic array, A049456.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 13, 20, 31, 48, 78, 118, 191, 300, 465, 734, 1175, 1850, 2926, 4597, 7296, 11552, 18278, 28863, 45832, 72356, 114742, 181721, 287926, 455748, 722458, 1144370, 1813975, 2873751, 4553643, 7213620, 11432169, 18120733, 28716294, 45491133
Offset: 0

Views

Author

N. J. A. Sloane, Oct 12 2017, answering a question raised by Barry Carter in an email message. Barry Carter worked out the first 26 terms

Keywords

Comments

Equivalently, a(n) is the number of distinct terms in row n of the Stern-Brocot sequence (A002487) when that sequence is divided into blocks of lengths 1, 2, 4, 8, 16, 32, ...
It would be nice to have a formula or recurrence, or even some bounds. Empirically, a(n) seems to be roughly 2^(2n/3) for the known values. Note that the first half of row n has about 2^(n-2) terms, and the maximal multiplicity is given by A293957(n), so 2^(n-2)/A293957(n) is a lower bound on a(n), which seems not too bad for the known values. - N. J. A. Sloane, Nov 04 2017
The multiset of terms in row n of Stern's diatomic array, with unique elements counted by a(n), is the same as the multiset of numerators of fractions in row n of Kepler's tree. Indeed, a fraction p/q is in row n-1 of Kepler's tree if and only if p/q and q/p are in row n of Calkin-Wilf tree. To form row n of Stern's diatomic array, one should take either numerators and denominators of fractions less than 1 or all numerators from Calkin-Wilf tree row n, in either case p/q and q/p contribute p and q. In Kepler's tree, a fraction p/q contributes p and q as numerators to the next row, i.e. row n. See A294442 for Kepler's tree and A294444 for the number of distinct denominators in it. - Andrey Zabolotskiy, Dec 05 2024

Examples

			Row 4 of A294442 contains eight fractions: 1/5, 4/5, 3/7, 4/7, 2/7, 2/7, 3/8, 5/8.
There are five distinct numerators, so a(4) = 5.
		

Crossrefs

See A135510 for the smallest positive missing number in each row.
Cf. A294442, A294444, A295783 (first differences).

Programs

  • Maple
    A049456 := proc(n, k)
        option remember;
        if n =1 then
            if k >= 0 and k <=1 then
                1;
            else
                0 ;
            end if;
        elif type(k, 'even') then
            procname(n-1, k/2) ;
        else
            procname(n-1, (k+1)/2)+procname(n-1, (k-1)/2) ;
        end if;
    end proc: # R. J. Mathar, Dec 12 2014
    # A293160. This is not especially fast, but it will easily calculate the first 26 terms and confirm Barry Carter's values.
    rho:=n->[seq(A049456(n,k),k=0..2^(n-1))];
    w:=n->nops(convert(rho(n),set));
    [seq(w(n),n=1..26)];
    # Alternative program:
    # S[n] is the list of fractions, written as pairs [i, j], in row n of Kepler's triangle; nc is the number of distinct numerators, and dc the number of distinct denominators
    S[0]:=[[1, 1]]; S[1]:=[[1, 2]];
    nc:=[1, 1]; dc:=[1, 1];
    for n from 2 to 18 do
    S[n]:=[];
    for k from 1 to nops(S[n-1]) do
    t1:=S[n-1][k];
    a:=[t1[1], t1[1]+t1[2]];
    b:=[t1[2], t1[1]+t1[2]];
    S[n]:=[op(S[n]), a, b];
    od:
    listn:={};
    for k from 1 to nops(S[n]) do listn:={op(listn), S[n][k][1]}; od:
    c:=nops(listn);  nc:=[op(nc), c];
    listd:={};
    for k from 1 to nops(S[n]) do listd:={op(listd), S[n][k][2]}; od:
    c:=nops(listd);  dc:=[op(dc), c];
    od:
    nc; # this sequence
    dc; # A294444
    # N. J. A. Sloane, Nov 20 2017
  • Mathematica
    Length[Union[#]]& /@ NestList[Riffle[#, Total /@ Partition[#, 2, 1]]&, {1, 1}, 26] (* Jean-François Alcover, Mar 25 2020, after Harvey P. Dale in A049456 *)
    Map[Length@ Union@ Numerator@ # &, #] &@ Nest[Append[#, Flatten@ Map[{#1/(#1 + #2), #2/(#1 + #2)} & @@ {Numerator@ #, Denominator@ #} &, Last@ #]] &, {{1/1}, {1/2}}, 21] (* Michael De Vlieger, Apr 18 2018 *)
  • Python
    from itertools import chain, product
    from functools import reduce
    def A293160(n): return n if n <= 1 else len({1}|set(sum(reduce(lambda x,y:(x[0],x[0]+x[1]) if y else (x[0]+x[1],x[1]),chain(k,(1,)),(1,0))) for k in product((False,True),repeat=n-2))) # Chai Wah Wu, Jun 20 2022

Extensions

a(28)-a(39) from Don Reble, Oct 16 2017
a(0) prepended and content related to Kepler's tree added from a duplicate entry (where the terms up to a(28) have been independently obtained by Michael De Vlieger) by Andrey Zabolotskiy, Dec 06 2024

A088208 Table read by rows where T(0,0)=1; n-th row has 2^n terms T(n,j),j=0 to 2^n-1. For j==0 mod 2, T(n+1,2j)=T(n,j) and T(n+1,2j+1)=T(n,j)+2^n. For j==1 mod 2, T(n+1,2j+1)=T(n,j) and T(n+1,2j)=T(n,j)+2^n.

Original entry on oeis.org

1, 1, 2, 1, 3, 4, 2, 1, 5, 7, 3, 4, 8, 6, 2, 1, 9, 13, 5, 7, 15, 11, 3, 4, 12, 16, 8, 6, 14, 10, 2, 1, 17, 25, 9, 13, 29, 21, 5, 7, 23, 31, 15, 11, 27, 19, 3, 4, 20, 28, 12, 16, 32, 24, 8, 6, 22, 30, 14, 10, 26, 18, 2, 1, 33, 49, 17, 25, 57, 41, 9, 13, 45, 61, 29, 21, 53, 37, 5, 7, 39, 55, 23
Offset: 1

Views

Author

Gary W. Adamson, Sep 23 2003

Keywords

Comments

Schroeder, p. 281 states "The ordering with which the iterates x_n fall into the 2^m different chaos bands [order as to magnitude] is also the same as the ordering of the iterates in a stable orbit of period length P = 2^m. For example, for both the period-4 orbit and the four chaos bands, the iterates, starting with the largest iterate x_1, are ordered as follows: x_1 > x_3 > x_4 > x_2."
From Andrey Zabolotskiy, Dec 06 2024: (Start)
For n>0, row n-1 is the permutation relating row n of the left half of Stern-Brocot tree with row n of Kepler's tree of fractions. Specifically, if K_n(k) [resp. SB_n(k)] is the k-th fraction in the n-th row of A294442 [resp. A057432], where 1/2 is in row 1 and k=1..2^(n-1), then SB_n(k) = K_n(T(n-1, k)).
The inverse permutation is row n of A131271.
Equals A362160+1. (End)

Examples

			1
1 2
1 3 4 2
1 5 7 3 4 8 6 2
1 9 13 5 7 15 11 3 4 12 16 8 6 14 10 2
		

References

  • Manfred R. Schroeder, "Fractals, Chaos, Power Laws", W.H. Freeman, 1991, p. 282.

Crossrefs

Programs

  • Haskell
    a088208 n k = a088208_tabf !! (n-1) !! (k-1)
    a088208_row n = a088208_tabf !! (n-1)
    a088208_tabf = iterate f [1] where
       f vs = (map (subtract 1) ws) ++ reverse ws where ws = map (* 2) vs
    -- Reinhard Zumkeller, Mar 14 2015
  • Mathematica
    nmax = 6;
    T[, 0] = 1; T[n, j_] /; j == 2^n = n;
    Do[Which[
      EvenQ[j], T[n+1, 2j] = T[n, j]; T[n+1, 2j+1] = T[n, j] + 2^n,
      OddQ[j], T[n+1, 2j+1] = T[n, j]; T[n+1, 2j] = T[n, j] + 2^n],
    {n, 0, nmax}, {j, 0, 2^n-1}];
    Table[T[n, j], {n, 0, nmax}, {j, 0, 2^n-1}] // Flatten (* Jean-François Alcover, Aug 03 2018 *)

Extensions

Edited by Ray Chandler and N. J. A. Sloane, Oct 08 2003

A086593 Bisection of A086592, denominators of the left-hand half of Kepler's tree of fractions.

Original entry on oeis.org

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, 19, 14, 19, 18, 21, 8, 13, 16, 17, 17, 22, 19, 23, 16, 23, 24, 27, 19, 26, 25, 29, 13, 20, 23, 25, 22, 29, 26, 31, 17, 25, 27, 30, 23, 31, 29, 34, 9, 15, 19, 20, 21, 27, 23, 28, 21, 30
Offset: 1

Views

Author

Antti Karttunen, Aug 28 2003

Keywords

Comments

Also denominator of alternate fractions in Kepler's tree as shown in A294442. - N. J. A. Sloane, Nov 20 2017

Crossrefs

Programs

  • Mathematica
    (* 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 *)
  • R
    maxlevel <- 15
    d <- c(1,2)
    for(m in 0:maxlevel)
     for(k in 1:2^m) {
       d[2^(m+1)    +k] <- d[k] + d[2^m+k]
       d[2^(m+1)+2^m+k] <- d[2^(m+1)+k]
    }
    a <- vector()
    for(m in 0:maxlevel) for(k in 0:(2^m-1)) a[2^m+k] <- d[2^(m+1)+k]
      a[1:63]
    # Yosu Yurramendi, May 16 2018

Formula

a(n) = A086592(2n-1) = A020650(4n-2).
a(n+1) = A071585(n) + A071766(n), n >= 0. - Yosu Yurramendi, Jun 30 2014
From Yosu Yurramendi, Jan 04 2016: (Start)
a(2^(m+1)+k+1) - a(2^m+k+1) = A071585(k), m >= 0, 0 <= k < 2^m.
a(2^(m+2)-k) = a(2^(m+1)-k) + a(2^m-k), m > 0, 0 <= k < 2^m-1.
(End)
a(2^n) = A000045(n+3). - Antti Karttunen, Jan 29 2016, based on above.
a(n) = A020651(4n-1), a(n+1) = A020651(4n+1), n > 0. - Yosu Yurramendi, May 08 2018
a(2^m+k) = A071585(2^(m+1)+k), m >= 0, 0 <= k < 2^m. - Yosu Yurramendi, May 16 2018

A131271 Triangular array T(n,k), n>=0, k=1..2^n, read by rows in bracketed pairs such that highest ranked element is bracketed with lowest ranked.

Original entry on oeis.org

1, 1, 2, 1, 4, 2, 3, 1, 8, 4, 5, 2, 7, 3, 6, 1, 16, 8, 9, 4, 13, 5, 12, 2, 15, 7, 10, 3, 14, 6, 11, 1, 32, 16, 17, 8, 25, 9, 24, 4, 29, 13, 20, 5, 28, 12, 21, 2, 31, 15, 18, 7, 26, 10, 23, 3, 30, 14, 19, 6, 27, 11, 22, 1, 64, 32, 33, 16, 49, 17, 48, 8, 57
Offset: 0

Views

Author

J. Demongeot (Jacques.Demongeot(AT)imag.fr), Jun 24 2007

Keywords

Comments

In a knockout competition with 2^n players, arranging the competition brackets (see Wikipedia) in T(n,k) order, where T(n,k) is the rank of the k-th player, ensures that highest ranked players cannot meet until the later stages of the competition. None of the top 2^p ranked players can meet earlier than the p-th from last round of the competition. At the same time the top ranked players in each match meet the lowest ranked player possible consistent with this rule. The sequence for the top ranked players meeting the highest ranked player possible is A049773. - Colin Hall, Feb 28 2012
Ranks in natural order of 2^n increasing real numbers appearing in limit cycles of interval iterations, or Median Spiral Order. [See the works by Demongeot & Waku]
From Andrey Zabolotskiy, Dec 06 2024 (Start):
For n>0, row n-1 is the permutation relating row n of Kepler's tree of fractions with row n of the left half of Stern-Brocot tree. Specifically, if K_n(k) [resp. SB_n(k)] is the k-th fraction in the n-th row of A294442 [resp. A057432], where 1/2 is in row 1 and k=1..2^(n-1), then K_n(k) = SB_n(T(n-1, k)).
The inverse permutation is row n of A088208.
When 1 is subtracted from each term, rows 3-5 become A240908, A240909, A240910. (End)

Examples

			Triangle begins:
1;
1,  2;
1,  4, 2, 3;
1,  8, 4, 5, 2,  7, 3,  6;
1, 16, 8, 9, 4, 13, 5, 12, 2, 15, 7, 10, 3, 14, 6, 11;
...
		

Crossrefs

Cf. A005578 (last elements in rows), A155944 (T(n,2^(n-1)) for n>0).

Programs

  • Maple
    T:= proc(n,k) option remember;
          `if`({n, k} = {1}, 1,
          `if`(irem(k, 2)=1, T(n-1, (k+1)/2), 2^(n-1)+1 -T(n-1, k/2)))
        end:
    seq(seq(T(n, k), k=1..2^(n-1)), n=1..7); # Alois P. Heinz, Feb 28 2012, with offset 1
  • Mathematica
    T[0, 1] = 1;
    T[n_, k_] := T[n, k] = If[Mod[k, 2] == 1, T[n, (k + 1)/2], 2^n + 1 - T[n, k/2]];
    Table[T[n, k], {n, 0, 6}, {k, 2^n}] // Flatten (* Jean-François Alcover, May 31 2018, after Alois P. Heinz *)

Formula

T(0,1) = 1, T(n,2k-1) = T(n-1,k), T(n,2k) = 2^n+1 - T(n-1,k).
T(n,1) = 1; for 1 < k <= 2^n, T(n,k) = 1 + (2^n)/m - T(n,k-m), where m = A006519(k-1). - Mathew Englander, Jun 20 2021

Extensions

Edited (with new name from Colin Hall) by Andrey Zabolotskiy, Dec 06 2024

A295511 The Schinzel-Sierpiński tree of fractions, read across levels.

Original entry on oeis.org

2, 2, 2, 3, 3, 2, 3, 7, 7, 5, 5, 7, 7, 3, 2, 5, 17, 13, 7, 11, 11, 5, 5, 11, 11, 7, 13, 17, 5, 2, 3, 11, 241, 193, 17, 29, 29, 13, 7, 17, 17, 11, 31, 43, 43, 13, 13, 43, 43, 31, 11, 17, 17, 7, 13, 29, 29, 17, 193, 241, 11, 3
Offset: 1

Views

Author

Peter Luschny, Nov 23 2017

Keywords

Comments

A conjecture of Schinzel and Sierpiński asserts that every positive rational number r can be represented as a quotient of shifted primes such that r = (p-1)/(q-1).
The function r -> [p, q] will be called the Schinzel-Sierpiński encoding of r if q is the smallest prime such that r = (p-1)/(q-1) for some prime p. In the case that no pair of such primes exists we set by convention p = q = 1.
The Schinzel-Sierpiński tree of fractions is the Euclid tree A295515 with root 1 and fractions represented by the Schinzel-Sierpiński encoding.

Examples

			The tree starts:
                                        2/2
                   2/3                                     3/2
         3/7                  7/5                 5/7                   7/3
   2/5       17/13     7/11        11/5     5/11       11/7     13/17        5/2
.
The numerators of the terms written as an array (the denominators are given by reversion of the arrays):
1: 2
2: 2, 3
3: 3, 7, 5, 7
4: 2, 17, 7, 11, 5, 11, 13, 5
5: 3, 241, 17, 29, 7, 17, 31, 43, 13, 43, 11, 17, 13, 29, 193, 11
		

References

  • E. Dijkstra, Selected Writings on Computing, Springer, 1982, p. 232.

Crossrefs

Programs

  • Sage
    def EuclidTree(n): # with root 1
        def DijkstraFusc(m):
            a, b, k = 1, 0, m
            while k > 0:
                if k % 2 == 1: b += a
                else: a += b
                k = k >> 1
            return b
        DF = [DijkstraFusc(k) for k in (2^(n-1)..2^n)]
        return [DF[j]/DF[j+1] for j in (0..2^(n-1)-1)]
    def SchinzelSierpinski(l):
        a, b = l.numerator(), l.denominator()
        p, q = 1, 2
        while q < 1000000000: # search limit
            r = a*(q - 1)
            if b.divides(r):
                p = r // b + 1
                if is_prime(p): return p/q
            q = next_prime(q)
        print("Search limit reached for ", l); return 0
    def SSETree(level):
        return [SchinzelSierpinski(l) for l in EuclidTree(level)]
    # With the imperfection that Sage reduces 2/2 automatically to 1.
    for level in (1..6): print(SSETree(level))

A295512 The Euclid tree with root 1 encoded by semiprimes, read across levels.

Original entry on oeis.org

4, -6, 6, -21, 35, -35, 21, -10, 221, -77, 55, -55, 77, -221, 10, -33, 46513, -493, 377, -119, 187, -1333, 559, -559, 1333, -187, 119, -377, 493, -46513, 33, -14, 143, -209, 629, -14527, 2881, -1189, 533, -161, 391, -15229, 2449, -2263, 3139, -1073, 95, -95
Offset: 1

Views

Author

Peter Luschny, Nov 23 2017

Keywords

Comments

The Euclid tree with root 1 is A295515 (sometimes called Calkin-Wilf tree).
For a positive rational r we use the Schinzel-Sierpiński encoding r -> [p, q] as described in A295511 and encode r as sgn*p*q where sgn is -1 if r < 1, else +1.
Apart from a(1) all terms are squarefree.

Examples

			The tree starts:
                     4
            -6                 6
       -21       35      -35       21
    -10  221  -77  55  -55  77  -221  10
		

References

  • E. Dijkstra, Selected Writings on Computing, Springer, 1982, p. 232.

Crossrefs

Programs

  • Maple
    EuclidTree := proc(n) local k, DijkstraFusc;
    DijkstraFusc := proc(m) option remember; local a, b, n; a := 1; b := 0; n := m;
    while n > 0 do if type(n, odd) then b := a+b else a := a+b fi; n := iquo(n,2) od; b end:
    seq(DijkstraFusc(k)/DijkstraFusc(k+1), k=2^(n-1)..2^n-1) end:
    SchinzelSierpinski := proc(l) local a, b, r, p, q, sgn;
    a := numer(l); b := denom(l); q := 2; sgn := `if`(a < b, -1, 1);
    while q < 1000000000 do r := a*(q - 1); if r mod b = 0 then p := r/b + 1;
    if isprime(p) then return(sgn*p*q) fi fi; q := nextprime(q); od;
    print("Search limit reached!", a, b) end:
    Tree := level -> seq(SchinzelSierpinski(l), l=EuclidTree(level)): seq(Tree(n), n=1..6);

A295515 The Euclid tree, read across levels.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 2, 1, 1, 3, 3, 2, 2, 3, 3, 1, 1, 4, 4, 3, 3, 5, 5, 2, 2, 5, 5, 3, 3, 4, 4, 1, 1, 5, 5, 4, 4, 7, 7, 3, 3, 8, 8, 5, 5, 7, 7, 2, 2, 7, 7, 5, 5, 8, 8, 3, 3, 7, 7, 4, 4, 5, 5, 1, 1, 6, 6, 5, 5, 9, 9, 4, 4, 11, 11, 7, 7, 10, 10, 3, 3, 11, 11, 8, 8
Offset: 1

Views

Author

Peter Luschny, Nov 25 2017

Keywords

Comments

Set N(x) = 1 + floor(x) - frac(x) and let '"' denote the ditto operator, referring to the previously computed expression. Assume the first expression is '0'. Then [0, repeat(N("))] will generate the natural numbers 0, 1, 2, 3, ... and [0, repeat(1/N("))] will generate the rational numbers 0/1, 1/1, 1/2, 2/1, 1/3, 3/2, ... Every reduced nonnegative rational number r appears exactly once in this list as a relatively prime pair [n, d] = r = n/d. We list numerator and denominator one after the other in the sequence.
The apt name 'Euclid tree' is taken from the exposition of Malter, Schleicher and Don Zagier. It is sometimes called the Calkin-Wilf tree. The enumeration is based on Stern's diatomic series (which is a subsequence) and computed by a modification of Dijkstra's 'fusc' function.
The tree listed has root 0, the variant with root 1 is more widely used. Seen as sequences the difference between the two trees is trivial: it is enough to leave out the first two terms; but as trees they are markedly different (see the example section).

Examples

			The tree with root 0 starts:
                                      [0/1]
                  [1/1,                                    1/2]
        [2/1,                1/3,                3/2,                2/3]
   [3/1,      1/4,      4/3,      3/5,      5/2,      2/5,      5/3,      3/4]
[4/1, 1/5, 5/4, 4/7, 7/3, 3/8, 8/5, 5/7, 7/2, 2/7, 7/5, 5/8, 8/3, 3/7, 7/4, 4/5]
.
The tree with root 1 starts:
                                      [1/1]
                  [1/2,                                    2/1]
        [1/3,                3/2,                2/3,                3/1]
   [1/4,      4/3,      3/5,      5/2,      2/5,      5/3,      3/4,      4/1]
[1/5, 5/4, 4/7, 7/3, 3/8, 8/5, 5/7, 7/2, 2/7, 7/5, 5/8, 8/3, 3/7, 7/4, 4/5, 5/1]
		

References

  • M. Aigner and G. M. Ziegler, Proofs from The Book, Springer-Verlag, Berlin, 3rd ed., 2004.

Crossrefs

Cf. A002487, A174981, A294446 (Stern-Brocot tree), A294442 (Kepler's tree), A295511 (Schinzel-Sierpiński tree), A295512 (encoded by semiprimes).

Programs

  • Maple
    # First implementation: use it only if you are not afraid of infinite loops.
    a := x -> 1/(1+floor(x)-frac(x)): 0; do a(%) od;
    # Second implementation:
    lusc := proc(m) local a, b, n; a := 0; b := 1; n := m; while n > 0 do
    if n mod 2 = 1 then b := a + b else a := a + b fi; n := iquo(n, 2) od; a end:
    R := n -> 3*2^(n-1)-1 .. 2^n: # The range of level n.
    EuclidTree_rat := n -> [seq(lusc(k+1)/lusc(k), k=R(n), -1)]:
    EuclidTree_num := n -> [seq(lusc(k+1), k=R(n), -1)]:
    EuclidTree_den := n -> [seq(lusc(k), k=R(n), -1)]:
    EuclidTree_pair := n -> ListTools:-Flatten([seq([lusc(k+1), lusc(k)], k=R(n), -1)]):
    seq(print(EuclidTree_pair(n)), n=1..5);
  • Sage
    def A295515(n):
        if n == 1: return 0
        M = [0, 1]
        for b in (n//2 - 1).bits():
            M[b] = M[0] + M[1]
        return M[1]
    print([A295515(n) for n in (1..85)])

Formula

Some characteristics in comparison to the tree with root 1, seen as a table with T(n,k) for n >= 1 and 1 <= k <= 2^(n-1). Here Tr(n,k), Tp(n,k), Tq(n,k) denotes the fraction r, the numerator of r and the denominator of r in row n and column k respectively.
With root 0: With root 1:
Root Tr(1,1) 0/1 1/1
Tp(n,1) 0,1,2,3,... 1,1,1,1,...
Tp(n,2^(n-1)) 0,1,2,3,... 1,2,3,4,...
Tq(n,1) 1,1,1,1,... 1,2,3,4,...
Tq(n,2^(n-1)) 1,2,3,4,... 1,1,1,1,...
Sum_k Tp(n,k) 0,2,8,26,... A024023 1,3,9,27,... A000244
Sum_k Tq(n,k) 1,3,9,27,... A000244 1,3,9,27,... A000244
Sum_k 2Tr(n,k) 0,3,9,21,... A068156 2,5,11,23,... A083329
Sum_k Tp(n,k)Tq(n,k) 0,3,17,81,... A052913-1 1,4,18,82,... A052913
----
a(n) = A002487(floor(n/2)). - Georg Fischer, Nov 29 2022
Showing 1-10 of 13 results. Next