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 21-30 of 47 results. Next

A244306 Table T(n,k), n>=1, k>=1, read by antidiagonals: T(n,k) = number of equivalence classes of ways of placing two 1 X 1 tiles in an n X k rectangle under all symmetry operations of the rectangle.

Original entry on oeis.org

0, 1, 1, 2, 3, 2, 4, 6, 6, 4, 6, 10, 13, 10, 6, 9, 15, 22, 22, 15, 9, 12, 21, 34, 36, 34, 21, 12, 16, 28, 48, 56, 56, 48, 28, 16, 20, 36, 65, 78, 88, 78, 65, 36, 20, 25, 45, 84, 106, 123, 123, 106, 84, 45, 25, 30, 55, 106, 136, 168, 171, 168, 136, 106, 55, 30
Offset: 1

Views

Author

Keywords

Examples

			T(n,k) for 1<=n<=11 and 1<=k<=11 is:
    k  1    2    3    4    5    6    7    8    9   10   11 ...
.n
.1     0    1    2    4    6    9   12   16   20   25   30
.2     1    3    6   10   15   21   28   36   45   55   66
.3     2    6   13   22   34   48   65   84  106  130  157
.4     4   10   22   36   56   78  106  136  172  210  254
.5     6   15   34   56   88  123  168  216  274  335  406
.6     9   21   48   78  123  171  234  300  381  465  564
.7    12   28   65  106  168  234  321  412  524  640  777
.8    16   36   84  136  216  300  412  528  672  820  996
.9    20   45  106  172  274  381  524  672  856 1045 1270
10    25   55  130  210  335  465  640  820 1045 1275 1550
11    30   66  157  254  406  564  777  996 1270 1550 1885
		

Crossrefs

Formula

Empirically,
T(n,k) = (4*k^2*n^2 + 2*k^2 + 8*k*n + 2*n^2 - 4*k - 4*n - 1 - (2*k^2 - 4*k - 1)*(-1)^n - (2*n^2 - 4*n - 1)*(-1)^k - (-1)^k*(-1)^n)/32.
T(1,k) = A002620(k) = floor(k^2/4).
T(2,k) = A000217(k) = k*(k+1)/2.
= T(1,k) + T(1,k+1) = floor(k^2/4) + floor((k+1)^2/4).
T(3,k) = 2*A000217(k) + A024206(k-2)
= k*(k+1) + floor((k-1)^2/4) - 1.

Extensions

Terms corrected and extended by Christopher Hunt Gribble, Apr 02 2015

A055084 Number of 6 X n binary matrices with no zero rows or columns, up to row and column permutation.

Original entry on oeis.org

1, 15, 180, 2001, 20755, 200082, 1781941, 14637962, 111011667, 779695050, 5093379110, 31092553357, 178203364143, 963217652830, 4930357535218, 23989343505296, 111335037107474, 494383391324356, 2106346854756098
Offset: 1

Views

Author

Vladeta Jovovic, Jun 13 2000

Keywords

Crossrefs

Column k=6 of A056152.

Programs

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

A055082 Number of 4 X n binary matrices with no zero rows or columns, up to row and column permutation.

Original entry on oeis.org

1, 8, 42, 179, 633, 2001, 5745, 15274, 38000, 89331, 199715, 427184, 878152, 1741964, 3345562, 6239390, 11327863, 20065972, 34747460, 58924066, 98002370, 160086580, 257148244, 406637336, 633669040, 973971441, 1477810227, 2215179768, 3282598034, 4811946882
Offset: 1

Views

Author

Vladeta Jovovic, Jun 13 2000

Keywords

Crossrefs

Column k=4 of A056152.

Programs

Extensions

Terms a(21) and beyond from Andrew Howroyd, Mar 25 2020

A055083 Number of 5 X n binary matrices with no zero rows or columns, up to row and column permutation.

Original entry on oeis.org

1, 11, 91, 633, 3835, 20755, 102089, 461272, 1930310, 7534742, 27602968, 95428291, 312864361, 976985630, 2917175450, 8357692894, 23046527311, 61337188725, 157950527167, 394427897066, 957058104818, 2260601179661, 5206447640059, 11709619965923, 25752660738209
Offset: 1

Views

Author

Vladeta Jovovic, Jun 13 2000

Keywords

Crossrefs

Column k=5 of A056152.

Programs

Extensions

Terms a(21) and beyond from Andrew Howroyd, Mar 25 2020

A078126 Negative determinant of n X n matrix M_{i,j}=1 if i=j or i+j=1 (mod 2).

Original entry on oeis.org

-1, -1, 0, 1, 3, 5, 8, 11, 15, 19, 24, 29, 35, 41, 48, 55, 63, 71, 80, 89, 99, 109, 120, 131, 143, 155, 168, 181, 195, 209, 224, 239, 255, 271, 288, 305, 323, 341, 360, 379, 399, 419, 440, 461, 483, 505, 528, 551, 575, 599, 624, 649, 675, 701, 728, 755, 783, 811
Offset: 0

Views

Author

Michael Somos, Nov 18 2002

Keywords

Comments

Apparently, also 6(n+3) times the Dedekind sum s(2,n+3). - Ralf Stephan, Sep 16 2013

Examples

			G.f. = -1 - x + x^3 + 3*x^4 + 5*x^5 + 8*x^6 + 11*x^7 + 15*x^8 + 19*x^9 + ...
		

Crossrefs

Programs

  • Maple
    A078126:=n->floor((n + 2)*(n - 2)/4); seq(A078126(n), n=0..100); # Wesley Ivan Hurt, Jan 30 2014
  • Mathematica
    Table[Floor[(n + 2)(n - 2)/4], {n, 0, 100}] (* Wesley Ivan Hurt, Jan 30 2014 *)
    LinearRecurrence[{2,0,-2,1},{-1,-1,0,1},60] (* Harvey P. Dale, Sep 10 2015 *)
    a[n_]:=-Det[Table[If[i==j ||Mod[i+j,2]==1,1,0],{i,n},{j,n}]]; Join[{-1},Array[a,57]] (* Stefano Spezia, Aug 06 2024 *)
  • PARI
    a(n)=-matdet(matrix(n,n,i,j,i==j||((i+j)%2))) /* Ralf Stephan, Sep 16 2013 */
    
  • PARI
    a(n)=sumdedekind(2,n+3)*6*(n+3) /* Ralf Stephan, Sep 16 2013 */

Formula

G.f.: (-1 + x + 2*x^2 - x^3) / ((1 - x^2) * (1 - x)^2).
a(n) = A002620(n) - 1.
a(n) = A002623(n-2) - A002623(n-3) - 1.
a(n) = A024206(n-1) for all n in Z.
a(n) = floor( (n+2)(n-2)/4 ). - Wesley Ivan Hurt, Jun 16 2013
A004526(n) = a(n) - a(n-1) for all n in Z. - Michael Somos, Aug 22 2016
a(n) = Sum_{i=1..n+2} floor((n-i+1)/2). - Wesley Ivan Hurt, Sep 12 2017
E.g.f.: ((x^2 + x - 4)*cosh(x) + (x^2 + x - 5)*sinh(x))/4. - Stefano Spezia, Aug 06 2024

Extensions

A-number twister corrected in cross-refs by R. J. Mathar, Feb 11 2010

A079524 Expansion of (x + b*x^2 - b*x^3)/((1 - x^2)*(1 - x)^2) with b=2.

Original entry on oeis.org

0, 1, 4, 6, 10, 13, 18, 22, 28, 33, 40, 46, 54, 61, 70, 78, 88, 97, 108, 118, 130, 141, 154, 166, 180, 193, 208, 222, 238, 253, 270, 286, 304, 321, 340, 358, 378, 397, 418, 438, 460, 481, 504, 526, 550, 573, 598, 622, 648, 673, 700, 726, 754, 781, 810, 838
Offset: 0

Views

Author

Carlos Alves, Jan 21 2003

Keywords

Examples

			G.f. = x + 4*x^2 + 6*x^3 + 10*x^4 + 13*x^5 + 18*x^6 + 22*x^7 + 28*x^8 + 33*x^9 + ...
Here b=2; a(0)=a(1)=1. a(2)= a(0)+2 (mod a(1)+2) = 3 (mod 3) =0 a(3)= a(1)+3 (mod a(2)+3) = 4 (mod 3) =1 a(4)= a(2)+4 (mod a(3)+4) = 4 (mod 5) =4 etc... we get 6, 13, 18, ...
		

Crossrefs

Cf. A024206 and A078126.

Programs

  • GAP
    a:=[0,1,4,6];; for n in [5..50] do a[n]:=2*a[n-1]-2*a[n-3]+a[n-4]; od; a; # G. C. Greubel, Jan 15 2019
  • Magma
    m:=50; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!( (x+2*x^2-2*x^3)/((1-x)^2*(1- x^2)) )); // G. C. Greubel, Jan 15 2019
    
  • Maple
    a:= n-> (Matrix([[4,1,0,-2]]). Matrix(4, (i,j)-> if (i=j-1) then 1 elif j=1 then [2,0,-2,1][i] else 0 fi)^n)[1,4]: seq(a(n), n=1..60); # Alois P. Heinz, Aug 06 2008
  • Mathematica
    b = 2; aa = {1, 1}; Do[AppendTo[aa, Mod[ aa[[ -2]] + n, aa[[ -1]] + n]], {n, b, 50}]; Drop[aa, 2]
    CoefficientList[Series[(x+2*x^2-2*x^3)/((1-x)^2*(1-x^2)),{x,0,60}],x]
    LinearRecurrence[{2,0,-2,1},{0,1,4,6},50] (* Harvey P. Dale, Apr 20 2015 *)
  • PARI
    my(x='x+O('x^50)); concat([0], Vec( (x+2*x^2-2*x^3)/((1-x)^2*(1- x^2)) )) \\ G. C. Greubel, Jan 15 2019
    
  • Sage
    ((x+2*x^2-2*x^3)/((1-x)^2*(1- x^2))).series(x, 50).coefficients(x, sparse=False) # G. C. Greubel, Jan 15 2019
    

Formula

G.f.: (x + 2*x^2 - 2*x^3) / ((1 - x)^2 * (1 - x^2)).
a(n) = a(n-2)+n (mod a(n-1)+n) with n>=2 and initial values (1, 1).
a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4), a(0)=0, a(1)=1, a(2)=4, a(3)=6. - Harvey P. Dale, Apr 20 2015
a(n) = (2*n*(n+6)-3*(1-(-1)^n))/8. - Luce ETIENNE, Jun 05 2015

A128223 a(n) = if n mod 2 = 0 then n*(n+1)/2 otherwise (n+1)^2/2-1.

Original entry on oeis.org

0, 1, 3, 7, 10, 17, 21, 31, 36, 49, 55, 71, 78, 97, 105, 127, 136, 161, 171, 199, 210, 241, 253, 287, 300, 337, 351, 391, 406, 449, 465, 511, 528, 577, 595, 647, 666, 721, 741, 799, 820, 881, 903, 967, 990, 1057, 1081, 1151, 1176, 1249, 1275, 1351, 1378, 1457, 1485
Offset: 0

Views

Author

Gary W. Adamson, Feb 19 2007

Keywords

Comments

a(n-1) is the length of the shortest path along the edges of the complete graph with n vertices. - Martin Fuller, Dec 06 2007
From Peter Kagey, Jan 25 2015: (Start)
For an irreflexive, non-transitive, symmetric relation, a(n) is the length of a relation chain required to demonstrate that a != b for all distinct elements a and b in S, where S contains n+1 elements.
For example, for the set {1,2,3} the chain requires a(2) = 3 relations (e.g., 1 != 2 != 3 != 1). For the set {1,2,3,4}, the chain requires a(3) = 7 relations (e.g., 1 != 2 != 3 != 4 != 1 != 3 != 2 != 4 -- noting the redundancy of 2!=3 and 3!=2). (End)
Given a set of n lots of n distinct items, it is possible to sort the items from fully collated (ABCABCABC) to fully sorted (AAABBBCCC), or vice versa, using a sorting algorithm whereby at each step a portion of the overall string is selected and its contents reversed. The minimum number of steps such an algorithm will take is a(n-1). For example, when n=3, a(n-1)=3: ABCABCABC -> ABBACBACC -> ABBAABCCC -> AAABBBCCC. - Elliott Line, Aug 02 2019

Examples

			a(5) = 17 = (5 + 1 + 5 + 1 + 5), row 5 of A128222.
		

Crossrefs

Row sums of A128222.
Cf. A024206, row sums of A128221 = A128174 * A127701.

Programs

  • Haskell
    a128223 n = if even n then n*(n + 1) `div` 2 else (n+1)^2 `div` 2 - 1 -- Peter Kagey, Jul 14 2015
    
  • Magma
    [(-1+(-1)^n-(-3+(-1)^n)*n+2*n^2)/4: n in [0..60]]; // Vincenzo Librandi, Mar 18 2015
    
  • Maple
    f:=n-> if n mod 2 = 0 then n*(n+1)/2 else (n+1)^2/2-1; fi;
  • Mathematica
    f[n_] := If[EvenQ@ n, n (n + 1)/2, (n + 1)^2/2 - 1]; Array[f, 54] (* Michael De Vlieger, Mar 17 2015 *)
    Table[(- 1 + (-1)^n - (- 3 + (-1)^n) n + 2 n^2) / 4, {n, 0, 60}] (* Vincenzo Librandi, Mar 18 2015 *)
    CoefficientList[ Series[(-x - 2x^2 - 2x^3 + x^4)/((-1 + x)^3 (1 + x)^2), {x, 0, 54}], x] (* Robert G. Wilson v, Nov 16 2016 *)
    LinearRecurrence[{1,2,-2,-1,1},{0,1,3,7,10},60] (* Harvey P. Dale, Mar 17 2020 *)
  • PARI
    main(size)={my(n,m,v=vector(size),i);for(i=0,size-1,v[i+1]=if(i%2==0,i*(i+1)/2,(i+1)^2/2-1));return(v);} /* Anders Hellström, Jul 14 2015 */

Formula

a(n) = (-1+(-1)^n-(-3+(-1)^n)*n+2*n^2)/4. a(n) = a(n-1)+2*a(n-2)-2*a(n-3)-a(n-4)+a(n-5). G.f.: x*(x^3-2*x^2-2*x-1) / ((x-1)^3*(x+1)^2). - Colin Barker, Oct 16 2013
a(n) = A053439(n) - 1. - Peter Kagey, Nov 16 2016

Extensions

Edited by N. J. A. Sloane, Dec 06 2007

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

A327615 Irregular triangle read by rows: T(n,k) is the number of unlabeled multigraphs with loops allowed and n edges covering k vertices, n >= 1, 1 <= k <= 2*n.

Original entry on oeis.org

1, 1, 1, 3, 2, 1, 1, 5, 8, 6, 2, 1, 1, 8, 19, 25, 16, 7, 2, 1, 1, 11, 40, 73, 73, 47, 19, 7, 2, 1, 1, 15, 77, 194, 263, 232, 133, 58, 20, 7, 2, 1, 1, 19, 132, 454, 835, 951, 719, 397, 164, 61, 20, 7, 2, 1, 1, 24, 217, 984, 2385, 3507, 3365, 2306, 1177, 490, 175, 62, 20, 7, 2, 1
Offset: 1

Views

Author

Andrew Howroyd, Oct 23 2019

Keywords

Comments

Covering k vertices means there are no vertices of degree zero.

Examples

			Triangle begins:
  1,  1;
  1,  3,   2,   1;
  1,  5,   8,   6,   2,   1;
  1,  8,  19,  25,  16,   7,   2,   1;
  1, 11,  40,  73,  73,  47,  19,   7,   2,  1;
  1, 15,  77, 194, 263, 232, 133,  58,  20,  7,  2, 1;
  1, 19, 132, 454, 835, 951, 719, 397, 164, 61, 20, 7, 2, 1;
  ...
		

Crossrefs

Row sums are A007717.
Columns k=2..3 are A024206, A327728.

Programs

  • PARI
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    edges(v, t) = {prod(i=2, #v, prod(j=1, i-1, my(g=gcd(v[i], v[j])); t(v[i]*v[j]/g)^g )) * prod(i=1, #v, my(c=v[i]); t(c)^((c+1)\2)*if(c%2, 1, t(c/2)))}
    C(n,m)={my(s=O(x*x^m)); forpart(p=n, s+=permcount(p)/edges(p, i->1-x^i+O(x*x^m))); Col(s/n!)}
    T(m) = {my(n=2*m, A=Mat(vector(n+1, n, C(n-1,m)))); A[2..m+1, 2..n+1]-A[2..m+1, 1..n]}
    { my(A=T(8)); for(n=1, matsize(A)[1], print(A[n, 1..2*n])) }

Formula

T(n,k) = A290428(n,k) - A290428(n,k-1).
Previous Showing 21-30 of 47 results. Next