|
|
A323780
|
|
a(n) = denominator of Sum_{d|n} (tau(d)/sigma(d)) where tau(k) = the number of the divisors of k (A000005) and sigma(k) = the sum of the divisors of k (A000203).
|
|
5
|
|
|
1, 3, 2, 21, 3, 2, 4, 105, 26, 9, 6, 7, 7, 12, 1, 3255, 9, 26, 10, 63, 8, 18, 12, 35, 93, 21, 65, 21, 15, 3, 16, 1085, 4, 27, 3, 91, 19, 6, 7, 315, 21, 8, 22, 9, 13, 36, 24, 2170, 76, 279, 3, 147, 27, 39, 9, 21, 20, 9, 30, 21, 31, 48, 104, 137795, 21, 12, 34
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Sum_{d|n} (tau(d)/sigma(d)) >= 1 for all n >= 1.
|
|
LINKS
|
|
|
FORMULA
|
a(p) = (p+1) / gcd(p+3, p+1) for p = primes p.
|
|
EXAMPLE
|
For n=4; Sum_{d|4} (tau(d)/sigma(d)) = (tau(1)/sigma(1))+(tau(2)/sigma(2))+(tau(4)/sigma(4)) = (1/1)+(2/3)+(3/7) = 44/21; a(4) = 21.
|
|
MATHEMATICA
|
Array[Denominator@ DivisorSum[#, Divide @@ DivisorSigma[{0, 1}, #] &] &, 67] (* Michael De Vlieger, Feb 15 2019 *)
|
|
PROG
|
(Magma) [Denominator(&+[NumberOfDivisors(d) / SumOfDivisors(d): d in Divisors(n)]): n in [1..100]]
(PARI) a(n) = denominator(sumdiv(n, d, numdiv(d)/sigma(d))); \\ Michel Marcus, Feb 13 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,frac
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|