A008307 -id:A008307 - cp's OEIS frontend

A001470 Number of degree-n permutations of order dividing 3.

1, 1, 1, 3, 9, 21, 81, 351, 1233, 5769, 31041, 142011, 776601, 4874013, 27027729, 168369111, 1191911841, 7678566801, 53474964993, 418199988339, 3044269834281, 23364756531621, 199008751634001, 1605461415071823, 13428028220072049, 123280859122040601

Offset: 0

N. J. A. Sloane, J. H. Conway and Simon Plouffe

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, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Seiichi Manyama, Table of n, a(n) for n = 0..631 (terms 0..100 from T. D. Noe)
Joerg Arndt, Generating Random Permutations, PhD thesis, Australian National University, Canberra, Australia, (2010).
Marcello Artioli, Giuseppe Dattoli, Silvia Licciardi, and Simonetta Pagnutti, Motzkin Numbers: an Operational Point of View, arXiv:1703.07262 [math.CO], 2017. See p. 7.
L. Moser and M. Wyman, On Solutions of x^d = 1 in Symmetric Groups, Canad. J. Math., 7 (1955), 159-168.

Cf. A000085, A001472.

Column k=3 of A008307.

R:=PowerSeriesRing(Rationals(), 40); Coefficients(R!(Laplace( Exp(x+x^3/3) ))); // G. C. Greubel, Sep 03 2023

spec := [S, {S=Set(Union(Cycle(Z, card=1), Cycle(Z, card=3)))}, labeled]: seq(combstruct[count](spec, size=n), n=0..25) # David Radcliffe, Aug 29 2025

a[n_] := HypergeometricPFQ[{(1-n)/3, (2-n)/3, -n/3}, {}, -9]; Table[a[n], {n, 0, 23}] (* Jean-François Alcover, Nov 03 2011 *)
With[{nn=30},CoefficientList[Series[Exp[x+x^3/3],{x,0,nn}],x] Range[ 0,nn]!] (* Harvey P. Dale, Aug 12 2016 *)

a(n):=n!*sum(if mod(n-k,2)=0 then binomial(k,(3*k-n)/2)*(1/3)^((n-k)/2)/k! else 0,k,floor(n/3),n); /* Vladimir Kruchinin, Sep 07 2010 */

def A001470_list(prec):
    P. = PowerSeriesRing(QQ, prec)
    return P( exp(x+x^3/3) ).egf_to_ogf().list()
A001470_list(40) # G. C. Greubel, Sep 03 2023

a(n) = Sum_{j=0..floor(n/3)} n!/(j!*(n-3j)!*(3^j)) (the latter formula from Roger Cuculière).
E.g.f.: exp(x + (1/3)*x^3).
D-finite with recurrence: a(n) = a(n-1) + (n-1)*(n-2)*a(n-3). - Geoffrey Critzer, Feb 03 2009
a(n) = n!*Sum_{k=floor(n/3)..n, n - k == 0 (mod 2)} binomial(k,(3*k-n)/2)*(1/3)^((n-k)/2)/k!. - Vladimir Kruchinin, Sep 07 2010
a(n) ~ n^(2*n/3)*exp(n^(1/3)-2*n/3)/sqrt(3) * (1 - 1/(6*n^(1/3)) + 25/(72*n^(2/3))). - Vaclav Kotesovec, Jul 28 2013

A001472 Number of degree-n permutations of order dividing 4.

Original entry on oeis.org

1, 1, 2, 4, 16, 56, 256, 1072, 6224, 33616, 218656, 1326656, 9893632, 70186624, 574017536, 4454046976, 40073925376, 347165733632, 3370414011904, 31426411211776, 328454079574016, 3331595921852416, 37125035407900672, 400800185285464064, 4744829049220673536

Offset: 0

N. J. A. Sloane and J. H. Conway

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, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Alois P. Heinz, Table of n, a(n) for n = 0..570 (first 201 terms from T. D. Noe)
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 25 (Dead link)
Vladimir Victorovich Kruchinin, Composition of ordinary generating functions, arXiv:1009.2565 [math.CO], 2010.
L. Moser and M. Wyman, On solutions of x^d = 1 in symmetric groups, Canad. J. Math., 7 (1955), 159-168.

Cf. A000085, A001470, A053495.

Column k=4 of A008307.

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x + x^2/2 +x^4/4) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, May 14 2019

spec := [S, {S = Set(Union(Cycle(Z, card = 1), Cycle(Z, card = 2), Cycle(Z, card = 4)))}, labeled]; seq(combstruct[count](spec, size = n), n = 0 .. 23); # David Radcliffe, Aug 29 2025

n = 23; CoefficientList[Series[Exp[x+x^2/2+x^4/4], {x, 0, n}], x] * Table[k!, {k, 0, n}] (* Jean-François Alcover, May 18 2011 *)

a(n):=n!*sum(sum(binomial(k,j)*binomial(j,n-4*k+3*j)*(1/2)^(n-4*k+3*j)*(1/4)^(k-j),j,floor((4*k-n)/3),k)/k!,k,1,n); /* Vladimir Kruchinin, Sep 07 2010 */

my(N=33, x='x+O('x^N)); egf=exp(x+x^2/2+x^4/4); Vec(serlaplace(egf)) /* Joerg Arndt, Sep 15 2012 */

m = 30; T = taylor(exp(x + x^2/2 + x^4/4), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, May 14 2019

E.g.f.: exp(x + x^2/2 + x^4/4).
D-finite with recurrence: a(0)=1, a(1)=1, a(2)=2, a(3)=4, a(n) = a(n-1) + (n-1)*a(n-2) + (n^3-6*n^2+11*n-6)*a(n-4) for n>3. - H. Palsdottir (hronn07(AT)ru.is), Sep 19 2008
a(n) = n!*Sum_{k=1..n} (1/k!)*Sum_{j=floor((4*k-n)/3)..k} binomial(k,j) * binomial(j,n-4*k+3*j) * (1/2)^(n-4*k+3*j)*(1/4)^(k-j), n>0. - Vladimir Kruchinin, Sep 07 2010
a(n) ~ n^(3*n/4)*exp(n^(1/4)-3*n/4+sqrt(n)/2-1/8)/2 * (1 - 1/(4*n^(1/4)) + 17/(96*sqrt(n)) + 47/(128*n^(3/4))). - Vaclav Kotesovec, Aug 09 2013

A053505 Number of degree-n permutations of order dividing 30.

Original entry on oeis.org

1, 1, 2, 6, 18, 90, 540, 3060, 20700, 145980, 1459800, 13854600, 140059800, 1514748600, 15869034000, 285268878000, 4109761962000, 59488383690000, 935767530036000, 13364309726748000, 240338216104020000, 4540941256642020000, 79739974380153240000

Offset: 0

N. J. A. Sloane, Jan 15 2000

nonn

R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Alois P. Heinz, Table of n, a(n) for n = 0..200
L. Moser and M. Wyman, On solutions of x^d = 1 in symmetric groups, Canad. J. Math., 7 (1955), 159-168.

Cf. A000085, A001470, A001472, A053495-A053505, A005388.

Column k=30 of A008307.

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x +x^2/2 +x^3/3 +x^5/5 +x^6/6 +x^10/10 +x^15/15 +x^30/30) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, May 15 2019

a:= proc(n) option remember; `if`(n<0, 0, `if`(n=0, 1,
       add(mul(n-i, i=1..j-1)*a(n-j), j=[1, 2, 3, 5, 6, 10, 15, 30])))
    end:
seq(a(n), n=0..25);  # Alois P. Heinz, Feb 14 2013

a[n_]:= a[n] = If[n<0, 0, If[n==0, 1, Sum[Product[n-i, {i, 1, j-1}]*a[n-j], {j, {1, 2, 3, 5, 6, 10, 15, 30}}]]]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Mar 03 2014, after Alois P. Heinz *)
With[{m = 30}, CoefficientList[Series[Exp[x +x^2/2 +x^3/3 +x^5/5 +x^6/6 + x^10/10 +x^15/15 +x^30/30], {x, 0, m}], x]*Range[0, m]!] (* G. C. Greubel, May 15 2019 *)

my(x='x+O('x^30)); Vec(serlaplace( exp(x +x^2/2 +x^3/3 +x^5/5 + x^6/6 +x^10/10 +x^15/15 +x^30/30) )) \\ G. C. Greubel, May 15 2019

m = 30; T = taylor(exp(x +x^2/2 +x^3/3 +x^5/5 +x^6/6 +x^10/10 +x^15/15 +x^30/30), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, May 15 2019

E.g.f.: exp(x + x^2/2 + x^3/3 + x^5/5 + x^6/6 + x^10/10 + x^15/15 + x^30/30).

A052501 Number of permutations sigma such that sigma^5=Id; degree-n permutations of order dividing 5.

Original entry on oeis.org

1, 1, 1, 1, 1, 25, 145, 505, 1345, 3025, 78625, 809425, 4809025, 20787625, 72696625, 1961583625, 28478346625, 238536558625, 1425925698625, 6764765838625, 189239120970625, 3500701266525625, 37764092547420625, 288099608198025625

Offset: 0

N. J. A. Sloane, Jan 15 2000; encyclopedia(AT)pommard.inria.fr, Jan 25 2000

The number of degree-n permutations of order exactly p (where p is prime) satisfies a(n) = a(n-1) + (1 + a(n-p))*(n-1)!/(n-p)! with a(n)=0 if p>n. Also a(n) = Sum_{j=1..floor(n/p)} (n!/(j!*(n-p*j)!*(p^j))).

These are the telephone numbers T^(5)n of [Artioli et al., p. 7]. - _Eric M. Schmidt, Oct 12 2017

R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Alois P. Heinz, Table of n, a(n) for n = 0..300
Marcello Artioli, Giuseppe Dattoli, Silvia Licciardi, Simonetta Pagnutti, Motzkin Numbers: an Operational Point of View, arXiv:1703.07262 [math.CO], 2017.
Tomislav Došlic, Darko Veljan, Logarithmic behavior of some combinatorial sequences, Discrete Math. 308 (2008), no. 11, 2182--2212. MR2404544 (2009j:05019). - _N. J. A. Sloane_, May 01 2012
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 26
M. B. Kutler, C. R. Vinroot, On q-Analogs of Recursions for the Number of Involutions and Prime Order Elements in Symmetric Groups, JIS 13 (2010) #10.3.6.
L. Moser and M. Wyman, On solutions of x^d = 1 in symmetric groups, Canad. J. Math., 7 (1955), 159-168.

Cf. A000085, A001470, A001472, A053495-A053505, A005388.

Column k=5 of A008307.

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x + x^5/5) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, May 14 2019

spec := [S,{S=Set(Union(Cycle(Z,card=1),Cycle(Z,card=5)))},labeled]: seq(combstruct[count](spec,size=n), n=0..20);

max = 30; CoefficientList[ Series[ Exp[x + x^5/5], {x, 0, max}], x]*Range[0, max]! (* Jean-François Alcover, Feb 15 2012, after e.g.f. *)

my(x='x+O('x^30)); Vec(serlaplace( exp(x + x^5/5) )) \\ G. C. Greubel, May 14 2019

m = 30; T = taylor(exp(x + x^5/5), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, May 14 2019

E.g.f.: exp(x + x^5/5).
a(n+5) = a(n+4) + (24 +50*n +35*n^2 +10*n^3 +n^4)*a(n), with a(0)= ... = a(4) = 1.
a(n) = a(n-1) + a(n-5)*(n-1)!/(n-5)!.
a(n) = Sum_{j = 0..floor(n/5)} n!/(5^j * j! * (n-5*j)!).
a(n) = A059593(n) + 1.

A053496 Number of degree-n permutations of order dividing 6.

Original entry on oeis.org

1, 1, 2, 6, 18, 66, 396, 2052, 12636, 91548, 625176, 4673736, 43575192, 377205336, 3624289488, 38829340656, 397695226896, 4338579616272, 54018173703456, 641634784488288, 8208962893594656, 113809776294348576, 1526808627197721792, 21533423236302943296

Offset: 0

N. J. A. Sloane, Jan 15 2000

nonn

R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Alois P. Heinz, Table of n, a(n) for n = 0..200
L. Moser and M. Wyman, On solutions of x^d = 1 in symmetric groups, Canad. J. Math., 7 (1955), 159-168.

Cf. A000085, A001470, A001472, A053495-A053505, A005388, A261317.

Column k=6 of A008307.

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x + x^2/2 +x^3/3 +x^6/6) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, May 14 2019

a:= proc(n) option remember; `if`(n<0, 0, `if`(n=0, 1,
       add(mul(n-i, i=1..j-1)*a(n-j), j=[1, 2, 3, 6])))
    end:
seq(a(n), n=0..25);  # Alois P. Heinz, Feb 14 2013

a[n_] := a[n] = If[n<0, 0, If[n == 0, 1, Sum[Product[n-i, {i, 1, j-1}]*a[n-j], {j, {1, 2, 3, 6}}]]]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Mar 03 2014, after Alois P. Heinz *)
With[{m = 30}, CoefficientList[Series[Exp[x +x^2/2 +x^3/3 +x^6/6], {x, 0, m}], x]*Range[0, m]!] (* G. C. Greubel, May 14 2019 *)

my(x='x+O('x^30)); Vec(serlaplace( exp(x+x^2/2+x^3/3+x^6/6) )) \\ G. C. Greubel, May 14 2019

m = 30; T = taylor(exp(x +x^2/2 +x^3/3 +x^6/6), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, May 14 2019

E.g.f.: exp(x +x^2/2 +x^3/3 +x^6/6).

D-finite with recurrence a(n) -a(n-1) +(-n+1)*a(n-2) -(n-1)*(n-2)*a(n-3) -(n-5)*(n-1)*(n-2)*(n-3)*(n-4)*a(n-6)=0. - R. J. Mathar, Jul 04 2023

A261430 Number A(n,k) of permutations p of [n] without fixed points such that p^k = Id; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 2, 3, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 9, 0, 15, 0, 0, 1, 0, 0, 2, 0, 0, 40, 0, 0, 0, 1, 0, 1, 0, 3, 24, 105, 0, 105, 0, 0, 1, 0, 0, 0, 0, 20, 0, 0, 0, 0, 0, 0, 1, 0, 1, 2, 9, 0, 175, 0, 2625, 2240, 945, 0, 0

Offset: 0

Alois P. Heinz, Aug 18 2015

			Square array A(n,k) begins:
  1, 1,   1,  1,    1,  1,    1,   1,    1, ...
  0, 0,   0,  0,    0,  0,    0,   0,    0, ...
  0, 0,   1,  0,    1,  0,    1,   0,    1, ...
  0, 0,   0,  2,    0,  0,    2,   0,    0, ...
  0, 0,   3,  0,    9,  0,    3,   0,    9, ...
  0, 0,   0,  0,    0, 24,   20,   0,    0, ...
  0, 0,  15, 40,  105,  0,  175,   0,  105, ...
  0, 0,   0,  0,    0,  0,  210, 720,    0, ...
  0, 0, 105,  0, 2625,  0, 4585,   0, 7665, ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Columns k=0+1,2-10 give: A000007, A001147, A052502, A052503, A052504, A261317, A261427, A261428, A261429, A261381.
Main diagonal gives A261431.
Cf. A008307.

with(numtheory):
A:= proc(n, k) option remember; `if`(n<0, 0, `if`(n=0, 1,
      add(mul(n-i, i=1..j-1)*A(n-j, k), j=divisors(k) minus {1})))
    end:
seq(seq(A(n, d-n), n=0..d), d=0..14);

A[0, 0] = A[0, 1] = 1; A[, 0|1] = 0; A[n, k_] := A[n, k] = If[n < 0, 0, If[n == 0, 1, Sum[Product[n - i, {i, 1, j - 1}]*A[n - j, k], {j, Rest @ Divisors[k]}]]]; Table[A[n, d - n], {d, 0, 14}, {n, 0, d}] // Flatten (* Jean-François Alcover, Jan 21 2017, after Alois P. Heinz *)

E.g.f. of column k: exp(Sum_{d|k, d>1} x^d/d).

A053497 Number of degree-n permutations of order dividing 7.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 721, 5761, 25921, 86401, 237601, 570241, 1235521, 892045441, 13348249201, 106757164801, 604924594561, 2722120577281, 10344007402561, 34479959558401, 24928970490633601, 546446134633639681, 6281586217487489041, 50248618811434961281

Offset: 0

N. J. A. Sloane, Jan 15 2000

nonn

R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Alois P. Heinz, Table of n, a(n) for n = 0..200
L. Moser and M. Wyman, On solutions of x^d = 1 in symmetric groups, Canad. J. Math., 7 (1955), 159-168.

Sequences with e.g.f. exp(x + x^m/m): A000079 (m=1), A000085 (m=2), A001470 (m=3), A118934 (m=4), A052501 (m=5), A293588 (m=6), this sequence (m=7).
Cf. A000085, A001470, A001472, A005388, A053495 - A053505, A261427.
Column k=7 of A008307.

R:=PowerSeriesRing(Rationals(), 31); Coefficients(R!(Laplace( Exp(x + x^7/7) ))); // G. C. Greubel, May 14 2019, Mar 07 2021

a:= proc(n) option remember; `if`(n<0, 0, `if`(n=0, 1,
       add(mul(n-i, i=1..j-1)*a(n-j), j=[1, 7])))
    end:
seq(a(n), n=0..25); # Alois P. Heinz, Feb 14 2013

CoefficientList[Series[Exp[x+x^7/7], {x, 0, 24}], x]*Range[0, 24]! (* Jean-François Alcover, Mar 24 2014 *)

my(x='x+O('x^30)); Vec(serlaplace( exp(x+x^7/7) )) \\ G. C. Greubel, May 14 2019

f=factorial; [sum(f(n)/(7^j*f(j)*f(n-7*j)) for j in (0..n/7)) for n in (0..30)] # G. C. Greubel, May 14 2019

E.g.f.: exp(x + x^7/7).

a(n) = Sum_{k=0..floor(n/7)} n!/(7^k*k!*(n-7*k)!). - G. C. Greubel, Mar 07 2021

A053499 Number of degree-n permutations of order dividing 9.

Original entry on oeis.org

1, 1, 1, 3, 9, 21, 81, 351, 1233, 46089, 434241, 2359611, 27387801, 264333213, 1722161169, 16514298711, 163094452641, 1216239520401, 50883607918593, 866931703203699, 8473720481213481, 166915156382509221, 2699805625227141201, 28818706120636531023, 439756550972215638129, 6766483260087819272601, 77096822666547068590401, 3568144263578808757678251

Offset: 0

N. J. A. Sloane, Jan 15 2000

nonn

Differs from A218003 first at n=27. - Alois P. Heinz, Jan 25 2014

R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Alois P. Heinz, Table of n, a(n) for n = 0..200
L. Moser and M. Wyman, On solutions of x^d = 1 in symmetric groups, Canad. J. Math., 7 (1955), 159-168.

Cf. A000085, A001470, A001472, A053495-A053505, A005388, A261429.

Column k=9 of A008307.

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x + x^3/3 + x^9/9) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, May 15 2019

a:= proc(n) option remember; `if`(n<0, 0, `if`(n=0, 1,
       add(mul(n-i, i=1..j-1)*a(n-j), j=[1, 3, 9])))
    end:
seq(a(n), n=0..25);  # Alois P. Heinz, Feb 14 2013

CoefficientList[Series[Exp[x+x^3/3+x^9/9], {x, 0, 30}], x]*Range[0, 30]! (* Jean-François Alcover, Mar 24 2014 *)

my(x='x+O('x^30)); Vec(serlaplace( exp(x + x^3/3 + x^9/9) )) \\ G. C. Greubel, May 15 2019

m = 30; T = taylor(exp(x + x^3/3 + x^9/9), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, May 15 2019

E.g.f.: exp(x + x^3/3 + x^9/9).

A053502 Number of degree-n permutations of order dividing 12.

Original entry on oeis.org

1, 1, 2, 6, 24, 96, 576, 3312, 21456, 152784, 1237536, 9984096, 133494912, 1412107776, 16369357824, 206123325696, 2866280276736, 36809077162752, 592066290710016, 8800038127378944, 136876273991755776, 2197453620220010496, 37915306084793106432

Offset: 0

N. J. A. Sloane, Jan 15 2000

nonn

R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Alois P. Heinz, Table of n, a(n) for n = 0..200
L. Moser and M. Wyman, On solutions of x^d = 1 in symmetric groups, Canad. J. Math., 7 (1955), 159-168.

Cf. A000085, A001470, A001472, A053495-A053505, A005388.

Column k=12 of A008307.

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x + x^2/2 + x^3/3 + x^4/4 + x^6/6 + x^12/12) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, May 15 2019

a:= proc(n) option remember; `if`(n<0, 0, `if`(n=0, 1,
       add(mul(n-i, i=1..j-1)*a(n-j), j=[1, 2, 3, 4, 6, 12])))
    end:
seq(a(n), n=0..25); # Alois P. Heinz, Feb 14 2013

a[n_]:= a[n] = If[n<0, 0, If[n==0, 1, Sum[Product[n-i, {i, 1, j-1}]*a[n-j], {j, {1, 2, 3, 4, 6, 12}}]]]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Apr 24 2014, after Alois P. Heinz *)
With[{m = 30}, CoefficientList[Series[Exp[x +x^2/2 +x^3/3 +x^4/4 +x^6/6 + x^12/12], {x, 0, m}], x]*Range[0, m]!] (* G. C. Greubel, May 15 2019 *)

my(x='x+O('x^30)); Vec(serlaplace( exp(x + x^2/2 + x^3/3 + x^4/4 + x^6/6 + x^12/12) )) \\ G. C. Greubel, May 15 2019

m = 30; T = taylor(exp(x + x^2/2 + x^3/3 + x^4/4 + x^6/6 + x^12/12), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, May 15 2019

E.g.f.: exp(x + x^2/2 + x^3/3 + x^4/4 + x^6/6 + x^12/12).

A074759 Number of degree-n permutations of order dividing n. Number of solutions to x^n = 1 in S_n.

Original entry on oeis.org

1, 1, 2, 3, 16, 25, 396, 721, 11264, 46089, 602200, 3628801, 133494912, 479001601, 7692266960, 95904273375, 1914926104576, 20922789888001, 628693317946656, 6402373705728001, 182635841123840000, 2496321046987530021, 55826951075231672512, 1124000727777607680001

Offset: 0

Vladeta Jovovic, Sep 28 2002

Alois P. Heinz, Table of n, a(n) for n = 0..450

Cf. A057731, A074351, A261431.

Main diagonal of A008307.

A:= proc(n,k) option remember; `if`(n<0, 0, `if`(n=0, 1,
       add(mul(n-i, i=1..j-1)*A(n-j,k), j=numtheory[divisors](k))))
    end:
a:= n-> A(n, n):
seq(a(n), n=0..25);  # Alois P. Heinz, Feb 14 2013

Table[a = Sum[x^i/i, {i, Divisors[n]}]; Part[Range[0, 20]! CoefficientList[Series[Exp[a], {x, 0, 20}], x],n + 1], {n, 0, 20}]  (* Geoffrey Critzer, Dec 04 2011 *)

a(n) = n! * [x^n] exp(Sum_{k divides n} x^k/k).

a(n) = Sum_{d|n} A057731(n,d) for n >= 1. - Alois P. Heinz, Jul 05 2021

cp's OEIS Frontend

A001470 Number of degree-n permutations of order dividing 3.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

SageMath

Formula

A001472 Number of degree-n permutations of order dividing 4.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

PARI

Sage

Formula

A053505 Number of degree-n permutations of order dividing 30.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A052501 Number of permutations sigma such that sigma^5=Id; degree-n permutations of order dividing 5.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A053496 Number of degree-n permutations of order dividing 6.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A261430 Number A(n,k) of permutations p of [n] without fixed points such that p^k = Id; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Examples

Links

Crossrefs