|
%I #22 Nov 21 2022 09:38:51
%S 1,3,7,73,21,21,43,39,757,63,111,73,157,129,147,151,273,2271,343,219,
%T 301,333,507,273,15751,471,511,3139,813,441,931,4161,777,819,903,
%U 55261,1333,1029,1099,819,1641,903,1807,8103,15897,1521,2163,1057,39331,47253
%N Numerator of sigma_3(n)/sigma(n).
%C If n is prime then a(n) = A002061(n). - _Robert Israel_, Jan 25 2018
%H Robert Israel, <a href="/A091259/b091259.txt">Table of n, a(n) for n = 1..10000</a>
%F Sum_{k=1..n} a(k)/A091258(k) ~ c * n^3, where c = (Pi^2/18)*zeta(3)^2 * Product_{p prime} (1 - 2/p^2 - 1/p^3 + 5/p^5 - 3/p^6) = 0.2382648075... . - _Amiram Eldar_, Nov 21 2022
%p seq(numer(numtheory:-sigma[3](n)/numtheory:-sigma(n)),n=1..100); # _Robert Israel_, Jan 25 2018
%t Array[Numerator[DivisorSigma[3,#]/DivisorSigma[1,#]]&,50] (* _Harvey P. Dale_, Feb 29 2016 *)
%o (PARI) a(n) = numerator(sigma(n, 3)/sigma(n)); \\ _Michel Marcus_, Jan 26 2018
%o (Magma) [Numerator(DivisorSigma(3,n)/DivisorSigma(1,n)): n in [1..50]]; // _Vincenzo Librandi_, Jan 26 2018
%Y Cf. A001158, A000203, A002061, A002117, A086463, A091258 (denominators), A091260.
%K easy,nonn,frac,look
%O 1,2
%A _Labos Elemer_, Feb 12 2004
|
