OFFSET
0,3
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Peter Luschny, Multifactorials
Eric Weisstein's World of Mathematics, Multifactorial.
FORMULA
a(n) = 1 for n < 1, otherwise a(n) = n*a(n-5).
Sum_{n>=0} 1/a(n) = A288092. - Amiram Eldar, Nov 10 2020
EXAMPLE
a(12) = 168 because 12*a(12-5) = 12*a(7) = 12*14 = 168.
MAPLE
a:= n-> `if`(n < 1, 1, n*a(n-5)) end proc; seq(a(n), n = 0..40); # G. C. Greubel, Aug 18 2019
MATHEMATICA
a[n_]:= If[n<1, 1, n*a[n-5]]; Table[a[n], {n, 0, 40}] (* G. C. Greubel, Aug 18 2019 *)
Table[Times@@Range[n, 1, -5], {n, 0, 40}] (* Harvey P. Dale, May 12 2020 *)
PROG
(PARI) a(n)=if(n<1, 1, n*a(n-5))
for(n=0, 50, print1(a(n), ", ")) \\ Herman Jamke (hermanjamke(AT)fastmail.fm), Oct 19 2006
(Magma)
b:= func< n | (n lt 6) select n else n*Self(n-5) >;
[1] cat [b(n): n in [1..40]]; // G. C. Greubel, Aug 18 2019
(Sage)
def a(n):
if (n<1): return 1
else: return n*a(n-5)
[a(n) for n in (0..40)] # G. C. Greubel, Aug 18 2019
(GAP)
a:= function(n)
if n<1 then return 1;
else return n*a(n-5);
fi;
end;
List([0..40], n-> a(n) ); # G. C. Greubel, Aug 18 2019
(Python)
def A085157(n):
if n <= 0:
return 1
else:
return n*A085157(n-5)
n = 0
while n <= 40:
print(n, A085157(n))
n = n+1 # A.H.M. Smeets, Aug 18 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Jun 21 2003
STATUS
approved