A215500 a(n) = ((sqrt(5) + 3)^n + (-sqrt(5) -1)^n + (-sqrt(5) + 3)^n + (sqrt(5) - 1)^n) / 2^n.

4, 2, 10, 14, 54, 112, 340, 814, 2254, 5702, 15250, 39404, 104004, 270922, 711490, 1859134, 4873054, 12748472, 33391060, 87394454, 228841254, 599050102, 1568437210, 4106054164, 10750060804, 28143585362, 73681573690, 192899714414, 505019869254, 1322156172352

Offset: 0

Author

Peter Luschny, Aug 13 2012

Keywords

• nonn
• easy

Comments

Inverse binomial transform is (-1)^n * a(n). - Michael Somos, Jun 02 2014

Examples

			G.f. = 4 + 2*x + 10*x^2 + 14*x^3 + 54*x^4 + 112*x^5 + 340*x^6 + ...
		

Links

• G. C. Greubel, Table of n, a(n) for n = 0..1000
• Index entries for linear recurrences with constant coefficients, signature (2,3,-4,1).

Crossrefs

Cf. A051927, A215502, A215503.

Cf. A124697.

Programs

Magma

I:=[4,2,10,14]; [n le 4 select I[n] else 2*Self(n-1) + 3*Self(n-2) - 4*Self(n-3) + Self(n-4): n in [1..30]]; // G. C. Greubel, Apr 23 2018

Maple

A215500 := x -> ((sqrt(5)+3)^x+(-sqrt(5)-1)^x+(-sqrt(5)+3)^x+(sqrt(5)-1)^x)/2^x;
seq(simplify(A215500(i)),i=0..29);

Mathematica

a[n_] := ((Sqrt[5] + 3)^n + (-Sqrt[5] - 1)^n + (-Sqrt[5] + 3)^n + (Sqrt[5] - 1)^n)/2^n; Table[a[n] // Simplify, {n, 0, 29}] (* Jean-François Alcover, Jul 02 2013 *)
LinearRecurrence[{2,3,-4,1}, {4, 2, 10, 14}, 50] (* G. C. Greubel, Apr 23 2018 *)

PARI

{a(n) = polsym( (1 + (-1)^(n>0)*x - x^2) * (1 - 3*x + x^2), abs(n))[1 + abs(n)]}; /* Michael Somos, Jun 02 2014 */

Sage

def A215500(x) :
    return ((sqrt(5)+3)^x+(-sqrt(5)-1)^x+(-sqrt(5)+3)^x+(sqrt(5)-1)^x)/2^x
[A215500(i).round() for i in (0..29)]

Formula

G.f.: 2*(2-x)*(1+x)*(1-2*x)/((1-3*x+x^2)*(1+x-x^2)). - Colin Barker, Aug 19 2012

a(n) = 2*a(n-1)+3*a(n-2)-4*a(n-3)+a(n-4). - Colin Barker, Aug 20 2012

a(-n) = A124697(n) if n>0. - Michael Somos, Jun 02 2014