%I #12 Oct 26 2022 03:05:07
%S 1,1,1,0,2,3,5,3,6,8,16,16,24,28,50,61,91,109,170,220,327,415,607,800,
%T 1164,1536,2192,2928,4172,5616,7921,10705,15049,20460,28638,39027,
%U 54453,74451,103662,141996,197288,270704,375632,516096,715258,983661,1362091
%N Expansion of 1/(1 - x + x^3 - 3*x^4 + x^5 - x^7 + x^8).
%H G. C. Greubel, <a href="/A147593/b147593.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_08">Index entries for linear recurrences with constant coefficients</a>, signature (1,0,-1,3,-1,0,1,-1).
%F G.f.: -1/(x^4*f(x)*f(1/x)), where f(x) = -1 - x^3 + x^4.
%F G.f.: 1/((1+x^3-x^4)*(1-x-x^4)). - _Colin Barker_, Nov 04 2012
%t f[x_]= x^4-x^3-1; CoefficientList[Series[-1/(x^4*f[x]*f[1/x]), {x,0,50}], x]
%o (Magma) R<x>:=PowerSeriesRing(Integers(), 50); Coefficients(R!( 1/((1+x^3-x^4)*(1-x-x^4)) )); // _G. C. Greubel_, Oct 25 2022
%o (SageMath)
%o def A147593_list(prec):
%o P.<x> = PowerSeriesRing(ZZ, prec)
%o return P( 1/((1+x^3-x^4)*(1-x-x^4)) ).list()
%o A147593_list(50) # _G. C. Greubel_, Oct 25 2022
%Y Cf. A147598, A147605, A147606, A147607, A147617, A147620.
%K nonn,easy
%O 0,5
%A _Roger L. Bagula_, Nov 08 2008
%E Edited by _Joerg Arndt_ and _Colin Barker_, Nov 04 2012.
|