A008817 Expansion of (1+x^10)/((1-x)^2*(1-x^10)).

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 13, 16, 19, 22, 25, 28, 31, 34, 37, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 97, 104, 111, 118, 125, 132, 139, 146, 153, 160, 169, 178, 187, 196, 205, 214, 223, 232, 241, 250, 261, 272, 283, 294, 305, 316, 327, 338, 349, 360

Offset

0,2

Links

G. C. Greubel, Table of n, a(n) for n = 0..1000

Index entries for linear recurrences with constant coefficients, signature (2,-1,0,0,0,0,0,0,0,1,-2,1).

Formula

G.f.: (1+x^10)/((1-x)^2*(1-x^10)).

a(0)=1, a(1)=2, a(2)=3, a(3)=4, a(4)=5, a(5)=6, a(6)=7, a(7)=8, a(8)=9, a(9)=10, a(10)=13, a(11)=16, a(n) = 2*a(n-1) - a(n-2) + a(n-10) - 2*a(n-11) + a(n-12). - Harvey P. Dale, Jul 31 2014

Maple

seq(coeff(series((1+x^10)/((1-x)^2*(1-x^10)), x, n+1), x, n), n = 0..80); # G. C. Greubel, Sep 12 2019

Mathematica

CoefficientList[Series[(1+x^10)/(1-x)^2/(1-x^10), {x, 0, 80}], x] (* or *) LinearRecurrence[{2, -1, 0, 0, 0, 0, 0, 0, 0, 1, -2, 1}, {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 13, 16}, 80] (* Harvey P. Dale, Jul 31 2014 *)

Prog

(PARI) my(x='x+O('x^80)); Vec((1+x^10)/((1-x)^2*(1-x^10))) \\ G. C. Greubel, Sep 12 2019
(Magma) R<x>:=PowerSeriesRing(Integers(), 80); Coefficients(R!( (1+x^10)/((1-x)^2*(1-x^10)) )); // G. C. Greubel, Sep 12 2019
(Sage)
def A008817_list(prec):
    P.<x> = PowerSeriesRing(ZZ, prec)
    return P((1+x^10)/((1-x)^2*(1-x^10))).list()
A008817_list(80) # G. C. Greubel, Sep 12 2019
(GAP) a:=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 13, 16];; for n in [13..80] do a[n]:=2*a[n-1]-a[n-2]+a[n-10]-2*a[n-11]+a[n-12]; od; a; # G. C. Greubel, Sep 12 2019

Crossrefs

Cf. Expansions of the form (1+x^m)/((1-x)^2*(1-x^m)): A000290 (m=1), A000982 (m=2), A008810 (m=3), A008811 (m=4), A008812 (m=5), A008813 (m=6), A008814 (m=7), A008815 (m=8), A008816 (m=9), this sequence (m=10).

Sequence in context: A092968 A177053 A151547 * A076356 A308493 A263443

Adjacent sequences: A008814 A008815 A008816 * A008818 A008819 A008820

Keyword

nonn

Author

N. J. A. Sloane

Status

approved