login
a(0)=1; a(n) = sigma_2(n) + sigma_3(n).
0

%I #9 Feb 05 2020 14:19:50

%S 1,2,14,38,94,152,302,394,670,848,1264,1454,2254,2368,3346,3788,5022,

%T 5204,7268,7222,9744,10132,12598,12698,17230,16402,20632,21260,26162,

%U 25232,33052,30754,38814,38516,45676,44644,57172,52024,63550,63244,75920,70604,89188

%N a(0)=1; a(n) = sigma_2(n) + sigma_3(n).

%F G.f.: 1 + Sum_{k>=1} k^2*(k + 1)*x^k/(1 - x^k). - _Ilya Gutkovskiy_, Mar 17 2017

%t Join[{1},Table[DivisorSigma[2,n]+DivisorSigma[3,n],{n,50}]] (* _Harvey P. Dale_, Feb 05 2020 *)

%Y Cf. A001157, A001158.

%K nonn

%O 0,2

%A _N. J. A. Sloane_, Mar 20 2004