login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A126116 a(n) = a(n-1) + a(n-3) + a(n-4), with a(0)=a(1)=a(2)=a(3)=1. 6
1, 1, 1, 1, 3, 5, 7, 11, 19, 31, 49, 79, 129, 209, 337, 545, 883, 1429, 2311, 3739, 6051, 9791, 15841, 25631, 41473, 67105, 108577, 175681, 284259, 459941, 744199, 1204139, 1948339, 3152479, 5100817, 8253295, 13354113, 21607409, 34961521 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,5

COMMENTS

This sequence has the same growth rate as the Fibonacci sequence, since x^4 - x^3 - x - 1 has the real roots phi and -1/phi.

The Ca1 sums, see A180662 for the definition of these sums, of triangle A035607 equal the terms of this sequence without the first term. - Johannes W. Meijer, Aug 05 2011

REFERENCES

Wolfram, S., A New Kind of Science. Champaign, IL: Wolfram Media, pp. 82-92, 2002

LINKS

Seiichi Manyama, Table of n, a(n) for n = 0..4786

K. T. Atanassov, D. R. Deford, A. G. Shannon, Pulsated Fibonacci recurrences, Fibonacci Quarterly, Vol. 52, No. 5, Dec. 2014, pp. 22-27.

Kelley L. Ross, The Golden Ratio and The Fibonacci Numbers

Eric Weisstein's World of Mathematics, MathWorld: Golden Ratio

Wikipedia, Golden Ratio

Index entries for linear recurrences with constant coefficients, signature (1,0,1,1). [From R. J. Mathar, Jul 22 2010]

FORMULA

From R. J. Mathar, Jul 22 2010: (Start)

G.f.: (1-x)*(1+x+x^2)/((1-x-x^2)*(1+x^2)).

a(n) = ( (-1)^floor(n/2) * A010684(n) + 2*A000032(n))/5.

a(2*n) = A061646(n). (End)

From Johannes W. Meijer, Aug 05 2011: (Start)

a(n) = F(n-1) + A070550(n-4) with F(n) = A000045(n).

a(n) = F(n-1) + F(floor((n-4)/2) + 1)*F(ceiling((n-4)/2) + 2). (End)

a(n) = 1/5*((sqrt(5)-1)*(1/2*(1+sqrt(5)))^n - (1+sqrt(5))*(1/2*(1-sqrt(5)))^n + sin((Pi*n)/2) - 3*cos((Pi*n)/2)). - Harvey P. Dale, Nov 08 2011

(-1)^n * a(-n) = a(n) = F(n) - A070550(n - 6). - Michael Somos, Feb 05 2012

EXAMPLE

G.f. = 1 + x + x^2 + x^3 + 3*x^4 + 5*x^5 + 7*x^6 + 11*x^7 + 19*x^8 + 31*x^9 + ...

MAPLE

# From R. J. Mathar, Jul 22 2010: (Start)

A010684 := proc(n) 1+2*(n mod 2) ; end proc:

A000032 := proc(n) coeftayl((2-x)/(1-x-x^2), x=0, n) ; end proc:

A126116 := proc(n) ((-1)^floor(n/2)*A010684(n)+2*A000032(n))/5 ; end proc: seq(A126116(n), n=0..80) ; # (End)

with(combinat): A126116 := proc(n): fibonacci(n-1) + fibonacci(floor((n-4)/2)+1)* fibonacci(ceil((n-4)/2)+2) end: seq(A126116(n), n=0..38); # Johannes W. Meijer, Aug 05 2011

MATHEMATICA

LinearRecurrence[{1, 0, 1, 1}, {1, 1, 1, 1}, 50] (* Harvey P. Dale, Nov 08 2011 *)

PROG

(PARI) Vec((x-1)*(1+x+x^2)/((x^2+x-1)*(x^2+1)) + O(x^50)) \\ Altug Alkan, Dec 25 2015

(MAGMA) [n le 4 select 1 else Self(n-1) + Self(n-3) + Self(n-4): n in [1..50]]; // Vincenzo Librandi, Dec 25 2015

(Sage) ((1-x)*(1+x+x^2)/((1-x-x^2)*(1+x^2))).series(x, 50).coefficients(x, sparse=False) # G. C. Greubel, Jul 15 2019

(GAP) a:=[1, 1, 1, 1];; for n in [5..50] do a[n]:=a[n-1]+a[n-3]+a[n-4]; od; a; # G. C. Greubel, Jul 15 2019

CROSSREFS

Cf. Fibonacci numbers A000045; Lucas numbers A000032; Tribonacci numbers A000213; Tetranacci numbers A000288; Pentanacci numbers A000322; Hexanacci numbers A000383; 7th-order Fibonacci numbers A060455; Octanacci numbers A079262; 9th-order Fibonacci sequence A127193; 10th-order Fibonacci sequence A127194; 11th-order Fibonacci sequence A127624, A128429.

Sequence in context: A175196 A077858 A327461 * A161423 A133846 A056208

Adjacent sequences:  A126113 A126114 A126115 * A126117 A126118 A126119

KEYWORD

nonn

AUTHOR

Luis A Restrepo (luisiii(AT)mac.com), Mar 05 2007

EXTENSIONS

Edited by Don Reble, Mar 09 2007

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 July 6 21:41 EDT 2020. Contains 335483 sequences. (Running on oeis4.)