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-3 of 3 results.

A232206 Triangle read by rows: T(n,k) is the number of non-equivalent regular polygons with n+1 edges, one of which is rooted, which are dissected by non-intersecting diagonals into k regions, such that two such polygons are identified up to reflection along the rooted edge and twisting along the diagonals that does not affect the root edge (for 1 <= k <= n-1 and n >= 2).

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 5, 8, 3, 1, 8, 22, 20, 6, 1, 11, 46, 73, 49, 11, 1, 15, 87, 206, 233, 119, 23, 1, 19, 147, 485, 807, 689, 288, 46, 1, 24, 236, 1021, 2320, 2891, 1988, 696, 98, 1, 29, 356, 1960, 5795, 9800, 9737, 5561, 1681, 207, 1, 35, 520, 3525, 13088, 28586, 38216, 31350, 15322, 4062, 451, 1, 41, 730, 5989, 27224, 74280, 127465, 139901, 97552, 41558, 9821, 983
Offset: 2

Views

Author

N. J. A. Sloane, Nov 21 2013

Keywords

Comments

From Petros Hadjicostas, Jan 20 2018: (Start)
According to Devadoss and Read (2001), the associahedron or Stasheff polytope K_n is a convex polytope of dim n-2 with co-dimension k faces corresponding to using k non-intersecting diagonals on a regular (n+1)-gon. If we redraw the picture of the dissected regular (n+1)-gon so that "every region of the dissection becomes a regular polygon without diagonals", then the resulting object is called a cluster, "where each regular polygon of the cluster is called a cell".
Two regular polygons G_1 and G_2 are of the same "dimension if their corresponding faces [on K_n] are of the same dimension"; of the same "type if their corresponding faces [on K_n] are the same polytope"; and of the same "class if they are identified under the actions of D_n [= dihedral group of order 2*n] and twisting [along the diagonals]".
In order to count how many faces of K_n belong to a given class (as defined above), Devadoss and Read (2001) follow the methodology of Read (1974, 1978) and first transform the regular (n+1)-gons (corresponding to the faces of K_n) into clusters, as mentioned above. Then they enumerate three kinds of clusters: those rooted at an outside edge (A-clusters); those rooted at a cell (B-clusters); and those rooted at an inside cell (C-clusters). In each one of these three kinds of clusters, two clusters are identified up to twisting along an edge (inside or outside depending on the case). Thus, the enumeration has to take twisting into account.
The g.f.s of these three kinds of clusters (A, B, and C) are derived, and these three g.f.s are used to derive the g.f. of the number of different classes of K_n of co-dimension k. See Section 3 in Devadoss and Read (2001).
The current triangular array, T(n,k), is the number of non-equivalent A-clusters (i.e., those rooted at an outside edge) having k cells and n outside edges, not counting the root edge, for 1 <= k <= n-1 and n>=2. Two such A-clusters are equivalent up to twisting the rooted edge or an inside edge separating two cells (as long as the inside twisting does not affect the rooted edge).
Read (1974, 1978) calls the A-clusters "out-rooted clusters". The difference between out-rooted clusters (= A-clusters) in Devadoss and Read (2001) and those in Read (1978) is that, in the first paper, two A-clusters are considered equivalent if one can be obtained from the other by twisting the rooted edge or an inside edge (separating two cells), while in the second paper twisting is not allowed.
If twisting is not allowed, then the number of (non-equivalent) A-clusters having k cells and n outside edges, not counting the root edge, is given by sequence A033282(n+1, k-1) = (1/k)*binomial(n-2,k-1)*binomial(n+k-1,k-1) for 1 <= k <= n-1 and n>=2. This is also the number of faces of K_n of co-dimension k-1. See Table 3 in Devadoss and Read (2001) and Table 1 in Read (1978). Unfortunately, in Table 1 in Read (1978), the label of each column is the number of outside edges including the rooted edge (i.e., n+1 rather than n).
Define T(n,k) = 0 for 0 <= n <= k, and T(n=1, k=0) = 1. Let also C(s,t) = Sum_{n>=0, k>=0} T(n,k)*s^n*t^k. (On p. 78 of their paper, Devadoss and Read (2001) use the notation a_{m,n} for T(n,m) and the notation A(x,y) for C(y,x), i.e., A(x,y) = Sum_{m>=0, n>=0} a_{m,n}*x^m*y^n in their paper.)
On p. 79, Eq. (3.2), of their paper, they prove that C(s,t) = s + (t/2)*(C(s,t)^2/(1-C(s,t)) + (1 + C(s,t))*C(s^2, t^2)/(1-C(s^2, t^2))). Unfortunately, the factor t 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).
We have C(s,t) = s+ t*s^2 + (t^2 + t)*s^3 + (2*t^3 + 3*t^2 + t)*s^4 + (3*t^4 + 8*t^3 + 5*t^2 + t)*s^5 + ...
Note that Sum_{0 <= k <= n-1} T(n,k) = A032132(n) for n>=1; T(n, k=1) = 1 for n>=2; T(k+1,k) = A001190(k+1) for k>=1, which are the Wedderburn-Etherington numbers; and T(n, k=2) = A024206(n-1) for n>=2.
According to Schöbel and Veselov (2014, 2015), for n>=2 and 0 <= p <= n-1, T(n+1, n-p) = A295380(n,p) is the number of inequivalent canonical forms for separation coordinates of the hypersphere S^n with p independent continuous parameters.
(End)

Examples

			From _Petros Hadjicostas_, Jan 20 2018: (Start)
According to Devadoss and Read (2015), triangle T(n,k) begins as follows:
(n=2)  1,
(n=3)  1,  1,
(n=4)  1,  3,   2,
(n=5)  1,  5,   8,    3,
(n=6)  1,  8,  22,   20,    6,
(n=7)  1, 11,  46,   73,   49,   11,
(n=8)  1, 15,  87,  206,  233,  119,   23,
(n=9)  1, 19, 147,  485,  807,  689,  288,   46,
(n=10) 1, 24, 236, 1021, 2320, 2891, 1988,  696,   98,
(n=11) 1, 29, 356, 1960, 5795, 9800, 9737, 5561, 1681, 207,
...
Geometrical example for n=5: If no twisting is allowed, the number of regular (n+1)-gons (= hexagons) with a rooted edge and dissected into k regions by non-intersecting diagonals is given by A033282(n+1, k-1) = A033282(6, k-1): 1, 9, 21, and 14 for k = 1, 2, 3, 4, respectively.
Recall that, according to Devadoss and Read (2001), two regular (unrooted) polygons G_1 and G_2 are of the same "class" if they are identified under the actions of D_n (= dihedral group of order 2*n) and twisting along the diagonals. To avoid confusion, we call two such unrooted equivalent polygons as being of the same DT-class (D for "dihedral" and T for "twisting"). According to Table 5 (p. 92) in Devadoss and Read (2001), the numbers of DT-classes for regular hexagons dissected into k regions (by k-1 non-intersecting diagonals) are 1, 2, 3, 2 for k = 1, 2, 3, 4, respectively. When an edge is rooted, each one of these DT-classes is subdivided into subclasses.
Call the regular hexagons ABCDEF with FA being the rooted edge (base).
For k=1, we have T(5,1) = 1 rooted regular hexagon with no intersecting diagonals.
For k=2, we have T(5,2) = 5 equivalent rooted regular hexagons with 1 diagonal. The equivalence classes are as follows according to the single dissecting diagonal:
Class 1a: AE, FB; Class 1b: AC, FD; Class 1c: BD, EC; (DT-class 1 has 6 hexagons)
Class 2a: AD, FC; Class 2b: BE; (DT-class 2 has 3 hexagons).
Note that 6 + 3 = 9 =  A033282(5+1, 2-1).
For k=3, we have T(5,3) = 8 equivalent rooted regular hexagons with 2 non-intersecting diagonals. The equivalence classes are as follows according to the two diagonals:
Class 1a: (AC, AE), (FB, FD), (BD, BF), (EA, EC);
Class 1b: (DB, DF), (CE, CA); (DT-class 1 has 6 hexagons)
Class 2a: (DB, EA), (CE, BF);
Class 2b: (DF, CA); (DT-class 2 has 3 hexagons)
Class 3a: (DA, EA), (FC, FB), (EA, EB), (BE, BF);
Class 3b: (AC, AD), (FC, FD), (DA, DB), (CE, CF);
Class 3c: (BD, BE), (EC, EB);
Class 3d: (CF, CA), (DA, DF); (DT-class 3 has 12 hexagons)
Note that 6 + 3 + 12 = 21 = A033282(5+1, 3-1).
For k=4, we have T(5,4) = 3 equivalent rooted regular hexagons with 3 non-intersecting diagonals. The equivalence classes are as follows according to the three diagonals:
Class 1: (EA, AC, CE), (BD, DF, FB); (DT-class 1 has 2 hexagons)
Class 2a: (DF, DA, AC), (DF, FC, CA), (CE, CF, CA), (DF, DA, DB);
Class 2b: (CE, BE, BF), (BD, BE, BF), (EC, EB, EA), (DB, EB, EA), (EC, CF, FB), (DB, DA, AE), (FD, FC, FB), (AE, AD, AC); (DT class 2 has 12 hexagons)
Note that 2 + 12 = 14 = A033282(5+1, 4-1).
Recall that two rooted hexagons are equivalent iff they are a reflection of each other along the rooted edge or one can be obtained from the other by twisting a diagonal as long as the twisting does not affect the rooted edge.
The case k = 4 = n-1 above is related to the triangulation of a convex polygon and the Wedderburn-Etherington commutative bracketing problem that appear in Comtet (1974, pp. 54-55). Devadoss and Read (2001, p. 89) claim that T(n,k) is the number of possible ways of using k-1 pairs of brackets on n commutative variables, but it is not clear how each one of the above hexagons (from k=1 to k=3) can be transformed into some kind of a generalized commutative bracketing problem.
(End)
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, pp. 54-55 and 74-75.

Crossrefs

Cf. A001190 (main diagonal), A024206 (second column), A032132 (row sums), A295380 (mirror image).

Programs

  • Mathematica
    rows = 12; A[, ] = 0;
    Do[A[s_, t_] = Series[s+(t/2)(A[s, t]^2/(1-A[s, t])+(1+A[s, t]) A[s^2, t^2] / (1-A[s^2, t^2])), {s, 0, rows+1}, {t, 0, rows+1}] // Normal, {rows+1}];
    Rest[CoefficientList[#, t]]& /@ Drop[CoefficientList[A[s, t], s], 2] // Flatten (* Jean-François Alcover, Sep 02 2018 *)
  • PARI
    BIKxy(p)={(1/(1-p) + (1+p)/subst(subst(1-p, x, x^2), y, y^2))/2}
    A(n)={my(p=x); for(i=2, n+1, p+=x^i*y*polcoeff(BIKxy(p + O(x*x^i)), i)); [Vecrev(q/y) | q<-Vecrev((p-x)/x^2)]}
    { my(T=A(10)); for(n=1, #T, print(T[n])) } \\ Andrew Howroyd, Aug 31 2018
  • SageMath
    #This is a crude Sage program on how to derive the g.f. of any column in the triangle given that you have correctly derived the g.f.s of the previous columns. Suppose you have the expansion for C(s,t) w.r.t. t up to the power of 3 and you want to get the coefficient of t^4.
    s,t=var('s,t')
    C(s,t) = s + t*s^2/(1-s) + t^2*s^3*(1+s-s^2)/((1+s)*(1-s)^3) + t^3*s^4*(2+2*s-2*s^3+s^4)/((1+s)^2*(1-s)^5)
    B(s,t)=s+(t/2)*(C(s,t)^2/(1-C(s,t))+(1+C(s,t))*C(s^2,t^2)/(1-C(s^2,t^2)))-C(s,t)
    factor(taylor(B(s,t),t,0,4)/t^4)
    # Petros Hadjicostas, Jan 21 2018
    

Formula

From Petros Hadjicostas, Jan 20 2018: (Start)
G.f.: If C(s,t) = Sum_{n>=0, k>=0} T(n,k)*s^n*t^k (with T(n=1,k=0) = 1), then C(s,t) = s + (t/2)*(C(s,t)^2/(1-C(s,t)) + (1 + C(s,t))*C(s^2, t^2)/(1-C(s^2, t^2))).
Note that t*C(s,t) = Sum_{k>=1} T(k,k-1)*(s*t)^k + Sum_{k>=2} (s*t)^k*Sum_{n >= k+1} T(n,k-1)*s^{n-k}. If we let w = s*t and D(s,w) = (w/s)*C(s,w/s), then the above functional equation implies that E(w): = lim_{s -> 0} D(s,w) = Sum_{k>=1} T(k,k-1)*w^k satisfies E(w) = w + (1/2)*(E(w)^2 + E(w^2)), which is the g.f. for the Wedderburn-Etherington numbers in sequence A001190. Thus, T(k,k-1) = A001190(k) for k>=1.
Expansion w.r.t to t: C(s,t) = s + t*s^2/(1-s) + t^2*s^3*(1+s-s^2)/((1+s)*(1-s)^3) + t^3*s^4*(2+2*s-2*s^3+s^4)/((1+s)^2*(1-s)^5) + t^4*s^5*(3+8*s+2*s^2-2*s^3-2*s^4+3*s^5-s^6)/((1+s)^3*(1-s)^7) + t^5*s^6*(6+19*s+30*s^2+15*s^3+23*s^4-4*s^5+2*s^6-4*s^7+6*s^8-4*s^9+s^10)/((1+s^2)*(1+s)^4*(1-s)^9) + ... (It can be proven using a symbolic computation package by manipulating the above functional equation by Devadoss and Read.)
This means that, in the above expansion for C(s,t), the coefficient of t is the g.f. of the first column of the triangle below, the coefficient of t^2 is the g.f. of the second column of the triangle below, and so on.
(End)

Extensions

Name edited by Petros Hadjicostas, Jan 20 2018
Offset corrected by Andrew Howroyd, Aug 31 2018

A032133 Number of series-reduced dyslexic planted planar trees with n leaves of 2 colors.

Original entry on oeis.org

2, 3, 12, 61, 340, 2070, 13178, 87115, 590720, 4089214, 28766442, 205072330, 1478149486, 10754615963, 78878136418, 582563752101, 4328922702536, 32341244773080, 242781826475456, 1830377795460438
Offset: 1

Views

Author

Keywords

Crossrefs

Cf. A032132.

Programs

  • PARI
    BIK(p)={(1/(1-p) + (1+p)/subst(1-p, x, x^2))/2}
    seq(n)={my(p=2*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.

A295380 Number of canonical forms for separation coordinates on hyperspheres S_n, ordered by increasing number of independent continuous parameters.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 3, 8, 5, 1, 6, 20, 22, 8, 1, 11, 49, 73, 46, 11, 1, 23, 119, 233, 206, 87, 15, 1, 46, 288, 689, 807, 485, 147, 19, 1, 98, 696, 1988, 2891, 2320, 1021, 236, 24, 1, 207, 1681, 5561, 9737, 9800, 5795, 1960, 356, 29, 1, 451, 4062, 15322, 31350, 38216, 28586, 13088, 3525, 520, 35, 1, 983, 9821, 41558, 97552, 139901, 127465, 74280, 27224, 5989, 730, 41, 1
Offset: 1

Views

Author

Tom Copeland, Nov 21 2017

Keywords

Comments

Table 1 of the Schöbel and Veselov paper with initial 1 added. Reverse of Table 2 of the Devadoss and Read paper.
Apparently A032132 contains the row sums.
From Petros Hadjicostas, Jan 28 2018: (Start)
In this triangle, which is read by rows, for 0 <= k <= n-1 and n>=1, let T(n,k) be the number of inequivalent canonical forms for separation coordinates of the hypersphere S^n with k independent continuous parameters. It is the mirror image of sequence A232206, that is, T(n, k) = A232206(n+1, n-k) for 0 <= k <= n-1 and n>=1. (Triangular array A232206(N, K) is defined for N >= 2 and 1 <= K <= N-1.)
If B(x,y) = Sum_{n,k>=0} T(n,k)*x^n*y^k (with T(0,0) = 1, T(0,k) = 0 for k>=1, and T(n,k) = 0 for 1 <= n <= k), then B(x,y) = 1 + (x/2)*(B(x,y)^2/(1-x*y*B(x,y)) + (1 + x*y*B(x,y))*B(x^2,y^2)/(1-x^2*y^2*B(x^2,y^2))). This can be derived from the bivariate g.f. of A232206. See the comments for that sequence.
Let S(n) := Sum_{k>=0} T(n,k). The g.f. of S(n) is B(x, y=1). If we let y=1 in the above functional equation, we get x*B(x,1) = x + (1/2)*((x*B(x,1))^2/(1-x*B(x,1)) + (1 + x*B(x,1))*x^2*B(x^2,1)/(1-x^2*B(x^2,1))). After some algebra, we get 2*x*B(x,1) = x + (1/2)(x*B(x,1)/(1-x*B(x,1)) + (x*B(x,1) + x^2*B(x^2,1))/(1-x^2*B(x,1))), i.e., 2*x*B(x,1) = x + BIK(x*B(x,1)), where we have the "BIK" (reversible, indistinct, unlabeled) transform of C. G. Bower. This proves that S(n) = A032132(n+1) for n>=0, which is Copeland's claim above.
Note that for the second column we have T(n,k=2) = A048739(n-2) for 2 <= n < = 10, but T(11,2) = 4062 <> 4059 = A048739(9). In any case, they have different g.f.s (see the formula section below).
(End)

Examples

			From _Petros Hadjicostas_, Jan 27 2018: (Start)
Triangle T(n,k) begins:
n\k      0     1     2     3     4     5     6    7   8  9
----------------------------------------------------------------
(S^1)    1,
(S^2)    1,    1,
(S^3)    2,    3,    1,
(S^4)    3,    8,    5,    1,
(S^5)    6,   20,   22,    8,    1,
(S^6)   11,   49,   73,   46,   11,    1,
(S^7)   23,  119,  233,  206,   87,   15,    1,
(S^8)   46,  288,  689,  807,  485,  147,   19,   1,
(S^9)   98,  696, 1988, 2891, 2320, 1021,  236,  24,  1,
(S^10) 207, 1681, 5561, 9737, 9800, 5795, 1960, 356, 29, 1,
...
(End)
		

Crossrefs

Formula

From Petros Hadjicostas, Jan 28 2018: (Start)
G.f.: If B(x,y) = Sum_{n,k>=0} T(n,k)*x^n*y^k (with T(0,0) = 1, T(0,k) = 0 for k>=1, and T(n,k) = 0 for 1 <= n <= k), then B(x,y) = 1 + (x/2)*(B(x,y)^2/(1-x*y*B(x,y)) + (1 + x*y*B(x,y))*B(x^2,y^2)/(1-x^2*y^2*B(x^2,y^2))).
If c(N,K) = A232206(N,K) and C(x,y) = Sum_{N,K>=0} c(N,K)*x^N*y^K (with c(1,0) = 1 and c(N,K) = 0 for 0 <= N <= K), then C(x,y) = x*B(x*y, 1/y) and B(x,y) = C(x*y, 1/y)/(x*y).
Setting y=0 in the above functional equation, we get x*B(x,0) = x + (1/2)*((x*B(x,0))^2 + x^2*B(x^2,0)), which is the functional equation for the g.f. of the first column. This proves that T(n,k=0) = A001190(n+1) for n>=0 (assuming T(0,0) = 1).
The g.f. of the second column is B_1(x,0) = Sum_{n>=0} T(n,2)*x^n = lim_{y->0} (B(x,y)-B(x,0))/y, where B(x,0) = 1 + x + x^2 + ... is the g.f. of the first column. We get B_1(x,0) = x*B(x,0)*(B(x,0) - 1)/(1 - x*B(x,0)).
(End)

Extensions

Typo for T(11,3)=15322 corrected by Petros Hadjicostas, Jan 28 2018
Showing 1-3 of 3 results.