OFFSET
0,2
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..500
Index entries for linear recurrences with constant coefficients, signature (31,31,31,31,31,31,31,31,31,31,-496).
FORMULA
G.f.: (t^11 + 2*t^10 + 2*t^9 + 2*t^8 + 2*t^7 + 2*t^6 + 2*t^5 + 2*t^4 + 2*t^3 + 2*t^2 + 2*t + 1)/(496*t^11 - 31*t^10 - 31*t^9 - 31*t^8 - 31*t^7 - 31*t^6 - 31*t^5 - 31*t^4 - 31*t^3 - 31*t^2 - 31*t + 1).
MATHEMATICA
With[{num=Total[2t^Range[10]]+t^11+1, den=Total[-31 t^Range[10]]+496t^11+ 1}, CoefficientList[Series[num/den, {t, 0, 30}], t]] (* Harvey P. Dale, Aug 16 2011 *)
With[{p=496, q=31}, CoefficientList[Series[(1+t)*(1-t^11)/(1-(q+1)*t + (p+q)*t^11-p*t^12), {t, 0, 40}], t]] (* G. C. Greubel, May 13 2016; Jul 25 2024 *)
coxG[{11, 496, -31, 30}] (* The coxG program is at A169452 *) (* G. C. Greubel, Jul 25 2024 *)
PROG
(Magma)
R<x>:=PowerSeriesRing(Integers(), 30);
Coefficients(R!( (1+x)*(1-x^11)/(1-32*x+527*x^11-496*x^12) )); // G. C. Greubel, Jul 25 2024
(SageMath)
def A166427_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P( (1+x)*(1-x^11)/(1-32*x+527*x^11-496*x^12) ).list()
A166427_list(30) # G. C. Greubel, Jul 25 2024
(PARI) a(n)=if(n, ([0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0; 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0; 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1; -496, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31]^(n-1)*[33; 1056; 33792; 1081344; 34603008; 1107296256; 35433480192; 1133871366144; 36283883716608; 1161084278931456; 37154696925806064])[1, 1], 1) \\ Charles R Greathouse IV, Jun 08 2026
(PARI) Vec((1+x+x^2+x^3+x^4+x^5+x^6+x^7+x^8+x^9+x^10)*(1+x)/(1-31*x-31*x^2-31*x^3-31*x^4-31*x^5-31*x^6-31*x^7-31*x^8-31*x^9-31*x^10+496*x^11)+O(x^99)) \\ Charles R Greathouse IV, Jun 08 2026
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
John Cannon and N. J. A. Sloane, Dec 03 2009
STATUS
approved
