OFFSET
1,2
LINKS
Daniel Suteu, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = Sum_{d|n, gcd(n/d, d) = 1} sigma(d), where sigma(d) is the sum of the divisors of d. - Daniel Suteu, Jun 27 2018
Sum_{k=1..n} a(k) ~ c * n^2, where c = (Pi^2/12) * Product_{p prime} (1 + 1/p^2 - 1/p^3) = 1.0741158... . - Amiram Eldar, Nov 01 2022
EXAMPLE
a(6) = (2+2)*(2+3) = 20.
MAPLE
A107758 := proc(n) local pf, p ; if n = 1 then 1; else pf := ifactors(n)[2] ; mul( 1+(op(1, p)^(op(2, p)+1)-1)/(op(1, p)-1), p=pf) ; end if; end proc:
seq(A107758(n), n=1..60) ; # R. J. Mathar, Jan 07 2011
MATHEMATICA
Table[DivisorSum[n, DivisorSigma[1, #] &, CoprimeQ[n/#, #] &], {n, 54}] (* Michael De Vlieger, Jun 27 2018 *)
f[p_, e_] := 1 + (p^(e + 1) - 1)/(p - 1); a[n_] := Times @@ f @@@ FactorInteger[n]; a[1] = 1; Array[a, 100] (* Amiram Eldar, Aug 26 2022 *)
PROG
(PARI) a(n) = sumdiv(n, d, if(gcd(n/d, d) == 1, sigma(d))); \\ Daniel Suteu, Jun 27 2018
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Yasutoshi Kohmoto, May 25 2005
STATUS
approved