%I #19 Jun 08 2022 03:27:36
%S 1,2,12,102,1176,16920,292320,5891760,135717120,3517032960,
%T 101268921600,3207514464000,110828037196800,4148515981209600,
%U 167232459621427200,7222900141416960000,332760193091149824000
%N Expansion of e.g.f. (1-x)/(1-3*x+x^3).
%H G. C. Greubel, <a href="/A052693/b052693.txt">Table of n, a(n) for n = 0..375</a>
%H INRIA Algorithms Project, <a href="http://ecs.inria.fr/services/structure?nbr=642">Encyclopedia of Combinatorial Structures 642</a>
%F E.g.f.: (1-x)/(1-3*x+x^3).
%F Recurrence: a(0)=1, a(1)=2, a(2)=12 a(n+3) = 3*(n+3)*a(n+2) - (n+1)*(n+2)*(n+3)*a(n).
%F a(n) = (n!/9)*Sum_{alpha=RootOf(1 -3*Z +Z^3)} (2 - alpha + alpha^2)*alpha^(-1-n).
%F a(n) = n! * A052536(n). - _G. C. Greubel_, Jun 01 2022
%p spec := [S,{S=Sequence(Union(Z,Prod(Z,Union(Z,Sequence(Z)))))},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
%t With[{nn=20},CoefficientList[Series[(1-x)/(1-3x+x^3),{x,0,nn}],x]Range[0,nn]!] (* _Harvey P. Dale_, Dec 17 2012 *)
%o (Magma) R<x>:=PowerSeriesRing(Rationals(), 40); Coefficients(R!(Laplace( (1-x)/(1-3*x+x^3) ))); // _G. C. Greubel_, Jun 01 2022
%o (SageMath)
%o @CachedFunction
%o def A052536(n):
%o if (n<3): return factorial(n+1)
%o else: return 3*A052536(n-1) - A052536(n-3)
%o def A052693(n): return factorial(n)*A052536(n)
%o [A052693(n) for n in (0..40)] # _G. C. Greubel_, Jun 01 2022
%Y Cf. A000142, A052536.
%K easy,nonn
%O 0,2
%A encyclopedia(AT)pommard.inria.fr, Jan 25 2000
|