%I #38 Sep 08 2022 08:44:51
%S 4,18,130,1134,10642,103158,1015690,10078254,100390882,1001953638,
%T 10009766650,100048830174,1000244144722,10001220711318,
%U 100006103532010,1000030517610894,10000152587956162,100000762939584198,1000003814697527770,10000019073486852414
%N Sum of n-th powers of divisors of 10.
%H T. D. Noe, <a href="/A034517/b034517.txt">Table of n, a(n) for n=0..200</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (18,-97,180,-100).
%F G.f.: -2*(90*x^3-97*x^2+27*x-2) / ((x-1)*(2*x-1)*(5*x-1)*(10*x-1)). - _Colin Barker_, Apr 20 2014
%t With[{dt = Divisors[10]}, Table[Total[dt^n], {n, 0, 20}]] (* _Harvey P. Dale_, Feb 02 2011 *)
%t Total[#^Range[0, 20]&/@Divisors[10]] (* _Vincenzo Librandi_, Apr 17 2014 *)
%t DivisorSigma[Range[0,20],10] (* _Harvey P. Dale_, May 26 2019 *)
%o (Sage) [sigma(10,n)for n in range(0,18)] # _Zerinvary Lajos_, Jun 04 2009
%o (PARI) a(n)=sigma(10,n) \\ _Charles R Greathouse IV_, May 16 2011
%o (Magma) [DivisorSigma(n, 10): n in [0..20]]; // _Vincenzo Librandi_, Apr 17 2014
%K nonn,easy
%O 0,1
%A _N. J. A. Sloane_
%E More terms from _Colin Barker_, Apr 20 2014
|