OFFSET
1,6
COMMENTS
k=4 case of the family of sequences whose k=1 case is the Fibonacci sequence A000045, k=2 case is the Padovan sequence A000931 (offset so as to begin 1,1,1) and k=3 case is A079398 (offset so as to begin 1,1,1,1).
The general case for integer k>1 is defined: a(1) = a(2) = ... = a(k+1) and for n>(k+1) a(n) = a(n-k) + a(n-[k+1]).
For this k=4 case, the ratio of successive terms a(n)/a(n-1) approaches the unique positive root of the irreducible characteristic polynomial: x^5 - x - 1 = 0, A160155.
REFERENCES
Zanten, A. J. van, The golden ratio in the arts of painting, building and mathematics, Nieuw Archief voor Wiskunde, 4 (17) (1999) 229-245.
LINKS
Indranil Ghosh, Table of n, a(n) for n = 1..14857
Jean-Paul Allouche and Tom Johnson, Narayana's cows and delayed morphisms, in G. Assayag, M. Chemillier, and C. Eloy, Troisièmes Journées d'Informatique Musicale, JIM '96, Île de Tatihou, France, 1996, pp. 2-7. [The hal link does not always work. - N. J. A. Sloane, Feb 19 2025]
Jean-Paul Allouche and Tom Johnson, Narayana's cows and delayed morphisms, in G. Assayag, M. Chemillier, and C. Eloy, Troisièmes Journées d'Informatique Musicale, JIM '96, Île de Tatihou, France, 1996, pp. 2-7. [Local copy with annotations and a correction from N. J. A. Sloane, Feb 19 2025]
Hùng Việt Chu and Zachary Louis Vasseur, Linear Recurrences of Generalized Schreier Sets Revisited, J. Int. Seq. 29 (2026), Article 26.2.2. See p. 4 (Table 3).
Richard Padovan, Dom Hans van der Laan and the Plastic Number.
Ernst S. Selmer, On the irreducibility of certain trinomials, Math. Scand., 4 (1956) 287-302.
Jeffrey Shallit, A generalization of automatic sequences, Theor. Comp. Sci. 61 (1988), 1-16.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,1,1).
FORMULA
G.f. -x*(1+x)*(1+x^2) / ( -1+x^4+x^5 ). - R. J. Mathar, Aug 26 2011
EXAMPLE
a(14) = 5 because a(14) = a(14-4) + a(14-5) = a(10) + a(9) = 3 + 2 = 5.
MATHEMATICA
k = 4; Do[a[n] = 1, {n, k + 1}]; a[n_] := a[n] = a[n - k] + a[n - k - 1]; Array[a, 61]
LinearRecurrence[{0, 0, 0, 1, 1}, {1, 1, 1, 1, 1}, 70] (* Harvey P. Dale, Apr 22 2015 *)
PROG
(PARI) a(n)=([0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1; 1, 1, 0, 0, 0]^(n-1)*[1; 1; 1; 1; 1])[1, 1] \\ Charles R Greathouse IV, Oct 03 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jonathan Vos Post, Feb 03 2005
EXTENSIONS
Edited by Ray Chandler and Robert G. Wilson v, Feb 06 2005
STATUS
approved
