A331321 a(n) = [x^n] ((x^2 - 1)*(x^2 + x - 1))/(x^2 + 2*x - 1)^2.

1, 3, 8, 23, 64, 175, 472, 1259, 3328, 8731, 22760, 59007, 152256, 391239, 1001656, 2556115, 6503936, 16505651, 41788616, 105571303, 266181440, 669923039, 1683255448, 4222878651, 10579130112, 26467818315, 66138242984, 165077936207, 411584855488, 1025162759287

Offset

0,2

Links

Colin Barker, Table of n, a(n) for n = 0..1000

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

Formula

a(n) = Sum_{k=0..n} A193737(n, k)*(1 + k).

Let h(k) = (1 + k)*exp((1 + k)*x)*(2*x + 10 - 5*k)/8 then

a(n) = n!*[x^n](h(sqrt(2)) + h(-sqrt(2)) + 1).

From Colin Barker, Jan 14 2020: (Start)

a(n) = 4*a(n-1) - 2*a(n-2) - 4*a(n-3) - a(n-4) for n>4.

a(n) = (-5*sqrt(2)*((1-sqrt(2))^n - (1+sqrt(2))^n) + 2*((1-sqrt(2))^n + (1+sqrt(2))^n)*n) / 8 for n>0.

(End)

Maple

gf := ((x^2 - 1)*(x^2 + x - 1))/(x^2 + 2*x - 1)^2:
ser := series(gf, x, 32): seq(coeff(ser, x, n), n=0..29);

Mathematica

LinearRecurrence[{4, -2, -4, -1}, {1, 3, 8, 23, 64}, 40] (* Harvey P. Dale, Feb 01 2022 *)

Prog

(PARI) Vec((1 - x)*(1 + x)*(1 - x - x^2) / (1 - 2*x - x^2)^2 + O(x^30)) \\ Colin Barker, Jan 14 2020
(Magma) R<x>:=PowerSeriesRing(Integers(), 33); Coefficients(R!( (1 - x)*(1 + x)*(1-x-x^2) / (1-2*x-x^2)^2 )); // Marius A. Burtea, Jan 15 2020

Crossrefs

Cf. A193737 (Fibonacci (with a(0)=1) triangle), A331319, A331320.

Sequence in context: A103819 A147484 A316980 * A017929 A017930 A305561

Adjacent sequences: A331318 A331319 A331320 * A331322 A331323 A331324

Keyword

nonn,easy

Author

Peter Luschny, Jan 14 2020

Status

approved