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.

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

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

Original entry on oeis.org

1, 3, 1, 6, 6, 1, 10, 21, 10, 1, 15, 55, 56, 15, 1, 21, 120, 220, 126, 21, 1, 28, 231, 680, 715, 252, 28, 1, 36, 406, 1771, 3060, 2002, 462, 36, 1, 45, 666, 4060, 10626, 11628, 5005, 792, 45, 1, 55, 1035, 8436, 31465, 53130, 38760, 11440, 1287, 55, 1
Offset: 1

Views

Author

Robert A. Russell, Mar 23 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. 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 orthoplex using up to k colors.

Examples

			Array begins with A(1,1):
1  3    6    10     15      21       28        36        45         55 ...
1  6   21    55    120     231      406       666      1035       1540 ...
1 10   56   220    680    1771     4060      8436     16215      29260 ...
1 15  126   715   3060   10626    31465     82251    194580     424270 ...
1 21  252  2002  11628   53130   201376    658008   1906884    5006386 ...
1 28  462  5005  38760  230230  1107568   4496388  15890700   50063860 ...
1 36  792 11440 116280  888030  5379616  26978328 115775100  436270780 ...
1 45 1287 24310 319770 3108105 23535820 145008513 752538150 3381098545 ...
For A(1,2) = 3, 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. A325004 (oriented), A325006 (chiral), A325007 (achiral), A325009 (exactly k colors).
Other n-dimensional polytopes: A325000 (simplex), A325013 (orthoplex).
Rows 1-3 are A000217, A002817, A198833.

Programs

  • Mathematica
    Table[Binomial[Binomial[d-n+2,2]+n-1,n],{d,1,11},{n,1,d}] // Flatten

Formula

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

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

Original entry on oeis.org

0, 1, 0, 3, 0, 0, 6, 3, 1, 0, 10, 15, 66, 94, 0, 15, 45, 920, 97974, 1047816, 0, 21, 105, 6350, 10700090, 481141220994, 400140831558512, 0, 28, 210, 29505, 390081800, 4802390808840576, 74515656021475803734579625, 527471421741473576372948457251328, 0
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. 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 up to k colors.

Examples

			Array begins with A(1,1):
0  1     3        6        10         15          21           28 ...
0  0     3       15        45        105         210          378 ...
0  1    66      920      6350      29505      106036       317856 ...
0 94 97974 10700090 390081800 7280687610 86121007714 730895668104 ...
For A(2,3)=3, each square has one of the three colors on two adjacent edges.
		

Crossrefs

Cf. A325012 (oriented), A325013 (unoriented), A325015 (achiral), A325018 (exactly k colors).
Other n-dimensional polytopes: A007318(k,n+1) (simplex), A325006 (orthotope).
Rows 1-2 are A161680, A050534.

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(k,n) = A325012(n,k) - A325013(n,k) = (A325012(n,k) - A325015(n,k)) / 2 = A325013(n,k) - A325015(n,k).
A(n,k) = Sum_{j=2..2^n} A325018(n,j) * binomial(k,j).

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

A128767 Number of inequivalent n-colorings of the 4D hypercube under the full orthogonal group of the cube (of order 2^4*4! = 384).

Original entry on oeis.org

1, 402, 132102, 11756666, 405385550, 7416923886, 86986719477, 735192450952, 4834517667381, 26073250910950, 119759687845446, 481750080584202, 1733588303252702, 5673534527793146, 17109303241791825, 48047227408513056
Offset: 1

Views

Author

Ricardo Perez-Aguila (ricardo.perez.aguila(AT)gmail.com), Apr 04 2007

Keywords

Comments

I assume this refers to colorings of the vertices of the cube. - N. J. A. Sloane, Apr 06 2007
Number of unoriented colorings of the 16 tetrahedral facets of a hyperoctahedron or of the 16 vertices of a tesseract (4-D cube) using up to n colors. Each chiral pair is counted as one when enumerating unoriented arrangements. 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. - Robert A. Russell, Oct 03 2020

Examples

			a(2)=402 because there are 402 inequivalent 2-colorings of the 4D hypercube.
		

References

  • Banks, D. C.; Linton, S. A. & Stockmeyer, P. K. Counting Cases in Substitope Algorithms. IEEE Transactions on Visualization and Computer Graphics, Vol. 10, No. 4, pp. 371-384. 2004.
  • Perez-Aguila, Ricardo. Enumerating the Configurations in the n-Dimensional Orthogonal Polytopes Through Polya's Counting and A Concise Representation. Proceedings of the 3rd International Conference on Electrical and Electronics Engineering and XII Conference on Electrical Engineering ICEEE and CIE 2006, pp. 63-66.
  • Polya, G. & Read R. C. Combinatorial Enumeration of Groups, Graphs and Chemical Compounds. Springer-Verlag, 1987.

Crossrefs

Cf. A337952 (oriented), A337954 (chiral), A337955 (achiral).
Other elements: A331359 (tesseract edges, hyperoctahedron faces), A331355 (tesseract faces, hyperoctahedron edges), A337957 (tesseract facets, hyperoctahedron vertices).
Other polychora: A000389(n+4) (4-simplex facets/vertices), A338949 (24-cell), A338965 (120-cell, 600-cell).
Row 4 of A325013 (orthoplex facets, orthotope vertices).

Programs

  • Mathematica
    Table[(1/384)*( 48*n^2 + 180*n^4 + 48*n^6 + 83*n^8 + 12*n^10 + 12*n^12 + n^16),{n,30}]

Formula

a(n) = (1/384)*(48*n^2 + 180*n^4 + 48*n^6 + 83*n^8 + 12*n^10 + 12*n^12 + n^16)
G.f.: -x*(x +1)*(x^14 +384*x^13 +125020*x^12 +9439904*x^11 +213777216*x^10 +1821620108*x^9 +6527222787*x^8 +10098845160*x^7 +6527222787*x^6 +1821620108*x^5 +213777216*x^4 +9439904*x^3 +125020*x^2 +384*x +1) / (x -1)^17. [Colin Barker, Dec 04 2012]
From Robert A. Russell, Oct 03 2020: (Start)
a(n) = 1*C(n,1) + 400*C(n,2) + 130899*C(n,3) + 11230666*C(n,4) + 347919225*C(n,5) + 5158324560*C(n,6) + 43174480650*C(n,7) + 225086553300*C(n,8) + 775894225050*C(n,9) + 1831178115900*C(n,10) + 3008073915000*C(n,11) + 3439243962000*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 unoriented colorings using exactly k colors.
a(n) = A337952(n) - A337954(n) = (A337952(n) + A337955(n)) / 2 = A337954(n) + A337955(n).
(End)

A325017 Triangle read by rows: T(n,k) is the number of unoriented 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

1, 1, 1, 4, 6, 3, 1, 20, 204, 1056, 2850, 4080, 2940, 840, 1, 400, 130899, 11230666, 347919225, 5158324560, 43174480650, 225086553300, 775894225050, 1831178115900, 3008073915000, 3439243962000, 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. 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 exactly k colors.

Examples

			Triangle begins with T(1,1):
1  1
1  4   6    3
1 20 204 1056 2850 4080 2940 840
For T(2,2)=4, two squares have three edges the same color, one has opposite edges the same color, and one has opposite edges different colors.
		

Crossrefs

Cf. A325016 (oriented), A325018 (chiral), A325019 (achiral), A325013 (up to k colors).
Other n-dimensional polytopes: A007318(n,k-1) (simplex), A325009 (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 (* A325013 *)
    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

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

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

Original entry on oeis.org

1, 2, 1, 3, 6, 1, 4, 21, 144, 1, 5, 55, 12111, 11251322, 1, 6, 120, 358120, 4825746875682, 314824456456819827136, 1, 7, 231, 5131650, 48038446526132256, 38491882660019692002988737797054040, 136221825854745676520058554256163406987047485113810944, 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 a cube with 12 edges. The number of edges is n*2^(n-1).
Also the number of unoriented colorings of the regular (n-2)-dimensional simplexes in a regular n-dimensional orthoplex.

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. A337407 (oriented), A337409 (chiral), A337410 (achiral).
Rows 1-4 are A000027, A002817, A199406, A331359.
Cf. A327084 (simplex edges), A337412 (orthoplex edges), A325013 (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 + 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, 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}]; 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) - A337409(n,k) = (A337407(n,k) - A337410(n,k)) / 2 = A337409(n,k) + A337410(n,k).

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

Original entry on oeis.org

1, 2, 1, 3, 10, 1, 4, 56, 49127, 1, 5, 220, 740360358, 314824532572147370464, 1, 6, 680, 733776248840, 38491882660671134164965704408524083, 38343035259947576596859948806931124970404417593861154473053467181056, 1
Offset: 2

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

Each chiral pair is counted as one when enumerating unoriented arrangements. 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 unoriented colorings of peaks of an n-dimensional orthoplex. A peak is an (n-3)-dimensional simplex.

Examples

			Array begins with T(2,1):
 1     2         3            4               5                 6 ...
 1    10        56          220             680              1771 ...
 1 49127 740360358 733776248840 155261523065875 12340612271439081 ...
		

Crossrefs

Cf. A337887 (oriented), A337889 (chiral), A337890 (achiral).
Other elements: A325013 (vertices), A337408 (edges).
Other polytopes: A337884 (simplex), A337892 (orthoplex).
Rows 2-4 are A000027, A198833, A331355.

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 + 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, 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}]; 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) = A337887(n,k) - A337889(n,k) = (A337887(n,k) + A337890(n,k)) / 2 = A337889(n,k) + A337890(n,k).
Showing 1-9 of 9 results.