%I #23 Dec 26 2022 15:31:55
%S 0,67,3079,65458,436705,3325420,21257887,137628082,852017725,
%T 5260500568,32028617995,194422680046,1174383558985,7081178928436,
%U 42616157629303,256244634375850,1539564650731285,9246057306575824,55510175964258211
%N Expansion of x*(67 +2476*x +38216*x^2 -124633*x^3 +129444*x^4)/((1-x)*(1+x)*(1-2*x)*(1+3*x)*(1-4*x)*(1-6*x)).
%H G. C. Greubel, <a href="/A120663/b120663.txt">Table of n, a(n) for n = 0..1000</a>
%H Roger L. Bagula, <a href="/A120663/a120663.txt">Mathematica program for A120663</a>
%H <a href="/index/Rec#order_06">Index entries for linear recurrences with constant coefficients</a>, signature (9,-7,-93,152,84,-144).
%F G.f.: x*(67 +2476*x +38216*x^2 -124633*x^3 +129444*x^4)/((1-x)*(1+x)*(1-2*x)*(1+3*x)*(1-4*x)*(1-6*x)). - _Colin Barker_, Nov 01 2012
%t See link for Mathematica program that uses matrices.
%t LinearRecurrence[{9,-7,-93,152,84,-144}, {0,67,3079,65458,436705, 3325420}, 31] (* _G. C. Greubel_, Dec 26 2022 *)
%o (Magma) R<x>:=PowerSeriesRing(Integers(), 30); [0] cat Coefficients(R!( x*(67+2476*x+38216*x^2-124633*x^3+129444*x^4)/(1-9*x+7*x^2+93*x^3 - 152*x^4-84*x^5+144*x^6) )); // _G. C. Greubel_, Dec 26 2022
%o (SageMath)
%o def f(x): return x*(67+2476*x+38216*x^2-124633*x^3+129444*x^4)/(1-9*x+7*x^2+93*x^3-152*x^4-84*x^5+144*x^6)
%o def A120663_list(prec):
%o P.<x> = PowerSeriesRing(ZZ, prec)
%o return P( f(x) ).list()
%o A120663_list(30) # _G. C. Greubel_, Dec 26 2022
%K nonn,easy,less
%O 0,2
%A _Roger L. Bagula_, Aug 10 2006
%E Edited by _N. J. A. Sloane_, Jul 13 2007
%E Meaningful name using g.f. from _Joerg Arndt_, Dec 26 2022
|