A386777 a(n) = n^2*sigma_6(n).
0, 1, 260, 6570, 66576, 390650, 1708200, 5764850, 17043520, 43105851, 101569000, 214359002, 437404320, 815730890, 1498861000, 2566570500, 4363141376, 6975757730, 11207521260, 16983563402, 26007914400, 37875064500, 55733340520, 78310985810, 111975926400, 152597656875
Offset: 0
Links
- Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Programs
-
Magma
[0] cat [n^2*DivisorSigma(6, n): n in [1..35]]; // Vincenzo Librandi, Aug 04 2025
-
Mathematica
Table[n^2*DivisorSigma[6, n], {n, 0, 30}] nmax = 30; CoefficientList[Series[Sum[k^8*x^k*(1 + x^k)/(1 - x^k)^3, {k, 1, nmax}], {x, 0, nmax}], x]
Formula
G.f.: Sum_{k>=1} k^8*x^k*(1 + x^k)/(1 - x^k)^3.
a(n) = n^2*A013954(n).
Dirichlet g.f.: zeta(s-2)*zeta(s-8). - R. J. Mathar, Aug 03 2025