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.

A305540 Triangle read by rows: T(n,k) is the number of achiral loops (necklaces or bracelets) of length n using exactly k different colors.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 4, 3, 1, 6, 6, 1, 10, 21, 12, 1, 14, 36, 24, 1, 22, 93, 132, 60, 1, 30, 150, 240, 120, 1, 46, 345, 900, 960, 360, 1, 62, 540, 1560, 1800, 720, 1, 94, 1173, 4980, 9300, 7920, 2520, 1, 126, 1806, 8400, 16800, 15120, 5040, 1, 190, 3801, 24612, 71400, 103320, 73080, 20160, 1, 254, 5796, 40824, 126000, 191520, 141120, 40320
Offset: 1

Views

Author

Robert A. Russell, Jun 04 2018

Keywords

Comments

The number of achiral necklaces is equivalent to the number of achiral bracelets.

Examples

			The triangle begins with T(1,1):
1;
1,   1;
1,   2;
1,   4,     3;
1,   6,     6;
1,  10,    21,     12;
1,  14,    36,     24;
1,  22,    93,    132,     60;
1,  30,   150,    240,    120;
1,  46,   345,    900,    960,     360;
1,  62,   540,   1560,   1800,     720;
1,  94,  1173,   4980,   9300,    7920,    2520;
1, 126,  1806,   8400,  16800,   15120,    5040;
1, 190,  3801,  24612,  71400,  103320,   73080,   20160;
1, 254,  5796,  40824, 126000,  191520,  141120,   40320;
1, 382, 11973, 113652, 480060, 1048320, 1234800,  745920, 181440;
1, 510, 18150, 186480, 834120, 1905120, 2328480, 1451520, 362880;
For a(4,2)=4, the achiral loops are AAAB, AABB, ABAB, and ABBB.

Crossrefs

Odd rows are A019538.
Even rows are A172106.
Columns 1-6 are A057427, A027383, A056489, A056490, A056491, and A056492.

Programs

  • Mathematica
    Table[(k!/2) (StirlingS2[Floor[(n + 1)/2], k] + StirlingS2[Ceiling[(n + 1)/2], k]), {n, 1, 15}, {k, 1, Ceiling[(n + 1)/2]}] // Flatten
  • PARI
    T(n, k) = (k!/2)*(stirling(floor((n+1)/2), k, 2)+stirling(ceil((n+1)/2), k, 2));
tabf(nn) = for(n=1, nn, for (k=1, ceil((n+1)/2), print1(T(n, k), ", ")); print); \\ Michel Marcus, Jul 02 2018

Formula

T(n,k) = (k!/2) * (S2(floor((n+1)/2),k) + S2(ceiling((n+1)/2),k)), where S2(n,k) is the Stirling subset number A008277.
T(n,k) = 2*A273891(n,k) - A087854(n,k).
G.f. for column k>1: (k!/2) * x^(2k-2) * (1+x)^2 / Product_{i=1..k} (1-i x^2). - Robert A. Russell, Sep 26 2018

A056345 Number of bracelets of length n using exactly five different colored beads.

Original entry on oeis.org

0, 0, 0, 0, 12, 150, 1200, 7905, 46400, 255636, 1346700, 6901725, 34663020, 171786450, 843130688, 4110958530, 19951305240, 96528492700, 466073976900, 2247627076731, 10832193571460, 52194109216950
Offset: 1

Views

Author

Marks R. Nester

Keywords

Comments

Turning over will not create a new bracelet.

Examples

			For a(5)=12, pair up the 24 permutations of BCDE, each with its reverse, such as BCDE-EDCB.  Precede the first of each pair with an A, such as ABCDE.  These are the 12 arrangements, all chiral.  If we precede the second of each pair with an A, such as AEDCB, 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 5 of A273891.
Equals (A056285 + A056491) / 2 = A056285 - A305544 = A305544 + A056491.

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, 5];
Array[a, 22] (* Jean-François Alcover, Nov 05 2017, after Andrew Howroyd *)
k=5; 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 *)

Formula

a(n) = A032276(n) - 5*A032275(n) + 10*A027671(n) - 10*A000029(n) + 5.
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=5 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=5 is the number of colors. (End)

A056501 Number of primitive (period n) periodic palindromes using exactly five different symbols.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 60, 120, 960, 1800, 9300, 16800, 71400, 126000, 480000, 834120, 2968440, 5103000, 17354340, 29607600, 97566000, 165528000, 533264700, 901020120, 2854995360, 4809004080, 15050445900, 25292030400, 78417321240, 131542866000, 404936052000
Offset: 1

Views

Author

Marks R. Nester

Keywords

Examples

			For example, aaabbb is not a (finite) palindrome but it is a periodic palindrome.

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]

Links

Crossrefs

Column 5 of A327878.
Cf. A056466, A056491.

Formula

a(n) = Sum_{d|n} mu(d)*A056491(n/d).

Extensions

Terms corrected by Andrew Howroyd, Sep 28 2019

A305544 Number of chiral pairs of color loops of length n with exactly 5 different colors.

Original entry on oeis.org

0, 0, 0, 0, 12, 150, 1200, 7845, 46280, 254676, 1344900, 6892425, 34646220, 171715050, 843004688, 4110478470, 19950471120, 96525524140, 466068873900, 2247609721431, 10832163963860, 52194011649150, 251522234238000, 1212501695554920, 5848043487355752, 28223528190496380, 136307124614215660, 658800774340433025, 3186621527711606940
Offset: 1

Views

Author

Robert A. Russell, Jun 04 2018

Keywords

Examples

			For a(5)=12, the chiral pairs of color loops are ABCDE-AEDCB, ABCED-ADECB, ABDCE-AECDB, ABDEC-ACEDB, ABECD-ADCEB, ABEDC-ACDEB, ACBDE-AEDBC, ACBED-ADEBC, ACDBE-AEBCD, ACEDB-ABDEC, ADBCE-AECBD, ADBEC-ACEBD, and ADCBE-AEBCD.

Crossrefs

Fifth column of A305541.

Programs

  • Mathematica
    k=5; Table[(k!/(2n)) DivisorSum[n, EulerPhi[#] StirlingS2[n/#, k] &] - (k!/4) (StirlingS2[Floor[(n+1)/2], k] + StirlingS2[Ceiling[(n+1)/2], k]), {n, 1, 40}]
  • PARI
    a(n) = my(k=5); -(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) = -(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), with k=5 different colors used and where S2(n,k) is the Stirling subset number A008277.
a(n) = (A052825(n) - A056491(n)) / 2.
a(n) = A305541(n,5).
G.f.: -30 * x^8 * (1+x)^2 / Product_{j=1..5} (1-j*x^2) - Sum_{d>0} (phi(d)/(2d)) * (log(1-5x^d) - 5*log(1-4x^d) + 10*log(1-3x^3) - 10*log(1-2x^d) + 5*log(1-x^d)).
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.