|
| |
|
|
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, 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; refs; listen; history; internal format)
|
|
|
|
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}]
RecurrenceTable[{a[0]==a[1]==a[2]==1, a[n]==(n-1)a[n-2]-(n-1)a[n-3]+ n-1}, a, {n, 30}] (* From Harvey P. Dale, Dec 08 2011 *)
|
|
|
CROSSREFS
| Sequence in context: A025562 A107401 A165153 * A010786 A005727 A118089
Adjacent sequences: A121986 A121987 A121988 * A121990 A121991 A121992
|
|
|
KEYWORD
| sign
|
|
|
AUTHOR
| Roger Bagula (rlbagulatftn(AT)yahoo.com), Sep 11 2006
|
|
|
EXTENSIONS
| Edited by N. J. A. Sloane (njas(AT)research.att.com), Sep 24 2006
|
| |
|
|
