|
%I #70 Feb 07 2023 07:23:55
%S 1,5,7,17,11,35,15,49,34,55,23,119,27,75,77,129,35,170,39,187,105,115,
%T 47,343,86,135,142,255,59,385,63,321,161,175,165,578,75,195,189,539,
%U 83,525,87,391,374,235,95,903,162,430,245,459,107,710,253,735,273,295,119
%N If n = Product p_i^e_i then a(n) = Product (1 + 2*p_i + 3*p_i^2 + ... + (e_i+1)*p_i^e_i).
%C Equals row sums of triangle A143313. - _Gary W. Adamson_, Aug 06 2008
%C Equals row sums of triangle A127099. - _Gary W. Adamson_, Jul 27 2008
%C Sum of the divisors d2 from the ordered pairs of divisors of n, (d1,d2) with d1<=d2, such that d1|d2. - _Wesley Ivan Hurt_, Mar 22 2022
%D D. M. Burton, Elementary Number Theory, Allyn and Bacon Inc., Boston, MA, 1976, p. 120.
%H Seiichi Manyama, <a href="/A060640/b060640.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Harry J. Smith)
%F a(n) = Sum_{d|n} d*tau(d), where tau(d) is the number of divisors of d, cf. A000005. a(n) = Sum_{d|n} d*sigma(n/d), where sigma(n)=sum of divisors of n, cf. A000203. - _Vladeta Jovovic_, Apr 23 2001
%F Multiplicative with a(p^e) = ((e+1)*p^{e+2} - (e+2)*p^{e+1} + 1) / (p-1)^2. Dirichlet g.f.: zeta(s)*zeta(s-1)^2. - _Franklin T. Adams-Watters_, Aug 03 2006
%F L.g.f.: Sum(A060640(n)*x^n/n) = -log( Product_{j>=1} P(x^j) ) where P(x) = Product_{k>=1} (1-x^k). - _Joerg Arndt_, May 03 2008
%F G.f.: Sum_{k>=1} k*tau(k)*x^k/(1 - x^k). - _Ilya Gutkovskiy_, Sep 06 2018
%F Sum_{k=1..n} a(k) ~ n^2/24 * ((4*gamma - 1)*Pi^2 + 2*Pi^2 * log(n) + 12*Zeta'(2)), where gamma is the Euler-Mascheroni constant A001620. - _Vaclav Kotesovec_, Feb 01 2019
%e a(4) = a(2^2) = 1 + (2)*(2) + (3)*(2^2) = 17;
%e a(6) = a(2)*a(3) = (1 + (2)*(2))*(1+(2)*(3)) = (5)*(7) = 35.
%e a(6) = tau(1) + 2*tau(2) + 3*tau(3) + 6*tau(6) = 1 + 2*2 + 3*2 + 6*4 = 35.
%p A060640 := proc(n) local ans, i, j; ans := 1: for i from 1 to nops(ifactors(n)[2]) do ans := ans*(1+sum((j+1)*ifactors(n)[2][i][1]^j,j=1..ifactors(n)[2][i][2])): od: RETURN(ans) end:
%t a[n_] := Total[#*DivisorSigma[1, n/#] & /@ Divisors[n]];
%t a /@ Range[59] (* _Jean-François Alcover_, May 19 2011, after _Vladeta Jovovic_ *)
%t f[p_, e_] := ((e + 1)*p^(e + 2) - (e + 2)*p^(e + 1) + 1)/(p - 1)^2; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Apr 10 2022 *)
%o (PARI) j=[]; for(n=1,200,j=concat(j,sumdiv(n,d,n/d*sigma(d)))); j
%o (PARI) a(n)=if(n<1,0,direuler(p=2,n,1/(1-X)/(1-p*X)^2)[n]) /* _Ralf Stephan_ */
%o (PARI) N=66; default(seriesprecision,N); x=z+O(z^(N+1))
%o c=sum(j=1,N,j*x^j); t=1/prod(j=1,N, eta(x^(j)));
%o t=log(t);t=serconvol(t,c);
%o Vec(t) /* _Joerg Arndt_, May 03 2008 */
%o (PARI) { for (n=1, 1000, write("b060640.txt", n, " ", direuler(p=2, n, 1/(1 - X)/(1 - p*X)^2)[n]); ) } /* _Harry J. Smith_, Jul 08 2009 */
%o (Haskell)
%o a060640 n = sum [d * a000005 d | d <- a027750_row n]
%o -- _Reinhard Zumkeller_, Feb 29 2012
%o (Sage)
%o def A060640(n) :
%o sigma = sloane.A000203
%o return add(sigma(k)*(n/k) for k in divisors(n))
%o [A060640(i) for i in (1..59)] # _Peter Luschny_, Sep 15 2012
%Y Cf. A000005, A000203, A001001, A006171, A038040 (Mobius trans), A049060, A057660, A057723.
%Y Cf. A143313. - _Gary W. Adamson_, Aug 06 2008
%Y Cf. A127099.
%Y Cf. A027750.
%K nonn,nice,easy,mult
%O 1,2
%A _N. J. A. Sloane_, Apr 17 2001
%E More terms from _James A. Sellers_, _Vladeta Jovovic_ and _Matthew Conroy_, Apr 17 2001
|
