A051168 -id:A051168 - cp's OEIS frontend

A050186 Triangular array T read by rows: T(h,k) = number of binary words of k 1's and h-k 0's which are not a juxtaposition of 2 or more identical subwords.

1, 1, 1, 0, 2, 0, 0, 3, 3, 0, 0, 4, 4, 4, 0, 0, 5, 10, 10, 5, 0, 0, 6, 12, 18, 12, 6, 0, 0, 7, 21, 35, 35, 21, 7, 0, 0, 8, 24, 56, 64, 56, 24, 8, 0, 0, 9, 36, 81, 126, 126, 81, 36, 9, 0, 0, 10, 40, 120, 200, 250, 200, 120, 40, 10, 0, 0, 11, 55, 165, 330, 462, 462, 330, 165, 55, 11

Offset: 0

Clark Kimberling

			For example, T(4,2) counts 1100,1001,0011,0110; T(2,1) counts 10, 01 (hence also counts 1010, 0101).
Rows:
  1;
  1,  1;
  0,  2,  0;
  0,  3,  3,  0;
  0,  4,  4,  4,  0;
  0,  5, 10, 10,  5,  0;

M. F. Hasler, A050186, rows 0..50, Sep 27 2018
M. F. Hasler, A050186, rows 0..100, Sep 27 2018
M. F. Hasler, A050186, rows 0..200, Sep 27 2018
N. J. A. Sloane, Transforms
Index entries for triangles and arrays related to Pascal's triangle

Same triangle as A053727 except this one includes column 0.

T(2n, n), T(2n+1, n) match A007727, A001700, respectively. Row sums match A027375.

T[n_, k_] := If[n == 0, 1, DivisorSum[GCD[k, n], MoebiusMu[#] Binomial[n/#, k/#]&]];
Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 16 2022 *)

A050186(n,k)=sumdiv(gcd(n+!n,k),d,moebius(d)*binomial(n/d,k/d)) \\ M. F. Hasler, Sep 27 2018

MOEBIUS transform of A007318 Pascal's Triangle.

If rows n > 1 are divided by n, this yields the triangle A051168, which equals A245558 surrounded by 0's (except for initial terms). This differs from A011847 from row n = 9 on. - M. F. Hasler, Sep 29 2018

A111492 Triangle read by rows: a(n,k) = (k-1)! * C(n,k).

Original entry on oeis.org

1, 2, 1, 3, 3, 2, 4, 6, 8, 6, 5, 10, 20, 30, 24, 6, 15, 40, 90, 144, 120, 7, 21, 70, 210, 504, 840, 720, 8, 28, 112, 420, 1344, 3360, 5760, 5040, 9, 36, 168, 756, 3024, 10080, 25920, 45360, 40320, 10, 45, 240, 1260, 6048, 25200, 86400, 226800, 403200, 362880

Offset: 1

Ross La Haye, Nov 15 2005

For k > 1, a(n,k) = the number of permutations of the symmetric group S_n that are pure k-cycles.
Reverse signed array is A238363. For a relation to (Cauchy-Euler) derivatives of the Vandermonde determinant, see Chervov link. - Tom Copeland, Apr 10 2014
Dividing the k-th column of T by (k-1)! for each column generates A135278 (the f-vectors, or face-vectors for the n-simplices). Then ignoring the first column gives A104712, so T acting on the column vector (-0,d,-d^2/2!,d^3/3!,...) gives the Euler classes for hypersurfaces of degree d in CP^n. Cf. A104712 and Dugger link therein. - Tom Copeland, Apr 11 2014
With initial i,j,n=1, given the n X n Vandermonde matrix V_n(x_1,...,x_n) with elements a(i=row,j=column)=(x_j)^(i-1), its determinant |V_n|, and the column vector of n ones C=(1,1,...,1), the n-th row of the lower triangular matrix T is given by the column vector determined by (1/|V_n|) * V_n(:x_1*d/dx_1:,...,:x_n*d/dx_n:)|V_n| * C, where :x_j*d/dx_j:^n = (x_j)^n*(d/dx_j)^n. - Tom Copeland, May 20 2014
For some other combinatorial interpretations of the first three columns of T, see A208535 and the link to necklace polynomials therein. Because of the simple relation of the array to the Pascal triangle, it can easily be related to many other arrays, e.g., T(p,k)/(p*(k-1)!) with p prime gives the prime rows of A185158 and A051168 when the non-integers are rounded to 0. - Tom Copeland, Oct 23 2014

			a(3,3) = 2 because (3-1)!C(3,3) = 2.
1;
2 1;
3 3 2;
4 6 8 6;
5 10 20 30 24;
6 15 40 90 144 120;
7 21 70 210 504 840 720;
8 28 112 420 1344 3360 5760 5040;
9 36 168 756 3024 10080 25920 45360 40320;

A. Chervov, A sum involving derivatives of Vandermonde

/* As triangle: */ [[Factorial(k-1)*Binomial(n,k): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Oct 21 2014

Flatten[Table[(k - 1)!Binomial[n, k], {n, 10}, {k, n}]]

a(n, k) = (k-1)!C(n, k) = P(n, k)/k.
E.g.f. (by columns) = exp(x)((x^k)/k).
a(n, 1) = A000027(n);
a(n, 2) = A000217(n-1);
a(n, 3) = A007290(n);
a(n, 4) = A033487(n-3).
a(n, n) = A000142(n-1);
a(n, n-1) = A001048(n-1) for n > 1.
Sum[a(n, k), {k, 1, n}] = A002104(n);
Sum[a(n, k), {k, 2, n}] = A006231(n).
a(n,k) = sum(j=k..n-1, j!/(j-k)!) (cf. Chervov link). - Tom Copeland, Apr 10 2014
From Tom Copeland, Apr 28 2014: (Start)
E.g.f. by row: [(1+t)^n-1]/t.
E.g.f. of row e.g.f.s: {exp[(1+t)*x]-exp(x)}/t.
O.g.f. of row e.g.f.s: {1/[1-(1+t)*x] - 1/(1-x)}/t.
E.g.f. of row o.g.f.s: -exp(x) * log(1-t*x). (End)

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

Robert A. Russell, Sep 28 2020

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.

			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.

E. M. Palmer and R. W. Robinson, Enumeration under two representations of the wreath product, Acta Math., 131 (1973), 123-143.

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.

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

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

Robert A. Russell, Sep 28 2020

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.

			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.

E. M. Palmer and R. W. Robinson, Enumeration under two representations of the wreath product, Acta Math., 131 (1973), 123-143.

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.

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

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

A011796 Number of irreducible alternating Euler sums of depth 6 and weight 6+2n.

Original entry on oeis.org

1, 3, 9, 20, 42, 75, 132, 212, 333, 497, 728, 1026, 1428, 1932, 2583, 3384, 4389, 5598, 7084, 8844, 10962, 13442, 16380, 19776, 23751, 28301, 33561, 39536, 46376, 54081, 62832, 72624, 83655, 95931, 109668, 124866, 141778, 160398

Offset: 1

N. J. A. Sloane, David Broadhurst

a(n-6) is the number of aperiodic necklaces (Lyndon words) with 6 black beads and n-6 white beads.

J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 147.

David J. Broadhurst, On the enumeration of irreducible k-fold Euler sums and their roles in knot theory and field theory, arXiv:hep-th/9604128, 1996.
David Broadhurst and Xavier Roulleau, Number of partitions of modular integers, arXiv:2502.19523 [math.NT], 2025. See pp. 3, 11, 19.
Index entries for sequences related to Lyndon words
Index entries for linear recurrences with constant coefficients, signature (2,1,-3,-1,1,4,-3,-3,4,1,-1,-3,1,2,-1).

Cf. A000031, A001037, A051168.

a:= n-> (Matrix([[42, 20, 9, 3, 1, 0$7, -1, -4, -9]]). Matrix(15, (i,j)-> if (i=j-1) then 1 elif j=1 then [2, 1, -3, -1, 1, 4, -3, -3, 4, 1, -1, -3, 1, 2, -1][i] else 0 fi)^(n-5))[1,1]: seq(a(n), n=1..50); # Alois P. Heinz, Aug 04 2008

a[n_] := Sum[Binomial[(n+6)/d, 6/d]*MoebiusMu[d],{d, Divisors[GCD[6, n+6]]}]/(n+6); Array[a, 40] (* Jean-François Alcover, Feb 02 2015 *)

G.f.: x*(1+x+2*x^2+2*x^3+3*x^4+2*x^6+x^7)/((1-x)^2*(1-x^2)^2*(1-x^3)*(1-x^6)).
G.f.: (1/(1-x)^6-1/(1-x^2)^3-1/(1-x^3)^2+1/(1-x^6))/6. - Herbert Kociemba, Oct 23 2016
a(n) = T(n,6), array T as in A051168.

A059861 a(n) = Product_{i=2..n} (prime(i) - 2).

Original entry on oeis.org

1, 1, 3, 15, 135, 1485, 22275, 378675, 7952175, 214708725, 6226553025, 217929355875, 8499244879125, 348469040044125, 15681106801985625, 799736446901266875, 45584977473372211875, 2689513670928960500625

Offset: 1

Labos Elemer, Feb 28 2001

nonn

Arises in Hardy-Littlewood k-tuple conjecture. Also a(n) is the exact number of d=2 and also d=4 differences in dRRS[modulus=n-th primorial]; see A049296 (dRRS[m]=set of first differences of reduced residue system modulo m).
For n>1 this is the determinant of the (n-1) X (n-1) matrix whose diagonal is A006093(n) = {1, 2, 4, 6, 10, 12, 16, 18..} = the first primes minus 1 and all other elements are 1's. The determinant begins: / (2-1) 1 1 1 1 1 1 ... / 1 (3-1) 1 1 1 1 1 ... / 1 1 (5-1) 1 1 1 1 ... / 1 1 1 (7-1) 1 1 1 ... / 1 1 1 1 (11-1) 1 1 ... / 1 1 1 1 1 (13-1) 1 ... - Alexander Adamchuk, May 21 2006
From Gary W. Adamson, Apr 21 2009: (Start)
Equals (-1)^n * (1, 1, 1, 3, 15, ...) dot (1, -2, 4, -6, 10, ...).
a(6) = 135 = (1, 1, 1, 3, 15) dot (1, -2, 4, -6, 10) = (1, -2, 4, -18, 150). (End)

			n=4, a(4) = 1*(3-2)*(5-2)*(7-2) = 15. 48 first terms of A049296 give one complete period of dRRS[210], in which 15 d=2, 15 d=4 and 18 larger differences occur. For n=1, 2, ..., 5 in the periods of length {1, 2, 8, 48, 480, ...} [see A005867] the number of d=2 and also d=4 differences is {1, 1, 3, 15, 135, ..}

Steven R. Finch, Mathematical Constants, Cambridge, 2003, pp. 84-94.
R. K. Guy, Unsolved Problems in Number Theory, Sections A8, A1.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979.
G. Polya, Mathematics and Plausible Reasoning, Vol. II, Appendix Princeton UP, 1954.

A.H.M. Smeets, Table of n, a(n) for n = 1..100
Steven Brown, Distance between consecutive elements of the multiplicative group of integers modulo n, arXiv:2311.06873 [math.NT], 2023. See Table 1 p. 25.
C. K. Caldwell, Prime k-tuple Conjecture
Steven R. Finch, Hardy-Littlewood Constants [Broken link]
Steven R. Finch, Hardy-Littlewood Constants [From the Wayback machine]
G. H. Hardy and J. E. Littlewood, Some problems of 'Partitio numerorum'; III: on the expression of a number as a sum of primes, Acta Mathematica, Vol. 44, pp. 1-70, 1923.
G. Niklasch, Some number theoretical constants: 1000-digit values [Cached copy]
G. Polya, Heuristic reasoning in the theory of numbers, Am. Math. Monthly, 66 (1959), 375-384.

Cf. A049296, A002110, A005867, A000847, A022008, A051160-A051168, A048298, A059861-A059865, A040976.

Cf. A067549, A006093.

Table[ Det[ DiagonalMatrix[ Table[ Prime[i-1] - 2, {i, 2, n} ] ] + 1 ], {n, 2, 20} ] (* Alexander Adamchuk, May 21 2006 *)
Table[Product[Prime@k - 2, {k, 2, n}], {n, 1, 18}] (* Harlan J. Brothers, Jul 02 2018 *)
a[1] = 1; a[n_] := a[n] = a[n - 1] (Prime[n] - 2);
Table[a[n], {n, 18}]  (* Harlan J. Brothers, Jul 02 2018 *)
Join[{1},FoldList[Times,Prime[Range[2,20]]-2]] (* Harvey P. Dale, Apr 19 2023 *)

a(n) = prod(i=2, n, prime(i)-2); \\ Michel Marcus, Apr 16 2017

a(n) = Det[ DiagonalMatrix[ Table[ Prime[i-1] - 2, {i, 2, n} ] ] + 1 ] for n>1. - Alexander Adamchuk, May 21 2006
a(n) = a(n-1) * (A000040(n) - 2) for n > 1. - A.H.M. Smeets, Dec 14 2019
a(n) = |{r | 0 <= r < primorial(n) and gcd(r, primorial(n)) = 1 and gcd(r + 2, primorial(n)) = 1}|. - Greg Tener, Oct 22 2021

Offset corrected by A.H.M. Smeets, Dec 14 2019

A245558 Square array read by antidiagonals: T(n,k) = number of n-tuples of nonnegative integers (u_0,...,u_{n-1}) satisfying Sum_{j=0..n-1} j*u_j == 1 mod n and Sum_{j=0..n-1} u_j = m.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 2, 3, 2, 1, 1, 3, 5, 5, 3, 1, 1, 3, 7, 8, 7, 3, 1, 1, 4, 9, 14, 14, 9, 4, 1, 1, 4, 12, 20, 25, 20, 12, 4, 1, 1, 5, 15, 30, 42, 42, 30, 15, 5, 1, 1, 5, 18, 40, 66, 75, 66, 40, 18, 5, 1, 1, 6, 22, 55, 99, 132, 132, 99, 55, 22, 6, 1

Offset: 1

N. J. A. Sloane, Aug 07 2014

The array is symmetric; for the entries on or below the diagonal see A245559.
If the congruence in the definition is changed from Sum_{j=0..n-1} j*u_j == 1 mod n  to Sum_{j=0..n-1} j*u_j == 0 mod n we get the array shown in A241926, A047996, and A037306.
Differs from A011847 from row n = 9, k = 4 on; if the rows are surrounded by 0's, this yields A051168 without its rows 0 and 1, i.e., a(1) is A051168(2,1). - M. F. Hasler, Sep 29 2018
This array was first studied by Fredman (1975). - Petros Hadjicostas, Jul 10 2019

			Square array begins:
  1, 1,  1,  1,   1,   1,    1,    1,    1,    1, ...
  1, 1,  2,  2,   3,   3,    4,    4,    5,    5, ...
  1, 2,  3,  5,   7,   9,   12,   15,   18,   22, ...
  1, 2,  5,  8,  14,  20,   30,   40,   55,   70, ...
  1, 3,  7, 14,  25,  42,   66,   99,  143,  200, ...
  1, 3,  9, 20,  42,  75,  132,  212,  333,  497, ...
  1, 4, 12, 30,  66, 132,  245,  429,  715, 1144, ...
  1, 4, 15, 40,  99, 212,  429,  800, 1430, 2424, ...
  1, 5, 18, 55, 143, 333,  715, 1430, 2700, 4862, ...
  1, 5, 22, 70, 200, 497, 1144, 2424, 4862, 9225, ...
  ...
Reading by antidiagonals, we get:
  1;
  1, 1;
  1, 1,  1;
  1, 2,  2,  1;
  1, 2,  3,  2,  1;
  1, 3,  5,  5,  3,   1;
  1, 3,  7,  8,  7,   3,   1;
  1, 4,  9, 14, 14,   9,   4,  1;
  1, 4, 12, 20, 25,  20,  12,  4,  1;
  1, 5, 15, 30, 42,  42,  30, 15,  5,  1;
  1, 5, 18, 40, 66,  75,  66, 40, 18,  5, 1;
  1, 6, 22, 55, 99, 132, 132, 99, 55, 22, 6, 1;
  ...

Taylor Brysiewicz, Necklaces count polynomial parametric osculants, arXiv:1807.03408 [math.AG], 2018.
A. Elashvili, M. Jibladze, Hermite reciprocity for the regular representations of cyclic groups, Indag. Math. (N.S.) 9 (1998), no. 2, 233-238. MR1691428 (2000c:13006).
A. Elashvili, M. Jibladze, D. Pataraia, Combinatorics of necklaces and "Hermite reciprocity", J. Algebraic Combin. 10 (1999), no. 2, 173-188. MR1719140 (2000j:05009). See p. 174.
M. L. Fredman, A symmetry relationship for a class of partitions, J. Combinatorial Theory Ser. A 18 (1975), 199-202.
I. M. Gessel and C. Reutenauer, Counting permutations with given cycle structure and descent set, J. Combin. Theory, Ser. A, 64, 1993, 189-215, Theorem 9.4.
J. E. Iglesias, Enumeration of closest-packings by the space group: a simple approach, Z. Krist. 221 (2006) 237-245, eq. (5).

This array is very similar to but different from A011847.
Cf. A051168, A092964, A241926, A047996, A037306, A245559.
Rows include A001840, A006918, A051170, A011796, A011797, A031164. Main diagonal is A022553.

# To produce the first 10 rows and columns (as on page 174 of the Elashvili et al. 1999 reference):
with(numtheory):
cnk:=(n,k) -> add(mobius(n/d)*d, d in divisors(gcd(n,k)));
anmk:=(n,m,k)->(1/(n+m))*add( cnk(d,k)*binomial((n+m)/d,n/d), d in divisors(gcd(n,m))); # anmk(n,m,k) is the value of a_k(n,m) as in Theorem 1, Equation (4), of the Elashvili et al. 1999 reference.
r2:=(n,k)->[seq(anmk(n,m,k),m=1..10)];
for n from 1 to 10 do lprint(r2(n,1)); od:

rows = 12;
cnk[n_, k_] := Sum[MoebiusMu[n/d] d, {d , Divisors[GCD[n, k]]}];
anmk[n_, m_, k_] := (1/(n+m)) Sum[cnk[d, k] Binomial[(n+m)/d, n/d], {d, Divisors[GCD[n, m]]}];
r2[n_, k_] := Table[anmk[n, m, k], {m, 1, rows}];
T = Table[r2[n, 1], {n, 1, rows}];
Table[T[[n-k+1, k]], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Nov 05 2018, from Maple *)

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

Robert A. Russell, Sep 28 2020

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.

			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.

E. M. Palmer and R. W. Robinson, Enumeration under two representations of the wreath product, Acta Math., 131 (1973), 123-143.

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.

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

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

Robert A. Russell, Sep 28 2020

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.

			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.

E. M. Palmer and R. W. Robinson, Enumeration under two representations of the wreath product, Acta Math., 131 (1973), 123-143.

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.

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

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

A011795 a(n) = floor(C(n,4)/5).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 3, 7, 14, 25, 42, 66, 99, 143, 200, 273, 364, 476, 612, 775, 969, 1197, 1463, 1771, 2125, 2530, 2990, 3510, 4095, 4750, 5481, 6293, 7192, 8184, 9275, 10472, 11781, 13209, 14763, 16450, 18278, 20254, 22386, 24682, 27150, 29799, 32637, 35673, 38916, 42375, 46060, 49980, 54145, 58565, 63250, 68211, 73458, 79002, 84854, 91025, 97527, 104371

Offset: 0

N. J. A. Sloane, David Broadhurst

a(n-1) = number of aperiodic necklaces (Lyndon words) with 5 black beads and n-5 white beads.

J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 147.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
David J. Broadhurst, On the enumeration of irreducible k-fold Euler sums and their roles in knot theory and field theory, arXiv:hep-th/9604128, 1996.
David Broadhurst and Xavier Roulleau, Number of partitions of modular integers, arXiv:2502.19523 [math.NT], 2025. See pp. 3, 11, 19.
Index entries for sequences related to Lyndon words
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1,1,-4,6,-4,1).

Cf. A000031, A001037, A011915, A051168.
Same as A051170(n+1).
A column of triangle A011847.

[Floor(Binomial(n+1,5)/(n+1)): n in [0..70]]; // Vincenzo Librandi Jun 19 2012

seq(floor(binomial(n,4)/5), n=0.. 70); # Zerinvary Lajos, Jan 12 2009

CoefficientList[Series[x^5(1+x^3)/((1-x)^3(1-x^2)(1-x^5)),{x,0,70}],x] (* Vincenzo Librandi, Jun 19 2012 *)
CoefficientList[Series[x^4/5 (1/(1-x)^5-1/(1- x^5)),{x,0,70}],x] (* Herbert Kociemba, Oct 16 2016 *)

a(n)=binomial(n,4)\5 \\ Charles R Greathouse IV, Oct 07 2015

[binomial(n,4)//5 for n in range(71)] # G. C. Greubel, Oct 20 2024

G.f.: x^5*(1+x^3)/((1-x)^3*(1-x^2)*(1-x^5)) = x^5*(1-x+x^2)/((1-x)^5*(1+x+x^2+x^3+x^4)).

a(n) = floor(binomial(n+1,5)/(n+1)). - Gary Detlefs, Nov 23 2011

cp's OEIS Frontend

A050186 Triangular array T read by rows: T(h,k) = number of binary words of k 1's and h-k 0's which are not a juxtaposition of 2 or more identical subwords.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A111492 Triangle read by rows: a(n,k) = (k-1)! * C(n,k).

Views

Author

Keywords

Comments

Examples

Links

Programs

Magma

Mathematica

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A011796 Number of irreducible alternating Euler sums of depth 6 and weight 6+2n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A059861 a(n) = Product_{i=2..n} (prime(i) - 2).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A245558 Square array read by antidiagonals: T(n,k) = number of n-tuples of nonnegative integers (u_0,...,u_{n-1}) satisfying Sum_{j=0..n-1} j*u_j == 1 mod n and Sum_{j=0..n-1} u_j = m.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs