|
|
A000322
|
|
Pentanacci numbers: a(n) = a(n-1) + a(n-2) + a(n-3) + a(n-4) + a(n-5) with a(0) = a(1) = a(2) = a(3) = a(4) = 1.
(Formerly M3786 N1542)
|
|
53
|
|
|
1, 1, 1, 1, 1, 5, 9, 17, 33, 65, 129, 253, 497, 977, 1921, 3777, 7425, 14597, 28697, 56417, 110913, 218049, 428673, 842749, 1656801, 3257185, 6403457, 12588865, 24749057, 48655365, 95653929, 188050673, 369697889, 726806913, 1428864769
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,6
|
|
COMMENTS
|
For n>=0: a(n+2) is the number of length-n words with letters {0,1,2,3,4} where the letter x is followed by at least x zeros, see fxtbook link below. - Joerg Arndt, Apr 08 2011
|
|
REFERENCES
|
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
B. G. Baumgart, Letter to the editor Part 1 Part 2 Part 3, Fib. Quart. 2 (1964), 260, 302.
|
|
MAPLE
|
A000322:=(-1+z**2+2*z**3+3*z**4)/(-1+z**2+z**3+z+z**4+z**5); # Simon Plouffe in his 1992 dissertation.
a:= n-> (Matrix([[1$5]]). Matrix(5, (i, j)-> if (i=j-1) or j=1 then 1 else 0 fi)^n)[1, 5]: seq (a(n), n=0..28); # Alois P. Heinz, Aug 26 2008
|
|
MATHEMATICA
|
LinearRecurrence[{1, 1, 1, 1, 1}, {1, 1, 1, 1, 1}, 50]
|
|
PROG
|
(Magma) [ n le 5 select 1 else Self(n-1)+Self(n-2)+Self(n-3)+Self(n-4)+Self(n-5): n in [1..40] ];
(J) (see www.jsoftware.com) First construct the generating matrix
(((+ +/), ]), :^:(1=#@$))/&.|.<:/~i.5
1 1 1 1 1
1 2 2 2 2
2 3 4 4 4
4 6 7 8 8
8 12 14 15 16
Given that matrix, one can produce the first 2000 numbers in almost 17 millisecs by
, ((((+ +/), ]), :^:(1=#@$))/&.|.<:/~i.5) (+/ . *)^:(i.400) 1 1 1 1 1x
|
|
CROSSREFS
|
Cf. A001591 (Pentanacci numbers starting 0, 0, 0, 0, 1).
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|