%I #47 Feb 21 2022 21:40:10
%S 1,1,1,1,2,1,1,1,4,6,6,5,2,1,1,1,8,28,60,94,116,114,94,69,44,26,14,5,
%T 2,1,1,1,16,120,568,1932,5096,10948,19788,30782,41944,50788,55308,
%U 54746,49700,41658,32398,23461,15864,10068,6036,3434,1860,958,470,221,100,42,14,5,2,1,1
%N Triangle by rows, related to the numbers of binary trees of height less than n, derived from the Mandelbrot set.
%C As shown on p. 74 [Diaconis & Graham], n-th row polynomials are cyclic with period n, given real roots, if the polynomials are divided through by n. For example, taking x^3 + 2x^2 + x + 1 = 0, the real root = -1.75487766... = c. Then using x^2 + c, we obtain the period three trajectory: -1.75487... -> 1.32471...-> 0.
%C The shuffling connection [p.75], resulting in a permutation that is the Gilbreath shuffle: "To make the connection with shuffling cards, write down a periodic sequence starting at zero. Write a one above the smallest point, a two above the next smallest point and so on. For example, if c = -1.75486...(a period three point), we have:
%C 2.............1.............3......
%C 0........-1.75487........1.32471... For a fixed value of c, the numbers written on top code up a permutation that is a Gilbreath shuffle".
%C Row sums = A003095: (1, 2, 5, 26, 677,...) relating to the number of binary trees of height less than n.
%C Let f(z) = z^2 + c, then row k lists the expansion of the n-fold composition f(f(...f(0)...) in falling powers of c (with the coefficients for c^0 omitted). The n initial terms of the reversed n-th row are the Catalan numbers (cf. A137560). - _Joerg Arndt_, Jun 04 2016
%D Persi Diaconis & R. L. Graham, "Magical Mathematics: The Mathematical Ideas That Animate Great Magic Tricks", Princeton University Press, 2012; pp. 73-83.
%H Alois P. Heinz, <a href="/A202019/b202019.txt">Rows n = 1..13, flattened</a>
%H Juan Carlos Nuño and Francisco J. Muñoz, <a href="https://arxiv.org/abs/2112.15563">Entropy-Variance curves of binary sequences generated by random substitutions of constant length</a>, arXiv:2112.15563 [math.PR], 2021.
%F Coefficients of x by rows such that given n-th row p(x), the next row is (p(x))^2 + x); starting x -> (x^2 + x) -> (x^4 + 2x^3 + x^2 + x)...
%F T(n,k) = A309049(2^(n-1)-1,k-1). - _Alois P. Heinz_, Jul 11 2019
%e Row 4 = (1, 4, 6, 6, 5, 2, 1, 1) since (x^4 + 2x^3 + x^2 + x)^2 + x = x^8 + 4x^7 + 6x^6 + 6x^5 + 5x^4 + 2x^3 + x^2 + x.
%e Triangle begins:
%e 1;
%e 1, 1;
%e 1, 2, 1, 1;
%e 1, 4, 6, 6, 5, 2, 1, 1;
%e 1, 8, 28, 60, 94, 116, 114, 94, 69, 44, 26, 14, 5, 2, 1, 1;
%e ...
%p b:= proc(n) option remember; `if`(n=0, 1, (g-> (f-> expand(
%p x^n+b(f)*b(n-1-f)))(min(g-1, n-g/2)))(2^ilog2(n)))
%p end:
%p T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(2^(n-1)-1)):
%p seq(T(n), n=1..7); # _Alois P. Heinz_, Jul 11 2019
%t b[n_] := b[n] = If[n == 0, 1, Function[g, Function[f, Expand[x^n + b[f]*
%t b[n-1-f]]][Min[g-1, n-g/2]]][2^(Length[IntegerDigits[n, 2]]-1)]];
%t T[n_] := CoefficientList[b[2^(n-1)-1], x];
%t Array[T, 7] // Flatten (* _Jean-François Alcover_, Feb 19 2021, after _Alois P. Heinz_ *)
%Y Row sums are A003095.
%Y Cf. A137560 (reversed rows).
%Y Cf. A309049.
%K nonn,tabf
%O 1,5
%A _Gary W. Adamson_, Dec 08 2011