|
|
A156626
|
|
a(0)=1; a(1)=2; a(2)=6; a(n+1) = 2*(n+1)*a(n) - n^2*a(n-1), n > 1.
|
|
1
|
|
|
1, 2, 6, 28, 170, 1252, 10774, 105764, 1164298, 14188468, 189461222, 2749300084, 43058394154, 723619035908, 12984464393398, 247704600763972, 5005042735932554, 106759075226130004, 2396869357456172038, 56491095210068416148, 1394373970361058540202
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
MAPLE
|
if n<=1 then
n+1 ;
elif n = 2 then
6 ;
else
2*n*procname(n-1)-(n-1)^2*procname(n-2) ;
end if;
end proc:
|
|
MATHEMATICA
|
Join[{1}, RecurrenceTable[{a[n] == 2*n*a[n-1] - (n-1)^2*a[n-2], a[1] == 2, a[2] == 6}, a, {n, 1, 50}]] (* G. C. Greubel, Sep 01 2018 *)
|
|
PROG
|
(PARI) m=30; v=concat([2, 6], vector(m-2)); for(n=3, m, v[n] = 2*n*v[n-1]-(n-1)^2*v[n-2]); concat([1], v) \\ G. C. Greubel, Sep 01 2018
(Magma) I:=[2, 6]; [1] cat [n le 2 select I[n] else 2*n*Self(n-1) - (n-1)^2*Self(n-2): n in [1..30]]; // G. C. Greubel, Sep 01 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|