%I M1295 N0496 #42 Aug 08 2021 11:34:54
%S 1,1,2,4,16,80,520,3640,29120,259840,2598400,28582400,343235200,
%T 4462057600,62468806400,936987251200,14991796019200,254860532326400,
%U 4587501779660800,87162533813555200,1743250676271104000,36608259566534656000,805381710463762432000
%N Expansion of e.g.f. exp((-x^3)/3)/(1-x).
%C a(n) is the number of permutations in the symmetric group S_n whose cycle decomposition contains no 3-cycle.
%D J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 85.
%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%D R. P. Stanley, Enumerative Combinatorics, Wadsworth, Vol. 1, 1986, page 93, problem 7.
%H Christian G. Bower, <a href="/A000090/b000090.txt">Table of n, a(n) for n = 0..100</a>
%H Simon Plouffe, <a href="http://www.plouffe.fr/simon/exact.htm">Exact formulas for integer sequences</a>
%H L. W. Shapiro & N. J. A. Sloane, <a href="/A006318/a006318_1.pdf">Correspondence, 1976</a>
%F a(n) = n! * Sum_{i=0..floor(n/3)} (-1)^i / (i! * 3^i); a(n)/n! ~ Sum_{i >= 0} (-1)^i / (i! * 3^i) = e^(-1/3); a(n) ~ e^(-1/3) * n!; a(n) ~ e^(-1/3) * (n/e)^n * sqrt(2 * Pi * n). - Avi Peretz (njk(AT)netvision.net.il), Apr 22 2001
%F 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=3, n>=0. - _Simon Plouffe_ from old notes, 1993
%F E.g.f.: E(x) = exp(-x^3/3)/(1-x)=G(0)/((1-x)^2); G(k) = 1 - x/(1 - x^2/(x^2 + 3*(k+1)/G(k+1))); (continued fraction). - _Sergei N. Gladkovskii_, Feb 11 2012
%e a(3) = 4 because the permutations in S_3 that contain no 3-cycles are the trivial permutation and the 3 transpositions.
%p seq(coeff(convert(series(exp((-x^3)/3)/(1-x),x,50),polynom),x,i)*i!,i=0..30);# series expansion A000090:=n->n!*add((-1)^i/(i!*3^i),i=0..floor(n/3));seq(A000090(n),n=0..30); # formula (Pab Ter)
%t nn=20;Range[0,nn]!CoefficientList[Series[Exp[-x^3/3]/(1-x),{x,0,nn}],x] (* _Geoffrey Critzer_, Oct 28 2012 *)
%o (PARI) {a(n) = if( n<0, 0, n! * polcoeff( exp( -(x^3 / 3) + x*O(x^n)) / (1 - x), n))} /* _Michael Somos_, Jul 28 2009 */
%Y Cf. A000142, A000138, A000266, A060725.
%K nonn,easy
%O 0,3
%A _N. J. A. Sloane_
%E More terms from Pab Ter (pabrlos2(AT)yahoo.com), Oct 22 2005
%E Entry improved by comments from _Michael Somos_, Jul 28 2009