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.

Previous Showing 11-20 of 21 results. Next

A359025 Number of inequivalent tilings of a 9 X n rectangle using integer-sided square tiles.

Original entry on oeis.org

1, 1, 30, 163, 2403, 32097, 459957, 6542578, 93604244
Offset: 0

Views

Author

John Mason, Dec 12 2022

Keywords

Crossrefs

Column k = 9 of A227690.
Sequences for fixed and free (inequivalent) tilings of m X n rectangles, for 2 <= m <= 10:

A359026 Number of inequivalent tilings of a 10 X n rectangle using integer-sided square tiles.

Original entry on oeis.org

1, 1, 51, 347, 7048, 130125, 2551794, 49828415
Offset: 0

Views

Author

John Mason, Dec 12 2022

Keywords

Crossrefs

Column k = 10 of A227690.
Sequences for fixed and free (inequivalent) tilings of m X n rectangles, for 2 <= m <= 10:

A068928 Number of incongruent ways to tile a 3 X 2n room with 1x2 Tatami mats. At most 3 Tatami mats may meet at a point.

Original entry on oeis.org

2, 2, 2, 4, 5, 9, 12, 21, 30, 51, 76, 127, 195, 322, 504, 826, 1309, 2135, 3410, 5545, 8900, 14445, 23256, 37701, 60813, 98514, 159094, 257608, 416325, 673933, 1089648, 1763581, 2852242, 4615823, 7466468, 12082291, 19546175, 31628466
Offset: 1

Views

Author

Dean Hickerson, Mar 11 2002

Keywords

Crossrefs

Cf. A068922 for total number of tilings, A068926 for more info.
Essentially the same as A001224.

Formula

For n >= 8, a(n) = a(n-1) + 2*a(n-2) - a(n-3) - a(n-5) - a(n-6).
O.g.f.: x(2-4x^2-x^4+x^6)/((1-x-x^2)(1-x^2-x^4)). a(n) = (A000045(n+1)+A053602(n+1))/2, n>1. [From R. J. Mathar, Aug 30 2008]

A032128 Number of dyslexic planted planar trees with n nodes.

Original entry on oeis.org

1, 1, 2, 4, 10, 25, 69, 193, 565, 1680, 5113, 15757, 49223, 155228, 493937, 1583002, 5106386, 16563542, 53995678, 176797966, 581196445, 1917446630, 6346554919, 21068877925, 70133571797, 234043258802, 782831380626, 2624022529690, 8813080348897, 29654400681966, 99953565213645, 337447946046906, 1140961171059563
Offset: 1

Views

Author

Keywords

Comments

From Petros Hadjicostas, Jan 14 2018: (Start)
For this sequence, if (b(n): n>=1) = BIK((a(n): n>=1)), then b(n) = a(n+1) for n>=1.
Let A(x) = Sum_{n>=1} a(n)*x^n be the g.f. for this sequence. For an explanation on how to derive the formula BIK(A(x)) = (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1 - A(x^2))) from Bower's formulae in the link below about transforms, see the comments for sequence A001224. (For that sequence, the roles of sequences (a(n): n>=1) and (b(n): n>=1) are reversed.)
According to Bower's theory in the link below, we have boxes of different sizes and colors. The size of a box is determined by the number of balls it can hold. Two boxes of the same size and color are considered identical (indistinct and unlabeled). We place the boxes on a line that can be read in either direction; i.e., we have a reversible line.
Here, a(n) = number of colors a box holding n balls can be, while b(n) = number of ways of placing boxes in a line that can be read in either direction when the total number of balls is n.
According to C. G. Bower in the weblink below, a "[d]yslexic planar tree is a planar tree where each sub-rooted tree extending from a node can be read from left to right or right to left." A dyslexic planar tree "can be thought of as viewed by an observer who does not know left from right or as sub-rooted trees that can be turned around independent of the rest of the tree."
Given the above definition, a "ball" is a "node", a "box" is a "sub-rooted tree" (without the single root of the whole planar tree), and a "set of colors" for a "box" with n "balls" (= "nodes") is "the set of all dyslexic trees with n nodes". Hence, a(n) = number of all dyslexic planar trees with n nodes. On the other hand, b(n) = number of "reversible" arrangements of having sub-rooted dyslexic planar trees (= boxes on a reversible line) with a total of n-1 nodes in all subtrees (the n-th node is the single root of the whole tree). This means that b(n-1) = a(n) for n >= 2.
(End)

Examples

			From _Petros Hadjicostas_, Jan 13 2018: (Start)
For this sequence, if (b(n): n>=1) = BIK((a(n): n>=1)), then b(n) = a(n+1) for n>=1.
Since a(1) = 1, a(2) = 1, a(3) = 2, a(4) = 4, etc., a box with 1 ball can be of 1 color only, a box with 2 balls can be of 1 color only, a box with 3 balls can be of 2 colors only, a box with 4 balls can be of 4 colors, and so on.
When we have n=2 balls, we have b(2) = a(3) = 2 because we either have two identical boxes on a line each holding 1 ball or a single box holding 2 balls (and all these boxes can be of 1 color only).
When we have n=3 balls, we have b(3) = a(4) = 4. Here, we consider three cases. In the first case, we have one box holding 3 balls and we have 2 possibilities. In the second case, we have a box with 2 balls and a box with 1 ball, and we have 1 possibility here because the line is reversible (i.e., 21 is considered the same as 12). In the third case, we have three identical boxes each holding 1 ball. Thus, b(3) = 2 + 1 + 1 = 4 = a(4).
When we have n=4 balls, we have b(4) = a(5) = 10. Here we consider 5 cases: a single box with 4 balls (a(4) = 4 possibilities); a box with 3 balls and a box with 1 ball (a(3) x a(1) = 2 x 1 = 2 possibilities); two identical boxes each with 2 balls (1 possibility because a(2) = 1); a box with 2 balls and two identical boxes each with 1 ball (2 possibilities because we have the cases 121 and 112); and 4 identical boxes each with 1 ball (1 possibility). Hence, b(4) = 4 + 2 + 1 + 2 + 1 = 10 = a(5).
Let us now switch the discussion to the counting of dyslexic planar trees. We explain why a(5) = 10. We have five nodes, but one of them is used for the single root of the whole tree. The other 4 nodes are used to create sub-rooted dyslexic planar trees. There are b(4) = 10 ways of doing that. As above, we consider 5 cases: a single subtree with 4 nodes (with a(4) = 4 possibilities); a subtree with 3 nodes and subtree with 1 node both connected to the single root (with a(3) x a(1) = 2 x 1 = 2 possibilities); two identical subtrees each with 2 nodes and connected to the single root (with a(2) = 1 possibilities); a subtree with 2 nodes and two identical subtrees each with 1 nodes, all connected to the single root (with 2 possibilities because we have the cases 121 and 112); and 4 identical sub-rooted trees each with 1 node (1 possibility). Hence, b(4) = 4 + 2 + 1 + 2 + 1 = 10 = a(5).
(End)
		

Crossrefs

When a(1) = 2, we get sequence A032130.

Programs

  • Mathematica
    m = 34; a[1] = 1; A[_] = 0;
    Do[A[x_] = x(a[1]+(1/2)(A[x]/(1-A[x])+(A[x]+A[x^2])/(1-A[x^2]))) + O[x]^m // Normal, {m}];
    CoefficientList[A[x], x] // Rest (* Jean-François Alcover, Sep 17 2019 *)
  • PARI
    BIK(p)={(1/(1-p) + (1+p)/subst(1-p, x, x^2))/2}
    seq(n)={my(p=O(1));for(i=1, n, p=BIK(x*p)); Vec(p)} \\ Andrew Howroyd, Aug 30 2018

Formula

Shifts left under "BIK" (reversible, indistinct, unlabeled) transform.
From Petros Hadjicostas, Jan 13 2018: (Start)
G.f.: If A(x) = Sum_{n>=1} a(n)*x^n, then -a(1) + A(x)/x = BIK(A(x)) = (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1-A(x^2))). Here a(1) = 1.
In general, if we let a(1) = c, we get:
a(2) = c,
a(3) = (1/2)*(c + 3)*c,
a(4) = (1/2)*(c + 3)*(c + 1)*c,
a(5) = (1/2)*(c^3 + 5*c^2 + 10*c + 4)*c,
a(6) = (1/4)*(2*c^4 + 15*c^3 + 38*c^2 + 37*c + 8)*c,
a(7) = (1/8)*(4*c^4 + 38*c^3 + 103*c^2 + 109*c + 22)*(c + 1)*c,
a(8) = (1/8)*(4*c^6 + 56*c^5 + 251*c^4 + 511*c^3 + 499*c^2 + 201*c + 22)*c,
and so on. No pattern is apparent in these formulae.
(End)

Extensions

a(28)-a(33) from Petros Hadjicostas, Jan 13 2018
Name edited by Petros Hadjicostas, Jan 14 2018

A060312 Number of distinct ways to tile a 2 X n rectangle with dominoes (solutions are identified if they are rotations or reflections of each other).

Original entry on oeis.org

1, 1, 2, 4, 5, 9, 12, 21, 30, 51, 76, 127, 195, 322, 504, 826, 1309, 2135, 3410, 5545, 8900, 14445, 23256, 37701, 60813, 98514, 159094, 257608, 416325, 673933, 1089648, 1763581, 2852242, 4615823, 7466468, 12082291, 19546175, 31628466
Offset: 1

Views

Author

Thomas Ward, Mar 27 2001

Keywords

Comments

Same as A001224 except that there a(2)=2 not 1. - N. J. A. Sloane, Mar 30 2015

Examples

			a(3)=2 because of the configurations |= and |||.
		

Crossrefs

Essentially the same as A001224, which is the main entry for this sequence. Other versions of the sequence can be found in A068928 and A102526.

Programs

  • Magma
    [n eq 1 select 1 else (1/2)*(Fibonacci(n+2)+Fibonacci(Floor((n-(-1)^n)/2)+2)): n in [0..40]]; // Vincenzo Librandi, Nov 22 2014
  • Maple
    # Maple code for A060312 and A001224 from N. J. A. Sloane, Mar 30 2015
    with(combinat); F:=fibonacci;
    f:=proc(n) option remember;
    if n=2 then 1 # change this to 2 to get A001224
    elif (n mod 2) = 0 then (F(n+1)+F(n/2+2))/2;
    else (F(n+1)+F((n+1)/2))/2; fi; end;
    [seq(f(n),n=1..50)];
  • Mathematica
    CoefficientList[Series[-(x^7 + x^6 + x^5 + 2 x^4 - x^3 + x^2 - 1) / ((x^2 + x - 1) (x^4 + x^2 - 1)), {x, 0, 40}], x] (* Vincenzo Librandi, Nov 22 2014 *)

Formula

If F(n) is the n-th Fibonacci number, then a(2n) = (F(2n+1) + F(n+2))/2 for n > 1 and a(2n-1) = (F(2n) + F(n))/2. [Corrected by Manfred Boergens, Aug 25 2025]
a(n) = (F(n+1)+F(floor((n+3+(-1)^n)/2)))/2 for n!=2. - Manfred Boergens, Aug 25 2025
G.f.: -x*(x^7 + x^6 + x^5 + 2*x^4 - x^3 + x^2 - 1) / ((x^2 + x - 1)*(x^4 + x^2 - 1)). - Colin Barker, Dec 15 2012

Extensions

Edited by N. J. A. Sloane, Mar 30 2015

A102526 Antidiagonal sums of Losanitsch's triangle (A034851).

Original entry on oeis.org

1, 1, 2, 2, 4, 5, 9, 12, 21, 30, 51, 76, 127, 195, 322, 504, 826, 1309, 2135, 3410, 5545, 8900, 14445, 23256, 37701, 60813, 98514, 159094, 257608, 416325, 673933, 1089648, 1763581, 2852242, 4615823, 7466468, 12082291, 19546175, 31628466
Offset: 0

Views

Author

Gerald McGarvey, Feb 24 2005

Keywords

Comments

This is an interleaving of A005207 and A051450. Thus a(2*m) = A005207(m) = (F(2*m-1) + F(m+1)) / 2, a(2*m - 1) = A051450(m) = (F(2*m) + F(m)) / 2 where F() are Fibonacci numbers (A000045). - Max Alekseyev, Jun 28 2006
The Kn11(n) and Kn21(n) sums, see A180662 for their definitions, of Losanitsch's triangle A034851 equal a(n), while the Kn12(n) and Kn22(n) sums equal (a(n+2)-A000012(n)) and the Kn13(n) and Kn23(n) sums equal (a(n+4)-A008619(n+4)). - Johannes W. Meijer, Jul 14 2011
a(n) is the number of homeomorphically irreducible caterpillars with n + 3 edges. - Christian Barrientos, Sep 12 2020

References

  • Jablan S. and Sazdanovic R., LinKnot: Knot Theory by Computer, World Scientific Press, 2007.

Crossrefs

Cf. A034851.
Essentially the same as A001224, A060312 and A068928.

Programs

  • Maple
    with(combinat): A102526 :=proc(n): if type(n, even) then (fibonacci(n+1)+fibonacci(n/2+2))/2 else (fibonacci(n+1)+fibonacci((n+1)/2))/2 fi: end: seq(A102526(n), n=0..38); # Johannes W. Meijer, Jul 14 2011
  • Mathematica
    LinearRecurrence[{1, 2, -1, 0, -1, -1}, {1, 1, 2, 2, 4, 5}, 40] (* Jean-François Alcover, Nov 17 2017 *)
  • PARI
    Vec((1+x)*(1-x-x^3)/(x^2+x-1)/(x^4+x^2-1)+O(x^99)) \\ Charles R Greathouse IV, Nov 17 2017
    
  • PARI
    a(n)=([0,1,0,0,0,0; 0,0,1,0,0,0; 0,0,0,1,0,0; 0,0,0,0,1,0; 0,0,0,0,0,1; -1,-1,0,-1,2,1]^n*[1;1;2;2;4;5])[1,1] \\ Charles R Greathouse IV, Nov 17 2017

Formula

G.f.: -(1+x)*(x^3+x-1) / ( (x^2+x-1)*(x^4+x^2-1) ). - R. J. Mathar, Nov 09 2013
a(n) = a(n-1) + 2*a(n-2) - a(n-3) - a(n-5) - a(n-6). - Wesley Ivan Hurt, Sep 17 2020

A032132 Number of series-reduced dyslexic planted planar trees with n leaves.

Original entry on oeis.org

1, 1, 2, 6, 17, 57, 191, 684, 2482, 9275, 35127, 135156, 525545, 2064329, 8173895, 32600082, 130823306, 527888023, 2140454687, 8716907165, 35638352814, 146221542191, 601870210193, 2484682879348, 10285116277096, 42679973961811, 177514171393035, 739881841810694, 3089914920914855, 12927860306782626
Offset: 1

Views

Author

Keywords

Comments

Apparently, beginning with a(3), number of non-equivalent canonical forms of separation coordinates on the hyperspheres. Cf. Schöbel and Veselov for this and other interpretations. - Tom Copeland, Nov 21 2017
From Petros Hadjicostas, Jan 17 2018: (Start)
Let A(x) = Sum_{n>=1} a(n)*x^n. For a derivation of the formula BIK(A(x)) = (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1-A(x^2))), see the comments for sequence A001224 and the weblink below containing Bower's theory of transforms.
We clarify the comment by T. Copeland above. Consider the material in Section 3 of Devadoss and Read (2001). According to their terminology, let b(m,n) be "the number of A-clusters having m cells and n outside edges not counting the root edge." Let B(x,y) = Sum_{m>=0, n>=0} b(m,n)*x^m*y^n. (See p. 78 in their paper, where they use the notations a_{m,n} and A(x,y) rather than b(m,n) and B(x,y), respectively, that we use here.)
On p. 79 (Eq. (3.1)) of their paper, they prove that B(x,y) = y + (x/2)*(B(x,y)^2/(1-B(x,y)) + (1 + B(x,y))*B(x^2, y^2)/(1-B(x^2,y^2))). Unfortunately, the factor x in the previous formula is left out (i.e., it is a typo), and the same typo is reproduced in Schöbel and Veselov (2014, 2015).
Using Table 2 (p. 92) from Devadoss and Read (2001) (and the material on p. 79), we get that B(x,y) = y+ x*y^2 + (x^2 + x)*y^3 + (2*x^3 + 3*x^2 + x)*y^4 + (3*x^4 + 8*x^3 + 5*x^2 + x)*y^5 + ...
We claim that a(n) = Sum_{m>=0} b(m,n) and A(y) = Sum_{n>=1} a(n)*y^n = B(x=1, y). To prove these claims, note that, for x=1, the above series becomes B(x=1,y) = y + y^2 + 2*y^3 + 6*y^4 + 17*y^5 + ..., while the functional equation above becomes B(1, y) = y + (1/2)*(B(1,y)^2/(1-B(1,y)) + (1 + B(1,y))*B(1,y^2)/(1-B(1,y^2))), which is equivalent to 2*B(1,y) = y + (1/2)*(B(1,y)/(1-B(1,y)) + (B(1,y) + B(1,y^2))/(1-B(1,y^2))). The latter formula is the one given in the formula section below (derived from Bower's theory) with x replaced with y and A(x) replaced with B(1,y). This proves that B(x=1, y) = A(y), from which we can easily get that a(n) = Sum_{m>=0} b(m,n).
Note that b(m=0, n) = 0 for n <> 1, but b(m=0, n=1) = 1; b(m,n) = 0 when m >= n >= 1; and b(m=1, n) = 1 for n>=2. Also, b(m,m+1) = A001190(m+1) for m>=1, which are the Wedderburn-Etherington numbers, and apparently b(m=2, n) = A024206(n-1) for n>=2 (conjecture).
In Section 6 of their paper, Schöbel and Veselov (2014, 2015) prove that b(m,n) is the "number of non-equivalent faces of [the Stasheff polytope] K_n of codimension m-1." Apparently then, for n>=2 and k>=0, b(n-k,n+1) is the "number of canonical forms for separation coordinates of [hypersphere] S^n" with k "independent continuous parameters". For k=0 and n>=2, b(n,n+1) = A001190(n+1) = "number of canonical forms for separation coordinates" of hypersphere S^n with 0 continuous parameters.
It turns out that for k, the number of continuous parameters of S^n, we have 0 <= k <= n-1 (see pp. 1269-1270 in Shobel and Veselov (2015)). Hence, for n>=2, Sum_{k=0..n-1} b(n-k, n+1) = Sum_{m=1..n} b(m, n+1) = Sum_{m=0..n} b(m, n+1) = a(n+1) (see above). As a result, for n>=2, a(n+1) is the "total number of [non-equivalent] canonical forms for separation coordinates on [hypersphere] S^n", which is the comment made by T. Copeland above.
(End)
For an explanation on the meaning of clusters of types A, B, and C see Section 3 (pp. 78-81) in Devadoos and Read (2001). See also the comments for sequence A232206. - Petros Hadjicostas, Mar 02 2018

Crossrefs

Programs

  • Mathematica
    BIK[p_] := (1/(1-p) + (1+p)/(1-p /. x -> x^2))/2;
    seq[n_] := Module[{p=x}, For[i=2, i <= n, i++, p += x^i Coefficient[BIK[p] + x O[x]^i // Normal, x, i]]; CoefficientList[p/x, x]];
    seq[30] (* Jean-François Alcover, Nov 22 2018, after Andrew Howroyd *)
  • PARI
    BIK(p)={(1/(1-p) + (1+p)/subst(1-p, x, x^2))/2}
    seq(n)={my(p=x); for(i=2, n, p+=x^i*polcoeff(BIK(p) + O(x*x^i), i)); Vecrev(p/x)} \\ Andrew Howroyd, Aug 30 2018

Formula

Doubles (index 2+) under "BIK" (reversible, indistinct, unlabeled) transform.
G.f.: If A(x) = Sum_{n>=1} a(n)*x^n, then 2*A(x) = x + BIK(A(x)) = x + (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1-A(x^2))). - Petros Hadjicostas, Jan 17 2018

Extensions

a(25)-a(30) from Petros Hadjicostas, Jan 17 2018

A089935 a(n) = order of recurrence generating row (or column) n of A089934.

Original entry on oeis.org

2, 2, 4, 5, 9, 11, 21, 30, 51, 76, 127, 195, 322, 504, 826
Offset: 1

Views

Author

Marc LeBrun, Nov 15 2003

Keywords

Comments

The only known value where a(n) is strictly less than the theoretical upper bound A001224(n+1) is a(6) = 11. - Andrew Howroyd, Dec 24 2019

Examples

			a(2)=2 because the recurrence relation for the second row/column is a(n) = 2 a(n-1) + a(n-2).
		

Crossrefs

Row/columns 1 through 7 of A089934 are A000045, A001333, A051736, A051737, A089936, A089937, A089938.
Cf. A001224.

Formula

a(n) <= A001224(n+1). - Andrew Howroyd, Dec 24 2019

Extensions

a(8)-a(10) from Andrew Howroyd, Dec 24 2019
a(11)-a(15) from Max Alekseyev, Dec 12 2024

A032130 Shifts left under the "BIK" (reversible, indistinct, unlabeled) transform with a(1) = 2.

Original entry on oeis.org

2, 2, 5, 15, 52, 193, 765, 3143, 13323, 57670, 254040, 1134249, 5122124, 23349966, 107310784, 496633774, 2312539465, 10826481544, 50929829953, 240616214596, 1141195080020, 5431477088428, 25933525825389, 124185539096075, 596268057962349, 2869992942831031, 13845453533124431, 66934180769445444, 324218809545624984
Offset: 1

Views

Author

Keywords

Comments

From Petros Hadjicostas, Jan 14 2018: (Start)
For this sequence, if (b(n): n>=1) = BIK((a(n): n>=1)), then b(n) = a(n+1) for n>=1.
Let A(x) = Sum_{n>=1} a(n)*x^n be the g.f. for this sequence. For an explanation on how to derive the formula BIK(A(x)) = (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1 - A(x^2))) from Bower's formulae in the link below about transforms, see the comments for sequence A001224. (For that sequence, the roles of sequences (a(n): n>=1) and (b(n): n>=1) are reversed.)
(End)

Examples

			From _Petros Hadjicostas_, Jan 14 2018: (Start)
According to Bower's theory in the link above, we have boxes of different sizes and colors. The size of a box is determined by the number of balls it can hold. Two boxes of the same size and color are considered identical (indistinct and unlabeled). We place the boxes on a line that can be read in either direction; i.e., we have a reversible line.
Here, a(n) = number of colors a box holding n balls can be, while b(n) = number of ways of placing boxes in a line that can be read in either direction when the total number of balls is n.
Since a(1) = 2, a(2) = 2, a(3) = 5, a(4) = 15, etc.,  a box with 1 ball can be of 2 colors only, a box with 2 balls can be of 2 colors only, a box with 3 balls can be of 5 colors only, a box with 4 balls can be of 15 colors, and so on.
When we have n=3 balls, we have b(3) = a(4) = 15. Here, we consider three cases. In the first case, we have one box holding 3 balls and we have 5 possibilities. In the second case, we have a box with 2 balls and a box with 1 ball, and we have 2 x 2 = 4 possibilities here because the line is reversible (i.e., 21 is considered the same as 12). In the third case, we have three identical boxes each holding 1 ball and we have 6 possibilities (if the colors are a and b, we have the possibilities aaa, aab, aba, bba, bab, and bbb). Thus, b(3) = 5 + 4 + 6 = 15 = a(4).
When we have n=4 balls, we have b(4) = a(5) = 52. Here we consider 5 cases: a single box with 4 balls (a(4) = 15 possibilities); a box with 3 balls and a box with 1 ball (a(3) x a(1) = 5 x 2 = 10 possibilities); two identical boxes each with 2 balls (3 possibilities, aa, ab, and bb); a box with 2 balls and two identical boxes each with 1 ball (14 possibilities, see below); and 4 identical boxes each with 1 ball (10 possibilities, aaaa, aaab, aaba, aabb, abba, baab, abab, abbb, babb, bbbb). Hence, b(4) = 15 + 10 + 3 + 14 + 10 = 52 = a(5).
We explain the fourth case above in more detail. Here, we have a box with 2 balls and two identical boxes each with 1 ball. Let a and b be the two colors for the 1-ball boxes and A and B be the colors for the 2-ball boxes. Then we have the following 14 cases: Aaa, Aab, Abb, Baa, Bab, Bbb, aAa, aAb, bAb, aBa, aBb, bBb, abA, and abB. Note that Aab = baA <> abA = Aba and abB = Bba <> Bab = baB.
(End)
		

Crossrefs

When a(1) = 1, we get sequence A032128.

Programs

  • Mathematica
    m = 30; a[1] = 2; A[_] = 0;
    Do[A[x_] = x(a[1]+(1/2)(A[x]/(1-A[x])+(A[x]+A[x^2])/(1-A[x^2]))) + O[x]^m // Normal, {m}];
    CoefficientList[A[x], x] // Rest (* Jean-François Alcover, Sep 17 2019 *)
  • PARI
    BIK(p)={(1/(1-p) + (1+p)/subst(1-p, x, x^2))/2}
    seq(n)={my(p=O(1));for(i=1, n, p=1+BIK(x*p)); Vec(p)} \\ Andrew Howroyd, Aug 30 2018

Formula

From Petros Hadjicostas, Jan 14 2018: (Start)
The sequence (a(n): n>=1) shifts left under the "BIK" (reversible, indistinct, unlabeled) transform with a(1) = 2.
G.f.: If A(x) = Sum_{n>=1} a(n)*x^n, then -a(1) + A(x)/x = BIK(A(x)) = (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1-A(x^2))). Here a(1) = 2.
In general, if we let a(1) = c, we get:
a(2) = c,
a(3) = (1/2)*(c + 3)*c,
a(4) = (1/2)*(c + 3)*(c + 1)*c,
a(5) = (1/2)*(c^3 + 5*c^2 + 10*c + 4)*c,
a(6) = (1/4)*(2*c^4 + 15*c^3 + 38*c^2 + 37*c + 8)*c,
a(7) = (1/8)*(4*c^4 + 38*c^3 + 103*c^2 + 109*c + 22)*(c + 1)*c,
a(8) = (1/8)*(4*c^6 + 56*c^5 + 251*c^4 + 511*c^3 + 499*c^2 + 201*c + 22)*c,
and so on. No pattern is apparent in these formulae.
(End)

Extensions

Name edited by Petros Hadjicostas, Jan 14 2018
a(24)-a(29) from Petros Hadjicostas, Jan 14 2018

A032131 Shifts left 2 places under the "BIK" (reversible, indistinct, unlabeled) transform with a(1) = a(2) = 1.

Original entry on oeis.org

1, 1, 1, 2, 3, 7, 13, 31, 66, 160, 369, 907, 2191, 5461, 13558, 34209, 86426, 220359, 563475, 1449282, 3739365, 9688104, 25173917, 65621067, 171498288, 449361649, 1180078602, 3105740797, 8189749105, 21636207962, 57257857968, 151771200002, 402899862910, 1071076810324, 2851165864937
Offset: 1

Views

Author

Keywords

Comments

From Petros Hadjicostas, Jan 14 2018: (Start)
For this sequence, if (b(n): n>=1) = BIK((a(n): n>=1)), then b(n) = a(n+2) for n>=1.
Let A(x) = Sum_{n>=1} a(n)*x^n be the g.f. for this sequence. For an explanation on how to derive the formula BIK(A(x)) = (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1 - A(x^2))) from Bower's formulae in the link below about transforms, see the comments for sequence A001224. (For that sequence, the roles of sequences (a(n): n>=1) and (b(n): n>=1) are reversed.)
(End)

Crossrefs

Programs

  • Mathematica
    m = 36; a[1] = a[2] = 1; A[_] = 0;
    Do[A[x_] = x^2 (a[1]/x + a[2] + (1/2)(A[x]/(1 - A[x]) + (A[x] + A[x^2])/(1 - A[x^2]))) + O[x]^m // Normal, {m}];
    CoefficientList[A[x], x] // Rest (* Jean-François Alcover, Sep 17 2019 *)
  • PARI
    BIK(p)={(1/(1-p) + (1+p)/subst(1-p, x, x^2))/2}
    seq(n)={my(p=1+O(x^(n%2)));for(i=1, n\2, p=1+x*BIK(x*p)); Vec(p)} \\ Andrew Howroyd, Aug 30 2018

Formula

From Petros Hadjicostas, Jan 14 2018: (Start)
G.f.: If A(x) = Sum_{n>=1} a(n)*x^n, then (A(x) - a(1)*x - a(2)*x^2)/x^2 = BIK(A(x)) = (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1-A(x^2))). Here, a(1) = a(2) = 1.
In general, we have:
a(3) = a(1),
a(4) = (1/2)*(a(1)^2 + a(1) + 2*a(2)),
a(5) = (1/2)*(a(1)^2 + a(1) + 2*a(2) + 2)*a(1),
a(6) = (1/2)*(a(1)^4 + 4*a(1)^2 + (3*a(1)^2 + a(1) + 3)*a(2) + a(2)^2 + a(1)),
a(7) = (1/2)*(a(1)^4 + 4*a(1)^2*a(2) + 6*a(1)^2 + 3*a(2)^2 + 3*a(1) + 7*a(2) + 2)*a(1),
and so on. No pattern is apparent here.
(End)

Extensions

Name edited by Petros Hadjicostas, Jan 14 2018
a(31)-a(35) from Petros Hadjicostas, Jan 14 2018
Previous Showing 11-20 of 21 results. Next