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-10 of 10 results.

A063843 Number of n-multigraphs on 5 nodes.

Original entry on oeis.org

0, 1, 34, 792, 10688, 90005, 533358, 2437848, 9156288, 29522961, 84293770, 217993600, 519341472, 1154658869, 2420188694, 4821091920, 9187076352, 16837177281, 29809183410, 51172613512, 85448030080, 139159855989, 221554769150, 345523218536, 528767663040
Offset: 0

Views

Author

N. J. A. Sloane, Aug 25 2001

Keywords

Comments

Equivalently, number of ways to color edges of complete graph on 5 nodes with n colors, under action of symmetric group S_5, of order 120, with cycle index on edges given by (1/120)*(24*x5^2 + 30*x2*x4^2 + 20*x3^3*x1 + 20*x3*x6*x1 + 15*x1^2*x2^4 + 10*x1^4*x2^3 + x1^10). Setting all x_i = n gives the sequence.
Number of vertex colorings of the Petersen graph. Marko Riedel, Mar 24 2016
Number of unoriented colorings of the 10 triangular edges or triangular faces of a pentachoron, Schläfli symbol {3,3,3}, using n or fewer colors. Also called a 5-cell or 4-simplex. - Robert A. Russell, Oct 17 2020

Crossrefs

Cf. A063842. A row of A063841.
Cf. A331350 (oriented), A331352 (chiral), A331353 (achiral), A000389(n+4) (vertices and facets)
Other polychora: A331359 (8-cell), A331355 (16-cell), A338953 (24-cell), A338965 (120-cell, 600-cell).
Row 4 of A327084 (simplex edges and ridges) and A337884 (simplex faces and peaks).

Programs

  • Maple
    f:=n-> 1/120*(24*n^2+50*n^3+20*n^4+15*n^6+10*n^7+n^10);
  • Mathematica
    Table[(24n^2+50n^3+20n^4+15n^6+10n^7+n^10)/120,{n,0,30}] (* or *) LinearRecurrence[{11,-55,165,-330,462,-462,330,-165,55,-11,1},{0,1,34,792,10688,90005,533358,2437848,9156288,29522961,84293770},30] (* Harvey P. Dale, Oct 20 2012 *)
  • PARI
    a(n)=n^2*(n^8+10*n^5+15*n^4+20*n^2+50*n+24)/120 \\ Charles R Greathouse IV, Jan 20 2012

Formula

a(n) = (1/120)*(24*n^2+50*n^3+20*n^4+15*n^6+10*n^7+n^10).
a(n+1) = (1/5!)*(n^10 + 10*n^9 + 45*n^8 + 130*n^7 + 295*n^6 + 552*n^5 + 805*n^4 + 900*n^3 + 774*n^2 + 448*n + 120).
G.f. = (1 + 23*x + 473*x^2 + 3681*x^3 + 10717*x^4 + 11221*x^5 + 3779*x^6 + 339*x^7 + 6*x^8)/(1-x)^11. - M. F. Hasler, Jan 19 2012
a(0)=0, a(1)=1, a(2)=34, a(3)=792, a(4)=10688, a(5)=90005, a(6)=533358, a(7)=2437848, a(8)=9156288, a(9)=29522961, a(10)=84293770, a(n)= 11*a(n-1)- 55*a(n-2)+165*a(n-3)-330*a(n-4)+462*a(n-5)-462*a(n-6)+ 330*a(n-7)- 165*a(n-8)+55*a(n-9)-11*a(n-10)+a(n-11). - Harvey P. Dale, Oct 20 2012
From Robert A. Russell, Oct 17 2020: (Start)
a(n) = A331350(n) - A331352(n) = (A331350(n) + A331353(n)) / 2 = A331352(n) + A331353(n).
a(n) = 1*C(n,1) + 32*C(n,2) + 693*C(n,3) + 7720*C(n,4) + 44150*C(n,5) + 138312*C(n,6) + 247380*C(n,7) + 252000*C(n,8) + 136080*C(n,9) + 30240*C(n,10), where the coefficient of C(n,k) is the number of unoriented colorings using exactly k colors. (End)

Extensions

More terms from Vladeta Jovovic, Sep 02 2001

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

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

Examples

			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.
		

Crossrefs

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

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[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 *)
  • PARI
    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
    
  • Python
    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

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

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

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

A337883 Array read by descending antidiagonals: T(n,k) is the number of oriented 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, 40, 1, 5, 36, 1197, 3504, 1, 6, 75, 18592, 9753615, 13724608, 1, 7, 141, 166885, 3056311808, 19854224207910, 3574466244480, 1, 8, 245, 1019880, 264940140875, 468488921670219776, 25959704193068472575379, 106607224611810055168, 1
Offset: 2

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

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

Examples

			The table begins with T(2,1):
 1    2       3          4            5              6               7 ...
 1    5      15         36           75            141             245 ...
 1   40    1197      18592       166885        1019880         4738153 ...
 1 3504 9753615 3056311808 264940140875 10156268150064 221646915632373 ...
For T(3,4)=36, 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. A337884 (unoriented), A337885 (chiral), A337886 (achiral), A051168 (binary Lyndon words).
Other elements: A324999 (vertices), A327083 (edges).
Other polytopes: A337887 (orthotope), A337891 (orthoplex).
Rows 2-4 are A000027, A006008, A331350.

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[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) + A337885(n,k) = 2*A337884(n,k) - A337886(n,k) = 2*A337885(n,k) + A337886(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).

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

A337892 Array read by descending antidiagonals: T(n,k) is the number of unoriented 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, 22, 1, 4, 267, 11251322, 1, 5, 1996, 4825746875682, 314824532572147370464, 1, 6, 10375, 48038446526132256, 38491882660671134164965704408524083, 31716615393638864931753532641338560302264320, 1
Offset: 2

Views

Author

Robert A. Russell, Sep 28 2020

Keywords

Comments

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

Examples

			Array begins with T(2,1):
 1        2             3                 4                    5 ...
 1       22           267              1996                10375 ...
 1 11251322 4825746875682 48038446526132256 60632984344185045000 ...
		

Crossrefs

Cf. A337891 (oriented), A337893 (chiral), A337894 (achiral).
Other elements: A325005 (vertices), A337412 (edges).
Other polytopes: A337884 (simplex), A337888 (orthotope).
Rows 2-4 are A000027, A128766, A331359

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; 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,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) = A337891(n,k) - A337893(n,k) = (A337891(n,k) + A337894(n,k)) / 2 = A337893(n,k) + A337894(n,k).

A338114 Triangle read by rows: T(n,k) is the number of unoriented 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, 1, 1, 32, 693, 7720, 44150, 138312, 247380, 252000, 136080, 30240, 1, 2134, 4971504, 1513872568, 124978335900, 4307090369304, 78010256156784, 849590196841344, 6053725780061400, 29824685516682000, 105382759395846240, 273441179492268480
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 unoriented colorings are the same if they are congruent; chiral pairs are counted as one.
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    1
  1 32 693 7720 44150 138312 247380 252000 136080 30240
  ...
For T(3,2)=3, the tetrahedron has one, two, or three faces (vertices) of one color. For T(3,4)=1, each of the four tetrahedron faces (vertices) is a different color.
		

Crossrefs

Cf. A338113 (oriented), A338115 (chiral), A338116 (achiral), A337884 (k or fewer colors), A007318(n,k-1) (vertices and facets), A327088 (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[pc[#] j^Total[CX[#, m+1]][[2]] & /@ IntegerPartitions[n+1]]/(n+1)!]
    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

A337884(n,k) = Sum_{j=1..C(n+1,3)} T(n,j) * binomial(k,j).
T(n,k) = A338113(n,k) - A338115(n,k) = (A338113(n,k) + A338116(n,k)) / 2 = A338115(n,k) + A338116(n,k).
T(3,k) = A007318(3,k-1); T(4,k) = A327088(4,k).

A053426 Pure 2-complexes on n nodes with at most 3 multiple 2-simplexes.

Original entry on oeis.org

1, 1, 4, 35, 10688, 1533771392, 234249157811872000, 128777207192421762546956433408, 1031041719465866666816675423468078224724119552
Offset: 1

Views

Author

Vladeta Jovovic, Jan 11 2000

Keywords

Crossrefs

Showing 1-10 of 10 results.