|
Search: id:A121989
|
|
|
| A121989 |
|
a(n) = (n - 1)*a(n - 2) - (n - 1)*a(n - 3) + n - 1; a(0) = a(1) = a(2). |
|
+0
1
|
|
| 1, 1, 1, 2, 3, 8, 10, 36, 21, 216, -126, 1960, -3751, 25044, -74230, 403144, -1489095, 7638000, -32168046, 164287728, -756314855, 3929115500, -19332654222, 103079467832, -535020703583, 2937890929320, -15952504285350, 90295702455504, -510040670796063, 2974949788743940
(list; graph; listen)
|
|
|
OFFSET
|
1,4
|
|
|
MATHEMATICA
|
a[n_] := a[n] = (n - 1)*a[n - 2] - (n - 1)*a[n - 3] + n - 1; a[0] = 1; a[1] = 1; a[2] = 1; Table[a[n], {n, 0, 30}] (* vector Matrix Markov*) M = {{0, 1, 0}, {0, 0, 1}, {-n - 1, n + 1, 0}}; v[0] = {1, 1, 1}; v[n_] := v[n] = M.v[n - 1] + {0, 0, 1 + n} a0 = Table[v[n][[1]], {n, 0, 30}]
|
|
CROSSREFS
|
Adjacent sequences: A121986 A121987 A121988 this_sequence A121990 A121991 A121992
Sequence in context: A063474 A025562 A107401 this_sequence A010786 A005727 A118089
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
Roger Bagula (rlbagulatftn(AT)yahoo.com), Sep 11 2006
|
|
EXTENSIONS
|
Edited by njas, Sep 24 2006
|
|
|
Search completed in 0.002 seconds
|
