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.
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
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.
Links
- G. Royle, Partitions and Permutations
Crossrefs
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
Comments