|
|
A176513
|
|
a(n+5) = a(n+3) + a(n+2) + a(n), with a(1) = a(2) = a(3) = a(4) = a(5) = 1.
|
|
2
|
|
|
1, 1, 1, 1, 1, 3, 3, 5, 7, 9, 15, 19, 29, 41, 57, 85, 117, 171, 243, 345, 499, 705, 1015, 1447, 2065, 2961, 4217, 6041, 8625, 12323, 17627, 25165, 35991, 51417, 73479, 105035, 150061, 214505, 306513, 438045, 626053, 894619, 1278603, 1827185
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,6
|
|
COMMENTS
|
Lim_{n -> infinity} a(n+1)/a(n) = s = 1.4291083198381..., where s is the root of the characteristic equation s^5 = s^3 + s^2 + 1.
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (x^4 + x^3 - x - 1)/(x^5 + x^3 + x^2 - 1).
|
|
EXAMPLE
|
a(7) = a(5) + a(4) + a(2) = 3.
|
|
MAPLE
|
seq(coeff(series((x^4+x^3-x-1)/(x^5+x^3+x^2-1), x, n+1), x, n), n = 0 .. 50); # Muniru A Asiru, Feb 17 2019
|
|
MATHEMATICA
|
LinearRecurrence[{0, 1, 1, 0, 1}, {1, 1, 1, 1, 1}, 50] (* G. C. Greubel, Feb 17 2019 *)
|
|
PROG
|
(PARI) a(n) = my(v=vector(n)); for(k=1, n, v[k]=if(k<=5, 1, v[k-2]+v[k-3]+v[k-5])); v[n] \\ Jianing Song, Feb 04 2019
(PARI) my(x='x+O('x^50)); Vec((x^4+x^3-x-1)/(x^5+x^3+x^2-1)) \\ G. C. Greubel, Feb 17 2019
(Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!( (x^4+x^3-x-1)/(x^5+x^3+x^2-1) )); // G. C. Greubel, Feb 17 2019
(Sage) ((x^4+x^3-x-1)/(x^5+x^3+x^2-1)).series(x, 50).coefficients(x, sparse=False) # G. C. Greubel, Feb 17 2019
(GAP) a:=[1, 1, 1, 1, 1];; for n in [6..30] do a[n]:=a[n-2]+a[n-3]+a[n-5]; od; a; # G. C. Greubel, Feb 17 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|