A337408 -id:A337408 - cp's OEIS frontend

A199406 The number of inequivalent ways to color the edges of a cube using at most n colors.

1, 144, 12111, 358120, 5131650, 45528756, 288936634, 1433251296, 5887880415, 20842168600, 65402344161, 185788177224, 485443851256, 1181242399260, 2703252560100, 5864398969216, 12138503871789, 24101498435616, 46112016365155, 85335258695400, 153249227870046

Offset: 1

Geoffrey Critzer, Nov 05 2011

Two edge colorings are equivalent if one is the mirror image of the other or the cube can be picked up and rotated in any manner to obtain the other.
The group here has order 48 (compare A060530). - N. J. A. Sloane, Aug 14 2012
Also the number of unoriented colorings of the 12 edges of a regular octahedron with n or fewer colors. The Schläfli symbols of the cube and octahedron are {4,3} and {3,4} respectively. They are mutually dual. For an unoriented coloring, chiral pairs are counted as one. - Robert A. Russell, Oct 17 2020

T. D. Noe, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (13,-78,286,-715,1287,-1716,1716,-1287,715,-286,78,-13,1).

Cf. A060530 (oriented), A337406 (chiral), A331351 (achiral), A128766 (cube vertices, octahedron faces), A198833 (cube faces, octahedron vertices), A063842(n-1) (tetrahedron), A337963 (dodecahedron, icosahedron).

Row 3 of A337408 (orthotope edges, orthoplex ridges) and A337412 (orthoplex edges, orthotope ridges).

Table[CycleIndex[KSubsetGroup[Automorphisms[CubicalGraph], Edges[CubicalGraph]],s] /. Table[s[i]->n, {i,1,6}], {n,1,15}]
Table[(8n^2+12n^3+8n^4+4n^6+12n^7+3n^8+n^12)/48, {n,20}] (* Robert A. Russell, Oct 17 2020 *)

a(n) = n^12/48 + n^8/16 + n^7/4 + n^6/12 + n^4/6 + n^3/4 + n^2/6.
Cycle index = 1/48(s_1^12+3s_1^4s_2^4+12s_1^2s_2^5+4s_2^6+8s_3^4+12s_4^3+8s_6^2).
G.f.: -x*(76*x^10 +10016*x^9 +212772*x^8 +1380453*x^7 +3384939*x^6 +3388593*x^5 +1380279*x^4 +211623*x^3 +10317*x^2 +131*x +1)/(x -1)^13. [Colin Barker, Aug 13 2012]
From Robert A. Russell, Oct 17 2020: (Start)
a(n) = A060530(n) - A337406(n) = (A060530(n) + A331351(n)) / 2 = A337406(n) + A331351(n).
a(n) = 1*C(n,1) + 142*C(n,2) + 11682*C(n,3) + 310536*C(n,4) + 3460725*C(n,5) + 19870590*C(n,6) + 65886660*C(n,7) + 133585200*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 unoriented colorings using exactly k colors. (End)

A327084 Array read by descending antidiagonals: A(n,k) is the number of unoriented 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, 10, 11, 1, 5, 20, 66, 34, 1, 6, 35, 276, 792, 156, 1, 7, 56, 900, 10688, 25506, 1044, 1, 8, 84, 2451, 90005, 1601952, 2302938, 12346, 1, 9, 120, 5831, 533358, 43571400, 892341888, 591901884, 274668

Offset: 1

Robert A. Russell, Aug 19 2019

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 unoriented colorings are the same if congruent; chiral pairs are counted as one.

A(n,k) is also the number of unoriented 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 unoriented colorings of the vertices (0-dimensional simplices) of an equilateral triangle.

			Array begins with A(1,1):
  1  2   3     4     5      6       7       8        9       10 ...
  1  4  10    20    35     56      84     120      165      220 ...
  1 11  66   276   900   2451    5831   12496    24651    45475 ...
  1 34 792 10688 90005 533358 2437848 9156288 29522961 84293770 ...
  ...
For A(2,3) = 10, the nine achiral colorings are AAA, AAB, AAC, ABB, ACC, BBB, BBC, BCC, and CCC. The chiral pair is ABC-ACB.

Chai Wah Wu, Table of n, a(n) for n = 1..1528 (terms 1..325 from Robert A. Russell)
Harald Fripertinger, The cycle type of the induced action on 2-subsets
E. M. Palmer and R. W. Robinson, Enumeration under two representations of the wreath product, Acta Math., 131 (1973), 123-143.

Cf. A327083 (oriented), A327085 (chiral), A327086 (achiral), A327088 (exactly k colors), A325000 (vertices, facets), A337884 (faces, peaks), A337408 (orthotope edges, orthoplex ridges), A337412 (orthoplex edges, orthotope ridges).
Rows 1-4 are A000027, A000292, A063842(n-1), A063843.
Cf. A063841 (k-multigraphs on n nodes).

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[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 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[pc[#]k^ex[#] &/@ IntegerPartitions[n+1]]/(n+1)!
Table[array[n,d-n+1], {d,10}, {n,d}] // Flatten
(* Another program (translated from Andrew Howroyd's PARI code): *)
permcount[v_] := Module[{m=1, s=0, k=0, t}, For[i=1, i <= Length[v], i++, t = v[[i]]; k = If[i>1 && t == v[[i-1]], k+1, 1]; m *= t*k; s += t]; s!/m];
edges[v_] := Sum[GCD[v[[i]], v[[j]]], {i, 2, Length[v]}, {j, 1, i-1}] + Total[Quotient[v, 2]];
T[n_, k_] := Module[{s = 0}, Do[s += permcount[p]*k^edges[p], {p, IntegerPartitions[n+1]}]; s/(n+1)!];
Table[T[n-k+1, k], {n, 1, 9}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Jan 08 2021 *)

permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
edges(v) = {sum(i=2, #v, sum(j=1, i-1, gcd(v[i], v[j]))) + sum(i=1, #v, v[i]\2)}
T(n, k) = {my(s=0); forpart(p=n+1, s+=permcount(p)*k^edges(p)); s/(n+1)!} \\ Andrew Howroyd, Sep 06 2019

from itertools import combinations
from math import prod, gcd, factorial
from fractions import Fraction
from sympy.utilities.iterables import partitions
def A327084_T(n,k): return int(sum(Fraction(k**(sum(p[r]*p[s]*gcd(r,s) for r,s in combinations(p.keys(),2))+sum((q>>1)*r+(q*r*(r-1)>>1) for q, r in p.items())),prod(q**r*factorial(r) for q, r in p.items())) for p in partitions(n+1))) # Chai Wah Wu, Jul 09 2024

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} A327088(n,j) * binomial(k,j).
A(n,k) = A327083(n,k) - A327085(n,k) = (A327083(n,k) + A327086(n,k)) / 2 = A327085(n,k) + A327086(n,k).
A(n,k) = A063841(n+1,k-1).

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

Robert A. Russell, Aug 26 2020

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.

			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.

K. Balasubramanian, Computational enumeration of colorings of hyperplanes of hypercubes for all irreducible representations and applications, J. Math. Sci. & Mod. 1 (2018), 158-180.

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

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

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

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

Robert A. Russell, Aug 26 2020

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.

			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.

K. Balasubramanian, Computational enumeration of colorings of hyperplanes of hypercubes for all irreducible representations and applications, J. Math. Sci. & Mod. 1 (2018), 158-180.

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

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

T(n,k) = A337411(n,k) - A337413(n,k) = (A337411(n,k) + A337414(n,k)) / 2 = A337413(n,k) + A337414(n,k).

A337407 Array read by descending antidiagonals: T(n,k) is the number of oriented 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, 24, 218, 1, 5, 70, 22815, 22409620, 1, 6, 165, 703760, 9651199594275, 629648865090036960064, 1, 7, 336, 10194250, 96076801068337216, 76983765319971869475595432431084156, 272443651709491352597039736725488834366101875164020736, 1

Offset: 1

Robert A. Russell, Aug 26 2020

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 a cube with 12 edges. The number of edges is n*2^(n-1).

Also the number of oriented colorings of the regular (n-2)-dimensional simplexes in a regular n-dimensional orthoplex.

			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.

K. Balasubramanian, Computational enumeration of colorings of hyperplanes of hypercubes for all irreducible representations and applications, J. Math. Sci. & Mod. 1 (2018), 158-180.

Cf. A337408 (unoriented), A337409 (chiral), A337410 (achiral).
Rows 1-4 are A000027, A006528, A060530, A331358.
Cf. A327083 (simplex edges), A337411 (orthoplex edges), A325012 (orthotope vertices).

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

T(n,k) = A337408(n,k) + A337409(n,k) = 2*A337408(n,k) - A337410(n,k) = 2*A337409(n,k) + A337410(n,k).

A337410 Array read by descending antidiagonals: T(n,k) is the number of achiral 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, 18, 70, 1, 5, 40, 1407, 93024, 1, 6, 75, 12480, 294157089, 47823602694208, 1, 7, 126, 69050, 91983927296, 67514530382043163023924, 443077371786837979607993095063601152, 1

Offset: 1

Robert A. Russell, Aug 26 2020

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 a cube with 12 edges. The number of edges is n*2^(n-1).

Also the number of achiral colorings of the regular (n-2)-dimensional simplexes in a regular n-dimensional orthoplex.

			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.

K. Balasubramanian, Computational enumeration of colorings of hyperplanes of hypercubes for all irreducible representations and applications, J. Math. Sci. & Mod. 1 (2018), 158-180.

Cf. A337407 (oriented), A337408 (unoriented), A337409 (chiral).
Rows 1-4 are A000027, A002411, A331351, A331361.
Cf. A327086 (simplex edges), A337414 (orthoplex edges), A325015 (orthotope vertices).

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

T(n,k) = 2*A337408(n,k) - A337407(n,k) = A337407(n,k) - 2*A337409(n,k) = A337408(n,k) - A337409(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

Robert A. Russell, Sep 28 2020

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.

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

K. Balasubramanian, Computational enumeration of colorings of hyperplanes of hypercubes for all irreducible representations and applications, J. Math. Sci. & Mod. 1 (2018), 158-180.

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.

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

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

A338143 Triangle read by rows: T(n,k) is the number of unoriented colorings of the edges of a regular n-D orthotope (or ridges of a regular n-D orthoplex) using exactly k colors. Row n has n*2^(n-1) columns.

Original entry on oeis.org

1, 1, 4, 6, 3, 1, 142, 11682, 310536, 3460725, 19870590, 65886660, 133585200, 168399000, 128898000, 54885600, 9979200, 1, 11251320, 4825713121719, 48019143606137456, 60392840368910627325

Offset: 1

Robert A. Russell, Oct 12 2020

Each chiral pair is counted as one when enumerating unoriented arrangements. 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 a cube (octahedron) with 12 edges. The number of edges (ridges) is n*2^(n-1). The Schläfli symbols for the n-D orthotope (hypercube) and the n-D orthoplex (hyperoctahedron, cross polytope) are {4,...,3,3} and {3,3,...,4} respectively, with n-2 3's in each case. The figures are mutually dual.

			Triangle begins with T(1,1):
  1
  1   4     6      3
  1 142 11682 310536 3460725 19870590 65886660 133585200 168399000
  ...

K. Balasubramanian, Computational enumeration of colorings of hyperplanes of hypercubes for all irreducible representations and applications, J. Math. Sci. & Mod. 1 (2018), 158-180.

Cf. A338142 (oriented), A338144 (chiral), A338145 (achiral), A337408 (k or fewer colors), A325017 (orthotope vertices, orthoplex facets).

Cf. A327088 (simplex), A338147 (orthoplex edges, orthotope ridges).

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

A337408(n,k) = Sum_{j=1..n*2^(n-1)} T(n,j) * binomial(k,j).
T(n,k) = A338142(n,k) - A338144(n,k) = (A338142(n,k) + A338145(n,k)) / 2 = A338144(n,k) + A338145(n,k).
T(2,k) = A338147(2,k) = A325017(2,k) = A325009(2,k); T(3,k) = A338147(3,k).

cp's OEIS Frontend

A199406 The number of inequivalent ways to color the edges of a cube using at most n colors.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A338143 Triangle read by rows: T(n,k) is the number of unoriented colorings of the edges of a regular n-D orthotope (or ridges of a regular n-D orthoplex) using exactly k colors. Row n has n*2^(n-1) columns.