OFFSET
1,4
LINKS
Robert Israel, Table of n, a(n) for n = 1..254
FORMULA
a(n) = a(n-2) + (n^2 - 3*n + 1)*a(n-1) with a(1)=0, a(2)=1, a(3)=1. - Robert Israel, Dec 18 2018
MAPLE
f:= gfun:-rectoproc({a(n) = a(n-2)+(n^2-3*n+1)*a(n-1), a(1)=0, a(2)=1}, a(n), remember):
map(f, [$1..30]); # Robert Israel, Dec 18 2018
MATHEMATICA
M[n_] := {{0, 1}, {1, n^2 - n - 1} }
v[1] = {0, 1}
v[n_] := v[n] = M[n].v[n - 1]
Table[v[n][[1]], {n, 1, 30}]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Roger L. Bagula, Sep 02 2006
EXTENSIONS
Edited by Robert Israel, Dec 18 2018
New name (using Robert Israel's formula) from Joerg Arndt, Jan 06 2022
STATUS
approved