%I M4973 N2135 #51 May 09 2022 03:51:40
%S 0,1,15,76,275,720,1666,3440,6129,11250,17545,28896,41405,65072,85950,
%T 128960,162996,238545,286995,404600,482160,662112,756470,1042560,
%U 1150625,1549730,1732590,2257920,2443105,3250800,3421160,4452096,4791600,6039522,6296500
%N a(n) = Sum_{k=1..n-1} k^2*sigma(k)*sigma(n-k).
%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%D J. Touchard, On prime numbers and perfect numbers, Scripta Math., 129 (1953), 35-39.
%H John Cerkan, <a href="/A000477/b000477.txt">Table of n, a(n) for n = 1..10000</a>
%H J. Touchard, <a href="/A000385/a000385.pdf">On prime numbers and perfect numbers</a>, Scripta Math., 129 (1953), 35-39. [Annotated scanned copy]
%F a(n) = Sum_{k=1..n-1} k^2*sigma(k)*sigma(n-k). - _Sean A. Irvine_, Nov 14 2010
%F G.f.: x*f(x)*g'(x), where f(x) = Sum_{k>=1} k*x^k/(1 - x^k) and g(x) = Sum_{k>=1} k^2*x^k/(1 - x^k)^2. - _Ilya Gutkovskiy_, May 02 2018
%F a(n) = (n^2/24 - n^3/6)*sigma_1(n) + (n^2/8)*sigma_3(n). - _Ridouane Oudra_, Sep 15 2020
%F Sum_{k=1..n} a(k) ~ Pi^4 * n^6 / 4320. - _Vaclav Kotesovec_, May 09 2022
%e G.f. = x^2 + 15*x^3 + 76*x^4 + 275*x^5 + 720*x^6 + 1666*x^7 + 3440*x^8 + ...
%p with(numtheory): S:=(n,e)->add(k^e*sigma(k)*sigma(n-k),k=1..n-1); f:=e->[seq(S(n,e),n=1..30)]; f(2); # _N. J. A. Sloane_, Jul 03 2015
%t a[n_] := Sum[k^2 DivisorSigma[1, k] DivisorSigma[1, n-k], {k, 1, n-1}]; Array[a, 35] (* _Jean-François Alcover_, Feb 08 2016 *)
%o (PARI) a(n) = sum(k=1, n-1, k^2*sigma(k)*sigma(n-k)); \\ _Michel Marcus_, Feb 02 2014
%Y Cf. A000385, A000441, A000499, A259692, A259693, A259694, A259695, A259696.
%Y Cf. A000203 (sigma_1), A001158 (sigma_3).
%K nonn
%O 1,3
%A _N. J. A. Sloane_
%E More terms from _Sean A. Irvine_, Nov 14 2010
%E a(1)=0 prepended by _Michel Marcus_, Feb 02 2014