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

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

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 4, 21, 6, 0, 0, 10, 140, 405, 28, 0, 0, 20, 575, 7904, 17154, 252, 0, 0, 35, 1785, 76880, 1415648, 1920375, 4726, 0, 0, 56, 4606, 486522, 41453650, 855834880, 547375212, 150324, 0
Offset: 1

Views

Author

Robert A. Russell, Aug 19 2019

Keywords

Comments

An n-dimensional simplex has n+1 vertices and (n+1)*n/2 edges. For n=1, the figure is a line segment with one edge. For n-2, the figure is a triangle with three edges. For n=3, the figure is a tetrahedron with six edges. The Schläfli symbol, {3,...,3}, of the regular n-dimensional simplex consists of n-1 threes. The chiral colorings of its edges come in pairs, each the reflection of the other.
A(n,k) is also the number of chiral pairs of colorings of (n-2)-dimensional regular simplices in an n-dimensional simplex using up to k colors. Thus, A(2,k) is also the number of chiral pairs of colorings of the vertices (0-dimensional simplices) of an equilateral triangle.

Examples

			Array begins with A(1,1):
  0 0   0    0     0      0       0       0        0        0         0 ...
  0 0   1    4    10     20      35      56       84      120       165 ...
  0 1  21  140   575   1785    4606   10416    21330    40425     71995 ...
  0 6 405 7904 76880 486522 2300305 8806336 28725192 82626270 214744629 ...
  ...
For A(2,3) = 1, the chiral pair is ABC-ACB.
		

Crossrefs

Cf. A327083 (oriented), A327084 (unoriented), A327086 (achiral), A327089 (exactly k colors), A325000(n,k-n) (vertices, facets), A337885 (faces, peaks), A337409 (orthotope edges, orthoplex ridges), A337413 (orthoplex edges, orthotope ridges).
Rows 1-4 are A000004, A000292(n-2), A337899, A331352.

Programs

  • Mathematica
    CycleX[{2}] = {{1,1}}; (* cycle index for permutation with given cycle structure *)
    CycleX[{n_Integer}] := CycleX[n] = If[EvenQ[n], {{n/2,1}, {n,(n-2)/2}}, {{n,(n-1)/2}}]
    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)
    CycleX[p_List] := CycleX[p] = compress[Join[CycleX[Drop[p, -1]], If[Last[p] > 1, CycleX[{Last[p]}], ## &[]], If[# == Last[p], {#, Last[p]}, {LCM[#, Last[p]], GCD[#, Last[p]]}] & /@ Drop[p, -1]]]
    pc[p_List] := 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[CycleX[#]][[2]] & /@ IntegerPartitions[n+1]]/(n+1)!]
    array[n_, k_] := row[n] /. j -> k
    Table[array[n,d-n+1], {d,1,10}, {n,1,d}] // Flatten
    (* Using Fripertinger's exponent per Andrew Howroyd's code in A063841: *)
    pc[p_] := Module[{ci, mb}, mb = DeleteDuplicates[p]; ci = Count[p, #] &/@ mb; Total[p]!/(Times @@ (ci!) Times @@ (mb^ci))]
    ex[v_] := Sum[GCD[v[[i]], v[[j]]], {i,2,Length[v]}, {j,i-1}] + Total[Quotient[v,2]]
    array[n_,k_] := Total[If[EvenQ[Total[1-Mod[#,2]]],1,-1] pc[#]k^ex[#] &/@ IntegerPartitions[n+1]]/(n+1)!
    Table[array[n,d-n+1], {d,10}, {n,d}] // Flatten

Formula

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.
A(n,k) = Sum_{j=1..(n+1)*n/2} A327089(n,j) * binomial(k,j).
A(n,k) = A327083(n,k) - A327084(n,k) = (A327083(n,k) - A327086(n,k)) / 2 = A327084(n,k) - A327086(n,k).

A331356 Number of chiral pairs of colorings of the edges of a regular 4-dimensional orthoplex with n available colors.

Original entry on oeis.org

0, 40927, 731279799, 732272925320, 155180061396500, 12338466190481025, 498892380429882397, 12297640855782563904, 207723543409061974215, 2604156223742219218875, 25650287482426463967550, 207022761847763612943192
Offset: 1

Views

Author

Robert A. Russell, Jan 14 2020

Keywords

Comments

A regular 4-dimensional orthoplex (also hyperoctahedron or cross polytope) has 8 vertices and 24 edges. Its Schläfli symbol is {3,3,4}. The chiral colorings of its edges come in pairs, each the reflection of the other. Also the number of chiral pairs of colorings of the square faces of a tesseract {4,3,3} with n available colors.

Crossrefs

Cf. A331354 (oriented), A331355 (unoriented), A331357 (achiral).
Other polychora: A331352 (5-cell), A331360 (8-cell), A338954 (24-cell), A338966 (120-cell, 600-cell).
Row 4 of A337413 (orthoplex edges, orthotope ridges) and A337889 (orthotope faces, orthoplex peaks).

Programs

  • Mathematica
    Table[(48n^3 - 20n^6 - 60n^7 + 8n^8 + 12n^9 - 3n^12 + 12n^13 + 18n^14 - 12n^15 - 4n^18 + n^24)/384, {n, 1, 25}]

Formula

a(n) = (48*n^3 - 20*n^6 - 60*n^7 + 8*n^8 + 12*n^9 - 3*n^12 + 12*n^13 + 18*n^14 - 12*n^15 - 4*n^18 + n^24) / 384.
a(n) = 40927*C(n,2) + 731157018*C(n,3) + 729348051686*C(n,4) + 151526009158620*C(n,5) + 11418355290999750*C(n,6) + 415756294427389020*C(n,7) + 8643340000393019040*C(n,8) + 113987930725267657695*C(n,9) + 1022999668724320645050*C(n,10) + 6559258733377155798300*C(n,11) + 31097930936416379343000*C(n,12) + 111710735118080165667600*C(n,13) + 309231158315533166512800*C(n,14) + 666846639586795403736000*C(n,15) + 1126625894182469352672000*C(n,16) + 1492173540716221595232000*C(n,17) + 1541987121926231652672000*C(n,18) + 1229356526029003532160000*C(n,19) + 741102367008078915840000*C(n,20) + 326583680209195368960000*C(n,21) + 99234043419574103040000*C(n,22) + 18581137031073576960000*C(n,23) + 1615751046180311040000*C(n,24), where the coefficient of C(n,k) is the number of colorings using exactly k colors.
a(n) = A331354(n) - A331355(n) = (A331354(n) - A331357(n)) / 2 = A331355(n) - A331357(n).

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).

A337411 Array read by descending antidiagonals: T(n,k) is the number of oriented colorings of the edges of a regular n-dimensional orthoplex (cross polytope) using k or fewer colors.

Original entry on oeis.org

1, 2, 1, 3, 6, 1, 4, 24, 218, 1, 5, 70, 2285, 90054, 1, 6, 165, 703760, 1471640157, 573439556, 1, 7, 336, 10194250, 1466049174160, 6332134720430727, 50043770249328, 1, 8, 616, 90775566, 310441584462375, 629648890639384572032, 1839894096099964270283469, 59966884221697869216, 1
Offset: 1

Views

Author

Robert A. Russell, Aug 26 2020

Keywords

Comments

Each chiral pair is counted as two when enumerating oriented arrangements. 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 an octahedron with 12 edges. The number of edges is 2n*(n-1) for n>1.
Also the number of oriented colorings of the regular (n-2)-dimensional orthotopes (hypercubes) in a regular n-dimensional orthotope.

Examples

			Table begins with T(1,1):
1   2     3      4        5        6         7          8           9 ...
1   6    24     70      165      336       616       1044        1665 ...
1 218 22815 703760 10194250 90775566 576941778 2863870080 11769161895 ...
For T(2,2)=6, the arrangements are AAAA, AAAB, AABB, ABAB, ABBB, and BBBB.
		

Crossrefs

Cf. A337412 (unoriented), A337413 (chiral), A337414 (achiral).
Rows 1-4 are A000027, A006528, A060530, A331354.
Cf. A327083 (simplex edges), A337407 (orthotope edges), A325004 (orthoplex 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 + 2 x^(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, m+1]];
    FiSum[] := (Do[Fi2[k2] = Fi1[k2], {k2, Divisors[per]}];DivisorSum[per, DivisorSum[d1 = #, MoebiusMu[d1/#] Fi2[#] &]/# &]);
    CCPol[r_List] := (r1 = r; r2 = cs - r1; If[EvenQ[Sum[If[EvenQ[j3], r1[[j3]], r2[[j3]]], {j3,n}]], (per = LCM @@ Table[If[cs[[j2]] == r1[[j2]], If[0 == cs[[j2]],1,j2], 2j2], {j2,n}]; Times @@ Binomial[cs, r1] 2^(n-Total[cs]) b^FiSum[]),0]);
    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[m]=b;
    row[n_Integer] := row[n] = Factor[(Total[(PartPol[#] pc[#])&/@ IntegerPartitions[n]])/(n! 2^(n-1))]
    array[n_, k_] := row[n] /. b -> k
    Table[array[n,d+m-n], {d,8}, {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) = A337412(n,k) + A337413(n,k) = 2*A337412(n,k) - A337414(n,k) = 2*A337413(n,k) + A337414(n,k).

A337412 Array read by descending antidiagonals: T(n,k) is the number of unoriented colorings of the edges of a regular n-dimensional orthoplex (cross polytope) using k or fewer colors.

Original entry on oeis.org

1, 2, 1, 3, 6, 1, 4, 21, 144, 1, 5, 55, 12111, 49127, 1, 6, 120, 358120, 740360358, 293122232, 1, 7, 231, 5131650, 733776248840, 3168520600399659, 25174334733080, 1, 8, 406, 45528756, 155261523065875, 314848558732420555904, 920040738175691418086226, 30035285091978202824, 1
Offset: 1

Views

Author

Robert A. Russell, Aug 26 2020

Keywords

Comments

Each chiral pair is counted as one when enumerating unoriented arrangements. 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 an octahedron with 12 edges. The number of edges is 2n*(n-1) for n>1.
Also the number of unoriented colorings of the regular (n-2)-dimensional orthotopes (hypercubes) in a regular n-dimensional orthotope.

Examples

			Table begins with T(1,1):
1   2     3      4       5        6         7          8          9 ...
1   6    21     55     120      231       406        666       1035 ...
1 144 12111 358120 5131650 45528756 288936634 1433251296 5887880415 ...
For T(2,2)=6, the arrangements are AAAA, AAAB, AABB, ABAB, ABBB, and BBBB.
		

Crossrefs

Cf. A337411 (oriented), A337413 (chiral), A337414 (achiral).
Rows 1-4 are A000027, A002817, A199406, A331355.
Cf. A327084 (simplex edges), A337408 (orthotope edges), A325005 (orthoplex 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 + 2 x^(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, m+1]];
    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}]; 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[m]=b;
    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,8}, {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) = A337411(n,k) - A337413(n,k) = (A337411(n,k) + A337414(n,k)) / 2 = A337413(n,k) + A337414(n,k).

A337414 Array read by descending antidiagonals: T(n,k) is the number of achiral colorings of the edges of a regular n-dimensional orthoplex (cross polytope) using k or fewer colors.

Original entry on oeis.org

1, 2, 1, 3, 6, 1, 4, 18, 70, 1, 5, 40, 1407, 8200, 1, 6, 75, 12480, 9080559, 12804908, 1, 7, 126, 69050, 1503323520, 4906480368591, 304899216832, 1, 8, 196, 281946, 81461669375, 48226825456539776, 187380251418565888983, 103685962258536432, 1
Offset: 1

Views

Author

Robert A. Russell, Aug 26 2020

Keywords

Comments

An achiral arrangement is identical to its reflection. 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 an octahedron with 12 edges. The number of edges is 2n*(n-1) for n>1.
Also the number of achiral colorings of the regular (n-2)-dimensional orthotopes (hypercubes) in a regular n-dimensional orthotope.

Examples

			Table begins with T(1,1):
1  2    3     4     5      6      7       8       9       10 ...
1  6   18    40    75    126    196     288     405      550 ...
1 70 1407 12480 69050 281946 931490 2632512 6598935 15041950 ...
For T(2,2)=6, the arrangements are AAAA, AAAB, AABB, ABAB, ABBB, and BBBB.
		

Crossrefs

Cf. A337411 (oriented), A337412 (unoriented), A337413 (chiral).
Rows 1-4 are A000027, A002411, A331351, A331357.
Cf. A327086 (simplex edges), A337410 (orthotope edges), A325007 (orthoplex 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 + 2 x^(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, m+1]];
    FiSum[] := (Do[Fi2[k2] = Fi1[k2], {k2, Divisors[per]}];DivisorSum[per, DivisorSum[d1 = #, MoebiusMu[d1/#] Fi2[#] &]/# &]);
    CCPol[r_List] := (r1 = r; r2 = cs - r1; If[EvenQ[Sum[If[EvenQ[j3], r1[[j3]], r2[[j3]]], {j3,n}]],0,(per = LCM @@ Table[If[cs[[j2]] == r1[[j2]], If[0 == cs[[j2]],1,j2], 2j2], {j2,n}]; 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[m]=b;
    row[n_Integer] := row[n] = Factor[(Total[(PartPol[#] pc[#])&/@ IntegerPartitions[n]])/(n! 2^(n-1))]
    array[n_, k_] := row[n] /. b -> k
    Table[array[n,d+m-n], {d,8}, {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) = 2*A337412(n,k) - A337411(n,k) = A337411(n,k) - 2*A337413(n,k) = A337412(n,k) - A337413(n,k).

A337406 Number of chiral pairs of colorings of the edges of a cube (or regular octahedron) using n or fewer colors.

Original entry on oeis.org

0, 74, 10704, 345640, 5062600, 45246810, 288005144, 1430618784, 5881281480, 20827126650, 65370603320, 185725346664, 485325996064, 1181031257770, 2702889008400, 5863794289280, 12137528310384, 24099966466794
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. Both the cube and the octahedron have 12 edges.

Crossrefs

Cf. A060530 (oriented), A199406 (unoriented), A331351 (achiral).
Row 3 of A337409 (orthotope edge colorings) and A337413 (orthoplex edge colorings).

Programs

  • Mathematica
    Table[(n-1)n^2(n+1)(n^8+n^6-2n^4+8)/48, {n,20}]
    LinearRecurrence[{13,-78,286,-715,1287,-1716,1716,-1287,715,-286,78,-13,1},{0,74,10704,345640,5062600,45246810,288005144,1430618784,5881281480,20827126650,65370603320,185725346664,485325996064},20] (* Harvey P. Dale, Jul 11 2025 *)

Formula

a(n) = (n-1) * n^2 * (n+1) * (n^8 + n^6 - 2n^4 + 8) / 48.
a(n) = 74*C(n,2) + 10482*C(n,3) + 303268*C(n,4) + 3440700*C(n,5) + 19842840*C(n,6) + 65867760*C(n,7) + 133580160*C(n,8) + 168399000*C(n,9) + 128898000*C(n,10) + 54885600*C(n,11) + 9979200*C(n,12), where the coefficient of C(n,k) is the number of colorings using exactly k colors.
a(n) = (A060530(n) - A331351(n)) / 2 = A060530(n) - A199406(n) = A199406(n) - A331351(n).
G.f.: 2 * (37*x^2 + 4871*x^3 + 106130*x^4 + 691514*x^5 + 1692248*x^6 + 1692248*x^7 + 691514*x^8 + 106130*x^9 + 4871*x^10 + 37*x^11) / (1-x)^13.

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

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 66, 11158298, 0, 0, 920, 4825452718593, 314824333015938998688, 0, 0, 6350, 48038354542204960, 38491882659300767730994725249684096, 31716615393292685397985382790580028572676096, 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. For n=2, the figure is a square with one square face. For n=3, the figure is an octahedron with 8 triangular faces. For higher n, the number of triangular faces is 8*C(n,3).
Also the number of chiral pairs of colorings of the peaks of an n-dimensional orthotope (hypercube). A peak is an (n-3)-dimensional orthotope.

Examples

			Table begins with T(2,1):
 0        0             0                 0                    0 ...
 0        1            66               920                 6350 ...
 0 11158298 4825452718593 48038354542204960 60632976384183154375 ...
		

Crossrefs

Cf. A337891 (oriented), A337892 (unoriented), A337894 (achiral).
Other elements: A325006 (vertices), A337413 (edges).
Other polytopes: A337885 (simplex), A337889 (orthotope).
Rows 2-4 are A000004, A337896, A331360.

Programs

  • Mathematica
    m=2; (* dimension of color element, here a face *)
    Fi1[p1_] := Module[{g, h}, Coefficient[Product[g = GCD[k1, p1]; h = GCD[2 k1, p1]; (1 + 2 x^(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, m+1]];
    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

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) = A337891(n,k) - A337892(n,k) = (A337891(n,k) - A337894(n,k)) / 2 = A337892(n,k) - A337894(n,k).

A338148 Triangle read by rows: T(n,k) is the number of chiral pairs of colorings of the edges of a regular n-D orthoplex (or ridges of a regular n-D orthotope) using exactly k colors. Row 1 has 1 column; row n>1 has 2*n*(n-1) columns.

Original entry on oeis.org

0, 0, 0, 3, 3, 0, 74, 10482, 303268, 3440700, 19842840, 65867760, 133580160, 168399000, 128898000, 54885600, 9979200, 0, 40927, 731157018, 729348051686, 151526009158620, 11418355290999750, 415756294427389020, 8643340000393019040
Offset: 1

Views

Author

Robert A. Russell, Oct 12 2020

Keywords

Comments

Chiral colorings come in pairs, each the reflection of the other. A ridge is an (n-2)-face of an n-D polytope. For n=1, the figure is a line segment with one edge. For n=2, the figure is a square with 4 edges (vertices). For n=3, the figure is an octahedron (cube) with 12 edges. For n>1, the number of edges (ridges) is 2*n*(n-1). The Schläfli symbols for the n-D orthotope (hypercube) and the n-D orthoplex (hyperoctahedron, cross polytope) are {4,3,...,3,3} and {3,3,...,3,4} respectively, with n-2 3's in each case. The figures are mutually dual.
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

			Triangle begins with T(1,1):
  0
  0  0     3      3
  0 74 10482 303268 3440700 19842840 65867760 133580160 168399000
  ...
For T(2,3)=3, the chiral pairs are AABC-AACB, ABBC-ACBB, and ABCC-ACCB. For T(2,4)=3, the chiral pairs are ABCD-ADCB, ACBD-ADBC, and ABDC-ACDB.
		

Crossrefs

Cf. A338146 (oriented), A338147 (unoriented), A338149 (achiral), A337413 (k or fewer colors), A325010 (orthoplex vertices, orthotope facets).
Cf. A327089 (simplex), A338144 (orthotope edges, orthoplex ridges).

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 + 2 x^(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, m+1]];
    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
    Join[{{0}},Table[LinearSolve[Table[Binomial[i,j],{i,2^(m+1)Binomial[n,m+1]},{j,2^(m+1)Binomial[n,m+1]}], Table[array[n,k],{k,2^(m+1)Binomial[n,m+1]}]], {n,m+1,m+4}]] // Flatten

Formula

For n>1, A337413(n,k) = Sum_{j=1..2*n*(n-1)} T(n,j) * binomial(k,j).
T(n,k) = A338146(n,k) - A338147(n,k) = (A338146(n,k) - A338149(n,k)) / 2 = A338147(n,k) - A338149(n,k).
T(2,k) = A338144(2,k) = A325018(2,k) = A325010(2,k); T(3,k) = A338144(3,k).
Showing 1-9 of 9 results.