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.

Previous Showing 11-12 of 12 results.

A327090 Triangle read by rows: T(n,k) is the number of achiral colorings of the edges of a regular n-dimensional simplex using exactly k colors. Row n has (n+1)*n/2 columns.

Original entry on oeis.org

1, 1, 2, 0, 1, 8, 18, 12, 0, 0, 1, 26, 306, 1400, 2800, 2520, 840, 0, 0, 0, 1, 126, 7971, 153660, 1268475, 5463990, 13534290, 20018880, 17478720, 8316000, 1663200, 0, 0, 0, 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. An achiral coloring is identical to its reflection.
T(n,k) is also the number of achiral colorings of (n-2)-dimensional regular simplices in an n-dimensional simplex using exactly k colors. Thus, T(2,k) is also the number of achiral colorings of the vertices (0-dimensional simplices) of an equilateral triangle.

Examples

			Triangle begins with T(1,1):
1
1  2   0
1  8  18   12    0    0
1 26 306 1400 2800 2520 840 0 0 0
For T(2,2) = 2, the two colorings of the triangle edges are AAB and ABB.
		

Crossrefs

Cf. A327087 (oriented), A327088 (unoriented), A327089 (chiral), A327086 (up to k colors), A325003 (vertices).

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]]], 0, pc[#] j^Total[CycleX[#]][[2]]] & /@ IntegerPartitions[n+1]]/((n+1)!/2)]
    array[n_, k_] := row[n] /. j -> k
    Table[LinearSolve[Table[Binomial[i,j], {i,1,(n+1)n/2}, {j,1,(n+1)n/2}], Table[array[n,k], {k,1,(n+1)n/2}]], {n,1,6}] // 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. The last n-1 entries in row n are zero.
A327086(n,k) = Sum_{j=1..(n+1)*n/2} T(n,j) * binomial(k,j).
A(n,k) = 2*A327084(n,k) - A327083(n,k) = A327083(n,k) - 2*A327085(n,k) = A327084(n,k) - A327085(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

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

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. A337408 (unoriented), A337409 (chiral), A337410 (achiral).
Rows 1-4 are A000027, A006528, A060530, A331358.
Cf. A327083 (simplex edges), A337411 (orthoplex edges), A325012 (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; 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

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) = A337408(n,k) + A337409(n,k) = 2*A337408(n,k) - A337410(n,k) = 2*A337409(n,k) + A337410(n,k).
Previous Showing 11-12 of 12 results.