OFFSET
1,2
COMMENTS
Let rad(n) be the radical of n, which equals the product of all prime factors of n (A007947). Let g(n) = 1/rad(n) and let f(n) = Sum_{d | n} g(d). This is a multiplicative function whose value on a prime power is f(p^k) = 1 + k/p. Hence f is a weighted divisor-counting function that weights divisors d higher when they have few and small prime divisors themselves. The sequence a(n) lists the numerators of the fractions f(n) in lowest terms.
If p is prime, then a(p^k) = p+k if p does not divide k, 1 + k/p if it does. In particular, a(p^p) = 2. - Robert Israel, Jan 25 2019
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)/A322966(k) = zeta(2)*zeta(3)/zeta(6) (A082695). - Amiram Eldar, Dec 09 2023
EXAMPLE
The divisors of 12 are 1,2,3,4,6,12, so f(12) = 1 + (1/2) + (1/3) + (1/2) + (1/6) + (1/6) = 8/3 and a(12) = 8. Alternately, since f is multiplicative, f(12) = f(4)*f(3) = (1+2/2)*(1+1/3) = 8/3.
MAPLE
rad:= n -> convert(numtheory:-factorset(n), `*`):
f:= proc(n) numer(add(1/rad(d), d=numtheory:-divisors(n))) end proc:
map(f, [$1..100]); # Robert Israel, Jan 25 2019
MATHEMATICA
Array[Numerator@ DivisorSum[#, 1/Apply[Times, FactorInteger[#][[All, 1]]] &] &, 71] (* Michael De Vlieger, Jan 19 2019 *)
PROG
(PARI) rad(n) = factorback(factor(n)[, 1]); \\ A007947
a(n) = numerator(sumdiv(n, d, 1/rad(d))); \\ Michel Marcus, Jan 10 2019
CROSSREFS
KEYWORD
frac,nonn
AUTHOR
David S. Metzler, Dec 31 2018
EXTENSIONS
More terms from Michel Marcus, Jan 19 2019
STATUS
approved