login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A017878
Expansion of 1/(1-x^9-x^10-x^11).
10
1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 2, 3, 2, 1, 0, 0, 0, 0, 1, 3, 6, 7, 6, 3, 1, 0, 0, 1, 4, 10, 16, 19, 16, 10, 4, 1, 1, 5, 15, 30, 45, 51, 45, 30, 15, 6, 7, 21, 50, 90, 126, 141, 126, 90, 51, 28, 34, 78, 161
OFFSET
0,20
COMMENTS
Number of compositions (ordered partitions) of n into parts 9, 10 and 11. - Ilya Gutkovskiy, May 27 2017
LINKS
FORMULA
a(n) = a(n-9) +a(n-10) +a(n-11) for n>10. - Vincenzo Librandi, Jul 01 2013
MATHEMATICA
CoefficientList[Series[1 / (1 - Total[x^Range[9, 11]]), {x, 0, 80}], x] (* Vincenzo Librandi, Jul 01 2013 *)
LinearRecurrence[{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1}, {1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1}, 70] (* Harvey P. Dale, May 25 2023 *)
PROG
(Magma)
m:=70; R<x>:=PowerSeriesRing(Integers(), m);
Coefficients(R!(1/(1-x^9-x^10-x^11))); // Vincenzo Librandi, Jul 01 2013
(SageMath)
def A017878_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P( (1-x)/(1-x-x^9+x^(12)) ).list()
A017878_list(80) # G. C. Greubel, Sep 25 2024
KEYWORD
nonn,easy
STATUS
approved