A325000 -id:A325000 - cp's OEIS frontend

A337885 Array read by descending antidiagonals: T(n,k) is the number of chiral pairs of colorings of the triangular faces of a regular n-dimensional simplex using k or fewer colors.

0, 0, 0, 0, 0, 0, 0, 0, 6, 0, 0, 1, 405, 1368, 0, 0, 5, 7904, 4775706, 6711288, 0, 0, 15, 76880, 1522540416, 9923557498416, 1785683627824, 0, 0, 35, 486522, 132342705750, 234239763858347776, 12979826761630383196344, 53302696800142157920, 0

Offset: 2

Robert A. Russell, Sep 28 2020

Each member of a chiral pair is a reflection, but not a rotation, of the other. An n-simplex has n+1 vertices. For n=2, the figure is a triangle with one triangular face. For n=3, the figure is a tetrahedron with 4 triangular faces. For higher n, the number of triangular faces is C(n+1,3).

Also the number of chiral pairs of colorings of the peaks of a regular n-dimensional simplex. A peak of an n-simplex is an (n-3)-dimensional simplex.

			Table begins with T(2,1):
 0    0       0          0            0             0               0 ...
 0    0       0          1            5            15              35 ...
 0    6     405       7904        76880        486522         2300305 ...
 0 1368 4775706 1522540416 132342705750 5076500214744 110809322249220 ...
For T(3,4)=1, the chiral pair is ABCD-ABDC.

E. M. Palmer and R. W. Robinson, Enumeration under two representations of the wreath product, Acta Math., 131 (1973), 123-143.

Cf. A337883 (oriented), A337884 (unoriented), A337886 (achiral), A051168 (binary Lyndon words).
Other elements: A325000(n,k-n) (vertices), A327085 (edges).
Other polytopes: A337889 (orthotope), A337893 (orthoplex).
Rows 2-4 are A000004, A000332, A331352.

m=2; (* dimension of color element, here a triangular face *)
lw[n_,k_]:=lw[n, k]=DivisorSum[GCD[n,k],MoebiusMu[#]Binomial[n/#,k/#]&]/n (*A051168*)
cxx[{a_, b_},{c_, d_}]:={LCM[a, c], GCD[a, c] b d}
compress[x:{{, } ...}] := (s=Sort[x];For[i=Length[s],i>1,i-=1,If[s[[i,1]]==s[[i-1,1]], s[[i-1,2]]+=s[[i,2]]; s=Delete[s,i], Null]]; s)
combine[a : {{, } ...}, b : {{, } ...}] := Outer[cxx, a, b, 1]
CX[p_List, 0] := {{1, 1}} (* cycle index for partition p, m vertices *)
CX[{n_Integer}, m_] := If[2m>n, CX[{n}, n-m], CX[{n},m] = Table[{n/k, lw[n/k, m/k]}, {k, Reverse[Divisors[GCD[n, m]]]}]]
CX[p_List, m_Integer] := CX[p, m] = Module[{v = Total[p], q, r}, If[2 m > v, CX[p, v - m], q = Drop[p, -1]; r = Last[p]; compress[Flatten[Join[{{CX[q, m]}}, Table[combine[CX[q, m - j], CX[{r}, j]], {j, Min[m, r]}]], 2]]]]
pc[p_] := Module[{ci, mb}, mb = DeleteDuplicates[p]; ci = Count[p, #] &/@ mb; Total[p]!/(Times @@ (ci!) Times @@ (mb^ci))] (* partition count *)
row[n_Integer] := row[n] = Factor[Total[If[EvenQ[Total[1-Mod[#, 2]]],1,-1] pc[#] j^Total[CX[#, m+1]][[2]] & /@ IntegerPartitions[n+1]]/(n+1)!]
array[n_, k_] := row[n] /. j -> k
Table[array[n,d+m-n], {d,8}, {n,m,d+m-1}] // Flatten

The algorithm used in the Mathematica program below assigns each permutation of the vertices to a partition of n+1. It then determines the number of permutations for each partition and the cycle index for each partition using a formula for binary Lyndon words. If the value of m is increased, one can enumerate colorings of higher-dimensional elements beginning with T(m,1).

T(n,k) = A337883(n,k) - A337884(n,k) = (A337883(n,k) - A337886(n,k)) / 2 = A337884(n,k) - A337886(n,k).

A243662 Triangle read by rows: the reversed x = 1+q Narayana triangle at m=2.

Original entry on oeis.org

1, 3, 1, 12, 8, 1, 55, 55, 15, 1, 273, 364, 156, 24, 1, 1428, 2380, 1400, 350, 35, 1, 7752, 15504, 11628, 4080, 680, 48, 1, 43263, 100947, 92169, 41895, 9975, 1197, 63, 1, 246675, 657800, 708400, 396704, 123970, 21560, 1960, 80, 1, 1430715, 4292145, 5328180, 3552120, 1381380, 318780, 42504, 3036, 99, 1

Offset: 1

N. J. A. Sloane, Jun 13 2014

See Novelli-Thibon (2014) for precise definition.
From Tom Copeland, Dec 13 2022: (Start)
The row polynomials are the nonvanishing numerator polynomials generated in the compositional, or Lagrange, inversion in x about the origin of the odd o.g.f. Od1(x,t) = x*(t*(1-x^2)-x^2) / (1-x^2) = t*x - x^3 - x^5 - x^7 - x^9 - ... .
For example, from the Lagrange inversion formula (LIF), the tenth derivative in x of (x/Od1(x,t))^11 / 11! = (1/((t*(1-x^2)-x^2) / (1-x^2)))^11 / 11! at x = 0 is (t^4 + 24*t^3 + 156*t^2 + 364*t + 273) / t^16. These polynomials are also generated by the iterated derivatives ((1/(D Od1(x,t)) D)^n g(x) evaluated at x = 0 where D = d/dx.
An explicit generating function for the polynomials can be obtained by finding the solution of the cubic equation y - t*x - y*x^2 + (1+t)*x^3 = 0 for x in terms of y and t that satisfies y(x=0;t) = 0 = x(y=0;t).
The row polynomials are also the polynomials generated in the compositional inverse of O(x,t) = x / (1+(1+t)x)*(1+x)^2) = x + (-t - 3)*x^2 + (t^2 + 4 t + 6)*x^3 + (-t^3 - 5*t^2 - 10*t - 10)*x^4 + ..., containing the truncated Pascal polynomials of A104712 / A325000.
For example, from the LIF, the third derivative of ((1 + (1+t)*x)*(1+x)^2)^4 / 4! at x = 0 is 55 + 55*t + 15*t^2 + t^3.
A natural refinement of this array was provided in a letter by Isaac Newton in 1676--a set of partition polynomials for generating the o.g.f. of the compositional inverse of the generic odd o.g.f. x + u_1 x^3 + u_2 x^5 + ... in the infinite set of indeterminates u_n. (End)
T(n,k) is the number of noncrossing cacti with n+1 nodes and n+1-k blocks. See A361242. - Andrew Howroyd, Apr 13 2023

			Triangle begins:
     1;
     3,    1;
    12,    8,    1;
    55,   55,   15,   1;
   273,  364,  156,  24,  1;
  1428, 2380, 1400, 350, 35, 1;
  ...

Michael De Vlieger, Table of n, a(n) for n = 1..11325 (rows 1 <= n <= 150, flattened)
Paul Barry, On the inversion of Riordan arrays, arXiv:2101.06713 [math.CO], 2021.
J.-C. Novelli and J.-Y. Thibon, Hopf Algebras of m-permutations,(m+1)-ary trees, and m-parking functions, arXiv preprint arXiv:1403.5962 [math.CO], 2014. See Fig. 10.

Cf. A001764, A001263, A243663 (m=3).
Row sums give A003168.
Row reversed triangle is A102537.

T[m_][n_, k_] := Binomial[(m + 1) n + 1 - k, n - k] Binomial[n, k - 1]/n;
Table[T[2][n, k], {n, 1, 10}, {k, 1, n}] // Flatten (* Jean-François Alcover, Feb 12 2019 *)

T(n)=[Vecrev(p) | p<-Vec(serreverse(x/((1+x+x*y)*(1+x)^2) + O(x*x^n)))]
{ my(A=T(10)); for(i=1, #A, print(A[i])) } \\ Andrew Howroyd, Apr 13 2023

T(n,k) = (binomial(3*n+1,n) * binomial(n,k-1) * binomial(n-1,k-1)) / (binomial(3*n,k-1) * (3*n+1)) = (A001764(n) * A001263(n,k) * k) / binomial(3*n,k-1) for 1 <= k <= n (conjectured). - Werner Schulte, Nov 22 2018
T(n,k) = binomial(3*n+1-k,n-k) * binomial(n,k-1) / n for 1 <= k <= n, more generally: T_m(n,k) = binomial((m+1)*n+1-k,n-k) * binomial(n,k-1) / n for 1 <= k <= n and some fixed integer m > 1. - Werner Schulte, Nov 22 2018
G.f.: A(x,y) is the series reversion of x/((1 + x + x*y)*(1 + x)^2). - Andrew Howroyd, Apr 13 2023

Data and Example (T(2,2) and T(5,3)) corrected and more terms added by Werner Schulte, Nov 22 2018

cp's OEIS Frontend

A337885 Array read by descending antidiagonals: T(n,k) is the number of chiral pairs of colorings of the triangular faces of a regular n-dimensional simplex using k or fewer colors.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula