cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A081405 a(n) = (n+1)*a(n-2) with a(0) = a(1) = 1.

This page as a plain text file.
%I A081405 #37 Sep 08 2022 08:45:09
%S A081405 1,1,3,4,15,24,105,192,945,1920,10395,23040,135135,322560,2027025,
%T A081405 5160960,34459425,92897280,654729075,1857945600,13749310575,
%U A081405 40874803200,316234143225,980995276800,7905853580625,25505877196800,213458046676875
%N A081405 a(n) = (n+1)*a(n-2) with a(0) = a(1) = 1.
%C A081405 A001147 and A002866 combined.
%H A081405 G. C. Greubel, <a href="/A081405/b081405.txt">Table of n, a(n) for n = 0..790</a>
%F A081405 a(0)=a(1)=1; a(2n) = A001147(2*n-2) odd terms, double factorial numbers; a(2n-1) = A002866(n) = 2^(n-1)*n!
%F A081405 0 = a(n)*(a(n+1) - a(n+3)) + a(n+1)*a(n+2) if n>=0. - _Michael Somos_, Jan 24 2014
%F A081405 a(n) = (n-1)-st term of column 1 of the array at A249159, for n >= 0. - _Clark Kimberling_, Oct 23 2014
%e A081405 G.f. = 1 + x + 3*x^2 + 4*x^3 + 15*x^4 + 24*x^5 + 105*x^6 + 192*x^7 + ...
%p A081405 a[0]:=1:a[1]:=1:for n from 2 to 50 do a[n]:=(a[n-2]*(n+1)^2) od: seq(sqrt(a[n]), n=0..26); # _Zerinvary Lajos_, Mar 04 2008
%t A081405 f[n_]:= (n+1)*f[n-2]; f[0] = 1; f[1] = 1; Table[f[n], {n, 1, 30}]
%t A081405 a[ n_]:= If[ n < 0, 0, If[OddQ[n], 2^((n-1)/2) ((n+1)/2)!, (n+1)!!]]; (* _Michael Somos_, Jan 24 2014 *)
%t A081405 RecurrenceTable[{a[0]==a[1]==1,a[n]==(n+1)a[n-2]},a,{n,30}] (* _Harvey P. Dale_, Nov 05 2021 *)
%o A081405 (PARI) {a(n) = if( n<2, n>=0, (n+1) * a(n-2))}; /* _Michael Somos_, Jan 24 2014 */
%o A081405 (PARI) {a(n) = if( n<0, 0, if( n%2, 2^(n\2) * (n\2 + 1)!, (n+1)! / (2^(n\2) * (n\2)!)))}; /* _Michael Somos_, Jan 24 2014 */
%o A081405 (Magma) [n le 1 select 1 else (n+1)*Self(n-1): n in [0..30]]; // _Vincenzo Librandi_, Oct 26 2014
%o A081405 (Sage)
%o A081405 def a(n):
%o A081405     if n<2: return 1
%o A081405     else: return (n+1)*a(n-2)
%o A081405 [a(n) for n in (0..30)] # _G. C. Greubel_, Aug 24 2019
%o A081405 (GAP)
%o A081405 a:= function(n)
%o A081405     if n<2 then return 1;
%o A081405     else return (n+1)*a(n-2);
%o A081405     fi;
%o A081405   end;
%o A081405 List([0..30], n-> a(n) ); # _G. C. Greubel_, Aug 24 2019
%Y A081405 Cf. A000142, A001147, A002866.
%K A081405 nonn
%O A081405 0,3
%A A081405 _Labos Elemer_, Apr 01 2003