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.

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

Original entry on oeis.org

1, 2, 1, 3, 4, 1, 4, 11, 12, 1, 5, 24, 87, 40, 1, 6, 45, 416, 1197, 184, 1, 7, 76, 1475, 18592, 42660, 1296, 1, 8, 119, 4236, 166885, 3017600, 4223313, 17072, 1, 9, 176, 10437, 1019880, 85025050, 1748176768, 1139277096, 424992
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. Two oriented colorings are the same if one is a rotation of the other; chiral pairs are counted as two.
A(n,k) is also the number of oriented 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 oriented colorings of the vertices (0-dimensional simplices) of an equilateral triangle.

Examples

			Array begins with A(1,1):
  1  2    3     4      5       6       7        8        9        10 ...
  1  4   11    24     45      76     119      176      249       340 ...
  1 12   87   416   1475    4236   10437    22912    45981     85900 ...
  1 40 1197 18592 166885 1019880 4738153 17962624 58248153 166920040 ...
  ...
For A(2,3) = 11, the nine achiral colorings are AAA, AAB, AAC, ABB, ACC, BBB, BBC, BCC, and CCC. The chiral pair is ABC-ACB.
		

Crossrefs

Cf. A327084 (unoriented), A327085 (chiral), A327086 (achiral), A327087 (exactly k colors), A324999 (vertices, facets), A337883 (faces, peaks), A337407 (orthotope edges, orthoplex ridges), A337411 (orthoplex edges, orthotope ridges).
Rows 1-4 are A000027, A006527, A046023, A331350.
Column 2 is A218144(n+1).

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]]], pc[#] j^Total[CycleX[#]][[2]], 0] & /@ IntegerPartitions[n+1]]/((n+1)!/2)]
    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]]], pc[#]k^ex[#], 0] &/@ IntegerPartitions[n+1]]/((n+1)!/2)
    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} A327087(n,j) * binomial(k,j).
A(n,k) = A327084(n,k) + A327085(n,k) = 2*A327084(n,k) - A327086(n,k) = 2*A327085(n,k) + A327086(n,k).

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

Original entry on oeis.org

1, 4, 1, 9, 4, 1, 16, 11, 5, 1, 25, 24, 15, 6, 1, 36, 45, 36, 21, 7, 1, 49, 76, 75, 56, 28, 8, 1, 64, 119, 141, 127, 84, 36, 9, 1, 81, 176, 245, 258, 210, 120, 45, 10, 1, 100, 249, 400, 483, 463, 330, 165, 55, 11, 1, 121, 340, 621, 848, 931, 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 oriented colorings are the same if one is a rotation of the other; chiral pairs are counted as two.

Examples

			The array begins with A(1,1):
  1  4  9  16  25   36   49    64    81   100   121    144    169    196 ...
  1  4 11  24  45   76  119   176   249   340   451    584    741    924 ...
  1  5 15  36  75  141  245   400   621   925  1331   1860   2535   3381 ...
  1  6 21  56 127  258  483   848  1413  2254  3465   5160   7475  10570 ...
  1  7 28  84 210  463  931  1744  3087  5215  8470  13300  20280  30135 ...
  1  8 36 120 330  792 1717  3440  6471 11560 19778  32616  52104  80952 ...
  1  9 45 165 495 1287 3003  6436 12879 24355 43923  76077 127257 206493 ...
  1 10 55 220 715 2002 5005 11440 24311 48630 92433 168180 294645 499422 ...
  ...
For A(1,2) = 4, the two achiral colorings use just one of the two colors for both vertices; the chiral pair uses two colors. For A(2,2)=4, the triangle may have 0, 1, 2, or 3 edges of one color.
		

Crossrefs

Cf. A325000 (unoriented), A325000(n,k-n) (chiral), A325001 (achiral), A325002 (exactly k colors), A327083 (edges, ridges), A337883 (faces, peaks), A325004 (orthotope facets, orthoplex vertices), A325012 (orthoplex facets, orthotope vertices).
Rows 1-4 are A000290, A006527, A006008, A337895.

Programs

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

Formula

A(n,k) = binomial(n+k,n+1) + binomial(k,n+1).
A(n,k) = Sum_{j=1..n+1} A325002(n,j) * binomial(k,j).
A(n,k) = A325000(n,k) + A325000(n,k-n) = 2*A325000(n,k) - A325001(n,k) = 2*A325000(n,k-n) + A325001(n,k).
G.f. for row n: (x + x^(n+1)) / (1-x)^(n+2).
Linear recurrence for row n: A(n,k) = Sum_{j=1..n+2} -binomial(j-n-3,j) * A(n,k-j).
G.f. for column k: (1 - 2*(1-x)^k + (1-x^2)^k) / (x*(1-x)^k) - 2*k.

A331350 Number of oriented colorings of the edges (or triangular faces) of a regular 4-dimensional simplex with n available colors.

Original entry on oeis.org

1, 40, 1197, 18592, 166885, 1019880, 4738153, 17962624, 58248153, 166920040, 432738229, 1032709536, 2298857821, 4822806184, 9613704465, 18329410048, 33605960689, 59516325288, 102196242685, 170682720160, 278019522837
Offset: 1

Views

Author

Robert A. Russell, Jan 14 2020

Keywords

Comments

A 4-dimensional simplex has 5 vertices and 10 edges. Its Schläfli symbol is {3,3,3}. Two oriented colorings are the same if one is a rotation of the other; chiral pairs are counted as two.
There are 60 elements in the rotation group of the 4-dimensional simplex. Each is an even permutation of the vertices and can be associated with a partition of 5 based on the conjugacy group of the permutation. The first formula is obtained by averaging their cycle indices after replacing x_i^j with n^j according to the Pólya enumeration theorem.
Partition Count Even Cycle Indices
5 24 x_5^2
311 20 x_1^1x_3^3
221 15 x_1^2x_2^4
11111 1 x_1^10

Crossrefs

Cf. A063843 (unoriented), A331352 (chiral), A331353 (achiral).
Other polychora: A331358 (8-cell), A331354 (16-cell), A338952 (24-cell), A338964 (120-cell, 600-cell).
Row 4 of A327083 (simplex edges and facets) and A337883 (simplex faces and peaks).

Programs

  • Mathematica
    Table[(24n^2 + 20n^4 + 15n^6 + n^10)/60, {n, 1, 25}]

Formula

a(n) = (24*n^2 + 20*n^4 + 15*n^6 + n^10) / 60.
a(n) = C(n,1) + 38*C(n,2) + 1080*C(n,3) + 14040*C(n,4) + 85500*C(n,5) + 274104*C(n,6) + 493920*C(n,7) + 504000*C(n,8) + 272160*C(n,9) + 60480*C(n,10), where the coefficient of C(n,k) is the number of colorings using exactly k colors.
a(n) = A063843(n) + A331352(n) = 2*A063843(n) - A331353(n) = 2*A331352(n) + A331353(n).

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

Original entry on oeis.org

1, 2, 1, 3, 5, 1, 4, 15, 34, 1, 5, 35, 792, 2136, 1, 6, 70, 10688, 4977909, 7013320, 1, 7, 126, 90005, 1533771392, 9930666709494, 1788782616656, 1, 8, 210, 533358, 132597435125, 234249157811872000, 12979877431438089379035, 53304527811667897248, 1
Offset: 2

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

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

Examples

			Table begins with T(2,1):
 1    2       3          4            5             6               7 ...
 1    5      15         35           70           126             210 ...
 1   34     792      10688        90005        533358         2437848 ...
 1 2136 4977909 1533771392 132597435125 5079767935320 110837593383153 ...
For T(3,4)=35, the 34 achiral arrangements are AAAA, AAAB, AAAC, AAAD, AABB, AABC, AABD, AACC, AACD, AADD, ABBB, ABBC, ABBD, ABCC, ABDD, ACCC, ACCD, ACDD, ADDD, BBBB, BBBC, BBBD, BBCC, BBCD, BBDD, BCCC, BCCD, BCDD, BDDD, CCCC, CCCD, CCDD, CDDD, and DDDD. The chiral pair is ABCD-ABDC.
		

Crossrefs

Cf. A337883 (oriented), A337885 (chiral), A337886 (achiral), A051168 (binary Lyndon words).
Other elements: A325000 (vertices), A327084 (edges).
Other polytopes: A337888 (orthotope), A337892 (orthoplex).
Rows 2-4 are A000027, A000332(n+3), A063843.

Programs

  • Mathematica
    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[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

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 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) - A337885(n,k) = (A337883(n,k) + A337886(n,k)) / 2 = A337885(n,k) + A337886(n,k).

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

Original entry on oeis.org

1, 2, 1, 3, 5, 1, 4, 15, 28, 1, 5, 34, 387, 768, 1, 6, 65, 2784, 202203, 302032, 1, 7, 111, 13125, 11230976, 7109211078, 3098988832, 1, 8, 175, 46836, 254729375, 9393953524224, 50669807706182691, 1831011525739328, 1
Offset: 2

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

An achiral arrangement is identical to its reflection. 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 achiral colorings of the peaks of a regular n-dimensional simplex. A peak of an n-simplex is an (n-3)-dimensional simplex.

Examples

			Table begins with T(2,1):
1   2      3        4         5          6           7            8 ...
1   5     15       34        65        111         175          260 ...
1  28    387     2784     13125      46836      137543       349952 ...
1 768 202203 11230976 254729375 3267720576 28271133933 183296831488 ...
For T(3,4)=34, the 34 achiral arrangements are AAAA, AAAB, AAAC, AAAD, AABB, AABC, AABD, AACC, AACD, AADD, ABBB, ABBC, ABBD, ABCC, ABDD, ACCC, ACCD, ACDD, ADDD, BBBB, BBBC, BBBD, BBCC, BBCD, BBDD, BCCC, BCCD, BCDD, BDDD, CCCC, CCCD, CCDD, CDDD, and DDDD.
		

Crossrefs

Cf. A337883 (oriented), A337884 (unoriented), A337885 (chiral), A051168 (binary Lyndon words).
Other elements: A325001 (vertices), A327086 (edges).
Other polytopes: A337890 (orthotope), A337894 (orthoplex).
Rows 2-4 are A000027, A006003, A331353.

Programs

  • Mathematica
    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[OddQ[Total[1-Mod[#, 2]]], pc[#] j^Total[CX[#, m+1]][[2]], 0] & /@ IntegerPartitions[n+1]]/((n+1)!/2)]
    array[n_, k_] := row[n] /. j -> 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 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) = A337884(n,k) - A337883(n,k) = A337883(n,k) - 2*A337885(n,k) = A337884(n,k) - A337885(n,k).

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.

Original entry on oeis.org

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

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

Examples

			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.
		

Crossrefs

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.

Programs

  • Mathematica
    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

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

A337887 Array read by descending antidiagonals: T(n,k) is the number of oriented 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, 57, 90054, 1, 5, 240, 1471640157, 629648865588086369152, 1, 6, 800, 1466049174160, 76983765319971901895960429658208179, 76686070519895153193719509580895099970955878067526648007224125292544, 1
Offset: 2

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

Each chiral pair is counted as two when enumerating oriented 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 oriented 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         57           240             800              2226 ...
 1 90054 1471640157 1466049174160 310441584462375 24679078461920106 ...
		

Crossrefs

Cf. A337888 (unoriented), A337889 (chiral), A337890 (achiral).
Other elements: A325012 (vertices), A337407 (edges).
Other polytopes: A337883 (simplex), A337891 (orthoplex).
Rows 2-4 are A000027, A047780, A331354.

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; 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[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,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) = A337888(n,k) + A337889(n,k) = 2*A337888(n,k) - A337890(n,k) = 2*A337889(n,k) + A337890(n,k).

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

Original entry on oeis.org

1, 2, 1, 3, 23, 1, 4, 333, 22409620, 1, 5, 2916, 9651199594275, 629648865588086369152, 1, 6, 16725, 96076801068337216, 76983765319971901895960429658208179, 63433230786931550329738915431918588874940416, 1
Offset: 2

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

Each chiral pair is counted as two when enumerating oriented arrangements. 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 oriented colorings of the peaks of an n-dimensional orthotope (hypercube). A peak is an (n-3)-dimensional orthotope.
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):
1        2             3                 4                     5 ...
1       23           333              2916                 16725 ...
1 22409620 9651199594275 96076801068337216 121265960728368199375 ...
		

Crossrefs

Cf. A337892 (unoriented), A337893 (chiral), A337894 (achiral).
Other elements: A325004 (vertices), A337411 (edges).
Other polytopes: A337883 (simplex), A337887 (orthotope).
Rows 2-4 are A000027, A000543, A331358

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; 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,6}, {n,m,d+m-1}] // Flatten

Formula

T(n,k) = A337892(n,k) + A337893(n,k) = 2*A337892(n,k) - A337894(n,k) = 2*A337893(n,k) + A337894(n,k).

A338113 Triangle read by rows: T(n,k) is the number of oriented colorings of the faces (and peaks) of a regular n-dimensional simplex using exactly k colors. Row n has C(n+1,3) columns.

Original entry on oeis.org

1, 1, 3, 3, 2, 1, 38, 1080, 14040, 85500, 274104, 493920, 504000, 272160, 60480, 1, 3502, 9743106, 3017318368, 249756082950, 8612276962188, 156010151929968, 1699145259725088, 12107373916276800, 59649257217110400
Offset: 2

Views

Author

Robert A. Russell, Oct 10 2020

Keywords

Comments

An n-dimensional simplex has n+1 vertices, C(n+1,3) faces, and C(n+1,3) peaks, which are (n-3)-dimensional simplexes. For n=2, the figure is a triangle with one face. For n=3, the figure is a tetrahedron with four triangular faces and four peaks (vertices). For n=4, the figure is a 4-simplex with ten triangular faces and ten peaks (edges). The Schläfli symbol {3,...,3}, of the regular n-dimensional simplex consists of n-1 3's. Two oriented colorings are the same if one is a rotation of the other; chiral pairs are counted as two.
The algorithm used in the Mathematica program below assigns each permutation of the vertices to a cycle-structure partition of n+1. It then determines the number of permutations for each partition and the cycle index for each partition. 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(2,1):
  1
  1  3    3     2
  1 38 1080 14040 85500 274104 493920 504000 272160 60480
  ...
For T(3,2)=3, the tetrahedron has one, two, or three faces (vertices) of one color.
		

Crossrefs

Cf. A338114 (unoriented), A338115 (chiral), A338116 (achiral), A337883 (k or fewer colors), A325002 (vertices and facets), A327087 (edges and ridges).

Programs

  • Mathematica
    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]]], pc[#] j^Total[CX[#, m+1]][[2]], 0] & /@ IntegerPartitions[n+1]]/((n+1)!/2)]
    array[n_, k_] := row[n] /. j -> k
    Table[LinearSolve[Table[Binomial[i,j],{i,Binomial[n+1,m+1]},{j,Binomial[n+1,m+1]}], Table[array[n,k],{k,Binomial[n+1,m+1]}]], {n,m,m+4}] // Flatten

Formula

A337883(n,k) = Sum_{j=1..C(n+1,3)} T(n,j) * binomial(k,j).
T(n,k) = A338114(n,k) + A338115(n,k) = 2*A338114(n,k) - A338116(n,k) = 2*A338115(n,k) + A338116(n,k).
T(3,k) = A325002(3,k); T(4,k) = A327087(4,k).
Showing 1-9 of 9 results.