login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A014292 a(n) = 2*a(n-1) - a(n-2) - a(n-4) with a(0) = a(1) = 0, a(2) = 1, a(3) = 2. 4
0, 0, 1, 2, 3, 4, 4, 2, -3, -12, -25, -40, -52, -52, -27, 38, 155, 324, 520, 678, 681, 360, -481, -2000, -4200, -6760, -8839, -8918, -4797, 6084, 25804, 54442, 87877, 115228, 116775, 63880, -76892, -332892, -705667, -1142322 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

COMMENTS

Imaginary part of the sequence of complex numbers defined by c(0) = 1, c(1) = 1, for n>1 c(n) = c(n-1) + i*c(n-2). - Gerald McGarvey, Apr 24 2005

a(n) = sqrt(3)*y where (x,y,y,y) is the quaternion b(n) of the sequence b of quaternions defined by b(0)=1,b(1)=1, b(n) = b(n-1) + b(n-2)*(0,s,s,s) where s = 1/sqrt(3). - Gerald McGarvey, Apr 25 2005

For n>=1, a(n) is equal to -1 times the imaginary part of the determinant of the n X n matrix with the sqrt(i)'s along the superdiagonal and the subdiagonal (i is the imaginary unit), 1's along the main diagonal, and 0's everywhere else (see Mathematica code below). - John M. Campbell, Jun 04 2011

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,-1).

FORMULA

a(n) = Sum_{k=0..floor((n+2)/2)} binomial(n-k+2, k)*sin(Pi*k/2). - Paul Barry, Apr 25 2005

G.f.: x^2/(1 - 2*x + x^2 + x^4). - R. J. Mathar, Oct 22 2008

MATHEMATICA

Table[-Im[Det[Array[KroneckerDelta[#1 + 1, #2]*Sqrt[I] &, {n, n}] + Array[KroneckerDelta[#1 - 1, #2]*Sqrt[I] &, {n, n}] + IdentityMatrix[n]]], {n, 1, 40}] (* John M. Campbell, Jun 04 2011 *)

LinearRecurrence[{2, -1, 0, -1}, {0, 0, 1, 2}, 40] (* G. C. Greubel, Jun 12 2019 *)

PROG

(PARI) my(x='x+O('x^40)); concat([0, 0], Vec(x^2/(1-2*x+x^2+x^4))) \\ G. C. Greubel, Jun 12 2019

(MAGMA) R<x>:=PowerSeriesRing(Integers(), 40); [0, 0] cat Coefficients(R!( x^2/(1-2*x+x^2+x^4) )); // G. C. Greubel, Jun 12 2019

(Sage) (x^2/(1-2*x+x^2+x^4)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Jun 12 2019

(GAP) a:=[0, 0, 1, 2];; for n in [5..40] do a[n]:=2*a[n-1]-a[n-2]-a[n-4]; od; a; # G. C. Greubel, Jun 12 2019

CROSSREFS

Cf. A104862, A104862.

Sequence in context: A160386 A307310 A174015 * A244445 A066078 A058339

Adjacent sequences:  A014289 A014290 A014291 * A014293 A014294 A014295

KEYWORD

sign

AUTHOR

N. J. A. Sloane

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 25 10:20 EDT 2019. Contains 326324 sequences. (Running on oeis4.)