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.

A052672 Expansion of e.g.f. (1-x)/(1-x-2*x^2+x^3).

This page as a plain text file.
%I A052672 #19 Jun 13 2022 02:38:12
%S A052672 1,0,4,6,120,600,10080,95760,1693440,23950080,475372800,8821612800,
%T A052672 199743667200,4533271142400,116906088499200,3112264995840000,
%U A052672 90679371374592000,2757644630028288000,89895729202126848000
%N A052672 Expansion of e.g.f. (1-x)/(1-x-2*x^2+x^3).
%H A052672 G. C. Greubel, <a href="/A052672/b052672.txt">Table of n, a(n) for n = 0..400</a>
%H A052672 INRIA Algorithms Project, <a href="http://ecs.inria.fr/services/structure?nbr=620">Encyclopedia of Combinatorial Structures 620</a>
%F A052672 E.g.f.: (1 - x)/(1 - x - 2*x^2 + x^3).
%F A052672 Recurrence: a(0)=1, a(1)=0, a(2)=4, a(n) = n*a(n-1) + 2*n*(n-1)*a(n-2) - n*(n-1)*(n-2)*a(n-3).
%F A052672 a(n) = (n!/7)*Sum_{alpha=RootOf(Z^3 -2*Z^2 -Z +1)} (3 - alpha)*alpha^(-n).
%F A052672 a(n) = n!*A052547(n). - _R. J. Mathar_, Nov 27 2011
%p A052672 spec := [S,{S=Sequence(Prod(Z,Union(Z,Prod(Z,Sequence(Z)))))},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
%t A052672 With[{nn=20},CoefficientList[Series[-(-1+x)/(x^3-2x^2-x+1),{x,0,nn}],x] Range[0,nn]!] (* _Harvey P. Dale_, Sep 25 2021 *)
%o A052672 (Magma) I:=[1,0,4]; [n le 3 select I[n] else (n-1)*(Self(n-1) +2*(n-2)*Self(n-2) -(n-2)*(n-3)*Self(n-3)): n in [1..31]]; // _G. C. Greubel_, Jun 13 2022
%o A052672 (SageMath)
%o A052672 def A052672_list(prec):
%o A052672     P.<x> = PowerSeriesRing(QQ, prec)
%o A052672     return P( (1-x)/(1-x-2*x^2+x^3) ).egf_to_ogf().list()
%o A052672 A052672_list(30) # _G. C. Greubel_, Jun 13 2022
%Y A052672 Cf. A000142, A052547.
%K A052672 easy,nonn
%O A052672 0,3
%A A052672 encyclopedia(AT)pommard.inria.fr, Jan 25 2000