A308401 Number of bracelets (turnover necklaces) of length n that have no reflection symmetry and consist of 6 white beads and n-6 black beads.
3, 6, 16, 30, 56, 91, 150, 224, 336, 477, 672, 912, 1233, 1617, 2112, 2700, 3432, 4290, 5340, 6552, 8008, 9678, 11648, 13888, 16503, 19448, 22848, 26658, 31008, 35853, 41346, 47424, 54264, 61803, 70224, 79464, 89733, 100947, 113344, 126840, 141680, 157780, 175416, 194480, 215280, 237708
Offset: 9
Examples
Using Frank Ruskey's website (listed above) to generate bracelets of fixed content (6, 3) with string length n = 9 and alphabet size 2, we get the following A005513(n = 9) = 7 bracelets: (1) WWWWWWBBB, (2) WWWWWBWBB, (3) WWWWBWWBB, (4) WWWWBWBWB, (5) WWWBWWWBB, (6) WWWBWWBWB, and (7) WWBWWBWWB. From these, bracelets 1, 4, 5, and 7 have reflection symmetry, while bracelets 2, 3 and 6 have no reflection symmetry (and thus, a(9) = 3). Starting with a black bead, we count that bead and how many white beads follow (in one direction), and continue this process until we count all beads around the circle. We thus use MacMahon's correspondence to get the following dihedral compositions of n = 9 into 3 parts: (1) 1 + 7 + 1, (2) 1 + 2 + 6, (3) 1 + 3 + 5, (4) 2 + 5 + 2, (5) 4 + 1 + 4, (6) 2 + 3 + 4, and (7) 3 + 3 + 3. Again, dihedral compositions 1, 4, 5, and 7 are symmetric (have reflection symmetry), while dihedral compositions 2, 3, and 6 are not symmetric (and thus, a(9) = 3). We may also start with a white bead and count that bead and how many black beads follow (in one direction), and continue this process until we count all beads around the circle. We thus use MacMahon's correspondence again to get the following (conjugate) dihedral compositions of n = 9 into 6 parts: (1) 1 + 1 + 1 + 1 + 1 + 4, (2) 1 + 1 + 1 + 1 + 2 + 3, (3) 1 + 1 + 1 + 2 + 1 + 3, (4) 1 + 1 + 1 + 2 + 2 + 2, (5) 1 + 1 + 2 + 1 + 1 + 3, (6) 1 + 1 + 2 + 1 + 2 + 2, and (7) 1 + 2 + 1 + 2 + 1 + 2. Again, dihedral compositions 1, 4, 5, and 7 have reflection symmetries, while dihedral compositions 2, 3, and 6 do not have reflection symmetries (and thus, a(9) = 3). For example, dihedral composition 1 is symmetric because we can draw an axis of symmetry through one of the 1s and 4. In addition, dihedral composition 5 is symmetric because we may draw an axis of symmetry through the numbers 2 and 3.
Links
- Colin Barker, Table of n, a(n) for n = 9..1000
- Hansraj Gupta, Enumeration of incongruent cyclic k-gons, Indian J. Pure and Appl. Math., 10 (1979), no. 8, 964-999.
- Petros Hadjicostas, The aperiodic version of Herbert Kociemba's formula for bracelets with no reflection symmetry, 2019.
- Arnold Knopfmacher and Neville Robbins, Some properties of dihedral compositions, Util. Math. 92 (2013), 207-220.
- Richard H. Reis, A formula for C(T) in Gupta's paper, Indian J. Pure and Appl. Math., 10 (1979), no. 8, 1000-1001.
- Frank Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc.
- Vladimir S. Shevelev, Necklaces and convex k-gons, Indian J. Pure and Appl. Math., 35 (2004), no. 5, 629-638.
- Vladimir S. Shevelev, Necklaces and convex k-gons, Indian J. Pure and Appl. Math., 35 (2004), no. 5, 629-638.
- Duncan M. Y. Sommerville, On certain periodic properties of cyclic compositions of numbers, Proc. London Math. Soc. S2-7(1) (1909), 263-313.
- Index entries for linear recurrences with constant coefficients, signature (2,1,-3,-1,1,4,-3,-3,4,1,-1,-3,1,2,-1).
Crossrefs
Programs
-
PARI
a(n) = (1/12)* (sumdiv(gcd(n, 6), d, eulerphi(d)*binomial((n/d) - 1, (6/d) - 1))) - (1/2)*binomial(floor(n/2), 3); \\ Michel Marcus, May 28 2019
-
PARI
Vec(x^9*(3 + x^2 + x^3 + x^4) / ((1 - x)^6*(1 + x)^3*(1 - x + x^2)*(1 + x + x^2)^2) + O(x^50)) \\ Colin Barker, Jun 02 2019
Formula
G.f.: (x^k/2) * (-(1 + x)/(1 - x^2)^floor((k/2) + 1) + (1/k) * Sum_{m|k} phi(m)/(1 - x^m)^(k/m)) with k = 6. (This formula is due to Herbert Kociemba.)
a(n) = -(1/2)*binomial(floor(n/2), 3) + (1/12)* Sum_{d|gcd(n, 6)} phi(d)*binomial((n/d) - 1, (6/d) - 1) for n >= 6. (This is a modification of formulas found in Gupta (1979) and Shevelev (2004).)
From Colin Barker, May 26 2019: (Start)
G.f.: x^9*(3 + x^2 + x^3 + x^4) / ((1 - x)^6*(1 + x)^3*(1 - x + x^2)*(1 + x + x^2)^2).
a(n) = 2*a(n-1) + a(n-2) - 3*a(n-3) - a(n-4) + a(n-5) + 4*a(n-6) - 3*a(n-7) - 3*a(n-8) + 4*a(n-9) + a(n-10) - a(n-11) - 3*a(n-12) + a(n-13) + 2*a(n-14) - a(n-15) for n > 23. (End)
Comments