OFFSET
1,3
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
For n = 4; pod(4) / tau(4) = 8 / 3; a(n) = 8.
MAPLE
f:= proc(n) local D; D:= numtheory:-divisors(n); numer(convert(D, `*`)/nops(D)) end proc:
map(f, [$1..100]); # Robert Israel, Sep 14 2017
MATHEMATICA
Table[Numerator[Apply[Times, Divisors@ n]/DivisorSigma[0, n]], {n, 51}] (* Michael De Vlieger, Sep 05 2017 *)
PROG
(Magma) [Numerator(&*[d: d in Divisors(n)] / #[d: d in Divisors(n)]): n in [1..1000]];
(PARI) a(n) = my(d=divisors(n)); numerator(prod(k=1, #d, d[k])/#d); \\ Michel Marcus, Sep 05 2017
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Jaroslav Krizek, Sep 05 2017
STATUS
approved