OFFSET
1,2
COMMENTS
Sum_{d|n} (pod(d)/sigma(d)) for n >= 1: 1, 5/3, 7/4, 59/21, 11/6, 65/12, 15/8, 743/105, ...
Sum_{d|n} (pod(d)/sigma(d)) > 1 for all n > 1.
a(p) = 2p + 1 for p = primes (A000040).
EXAMPLE
For n=4: Sum_{d|4} (pod(d)/sigma(d)) = pod(1)/sigma(1) + pod(2)/sigma(2) + pod(4)/sigma(4) = 1/1 + 2/3 + 8/7 = 59/21; a(4) = 59.
MATHEMATICA
Array[Numerator@ DivisorSum[#, Apply[Times, Divisors@ #]/DivisorSigma[1, #] &] &, 46] (* Michael De Vlieger, Mar 24 2019 *)
PROG
(Magma) [Numerator(&+[&*[c: c in Divisors(d)] / SumOfDivisors(d): d in Divisors(n)]): n in [1..100]]
(PARI) a(n) = numerator(sumdiv(n, d, vecprod(divisors(d))/sigma(d))); \\ Michel Marcus, Mar 23 2019
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Jaroslav Krizek, Mar 22 2019
STATUS
approved