A271391 Expansion of (1 + x + 2*x^2 + 6*x^3 + x^4 + x^5)/(1 - x^2)^3.

1, 1, 5, 9, 13, 25, 25, 49, 41, 81, 61, 121, 85, 169, 113, 225, 145, 289, 181, 361, 221, 441, 265, 529, 313, 625, 365, 729, 421, 841, 481, 961, 545, 1089, 613, 1225, 685, 1369, 761, 1521, 841, 1681, 925, 1849, 1013, 2025, 1105, 2209, 1201, 2401, 1301, 2601, 1405

Offset

0,3

Comments

Centered square numbers alternating with odd squares.

Links

Jinyuan Wang, Table of n, a(n) for n = 0..1000

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

Formula

G.f.: (x + 5*x^2 + 6*x^3 - 2*x^4 + x^5 + x^6)/(1 - x^2)^3.

E.g.f.: ((2 + x*(2 + x))*cosh(x) + x*(3 + 2*x)*sinh(x))/2.

a(n) = 3*a(n-2) - 3*a(n-4) + a(n-6).

a(n) = (3*n^2 + 2*n + 2 + (-1)^n*(-n^2 + 2*n + 2))/4.

Example

Illustration of initial terms:
                                                      o
                            o       o o o o o       o o o
          o     o o o     o o o     o o o o o     o o o o o
o   o   o o o   o o o   o o o o o   o o o o o   o o o o o o o
          o     o o o     o o o     o o o o o     o o o o o
                            o       o o o o o       o o o
                                                      o
0   1     2       3         4           5             6

Maple

a:=series((1+x+2*x^2+6*x^3+x^4+x^5)/(1-x^2)^3, x=0, 55): seq(coeff(a, x, n), n=0..54); # Paolo P. Lava, Mar 27 2019

Mathematica

LinearRecurrence[{0, 3, 0, -3, 0, 1}, {1, 1, 5, 9, 13, 25}, 55]
Table[(3 n^2 + 2 n + 2 + (-1)^n (-n^2 + 2 n + 2))/4, {n, 0, 54}]

Prog

(PARI) x='x+O('x^99); Vec((1+x+2*x^2+6*x^3+x^4+x^5)/(1-x^2)^3) \\ Altug Alkan, Apr 06 2016

Crossrefs

Cf. A001844, A016754.

Sequence in context: A314798 A208718 A208774 * A151907 A151895 A267190

Adjacent sequences: A271388 A271389 A271390 * A271392 A271393 A271394

Keyword

nonn,easy

Author

Ilya Gutkovskiy, Apr 06 2016

Status

approved