A163394 The odd part of Minkowski(n)/n!

1, 1, 1, 1, 1, 3, 1, 9, 9, 15, 3, 9, 3, 945, 135, 27, 27, 405, 45, 8505, 1701, 66825, 6075, 18225, 6075, 995085, 76545, 8505, 1215, 18225, 1215, 841995, 841995, 6506325, 382725, 32805, 3645, 850797675, 44778825, 3444525

Offset: 0

Author

Peter Luschny, Jul 26 2009

Keywords

• nonn

Crossrefs

Cf. A000265, A053657, A163176.

Programs

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 *)

Formula

a(n) = A000265(A163176(n)). - Michel Marcus, Jun 17 2019