cp's OEIS Frontend

A075374 a(n+2) = n*a(n+1) - a(n), with a(1)=1, a(2)=2.

%I A075374 #23 Mar 05 2022 01:44:43
%S A075374 1,2,1,0,-1,-4,-19,-110,-751,-5898,-52331,-517412,-5639201,-67153000,
%T A075374 -867349799,-12075744186,-180268812991,-2872225263670,-48647560669399,
%U A075374 -872783866785512,-16534245908255329,-329812134298321068,-6909520574356487099,-151679640501544395110
%N A075374 a(n+2) = n*a(n+1) - a(n), with a(1)=1, a(2)=2.
%C A075374 Starting with offset 5 unsigned: (1, 4, 19, 110, 751, ...) = eigensequence of triangle A003991. - _Gary W. Adamson_, May 17 2010
%H A075374 Seiichi Manyama, <a href="/A075374/b075374.txt">Table of n, a(n) for n = 1..451</a>
%F A075374 a(n+1) = (a(n) + a(n+2))/n with a(1) = 1, a(2) = 2.
%p A075374 a[1] := 1:a[2] := 2:for n from 1 to 45 do a[n+2] := n*a[n+1]-a[n]:od:seq(a[i],i=1..45);
%t A075374 a[n_]:= a[n]= If[n<3, n, (n-2)*a[n-1] -a[n-2]];
%t A075374 Table[a[n], {n,50}] (* _G. C. Greubel_, Mar 04 2022 *)
%o A075374 (Magma) [n le 2 select n else (n-2)*Self(n-1) - Self(n-2): n in [1..50]]; // _G. C. Greubel_, Mar 04 2022
%o A075374 (Sage)
%o A075374 @CachedFunction
%o A075374 def a(n): return n if (n<3) else (n-2)*a(n-1) - a(n-2) # A075374
%o A075374 [a(n) for n in (1..50)] # _G. C. Greubel_, Mar 04 2022
%Y A075374 Cf. A003991. - _Gary W. Adamson_, May 17 2010
%K A075374 sign
%O A075374 1,2
%A A075374 _Amarnath Murthy_, Sep 20 2002
%E A075374 More terms from _Sascha Kurz_, Jan 30 2003