OFFSET
1,8
COMMENTS
k=6 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), k=3 case is A079398 (offset so as to begin 1,1,1,1), k=4 case is A103372 and k=5 case is A103373.
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=6 case, the ratio of successive terms a(n)/a(n-1) approaches the unique positive root of the characteristic polynomial: x^7 - x - 1 = 0. This is the real constant 1.1127756842787... (see A230160).
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
G. C. Greubel, Table of n, a(n) for n = 1..1000
Richard Padovan, Dom Hans van der Laan and the Plastic Number.
J.-P. Allouche and T. Johnson, Narayana's Cows and Delayed Morphisms
E. S. Selmer, On the irreducibility of certain trinomials, Math. Scand., 4 (1956) 287-302.
J. Shallit, A generalization of automatic sequences, Theoretical Computer Science, 61 (1988) 1-16.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,1,1).
FORMULA
G.f.: x*(1+x)*(1+x+x^2)*(x^2-x+1) / ( 1-x^6-x^7 ). - R. J. Mathar, Aug 26 2011
EXAMPLE
a(32) = 17 because a(32) = a(32-6) + a(32-7) = a(26) + a(25) = 9 + 8 = 17.
MATHEMATICA
k = 6; Do[a[n] = 1, {n, k + 1}]; a[n_] := a[n] = a[n - k] + a[n - k - 1]; Array[a, 70]
RecurrenceTable[{a[n] == a[n - 6] + a[n - 7], a[1] == a[2] == a[3] == a[4] == a[5] == a[6] == a[7] == 1}, a, {n, 70}] (* or *)
Rest@ CoefficientList[Series[-x (1 + x) (1 + x + x^2) (x^2 - x + 1)/(-1 + x^6 + x^7), {x, 0, 70}], x] (* Michael De Vlieger, Oct 03 2016 *)
LinearRecurrence[{0, 0, 0, 0, 0, 1, 1}, {1, 1, 1, 1, 1, 1, 1}, 80] (* Harvey P. Dale, Sep 02 2024 *)
PROG
(PARI) a(n)=([0, 1, 0, 0, 0, 0, 0; 0, 0, 1, 0, 0, 0, 0; 0, 0, 0, 1, 0, 0, 0; 0, 0, 0, 0, 1, 0, 0; 0, 0, 0, 0, 0, 1, 0; 0, 0, 0, 0, 0, 0, 1; 1, 1, 0, 0, 0, 0, 0]^(n-1)*[1; 1; 1; 1; 1; 1; 1])[1, 1] \\ Charles R Greathouse IV, Oct 03 2016
(PARI) x='x+O('x^50); Vec(x*(1+x)*(1+x+x^2)*(x^2-x+1)/(1-x^6-x^7)) \\ G. C. Greubel, May 01 2017
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