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.

Previous Showing 11-13 of 13 results.

A330341 Triangle read by rows: T(n,k) is the number of n-bead bracelets using exactly k colors with no adjacent beads having the same color.

Original entry on oeis.org

0, 0, 1, 0, 0, 1, 0, 1, 3, 3, 0, 0, 3, 12, 12, 0, 1, 10, 46, 90, 60, 0, 0, 9, 120, 480, 720, 360, 0, 1, 27, 384, 2235, 5670, 6300, 2520, 0, 0, 29, 980, 9380, 36960, 68880, 60480, 20160, 0, 1, 75, 2904, 38484, 217152, 604800, 876960, 635040, 181440
Offset: 1

Views

Author

Andrew Howroyd, Dec 20 2019

Keywords

Comments

In the case of n = 1, the single bead is considered to be cyclically adjacent to itself giving T(1,1) = 0. If compatibility with A208544 is wanted then T(1,1) should be 1.

Examples

			Triangle begins:
  0;
  0, 1;
  0, 0,  1;
  0, 1,  3,   3;
  0, 0,  3,  12,   12;
  0, 1, 10,  46,   90,    60;
  0, 0,  9, 120,  480,   720,   360;
  0, 1, 27, 384, 2235,  5670,  6300,  2520;
  0, 0, 29, 980, 9380, 36960, 68880, 60480, 20160;
  ...
		

Crossrefs

Column 3 is A330632.
Row sums are A330621.

Programs

  • PARI
    \\ here U(n, k) is A208544(n, k) for n > 1.
    U(n, k) = (sumdiv(n, d, eulerphi(n/d)*(k-1)^d)/n + if(n%2, 1-k, k*(k-1)^(n/2)/2))/2;
    T(n, k)={sum(j=1, k, (-1)^(k-j)*binomial(k, j)*U(n, j))}

Formula

T(n,k) = Sum_{j=1..k} (-1)^(k-j)*binomial(k,j)*A208544(n,j) for n > 1.

A056346 Number of bracelets of length n using exactly six different colored beads.

Original entry on oeis.org

0, 0, 0, 0, 0, 60, 1080, 11970, 105840, 821952, 5874480, 39713550, 258136200, 1631273220, 10096734312, 61536377700, 370710950400, 2213749658880, 13132080672480, 77509456944318, 455754569692680
Offset: 1

Views

Author

Keywords

Comments

Turning over will not create a new bracelet.

Examples

			For a(6)=60, pair up the 120 permutations of BCDEF, each with its reverse, such as BCDEF-FEDCB.  Precede the first of each pair with an A, such as ABCDEF.  These are the 60 arrangements, all chiral.  If we precede the second of each pair with an A, such as AFEDCB, we get the chiral partner of each. - _Robert A. Russell_, Sep 27 2018
		

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 6 of A273891.
Equals (A056286 + A056492) / 2 = A056286 - A305545 = A305545 + A056492.
Cf. A008277.

Programs

  • Mathematica
    t[n_, k_] := (For[t1 = 0; d = 1, d <= n, d++, If[Mod[n, d] == 0, t1 = t1 + EulerPhi[d]*k^(n/d)]]; If[EvenQ[n], (t1 + (n/2)*(1 + k)*k^(n/2))/(2*n), (t1 + n*k^((n + 1)/2))/(2*n)]);
    T[n_, k_] := Sum[(-1)^i*Binomial[k, i]*t[n, k - i], {i, 0, k - 1}];
    a[n_] := T[n, 6];
    Array[a, 21] (* Jean-François Alcover, Nov 05 2017, after Andrew Howroyd *)
    k=6; Table[k! DivisorSum[n, EulerPhi[#] StirlingS2[n/#,k]&]/(2n) + k!(StirlingS2[Floor[(n+1)/2], k] + StirlingS2[Ceiling[(n+1)/2], k])/4, {n,1,30}] (* Robert A. Russell, Sep 27 2018 *)
  • PARI
    a(n) = my(k=6); (k!/4) * (stirling(floor((n+1)/2),k,2) + stirling(ceil((n+1)/2),k,2)) + (k!/(2*n))*sumdiv(n, d, eulerphi(d)*stirling(n/d,k,2)); \\ Michel Marcus, Sep 29 2018

Formula

a(n) = A056341(n) - 6*A032276(n) + 15*A032275(n) - 20*A027671(n) + 15*A000029(n) - 6.
From Robert A. Russell, Sep 27 2018: (Start)
a(n) = (k!/4) * (S2(floor((n+1)/2),k) + S2(ceiling((n+1)/2),k)) + (k!/2n) * Sum_{d|n} phi(d) * S2(n/d,k), where k=6 is the number of colors and S2 is the Stirling subset number A008277.
G.f.: (k!/4) * x^(2k-2) * (1+x)^2 / Product_{i=1..k} (1-i x^2) - Sum_{d>0} (phi(d)/2d) * Sum_{j} (-1)^(k-j) * C(k,j) * log(1-j x^d), where k=6 is the number of colors. (End)

A019537 Number of special orbits for dihedral group of degree n.

Original entry on oeis.org

1, 2, 4, 14, 61, 414, 3416, 34274, 394009, 5113712, 73758368, 1170495180, 20263806277, 380048113202, 7676106638884, 166114210737254, 3834434327929981, 94042629562443206, 2442147034770292496, 66942194906543381336, 1931543452346146410965, 58519191359170883258606
Offset: 1

Views

Author

Manfred Goebel (goebel(AT)informatik.uni-tuebingen.de)

Keywords

Comments

a(n) is the number of ways to color a necklace of n beads using at most n colors. Turning the necklace over does not count as different. - Robert A. Russell, May 31 2018

Examples

			For a(3) = 4, the necklaces are AAA, AAB, ABB, and ABC. Last one is chiral. For a(4) = 14, the necklacess are AAAA, AAAB, AABB, ABAB, ABBB, ABAC, ABCB, ACBC, AABC, ABBC, ABCC, ABCD, ABDC, and ACBD. Last six are chiral. - _Robert A. Russell_, May 31 2018
		

Crossrefs

Cf. A019536.
Row sums of A273891.

Programs

  • Mathematica
    Table[Sum[(k!/(2n)) DivisorSum[n, EulerPhi[#] StirlingS2[n/#,k] &] + (k!/4) (StirlingS2[Floor[(n+1)/2],k] + StirlingS2[Ceiling[(n+1)/2],k]), {k, 1, n}], {n, 1, 40}] (* Robert A. Russell, May 31 2018 *)
  • PARI
    a(n) = sum(k=1, n, (k!/4)*(stirling(floor((n+1)/2),k,2) + stirling(ceil((n+1)/2),k,2)) + (k!/(2*n))*sumdiv(n, d, eulerphi(d)*stirling(n/d,k,2))); \\ Michel Marcus, Jun 06 2018

Formula

a(n) = Sum_{k=1..n} ((k!/4)*(S2(floor((n+1)/2),k) + S2(ceiling((n+1)/2),k)) + (k!/(2 n))*Sum_{d|n} phi(d)*S2(n/d,k)), where S2(n,k) is the Stirling subset number A008277. - Robert A. Russell, May 31 2018
a(n) ~ (n-1)! / (4 * log(2)^(n+1)). - Vaclav Kotesovec, Jul 21 2019

Extensions

More terms (using A273891) from Alois P. Heinz, Jun 02 2016
Previous Showing 11-13 of 13 results.