OFFSET
0,6
FORMULA
MAPLE
a := proc(n) local L, p;
L := proc(n, p, r) local q, s, m; m:=n-r; q:=p-r; s:=0;
do if q > n then break fi; s := s+iquo(m, q);
q := q*p od; s end; mul(p^(L(n, p, 1)-L(n, p, 0)),
p = select(isprime, [$3..n])); end
MATHEMATICA
L[n_, p_, r_] := Module[{q, s, m}, m = n-r; q = p-r; s = 0; While[True, If[ q > n, Break[]]; s = s + Quotient[m, q]; q = q*p]; s];
a[n_] := Product[p^(L[n, p, 1]-L[n, p, 0]), {p, Select[Range[3, n], PrimeQ] }];
Table[a[n], {n, 0, 39}] (* Jean-François Alcover, Jun 17 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Jul 26 2009
STATUS
approved