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.

A000138 Expansion of e.g.f. exp(-x^4/4)/(1-x).

Original entry on oeis.org

1, 1, 2, 6, 18, 90, 540, 3780, 31500, 283500, 2835000, 31185000, 372972600, 4848643800, 67881013200, 1018215198000, 16294848570000, 277012425690000, 4986223662420000, 94738249585980000, 1894745192712372000, 39789649046959812000, 875372279033115864000
Offset: 0

Views

Author

Keywords

Comments

a(n) is the number of permutations in the symmetric group S_n whose cycle decomposition contains no 4-cycle.

Examples

			a(4) = 18 because in S_4 the permutations with no 4-cycle are the complement of the six 4-cycles so a(4) = 4! - 6 = 18.
		

References

  • J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 85.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Wadsworth, Vol. 1, 1986, page 93, problem 7.

Crossrefs

Programs

  • Mathematica
    nn=20;Range[0,nn]!CoefficientList[Series[Exp[-x^4/4]/(1-x),{x,0,nn}],x]  (* Geoffrey Critzer, Oct 28 2012 *)
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( exp( -(x^4/4) + x*O(x^n)) / (1 - x), n))} /* Michael Somos, Jul 28 2009 */

Formula

a(n) = n! * Sum_{i=0..floor(n/4)} (-1)^i / (i! * 4^i); a(n)/n! ~ Sum_{i >= 0} (-1)^i / (i! * 4^i) = e^(-1/4); a(n) ~ e^(-1/4) * n!; a(n) ~ e^(-1/4) * (n/e)^n * sqrt(2*Pi*n). - Avi Peretz (njk(AT)netvision.net.il), Apr 22 2001
a(n,k) = n!*floor(floor(n/k)!*k^floor(n/k)/exp(1/k) + 1/2)/(floor(n/k)!*k^floor(n/k)), here k=4, n>=0. Simon Plouffe, from old notes, 1993
E.g.f.: exp(-x^4/4)/(1-x) = 1/G(0); G(k) = 1 - x/(1 - (x^3)/(x^3 - 4*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Feb 28 2012

Extensions

Entry improved by comments from Michael Somos, Jul 28 2009
Name corrected by Joerg Arndt, May 27 2011

A060725 E.g.f.: exp(-(x^5/5))/(1-x).

Original entry on oeis.org

1, 1, 2, 6, 24, 96, 576, 4032, 32256, 290304, 2975616, 32731776, 392781312, 5106157056, 71486198784, 1070549415936, 17128790654976, 291189441134592, 5241409940422656, 99586788868030464, 1991897970827821056, 41829857387384242176, 920256862522453327872
Offset: 0

Views

Author

Avi Peretz (njk(AT)netvision.net.il), Apr 22 2001

Keywords

Comments

a(n) is the number of permutations in the symmetric group S_n whose cycle decomposition contains no 5-cycle.

Examples

			a(5) = 96 because in S_5 the permutations with no 5-cycle are the complement of the 24 5-cycles so a(5) = 5! - 24 = 96.
		

References

  • R. P. Stanley, Enumerative Combinatorics, Wadsworth, Vol. 1, 1986, page 93, problem 7.

Crossrefs

Programs

  • Maple
    for n from 0 to 30 do printf(`%d,`, n! * sum(( (-1)^i /(i! * 5^i)), i=0..floor(n/5))) od:
  • Mathematica
    With[{nn=30},CoefficientList[Series[Exp[-(x^5/5)]/(1-x),{x,0,nn}],x] Range[ 0,nn]!] (* Harvey P. Dale, Aug 24 2019 *)
  • PARI
    { for (n=0, 100, write("b060725.txt", n, " ", n! * sum(i=0, n\5, (-1)^i / (i! * 5^i))); ) } \\ Harry J. Smith, Jul 10 2009
    
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( exp( -(x^5 / 5) + x*O(x^n)) / (1 - x), n))} /* Michael Somos, Jul 28 2009 */
    
  • PARI
    { A060725_list(numterms) = Vec(serlaplace(exp(-x^5/5 + O(x^numterms))/(1-x))); } /* Eric M. Schmidt, Aug 22 2012 */

Formula

The formula for a(n) is: a(n) = n! * sum i=0 ... [ n/5 ]( (-1)^i /(i! * 5^i)) by this formula we have as n -> infinity: a(n)/n! ~ sum i >= 0 (-1)^i /(i! * 5^i) = e^(-1/5) or a(n) ~ e^(-1/5) * n! and using Stirling's formula in A000142: a(n) ~ e^(-1/5) * (n/e)^n * sqrt(2 * Pi * n).
a(n,k) = n!*floor(floor(n/k)!*k^floor(n/k)/exp(1/k) + 1/2)/(floor(n/k)! * k^floor(n/k)), k=5, n>=0. - Simon Plouffe, Feb 18 2011

Extensions

More terms from James Sellers, Apr 24 2001
Entry improved by comments from Michael Somos, Jul 28 2009

A060726 For n >= 1, a(n) is the number of permutations in the symmetric group S_n such that their cycle decomposition contains no 6-cycle.

Original entry on oeis.org

1, 1, 2, 6, 24, 120, 600, 4200, 33600, 302400, 3024000, 33264000, 405820800, 5275670400, 73859385600, 1107890784000, 17726252544000, 301346293248000, 5419293175296000, 102966570330624000, 2059331406612480000
Offset: 0

Views

Author

Avi Peretz (njk(AT)netvision.net.il), Apr 22 2001

Keywords

Comments

This is the expansion of exp ((-x^6)/6) /(1-x).

Examples

			a(6) = 600 because in S_6 the permutations with no 6-cycle are the complement of the 120 6-cycles so a(6) = 6! - 120 = 600.
		

References

  • R. P. Stanley, Enumerative Combinatorics, Wadsworth, Vol. 1, 1986, page 93, problem 7.

Crossrefs

Programs

  • Maple
    for n from 0 to 30 do printf(`%d,`, n! * sum(( (-1)^i /(i! * 6^i)), i=0..floor(n/6))) od:
  • PARI
    a(n)={n! * sum(i=0, n\6, (-1)^i / (i! * 6^i))} \\ Harry J. Smith, Jul 10 2009

Formula

The formula for a(n) is: a(n) = n! * Sum_{i=0..floor(n/6)} ((-1)^i /(i! * 6^i)) by this formula we have as n -> infinity: a(n)/n! ~ Sum_{i>= 0} (-1)^i /(i! * 6^i) = e^(-1/6) or a(n) ~ e^(-1/6) * n! and using Stirling's formula in A000142: a(n) ~ e^(-1/6) * (n/e)^n * sqrt(2 * Pi * n)
a(n,k) = n!*floor(floor(n/k)!*k^floor(n/k)/exp(1/k) + 1/2)/(floor(n/k)!*k^floor(n/k)), k=6, n >= 0. - Simon Plouffe, Feb 18 2011

Extensions

More terms from James Sellers, Apr 24 2001

A122974 Triangle T(n,k), the number of permutations on n elements that have no cycles of length k.

Original entry on oeis.org

0, 1, 1, 2, 3, 4, 9, 15, 16, 18, 44, 75, 80, 90, 96, 265, 435, 520, 540, 576, 600, 1854, 3045, 3640, 3780, 4032, 4200, 4320, 14833, 24465, 29120, 31500, 32256, 33600, 34560, 35280, 133496, 220185, 259840, 283500, 290304, 302400, 311040, 317520, 322560
Offset: 1

Views

Author

Dennis P. Walsh, Oct 27 2006

Keywords

Comments

Read as sequence, a(n) is the number of permutations on j elements with no cycles of length i where j=round((2*n)^.5) and i=n-C(j,2).
T(n,k) generalizes several sequences already in the On-Line Encyclopedia, such as A000166, the number of permutations on n elements with no fixed points and A000266, the number of permutations on n elements with no transpositions (i.e., no 2-cycles). See the cross references for further examples.

Examples

			T(3,2)=3 since there are exactly 3 permutations of 1,2,3 that have no cycles of length 2, namely, (1)(2)(3),(1 2 3) and (2 1 3).
Triangle T(n,k) begins:
      0;
      1,     1;
      2,     3,     4;
      9,    15,    16,    18;
     44,    75,    80,    90,    96;
    265,   435,   520,   540,   576,   600;
   1854,  3045,  3640,  3780,  4032,  4200,  4320;
  14833, 24465, 29120, 31500, 32256, 33600, 34560, 35280;
  ...
		

Crossrefs

Cf. T(n, 1)=A000166 for n=>1 T(n, 2)=A000266 for n=>2 T(n, 3)=A000090 for n=>3 T(n, 4)=A000138 for n=>4 T(n, 5)=A060725 for n=>5 T(n, 6)=A060726 for n=>6 T(n, 7)=A060727 for n=>7.
T(n,n) gives A094304(n+1).

Programs

  • Maple
    seq((round((2*n)^.5))!*sum((-1/(n-binomial(round((2*n)^.5),2)))^r/r!,r=0..floor(round((2*n)^.5)/(n-binomial(round((2*n)^.5),2)))),n=1..66);
    # second Maple program:
    T:= proc(n, k) option remember; `if`(n=0, 1, add(`if`(j=k, 0,
          T(n-j, k)*binomial(n-1, j-1)*(j-1)!), j=1..n))
        end:
    seq(seq(T(n, k), k=1..n), n=1..12);  # Alois P. Heinz, Nov 24 2019
  • Mathematica
    T[n_, k_] := T[n, k] = If[n==0, 1, Sum[If[j==k, 0, T[n - j, k] Binomial[n - 1, j - 1] (j - 1)!], {j, 1, n}]];
    Table[Table[T[n, k], {k, 1, n}], {n, 1, 12}] // Flatten (* Jean-François Alcover, Dec 08 2019, after Alois P. Heinz *)

Formula

T(n,k)=n!*sum r=0..floor(n/k)((-1/k)^r/r!) E.G.F: exp(-x^k/k)/(1-x) a(n)=(round((2*n)^.5))!*sum((-1/(n-binomial(round((2*n)^.5),2)))^r/r!,r=0..floor(round((2*n)^.5)/(n-binomial(round((2*n)^.5),2)))).
T(n,k) = n! - A293211(n,k). - Alois P. Heinz, Nov 24 2019
Showing 1-4 of 4 results.