%I #17 Sep 08 2022 08:46:04
%S 0,0,1,2,3,6,13,32,85,246,757,2480,8537,30858,116229,455668,1850417,
%T 7774102,33679941,150291472,689170529,3244125554,15649195077,
%U 77287580604,390271482145,2013310674830,10599283282021
%N a(0)=0, a(1)=0; for n>1, a(n) = a(n-1) + (n-4)*a(n-2) + 1
%H Vincenzo Librandi, <a href="/A220699/b220699.txt">Table of n, a(n) for n = 0..800</a>
%F a(0)=a(1)=0, a(2)=1, a(n) = 2*a(n-1)+(n-5)*a(n-2)-(n-5)*a(n-3).
%F a(n) ~ (sqrt(Pi)/2+sqrt(2)) * n^(n/2-3/2)*exp(sqrt(n)-n/2-1/4) * (1-29/(24*sqrt(n))). - _Vaclav Kotesovec_, Dec 27 2012
%t RecurrenceTable[{a[1]==0, a[2]==0, a[n]==a[n-1] + (n-5) a[n-2] + 1}, a, {n, 40}]
%o (Magma) [n le 2 select 0 else Self(n-1)+(n-5)*Self(n-2) + 1: n in [1..30]];
%Y Cf. A185108, A185109, A185308, A185309, A186738, A186739, A213720, A187830.
%K nonn,easy
%O 0,4
%A _Vincenzo Librandi_, Dec 25 2012