A017883 Expansion of 1/(1-x^9-x^10-x^11-x^12-x^13-x^14-x^15-x^16).
1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 7, 7, 8, 10, 13, 17, 22, 28, 36, 42, 47, 52, 58, 66, 77, 92, 112, 141, 176, 215, 257, 302, 351, 406, 470, 546, 645, 774, 937, 1136, 1372, 1646
Offset: 0
Links
- Vincenzo Librandi, Table of n, a(n) for n = 0..1000
- Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1).
Programs
-
Magma
m:=70; R
:=PowerSeriesRing(Integers(), m); Coefficients(R!(1/(1-x^9-x^10-x^11-x^12-x^13-x^14-x^15-x^16))); // Vincenzo Librandi, Jul 01 2013 -
Mathematica
CoefficientList[Series[1 / (1 - Total[x^Range[9, 16]]), {x, 0, 70}], x] (* Vincenzo Librandi, Jul 01 2013 *)
-
SageMath
def A017883_list(prec): P.
= PowerSeriesRing(ZZ, prec) return P( (1-x)/(1-x-x^9+x^(17)) ).list() A017883_list(65) # G. C. Greubel, Sep 25 2024
Formula
a(n) = a(n-9) +a(n-10) +a(n-11) +a(n-12) +a(n-13) +a(n-14) +a(n-15) +a(n-16) for n>15. - Vincenzo Librandi, Jul 01 2013
Comments