|
%I #31 Jul 06 2023 01:51:07
%S 1,6,10,28,21,78,36,120,91,171,78,406,105,300,300,496,171,780,210,903,
%T 528,666,300,1830,496,903,820,1596,465,2628,528,2016,1176,1485,1176,
%U 4186,741,1830,1596,4095,903,4656,990,3570,3081,2628,1176,7750,1653,4371
%N a(n) = sum of numbers from 1 to sigma(n), where sigma(n) = A000203(n).
%H Antti Karttunen, <a href="/A184387/b184387.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>.
%F a(n) = Sum_{i = 1..sigma(n)} i = A000217(A000203(n)) = A000203(n)*(A000203(n) + 1)/2.
%F Sum_{k=1..n} a(k) = (5*zeta(3)/12) * n^3 + O(n^2*log(n)^2). - _Amiram Eldar_, Dec 08 2022
%e For n = 4; sigma(4) = 7; a(4) = 1+2+3+4+5+6+7 = 28.
%t Array[PolygonalNumber@ DivisorSigma[1, #] &, 50] (* _Michael De Vlieger_, Nov 16 2017 *)
%o (PARI) a(n)=binomial(sigma(n)+1,2) \\ _Charles R Greathouse IV_, Feb 14 2013
%o (Python)
%o from sympy import divisor_sigma
%o def A184387(n): return (m:=divisor_sigma(n))*(m+1)>>1 # _Chai Wah Wu_, Jul 05 2023
%Y Cf. A000203, A000217, A002117, A184388, A184389, A184390, A184391, A130674.
%K nonn
%O 1,2
%A _Jaroslav Krizek_, Jan 12 2011
|
