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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A107786 a(n) = |b(n)| where b(n) = -b(n-1) + 6*b(n-3) with b(0)=0, b(1)=1, b(2)=1. 1
0, 1, 1, 1, 7, 1, 5, 47, 53, 23, 259, 577, 715, 839, 4301, 8591, 3557, 22249, 73795, 95137, 38357, 481127, 1051949, 821807, 2064955, 8376649, 13307491, 917761, 49342133, 129187079, 134693645, 161359153, 936481627, 1744643497, 776488579 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,5

COMMENTS

Original name was: Sequence obtained using characteristic polynomial that is Laplace transform of the minimal Pisot characteristic polynomial: -s^4*L(t^3-t-1) = s^3 + s^2-6.

LINKS

Robert Israel, Table of n, a(n) for n = 0..3380

FORMULA

G.f. of b(n): (x+2*x^2)/(1+x-6*x^3). - Robert Israel, Jun 27 2018

MAPLE

f:= gfun:-rectoproc({a(n) = -a(n-1)+6*a(n-3), a(0)=0, a(1)=1, a(2)=1}, a(n), remember):

map(abs@f, [$0..50]); # Robert Israel, Jun 27 2018

MATHEMATICA

(* first method *) F[1] = 0; F[2] = 1; F[3] = 1; F[n__] := F[n] = -F[n - 1] + 6*F[n - 3] a = Table[Abs[F[n]], {n, 1, 50}] (*second method*) M = {{0, 1, 0}, {0, 0, 1}, {6, 0, -1}} v[1] = {0, 1, 1} v[n_] := v[n] = M.v[n - 1] a = Table[Abs[v[n][[1]]], {n, 1, 50}] Det[M - x*IdentityMatrix[3]]

Abs[LinearRecurrence[{-1, 0, 6}, {0, 1, 1}, 40]] (* Vincenzo Librandi, Jun 28 2016 *)

CROSSREFS

Sequence in context: A047875 A064467 A089204 * A246949 A154932 A026497

Adjacent sequences:  A107783 A107784 A107785 * A107787 A107788 A107789

KEYWORD

nonn

AUTHOR

Roger L. Bagula, Jun 11 2005

EXTENSIONS

Edited by Robert Israel, Jun 27 2018

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 November 12 17:06 EST 2019. Contains 329058 sequences. (Running on oeis4.)