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

A213940 Triangle with entry a(n,m) giving the number of bracelets of n beads (dihedral D_n symmetry) with n colors available for each bead, but only m distinct fixed colors, say c[1],...,c[m], are present, with m from {1,...,n} and n>=1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 3, 2, 3, 1, 3, 6, 6, 12, 1, 7, 20, 26, 30, 60, 1, 8, 40, 93, 150, 180, 360, 1, 18, 106, 424, 633, 1050, 1260, 2520, 1, 22, 304, 1180, 3260, 5040, 8400, 10080, 20160, 1, 46, 731, 4844, 16212, 29244, 45360, 75600, 90720, 181440
Offset: 1

Views

Author

Wolfdieter Lang, Jul 20 2012

Keywords

Comments

This triangle is obtained from the partition array A213939 by summing in row n, for n>=1, all entries related to partitions of n with the same number of parts m.
a(n,m) is the number of bracelets of n beads (dihedral D_n symmetry) corresponding to the representative color multinomials obtained from all partitions of n with m parts by 'exponentiation', hence only m from the available n colors are present. As a representative multinomial of each of the p(n,m)=A008284(n,m) such m-color classes we take the one where the considered m part partition of n, [p[1],...,p[m]], written in nonincreasing order, is distributed as exponents on the color indices like c[1]^p[1]*...*c[m]^p[m]. That is only the first m colors from the n available ones are involved.
See the comments on A212359 for the Abramowitz-Stegun (A-St) order of partitions, and the 'exponentiation' to obtain multisets, used to encode color multinomials, from partitions.
The row sums of this triangle coincide with the ones of array A213939, and they are given by A213943.
Number of n-length bracelets w over a k-ary alphabet {a1,a2,...,ak} such that #(w,a1) >= #(w,a2) >= ... >= #(w,ak) >= 1, where #(w,x) counts the letters x in word w (bracelet analog of A226874). - Andrew Howroyd, Sep 26 2017

Examples

			n\m  1  2   3    4     5     6     7     8     9     10 ...
1    1
2    1  1
3    1  1   1
4    1  3   2    3
5    1  3   6    6    12
6    1  7  20   26    30    60
7    1  8  40   93   150   180   360
8    1 18 106  424   633  1050  1260  2520
9    1 22 304 1180  3260  5040  8400 10080 20160
10   1 46 731 4844 16212 29244 45360 75600 90720 181440
...
a(5,3) = 2 + 4 = 6, from A213939(5,4) + A213939(5,5), because k(5,3,1) = 4 and p(5,3) = 2.
a(2,1) = 1 because the partition [2] of n=2 with part number m=1 corresponds to the representative color multinomial (here monomial) c[1]^2 = c[1]*c[1], and there is one such representative bracelet. There is another bracelet color monomial in this class of n=2 colors where only m=1 color is active: c[2]*c[2]. See the triangle entry A213941(2,1)=2. The same holds for the necklace case.
a(3,1) = 1 from the color monomial representative c[1]^3. This class has 2 other members: c[2]^3 and c[3]^3. See A213941(3,1)=3. The same holds for the necklace case.
Like in the necklace case one has in general a(n,1)=1 and A213941(n,1) = n from the partition [n] providing the color signature and a representative c[1]^n.
a(3,2) = 1 from the representative color multinomial c[1]^2*c[2] (from the m=2 partition [2,1] of n=3) leading to just one representative bracelet (and necklace) cyclic(112) (when one uses j for color c[j]). The whole class consists of A213941(3,2)=6 bracelets (or necklaces): cyclic(112), cyclic(113), cyclic(221), cyclic(223), cyclic(331) and cyclic(332).
a(3,3) = 1. The representative color multinomial is c[1]*c[2]*c[3] (from the m=3 partition [1,1,1]). There is only one bracelet cyclic(1,2,3) which constitutes already the whole class (A213941(3,3)=1). The necklace cyclic(1,3,2) becomes equivalent under D_3.
a(4,2) = 3 from two representative color multinomials c[1]^3*c[2] and c[1]^2*c[2]^2 (from the two m=2 partitions of n=4: [3,1] and [2,2]). The first one has one representative bracelet, namely cyclic(1112), the second one leads to the two representative bracelets: cyclic(1122) and cyclic(1212). Together these are the 3 bracelets counted by a(4,2). The first color class c[.]^3*c[.] consists of 4*3=12 bracelets, when all 4 colors are used. The second one consists of 2*6=12 bracelets. Together they sum up to the 24 bracelets counted by A213941(4,2). In this example the necklace case does not differ from the bracelet one.
		

Crossrefs

Columns k=2..5 are A213942, A214307, A214309, A214311.
Cf. A213934 (cyclic symmetry).

Programs

  • PARI
    Cyc(v)={my(s=vecsum(v)); sumdiv(gcd(v), d, eulerphi(d)*(s/d)!/prod(i=1, #v, (v[i]/d)!))/s}
    CPal(v)={my(odds=#select(t->t%2,v), s=vecsum(v));  if(odds>2, 0, ((s-odds)/2)!/prod(i=1, #v, (v[i]\2)!))}
    T(n,k)={my(t=0); forpart(p=n, t+=Cyc(Vec(p))+CPal(Vec(p)), [1,n], [k,k]); t/2}
    for(n=1, 10, for(k=1,n, print1(T(n,k), ", ")); print); \\ Andrew Howroyd, Sep 26 2017
    
  • PARI
    \\ faster version; here U is A226874 as vector of polynomials.
    U(n)={Vec(serlaplace(prod(k=1, n, 1/(1-y*x^k/k!) + O(x*x^n))))}
    T(n)={my(t=U(n)); vector(n, n, vector(n, k, ((1/n)*sumdiv(n, d, eulerphi(n/d) * polcoeff(t[d+1], k)) + if(n%2, sum(d=0, (n-1)/2, binomial((n-1)/2, d)*polcoeff(t[d+1], (k-1))), polcoeff(t[n/2+1], k) + sum(d=0, n/2-1, binomial(n/2-1, d)*(2^d + if(d%2, 0, binomial(d, d/2)))*polcoeff(t[n/2-d], k-2))/2))/2))}
    { my(t=T(10)); for(n=1, #t, print(t[n])) } \\ Andrew Howroyd, Dec 22 2017

Formula

a(n,m) = Sum_{j=1..p(n,m)}A213939(n,k(n,m,1)+j-1), with k(n,m,1) the position where in the list of partitions of n in A-St order the first with m parts appears, and p(n,m) is the number of partitions of n with m parts shown in the array A008284. E.g., n=5, m=3: k(5,3,1)=4, p(5,3)=2.

A056358 Number of bracelet structures using exactly three different colored beads.

Original entry on oeis.org

0, 0, 1, 2, 5, 14, 31, 82, 202, 538, 1401, 3838, 10395, 28890, 80207, 225368, 634265, 1796648, 5100325, 14535298, 41513434, 118880650, 341094843, 980665898, 2824223495, 8146908210, 23535345372, 68084937912, 197211483155, 571915789978, 1660402195255, 4825554617686
Offset: 1

Views

Author

Keywords

Comments

Turning over will not create a new bracelet. Permuting the colors of the beads will not change the structure.

References

  • M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]

Crossrefs

Column 3 of A152176.

Formula

a(n) = A056353(n) - A000011(n).

Extensions

Terms a(28) and beyond from Andrew Howroyd, Oct 24 2019

A213942 a(n) is the number of representative two-color bracelets (necklaces with turnover allowed) with n beads for n >= 2.

Original entry on oeis.org

1, 1, 3, 3, 7, 8, 18, 22, 46, 62, 136, 189, 409, 611, 1344, 2055, 4535, 7154, 15881, 25481, 56533, 92204, 204759, 337593, 748665, 1246862, 2762111, 4636389, 10253938, 17334800, 38278784, 65108061, 143534770, 245492243, 540353057, 928772649, 2041154125
Offset: 2

Views

Author

Wolfdieter Lang, Jul 31 2012

Keywords

Comments

This is the second column (m=2) of triangle A213940.
The relevant floor(n/2) representative color multinomials are c[1]^(n-1)*c[2], c[1]^(n-2)*c[2]^2, ..., c[1]^(n-floor(n/2))* c[2]^(floor(n/2)). For such representative bracelets the color c[1] is therefore preferred. Only for even n can c[2] appear as often as c[1], namely, n/2 times.
Note that beads with different colors are always present. This is in contrast to, e.g., A000029, where not only representatives but also one-color bracelets are counted. This sequences gives the number of binary bracelets with at least as many 0's as 1's and at least one 1 (bracelet analog of A226881). The number of two-color bracelets up to permutations of colors is given by A056357. For odd n these two sequences are equal. For a(8), the bracelets 00011011 and 11100100 are equivalent in A056357 but distinct in this sequence. - Andrew Howroyd and Wolfdieter Lang, Sep 25 2017

Examples

			a(5) = A213939(5,2) + A213939(5,3) = 1 + 2 = 3 from the representative bracelets (with colors j for c[j], j=1,2) cyclic(11112), cyclic(11122) and cyclic(11212). The first one has color signature (exponents) [4,1] and the two others have signature [3,2]. For the number of all two-color 5-bracelets with beads of five colors available see A214308(5) = 60.
a(8) = 18 =  1 + 4 + 5 + 8 for the partitions of 8 with 2 parts (7,1), (6, 2), (5,3), (4,4), respectively. see A213939(5, k), k = 2..5). The 8 representative bracelets for the exponents (signature) from partition (4,4) are B1 = (11112222), B2 = (11121222), B3 = (11212122), B4 = (11212212), B5 = (11221122), B6 = (12121212), B7 = (11122122) and B8 = (11211222). B1 to B6 are color exchange (1 <-> 2) invariant (modulo D_8 symmetry, i.e., cyclic or anti-cyclic operations). B7 is equivalent to B8 under color exchange.
This explains why A056357(8) = 17. The difference between the present sequence and A056357 is that there, besides D_n symmetry, also color exchange is allowed. Here only color exchange compatible with D_n symmetry is allowed. - _Wolfdieter Lang_, Sep 28 2017
		

Crossrefs

Cf. A213939, A213940, A214307 (m=3), A214308 (m=2, all bracelets).

Programs

  • Mathematica
    a29[n_] := (1/4)*(Mod[n, 2] + 3)*2^Quotient[n, 2] + DivisorSum[n, EulerPhi[#]*2^(n/#)&]/(2*n);
    a5648[n_] := 1/2*(Binomial[2*Quotient[n, 2], Quotient[n, 2]] + DivisorSum[n, EulerPhi[#]*Binomial[2*n/#, n/#]&]/(2*n));
    a[n_] := a29[n]/2 - 1 + If[EvenQ[n], a5648[n/2]/2, 0];
    Array[a, 37, 2] (* Jean-François Alcover, Nov 05 2017, after Andrew Howroyd *)

Formula

a(n) = A213940(n,2), n >= 2.
a(n) = Sum_{k=2..A008284(n,2)+1} A213939(n,k), n >= 2, with A008284(n,2) = floor(n/2).
a(2n) = (A000029(2n) + A005648(n)) / 2 - 1, a(2n+1) = A000029(2n+1) / 2 - 1. - Andrew Howroyd, Sep 25 2017

Extensions

Terms a(26) and beyond from Andrew Howroyd, Sep 25 2017

A214310 a(n) is the number of all three-color bracelets (necklaces with turning over allowed) with n beads and the three colors are from a repertoire of n distinct colors, for n >= 3.

Original entry on oeis.org

1, 24, 180, 1120, 5145, 23016, 91056, 357480, 1327095, 4893680, 17525508, 62254920, 217457695, 753332160, 2581110000, 8779264032, 29624681763, 99350001360, 331159123260, 1098168382080, 3624003213369, 11908069219816, 38972450763000, 127087400895000
Offset: 3

Views

Author

Wolfdieter Lang, Jul 31 2012

Keywords

Comments

This is the third column (m=3) of triangle A214306.
Each 3 part partition of n, with the parts written in nonincreasing order, defines a color signature. For a given color signature, say [p[1], p[2], p[3]], with p[1] >= p[2] >= p[3] >= 1, there are A213941(n,k)= A035206(n,k)* A213939(n,k) bracelets if this signature corresponds (with the order of the parts reversed) to the k-th partition of n in Abramowitz-Stegun (A-St) order. See A213941 for more details. Here all p(n,3)= A008284(n,3) partitions of n with 3 parts are considered. The color repertoire for a bracelet with n beads is [c[1], ..., c[n]].
Compare this with A027671 where also single color bracelets are included, and the color repertoire is only [c[1], c[2], c[3]] for all n.

Examples

			a(5) = A213941(5,4) + A213941(5,5) = 60 + 120 = 180 from the bracelet (with colors j for c[j], j=1, 2, ..., 5) 11123 and 11213, both taken cyclically, each representing a class of order A035206(5,4)= 30 (if all 5 colors are used), and 11223, 11232, 12123 and 12213, all taken cyclically, each representing a class of order A035206(5,5)= 30. For example, cyclic(11322) becomes equivalent to cyclic(11223) by turning over or reflection. The multiplicity A035206 depends only on the color signature.
		

Crossrefs

Cf. A213941, A214306, A214307 (m=3, representative bracelets), A214312 (m=4).

Formula

a(n) = A214306(n,3), n >= 3.
a(n) = sum(A213941(n,k), k = A214314(n,3).. (A214314(n,3) - 1 + A008284(n,3))), n >= 3.
a(n) = binomial(n,3) * A056343(n). - Andrew Howroyd, Mar 25 2017

Extensions

a(26) from Andrew Howroyd, Mar 25 2017
Showing 1-4 of 4 results.