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

A325000 Array read by descending antidiagonals: T(n,k) is the number of unoriented colorings of the facets (or vertices) of a regular n-dimensional simplex using up to k colors.

Original entry on oeis.org

1, 3, 1, 6, 4, 1, 10, 10, 5, 1, 15, 20, 15, 6, 1, 21, 35, 35, 21, 7, 1, 28, 56, 70, 56, 28, 8, 1, 36, 84, 126, 126, 84, 36, 9, 1, 45, 120, 210, 252, 210, 120, 45, 10, 1, 55, 165, 330, 462, 462, 330, 165, 55, 11, 1, 66, 220, 495, 792, 924, 792, 495, 220, 66, 12, 1
Offset: 1

Views

Author

Robert A. Russell, Mar 23 2019

Keywords

Comments

For n=1, the figure is a line segment with two vertices. For n=2, the figure is a triangle with three edges. For n=3, the figure is a tetrahedron with four triangular faces. The Schläfli symbol, {3,...,3}, of the regular n-dimensional simplex consists of n-1 threes. Each of its n+1 facets is a regular (n-1)-dimensional simplex. Two unoriented colorings are the same if congruent; chiral pairs are counted as one.
Note that antidiagonals are part of rows of the Pascal triangle.
T(n,k-n) is the number of chiral pairs of colorings of the facets (or vertices) of a regular n-dimensional simplex using k or fewer colors. - Robert A. Russell, Sep 28 2020

Examples

			The array begins with T(1,1):
  1  3  6  10  15   21   28    36    45    55    66     78     91    105 ...
  1  4 10  20  35   56   84   120   165   220   286    364    455    560 ...
  1  5 15  35  70  126  210   330   495   715  1001   1365   1820   2380 ...
  1  6 21  56 126  252  462   792  1287  2002  3003   4368   6188   8568 ...
  1  7 28  84 210  462  924  1716  3003  5005  8008  12376  18564  27132 ...
  1  8 36 120 330  792 1716  3432  6435 11440 19448  31824  50388  77520 ...
  1  9 45 165 495 1287 3003  6435 12870 24310 43758  75582 125970 203490 ...
  1 10 55 220 715 2002 5005 11440 24310 48620 92378 167960 293930 497420 ...
  ...
For T(1,2) = 3, the two achiral colorings use just one of the two colors for both vertices; the chiral pair uses two colors. For T(2,2)=4, the triangle may have 0, 1, 2, or 3 edges of one color.
		

Crossrefs

Cf. A324999 (oriented), A325001 (achiral).
Unoriented: A007318(n,k-1) (exactly k colors), A327084 (edges, ridges), A337884 (faces, peaks), A325005 (orthotope facets, orthoplex vertices), A325013 (orthoplex facets, orthotope vertices).
Chiral: A327085 (edges, ridges), A337885 (faces, peaks), A325006 (orthotope facets, orthoplex vertices), A325014 (orthoplex facets, orthotope vertices).
Cf. A104712 (same sequence for a triangle; same sequence apart from offset).
Rows 1-4 are A000217, A000292, A000332(n+3), A000389(n+4). - Robert A. Russell, Sep 28 2020

Programs

  • Mathematica
    Table[Binomial[d+1,n+1], {d,1,15}, {n,1,d}] // Flatten

Formula

T(n,k) = binomial(n+k,n+1) = A007318(n+k,n+1).
T(n,k) = Sum_{j=1..n+1} A007318(n,j-1) * binomial(k,j).
T(n,k) = A324999(n,k) + T(n,k-n) = (A324999(n,k) - A325001(n,k)) / 2 = T(n,k-n) + A325001(n,k). - Robert A. Russell, Sep 28 2020
G.f. for row n: x / (1-x)^(n+2).
Linear recurrence for row n: T(n,k) = Sum_{j=1..n+2} -binomial(j-n-3,j) * T(n,k-j).
G.f. for column k: (1 - (1-x)^k) / (x * (1-x)^k) - k.
T(n,k-n) = A324999(n,k) - T(n,k) = (A324999(n,k) - A325001(n,k)) / 2 = T(n,k) - A325001(n,k). - Robert A. Russell, Oct 10 2020

A325012 Array read by descending antidiagonals: A(n,k) is the number of oriented colorings of the facets of a regular n-dimensional orthoplex using up to k colors.

Original entry on oeis.org

1, 4, 1, 9, 6, 1, 16, 24, 23, 1, 25, 70, 333, 496, 1, 36, 165, 2916, 230076, 2275974, 1, 49, 336, 16725, 22456756, 965227578201, 800648638402240, 1, 64, 616, 70911, 795467350, 9607713956430560, 149031415906337877339236058, 1054942853799126580390222487977120, 1
Offset: 1

Views

Author

Robert A. Russell, May 27 2019

Keywords

Comments

Also called cross polytope and hyperoctahedron. For n=1, the figure is a line segment with two vertices. For n=2 the figure is a square with four edges. For n=3 the figure is an octahedron with eight triangular faces. For n=4, the figure is a 16-cell with sixteen tetrahedral facets. The Schläfli symbol, {3,...,3,4}, of the regular n-dimensional orthoplex (n>1) consists of n-2 threes followed by a four. Each of its 2^n facets is an (n-1)-dimensional simplex. Two oriented colorings are the same if one is a rotation of the other; chiral pairs are counted as two.
Also the number of oriented colorings of the vertices of a regular n-dimensional orthotope (cube) using up to k colors.

Examples

			Array begins with A(1,1):
1   4      9       16        25          36           49            64 ...
1   6     24       70       165         336          616          1044 ...
1  23    333     2916     16725       70911       241913        701968 ...
1 496 230076 22456756 795467350 14697611496 173107727191 1466088119056 ...
For A(1,2) = 4, the two achiral colorings use just one of the two colors for both vertices; the chiral pair uses one color for each vertex.
		

Crossrefs

Cf. A325013 (unoriented), A325014 (chiral), A325015 (achiral), A325016 (exactly k colors).
Other n-dimensional polytopes: A324999 (simplex), A325004 (orthotope).
Rows 1-3 are A000290, A006528, A000543; column 2 is A237748.

Programs

  • Mathematica
    a48[n_] := a48[n] = DivisorSum[NestWhile[#/2&,n,EvenQ], MoebiusMu[#]2^(n/#)&]/(2n); (* A000048 *)
    a37[n_] := a37[n] = DivisorSum[n,MoebiusMu[n/#]2^#&]/n; (* A001037 *)
    CI0[{n_Integer}] := CI0[{n}] = CI[Transpose[If[EvenQ[n], p2 = IntegerExponent[n, 2]; sub = Divisors[n/2^p2]; {2^(p2+1) sub, a48 /@ (2^p2 sub) }, sub = Divisors[n]; {sub, a37 /@ sub}]]] 2^(n-1);(* even perm. *)
    CI1[{n_Integer}] := CI1[{n}] = CI[sub = Divisors[n]; Transpose[If[EvenQ[n], {sub, a37 /@ sub}, {2 sub, a48 /@ sub}]]] 2^(n-1); (* odd perm. *)
    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)
    cix[{a_, b_}, {c_, d_}] := {LCM[a, c], (a b c d)/LCM[a, c]};
    Unprotect[Times]; Times[CI[a_List], CI[b_List]] :=  (* combine *) CI[compress[Flatten[Outer[cix, a, b, 1], 1]]]; Protect[Times];
    CI0[p_List] := CI0[p] = Expand[CI0[Drop[p, -1]] CI0[{Last[p]}] + CI1[Drop[p, -1]] CI1[{Last[p]}]]
    CI1[p_List] := CI1[p] = Expand[CI0[Drop[p, -1]] CI1[{Last[p]}] + CI1[Drop[p, -1]] CI0[{Last[p]}]]
    pc[p_List] := Module[{ci,mb},mb = DeleteDuplicates[p]; ci = Count[p, #] & /@ mb; n!/(Times @@ (ci!) Times @@ (mb^ci))] (* partition count *)
    row[n_Integer] := row[n] = Factor[(Total[(CI0[#] pc[#]) & /@ IntegerPartitions[n]])/(n! 2^(n - 1))] /. CI[l_List] :> j^(Total[l][[2]])
    array[n_, k_] := row[n] /. j -> k
    Table[array[n, d-n+1], {d, 1, 10}, {n, 1, d}] // Flatten

Formula

The algorithm used in the Mathematica program below assigns each permutation of the axes to a partition of n. It then determines the number of permutations for each partition and the cycle index for each partition.
A(n,k) = A325013(n,k) + A325014(n,k) = 2*A325013(n,k) - A325015(n,k) = 2*A325014(n,k) + A325015(n,k).
A(n,k) = Sum_{j=1..2^n} A325016(n,j) * binomial(k,j).

A325006 Array read by descending antidiagonals: A(n,k) is the number of chiral pairs of colorings of the facets of a regular n-dimensional orthotope using up to k colors.

Original entry on oeis.org

0, 1, 0, 3, 0, 0, 6, 3, 0, 0, 10, 15, 1, 0, 0, 15, 45, 20, 0, 0, 0, 21, 105, 120, 15, 0, 0, 0, 28, 210, 455, 210, 6, 0, 0, 0, 36, 378, 1330, 1365, 252, 1, 0, 0, 0, 45, 630, 3276, 5985, 3003, 210, 0, 0, 0, 0, 55, 990, 7140, 20475, 20349, 5005, 120, 0, 0, 0, 0, 66, 1485, 14190, 58905, 98280, 54264, 6435, 45, 0, 0, 0, 0
Offset: 1

Views

Author

Robert A. Russell, May 27 2019

Keywords

Comments

Also called hypercube, n-dimensional cube, and measure polytope. For n=1, the figure is a line segment with two vertices. For n=2 the figure is a square with four edges. For n=3 the figure is a cube with six square faces. For n=4, the figure is a tesseract with eight cubic facets. The Schläfli symbol, {4,3,...,3}, of the regular n-dimensional orthotope (n>1) consists of a four followed by n-2 threes. Each of its 2n facets is an (n-1)-dimensional orthotope. The chiral colorings of its facets come in pairs, each the reflection of the other.
Also the number of chiral pairs of colorings of the vertices of a regular n-dimensional orthoplex using up to k colors.

Examples

			Array begins with A(1,1):
0 1 3  6  10   15     21       28        36         45          55 ...
0 0 3 15  45  105    210      378       630        990        1485 ...
0 0 1 20 120  455   1330     3276      7140      14190       26235 ...
0 0 0 15 210 1365   5985    20475     58905     148995      341055 ...
0 0 0  6 252 3003  20349    98280    376992    1221759     3478761 ...
0 0 0  1 210 5005  54264   376740   1947792    8145060    28989675 ...
0 0 0  0 120 6435 116280  1184040   8347680   45379620   202927725 ...
0 0 0  0  45 6435 203490  3108105  30260340  215553195  1217566350 ...
0 0 0  0  10 5005 293930  6906900  94143280  886163135  6358402050 ...
0 0 0  0   1 3003 352716 13123110 254186856 3190187286 29248649430 ...
For a(2,3)=3, each chiral pair consists of two adjacent edges of the square with one of the three colors.
		

Crossrefs

Cf. A325004 (oriented), A325005 (unoriented), A325007 (achiral), A325010 (exactly k colors)
Other n-dimensional polytopes: A007318(k,n+1) (simplex), A325014 (orthoplex)
Rows 1-3 are A161680, A050534, A093566(n+1), A234249(n-1)

Programs

  • Mathematica
    Table[Binomial[Binomial[d-n+1,2],n],{d,1,12},{n,1,d}] // Flatten
  • PARI
    a(n, k) = binomial(binomial(k, 2), n)
    array(rows, cols) = for(x=1, rows, for(y=1, cols, print1(a(x, y), ", ")); print(""))
    /* Print initial 10 rows and 11 columns of array as follows: */
    array(10, 11) \\ Felix Fröhlich, May 30 2019

Formula

A(n,k) = binomial(binomial(k,2),n).
A(n,k) = Sum_{j=1..2*n} A325010(n,j) * binomial(k,j).
A(n,k) = A325004(n,k) - A325005(n,k) = (A325004(n,k) - A325007(n,k)) / 2 = A325005(n,k) - A325007(n,k).
G.f. for row n: Sum{j=1..2*n} A325010(n,j) * x^j / (1-x)^(j+1).
Linear recurrence for row n: T(n,k) = Sum_{j=0..2*n} binomial(-2-j,2*n-j) * T(n,k-1-j).
G.f. for column k: (1+x)^binomial(k,2) - 1.

A325015 Array read by descending antidiagonals: T(n,k) is the number of achiral colorings of the facets of a regular n-dimensional orthoplex using up to k colors.

Original entry on oeis.org

1, 2, 1, 3, 6, 1, 4, 18, 21, 1, 5, 40, 201, 308, 1, 6, 75, 1076, 34128, 180342, 1, 7, 126, 4025, 1056576, 2945136213, 366975285216, 1, 8, 196, 11901, 15303750, 2932338749408, 103863386269870076808, 10316179427644325573474464, 1
Offset: 1

Views

Author

Robert A. Russell, May 27 2019

Keywords

Comments

Also called cross polytope and hyperoctahedron. For n=1, the figure is a line segment with two vertices. For n=2 the figure is a square with four edges. For n=3 the figure is an octahedron with eight triangular faces. For n=4, the figure is a 16-cell with sixteen tetrahedral facets. The Schläfli symbol, {3,...,3,4}, of the regular n-dimensional orthoplex (n>1) consists of n-2 threes followed by a four. Each of its 2^n facets is an (n-1)-dimensional simplex. An achiral coloring is identical to its reflection.
Also the number of achiral colorings of the vertices of a regular n-dimensional orthotope (cube) using up to k colors.

Examples

			Array begins with T(1,1):
1   2     3       4        5         6         7          8 ...
1   6    18      40       75       126       196        288 ...
1  21   201    1076     4025     11901     29841      66256 ...
1 308 34128 1056576 15303750 136236276 865711763 4296782848 ...
...
For T(2,2)=6, two squares have all edges the same color, two have three edges the same color, one has opposite edges the same color, and one has opposite edges different colors.
		

Crossrefs

Cf. A325012 (oriented), A325013 (unoriented), A325014 (chiral), A325019 (exactly k colors).
Other n-dimensional polytopes: A325001 (simplex), A325007 (orthotope).
Rows 1-2 are A000027, A002411.

Programs

  • Mathematica
    a48[n_] := a48[n] = DivisorSum[NestWhile[#/2&, n, EvenQ], MoebiusMu[#]2^(n/#)&]/(2n); (* A000048 *)
    a37[n_] := a37[n] = DivisorSum[n, MoebiusMu[n/#]2^#&]/n; (* A001037 *)
    CI0[{n_Integer}] := CI0[{n}] = CI[Transpose[If[EvenQ[n], p2 = IntegerExponent[n, 2]; sub = Divisors[n/2^p2]; {2^(p2+1) sub, a48 /@ (2^p2 sub) }, sub = Divisors[n]; {sub, a37 /@ sub}]]] 2^(n-1); (* even perm. *)
    CI1[{n_Integer}] := CI1[{n}] = CI[sub = Divisors[n]; Transpose[If[EvenQ[n], {sub, a37 /@ sub}, {2 sub, a48 /@ sub}]]] 2^(n-1); (* odd perm. *)
    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)
    cix[{a_, b_}, {c_, d_}] := {LCM[a, c], (a b c d)/LCM[a, c]};
    Unprotect[Times]; Times[CI[a_List], CI[b_List]] :=  (* combine *) CI[compress[Flatten[Outer[cix, a, b, 1], 1]]]; Protect[Times];
    CI0[p_List] := CI0[p] = Expand[CI0[Drop[p, -1]] CI0[{Last[p]}] + CI1[Drop[p, -1]] CI1[{Last[p]}]]
    CI1[p_List] := CI1[p] = Expand[CI0[Drop[p, -1]] CI1[{Last[p]}] + CI1[Drop[p, -1]] CI0[{Last[p]}]]
    pc[p_List] := Module[{ci,mb},mb = DeleteDuplicates[p]; ci = Count[p, #] & /@ mb; n!/(Times @@ (ci!) Times @@ (mb^ci))] (* partition count *)
    row[n_Integer] := row[n] = Factor[(Total[(CI1[#] pc[#]) & /@ IntegerPartitions[n]])/(n! 2^(n - 1))] /. CI[l_List] :> j^(Total[l][[2]])
    array[n_, k_] := row[n] /. j -> k
    Table[array[n, d-n+1], {d, 1, 10}, {n, 1, d}] // Flatten

Formula

The algorithm used in the Mathematica program below assigns each permutation of the axes to a partition of n. It then determines the number of permutations for each partition and the cycle index for each partition.
T(n,k) = 2*A325013(n,k) - A325012(n,k) = A325012(n,k) - 2*A325014(n,k) = A325013(n,k) - A325014(n,k).
T(n,k) = Sum_{j=1..3*2^(n-2)} A325019(n,j) * binomial(k,j).

A325013 Array read by descending antidiagonals: A(n,k) is the number of unoriented colorings of the facets of a regular n-dimensional orthoplex using up to k colors.

Original entry on oeis.org

1, 3, 1, 6, 6, 1, 10, 21, 22, 1, 15, 55, 267, 402, 1, 21, 120, 1996, 132102, 1228158, 1, 28, 231, 10375, 11756666, 484086357207, 400507806843728, 1, 36, 406, 41406, 405385550, 4805323147589984, 74515759884862073604656433, 527471432057653004017274030725792, 1
Offset: 1

Views

Author

Robert A. Russell, May 27 2019

Keywords

Comments

Also called cross polytope and hyperoctahedron. For n=1, the figure is a line segment with two vertices. For n=2 the figure is a square with four edges. For n=3 the figure is an octahedron with eight triangular faces. For n=4, the figure is a 16-cell with sixteen tetrahedral facets. The Schläfli symbol, {3,...,3,4}, of the regular n-dimensional orthoplex (n>1) consists of n-2 threes followed by a four. Each of its 2^n facets is an (n-1)-dimensional simplex. Two unoriented colorings are the same if congruent; chiral pairs are counted as one.
Also the number of unoriented colorings of the vertices of a regular n-dimensional orthotope (cube) using up to k colors.

Examples

			Array begins with A(1,1):
1   3      6       10        15         21          28           36 ...
1   6     21       55       120        231         406          666 ...
1  22    267     1996     10375      41406      135877       384112 ...
1 402 132102 11756666 405385550 7416923886 86986719477 735192450952 ...
For A(2,2)=6, two squares have all edges the same color, two have three edges the same color, one has opposite edges the same color, and one has opposite edges different colors.
		

Crossrefs

Cf. A325012 (oriented), A325014 (chiral), A325015 (achiral), A325017 (exactly k colors).
Other n-dimensional polytopes: A325000 (simplex), A325005 (orthotope).
Rows 1-4 are A000217, A002817, A128766, A128767; column 2 is A000616.

Programs

  • Mathematica
    a48[n_] := a48[n] = DivisorSum[NestWhile[#/2&, n, EvenQ], MoebiusMu[#]2^(n/#)&]/(2n); (* A000048 *)
    a37[n_] := a37[n] = DivisorSum[n, MoebiusMu[n/#]2^#&]/n; (* A001037 *)
    CI0[{n_Integer}] := CI0[{n}] = CI[Transpose[If[EvenQ[n], p2 = IntegerExponent[n, 2]; sub = Divisors[n/2^p2]; {2^(p2+1) sub, a48 /@ (2^p2 sub) }, sub = Divisors[n]; {sub, a37 /@ sub}]]] 2^(n-1); (* even perm. *)
    CI1[{n_Integer}] := CI1[{n}] = CI[sub = Divisors[n]; Transpose[If[EvenQ[n], {sub, a37 /@ sub}, {2 sub, a48 /@ sub}]]] 2^(n-1); (* odd perm. *)
    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)
    cix[{a_, b_}, {c_, d_}] := {LCM[a, c], (a b c d)/LCM[a, c]};
    Unprotect[Times]; Times[CI[a_List], CI[b_List]] :=  (* combine *) CI[compress[Flatten[Outer[cix, a, b, 1], 1]]]; Protect[Times];
    CI0[p_List] := CI0[p] = Expand[CI0[Drop[p, -1]] CI0[{Last[p]}] + CI1[Drop[p, -1]] CI1[{Last[p]}]]
    CI1[p_List] := CI1[p] = Expand[CI0[Drop[p, -1]] CI1[{Last[p]}] + CI1[Drop[p, -1]] CI0[{Last[p]}]]
    pc[p_List] := Module[{ci,mb},mb = DeleteDuplicates[p]; ci = Count[p, #] & /@ mb; n!/(Times @@ (ci!) Times @@ (mb^ci))] (* partition count *)
    row[n_Integer] := row[n] = Factor[(Total[((CI0[#] + CI1[#]) pc[#]) & /@ IntegerPartitions[n]])/(n! 2^n)] /. CI[l_List] :> j^(Total[l][[2]])
    array[n_, k_] := row[n] /. j -> k
    Table[array[n, d-n+1], {d, 1, 10}, {n, 1, d}] // Flatten

Formula

The algorithm used in the Mathematica program below assigns each permutation of the axes to a partition of n. It then determines the number of permutations for each partition and the cycle index for each partition.
A(n,k) = A325012(n,k) - A325014(n,k) = (A325012(n,k) + A325015(n,k)) / 2 = A325014(n,k) + A325015(n,k).
A(n,k) = Sum_{j=1..2^n} A325017(n,j) * binomial(k,j).

A337409 Array read by descending antidiagonals: T(n,k) is the number of chiral pairs of colorings of the edges of a regular n-dimensional orthotope (hypercube) using k or fewer colors.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 3, 74, 0, 0, 15, 10704, 11158298, 0, 0, 45, 345640, 4825452718593, 314824408633217132928, 0, 0, 105, 5062600, 48038354542204960, 38491882659952177472606694634030116, 136221825854745676076981182469325427379054390050209792, 0
Offset: 1

Views

Author

Robert A. Russell, Aug 26 2020

Keywords

Comments

Each member of a chiral pair is a reflection, but not a rotation, of the other. For n=1, the figure is a line segment with one edge. For n=2, the figure is a square with 4 edges. For n=3, the figure is a cube with 12 edges. The number of edges is n*2^(n-1).
Also the number of chiral pairs of colorings of the regular (n-2)-dimensional simplexes in a regular n-dimensional orthoplex.

Examples

			Table begins with T(1,1):
0  0     0      0       0        0         0          0          0 ...
0  0     3     15      45      105       210        378        630 ...
0 74 10704 345640 5062600 45246810 288005144 1430618784 5881281480 ...
For T(2,3)=3, the chiral arrangements are AABC-AACB, ABBC-ACBB, and ABCC-ACCB.
		

Crossrefs

Cf. A337407 (oriented), A337408 (unoriented), A337410 (achiral).
Rows 2-4 are A050534, A337406, A331360.
Cf. A327085 (simplex edges), A337413 (orthoplex edges), A325014 (orthotope vertices).

Programs

  • Mathematica
    m=1; (* dimension of color element, here an edge *)
    Fi1[p1_] := Module[{g, h}, Coefficient[Product[g = GCD[k1, p1]; h = GCD[2 k1, p1]; (1+2x^(k1/g))^(r1[[k1]] g) If[Divisible[k1, h], 1, (1+2x^(2 k1/h))^(r2[[k1]] h/2)], {k1, Flatten[Position[cs, n1_ /; n1 > 0]]}], x, n-m]];
    FiSum[] := (Do[Fi2[k2] = Fi1[k2], {k2, Divisors[per]}]; DivisorSum[per, DivisorSum[d1 = #, MoebiusMu[d1/#] Fi2[#] &]/# &]);
    CCPol[r_List] := (r1 = r; r2 = cs - r1; per = LCM @@ Table[If[cs[[j2]] == r1[[j2]], If[0 == cs[[j2]],1,j2], 2j2], {j2,n}]; If[EvenQ[Sum[If[EvenQ[j3], r1[[j3]], r2[[j3]]], {j3,n}]],1,-1]Times @@ Binomial[cs, r1] 2^(n-Total[cs]) b^FiSum[]);
    PartPol[p_List] := (cs = Count[p, #]&/@ Range[n]; Total[CCPol[#]&/@ Tuples[Range[0,cs]]]);
    pc[p_List] := Module[{ci, mb}, mb = DeleteDuplicates[p]; ci = Count[p, #]&/@ mb; n!/(Times@@(ci!) Times@@(mb^ci))] (*partition count*)
    row[n_Integer] := row[n] = Factor[(Total[(PartPol[#] pc[#])&/@ IntegerPartitions[n]])/(n! 2^n)]
    array[n_, k_] := row[n] /. b -> k
    Table[array[n,d+m-n], {d,7}, {n,m,d+m-1}] // Flatten

Formula

The algorithm used in the Mathematica program below assigns each permutation of the axes to a partition of n and then considers separate conjugacy classes for axis reversals. It uses the formulas in Balasubramanian's paper. If the value of m is increased, one can enumerate colorings of higher-dimensional elements beginning with T(m,1).
T(n,k) = A337407(n,k) - A337408(n,k) = (A337407(n,k) - A337410(n,k)) / 2 = A337408(n,k) - A337410(n,k).

A325018 Triangle read by rows: T(n,k) is the number of chiral pairs of colorings of the facets of a regular n-dimensional orthoplex using exactly k colors. Row n has 2^n columns.

Original entry on oeis.org

0, 1, 0, 0, 3, 3, 0, 1, 63, 662, 2400, 3900, 2940, 840, 0, 94, 97692, 10308758, 337560150, 5098740090, 42976836210, 224685801060, 775389028050, 1830791421900, 3007909258200, 3439214024400, 2685727044000, 1366701336000, 408648240000, 54486432000
Offset: 1

Views

Author

Robert A. Russell, Jun 09 2019

Keywords

Comments

Also called cross polytope and hyperoctahedron. For n=1, the figure is a line segment with two vertices. For n=2 the figure is a square with four edges. For n=3 the figure is an octahedron with eight triangular faces. For n=4, the figure is a 16-cell with sixteen tetrahedral facets. The Schläfli symbol, {3,...,3,4}, of the regular n-dimensional orthoplex (n>1) consists of n-2 threes followed by a four. Each of its 2^n facets is an (n-1)-dimensional simplex. The chiral colorings of its facets come in pairs, each the reflection of the other.
Also the number of chiral pairs of colorings of the vertices of a regular n-dimensional orthotope (cube) using exactly k colors.

Examples

			Triangle begins with T(1,1):
0 1
0 0  3   3
0 1 63 662 2400 3900 2940 840
For T(2,3)=3, each square has one of the three colors on two adjacent edges.
		

Crossrefs

Cf. A325016 (oriented), A325017 (unoriented), A325019 (achiral), A325014 (up to k colors).
Other n-dimensional polytopes: A325010 (orthotope).

Programs

  • Mathematica
    a48[n_] := a48[n] = DivisorSum[NestWhile[#/2&,n,EvenQ],MoebiusMu[#]2^(n/#)&]/(2n); (* A000048 *)
    a37[n_] := a37[n] = DivisorSum[n, MoebiusMu[n/#]2^#&]/n; (* A001037 *)
    CI0[{n_Integer}] := CI0[{n}] = CI[Transpose[If[EvenQ[n], p2 = IntegerExponent[n, 2]; sub = Divisors[n/2^p2]; {2^(p2+1) sub, a48 /@ (2^p2 sub) }, sub = Divisors[n]; {sub, a37 /@ sub}]]] 2^(n-1); (* even perm. *)
    CI1[{n_Integer}] := CI1[{n}] = CI[sub = Divisors[n]; Transpose[If[EvenQ[n], {sub, a37 /@ sub}, {2 sub, (a37 /@ sub)/2}]]] 2^(n-1); (* odd perm. *)
    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)
    cix[{a_, b_}, {c_, d_}] := {LCM[a, c], (a b c d)/LCM[a, c]};
    Unprotect[Times]; Times[CI[a_List], CI[b_List]] :=  (* combine *) CI[compress[Flatten[Outer[cix, a, b, 1], 1]]]; Protect[Times];
    CI0[p_List] := CI0[p] = Expand[CI0[Drop[p, -1]] CI0[{Last[p]}] + CI1[Drop[p, -1]] CI1[{Last[p]}]]
    CI1[p_List] := CI1[p] = Expand[CI0[Drop[p, -1]] CI1[{Last[p]}] + CI1[Drop[p, -1]] CI0[{Last[p]}]]
    pc[p_List] := Module[{ci,mb},mb = DeleteDuplicates[p]; ci = Count[p, #] & /@ mb; n!/(Times @@ (ci!) Times @@ (mb^ci))] (* partition count *)
    row[n_Integer] := row[n] = Factor[(Total[((CI0[#] - CI1[#]) pc[#]) & /@ IntegerPartitions[n]])/(n! 2^n)] /. CI[l_List] :> j^(Total[l][[2]])
    array[n_, k_] := row[n] /. j -> k (* A325014 *)
    Table[LinearSolve[Table[Binomial[i,j],{i,1,2^n},{j,1,2^n}],Table[array[n,k],{k,1,2^n}]],{n,1,6}] // Flatten

Formula

A325014(n,k) = Sum_{j=1..2^n} T(n,j) * binomial(k,j).
T(n,k) = A325016(n,k) - A325017(n,k) = (A325016(n,k) - A325019(n,k)) / 2 = A325017(n,k) - A325019(n,k).

A337954 Number of chiral pairs of colorings of the 16 tetrahedral facets of a hyperoctahedron or of the 16 vertices of a tesseract.

Original entry on oeis.org

0, 94, 97974, 10700090, 390081800, 7280687610, 86121007714, 730895668104, 4816861200630, 26010740238450, 119563513291420, 481192778757834, 1732132086737234, 5669991002636870, 17101193825828700, 48029634770843680
Offset: 1

Views

Author

Robert A. Russell, Oct 03 2020

Keywords

Comments

Each member of a chiral pair is a reflection, but not a rotation, of the other. The Schläfli symbols for the tesseract and the hyperoctahedron are {4,3,3} and {3,3,4} respectively. Both figures are regular 4-D polyhedra and they are mutually dual.

Crossrefs

Cf. A337952 (oriented), A128767 (unoriented), A337955 (achiral).
Other elements: A331360 (tesseract edges, hyperoctahedron faces), A331356 (tesseract faces, hyperoctahedron edges), A234249(n+1) (tesseract facets, hyperoctahedron vertices).
Other polychora: A000389 (4-simplex facets/vertices), A338950 (24-cell), A338966 (120-cell, 600-cell).
Row 4 of A325014 (orthoplex facets, orthotope vertices).

Programs

  • Mathematica
    Table[(n^16-12n^12+12n^10+43n^8-48n^6-44n^4+48n^2)/384,{n, 30}]

Formula

a(n) = (n-1) * n^2 * (n+1) * (n^12 + n^10 - 11*n^8 + n^6 + 44 n^4 - 4 n^2 - 48) / 384.
a(n) = 94*C(n,2) + 97692*C(n,3) + 10308758*C(n,4) + 337560150*C(n,5) + 5098740090*C(n,6) + 42976836210*C(n,7) + 224685801060*C(n,8) + 775389028050*C(n,9) + 1830791421900*C(n,10) + 3007909258200*C(n,11) + 3439214024400*C(n,12) + 2685727044000*C(n,13) + 1366701336000*C(n,14) + 408648240000*C(n,15) + 54486432000*C(n,16), where the coefficient of C(n,k) is the number of chiral pairs of colorings using exactly k colors.
a(n) = A337952(n) - A128767(n) = (A337952(n) - A337955(n)) / 2 = A128767(n) - A337955(n).

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

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 40927, 0, 0, 20, 731279799, 314824333015938998688, 0, 0, 120, 732272925320, 38491882659300767730994725249684096, 38343035259947576596859560773963975000551460473665493534170658111488, 0
Offset: 2

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

Each member of a chiral pair is a reflection, but not a rotation, of the other. Each face is a square bounded by four edges. For n=2, the figure is a square with one face. For n=3, the figure is a cube with 6 faces. For n=4, the figure is a tesseract with 24 faces. The number of faces is 2^(n-2)*C(n,2).
Also the number of chiral pairs of colorings of peaks of an n-dimensional orthoplex. A peak is an (n-3)-dimensional simplex.
The algorithm used in the Mathematica program below assigns each permutation of the axes to a partition of n and then considers separate conjugacy classes for axis reversals. It uses the formulas in Balasubramanian's paper. If the value of m is increased, one can enumerate colorings of higher-dimensional elements beginning with T(m,1).

Examples

			Array begins with T(2,1):
0     0         0            0               0                 0 ...
0     0         1           20             120               455 ...
0 40927 731279799 732272925320 155180061396500 12338466190481025 ...
		

Crossrefs

Cf. A337887 (oriented), A337888 (unoriented), A337890 (achiral).
Other elements: A325014 (vertices), A337409 (edges).
Other polytopes: A337885 (simplex), A337893 (orthoplex).
Rows 2-4 are A000004, A093566(n+1), A331356.

Programs

  • Mathematica
    m=2; (* dimension of color element, here a square face *)
    Fi1[p1_] := Module[{g, h}, Coefficient[Product[g = GCD[k1, p1]; h = GCD[2 k1, p1]; (1+2x^(k1/g))^(r1[[k1]] g) If[Divisible[k1, h], 1, (1+2x^(2 k1/h))^(r2[[k1]] h/2)], {k1, Flatten[Position[cs, n1_ /; n1 > 0]]}], x, n-m]];
    FiSum[] := (Do[Fi2[k2] = Fi1[k2], {k2, Divisors[per]}]; DivisorSum[per, DivisorSum[d1 = #, MoebiusMu[d1/#] Fi2[#] &]/# &]);
    CCPol[r_List] := (r1 = r; r2 = cs - r1; per = LCM @@ Table[If[cs[[j2]] == r1[[j2]], If[0 == cs[[j2]],1,j2], 2j2], {j2,n}]; If[EvenQ[Sum[If[EvenQ[j3], r1[[j3]], r2[[j3]]], {j3,n}]],1,-1]Times @@ Binomial[cs, r1] 2^(n-Total[cs]) b^FiSum[]);
    PartPol[p_List] := (cs = Count[p, #]&/@ Range[n]; Total[CCPol[#]&/@ Tuples[Range[0,cs]]]);
    pc[p_List] := Module[{ci, mb}, mb = DeleteDuplicates[p]; ci = Count[p, #]&/@ mb; n!/(Times@@(ci!) Times@@(mb^ci))] (*partition count*)
    row[n_Integer] := row[n] = Factor[(Total[(PartPol[#] pc[#])&/@ IntegerPartitions[n]])/(n! 2^n)]
    array[n_, k_] := row[n] /. b -> k
    Table[array[n,d+m-n], {d,6}, {n,m,d+m-1}] // Flatten

Formula

T(n,k) = A337887(n,k) - A337888(n,k) = (A337887(n,k) - A337890(n,k)) / 2 = A337888(n,k) - A337890(n,k).

A337896 Number of chiral pairs of colorings of the 8 triangular faces of a regular octahedron or the 8 vertices of a cube using n or fewer colors.

Original entry on oeis.org

0, 1, 66, 920, 6350, 29505, 106036, 317856, 832140, 1961025, 4248310, 8590296, 16398746, 29814785, 51983400, 87399040, 142333656, 225359361, 347978730, 525376600, 777308070, 1129138241, 1613050076, 2269437600
Offset: 1

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

Each member of a chiral pair is a reflection, but not a rotation, of the other.

Examples

			For a(2)=1, centering the octahedron (cube) at the origin and aligning the diagonals (edges) with the axes, color the faces (vertices) in the octants ---, --+, -++, and +++ with one color and the other 4 elements with the other color.
		

Crossrefs

Cf. A000543 (oriented), A128766(unoriented), A337897 (achiral).
Other elements: A337406 (edges), A093566(n+1) (cube faces, octahedron vertices).
Other polyhedra: A000332 (simplex), A093566(n+1) (cube/octahedron).
Row 3 of A325014 (chiral pairs of colorings of orthoplex facets or orthotope vertices).
Row 3 of A337893 (chiral pairs of colorings of orthoplex faces or orthotope peaks).

Programs

  • Mathematica
    Table[(n-1)n^2(n+1)(8-5n^2+n^4)/48, {n,30}]

Formula

a(n) = (n-1) * n^2 * (n+1) * (8 - 5*n^2 + n^4) / 48.
a(n) = 1*C(n,2) + 63*C(n,3) + 662*C(n,4) + 2400*C(n,5) + 3900*C(n,6) + 2940*C(n,7) + 840*C(n,8), where the coefficient of C(n,k) is the number of chiral pairs of colorings using exactly k colors.
G.f.: x^2 * (1+x) * (1+56*x+306*x^2+56*x^3+x^4) / (1-x)^9.
a(n) = A000543(n) - A128766(n) = (A000543(n) - A337897(n)) / 2 = A128766(n) - A337897(n).
Showing 1-10 of 10 results.