login
A339353
G.f.: Sum_{k>=1} k^2 * x^(k*(k + 1)) / (1 - x^k).
5
0, 1, 1, 1, 1, 5, 1, 5, 1, 5, 1, 14, 1, 5, 10, 5, 1, 14, 1, 21, 10, 5, 1, 30, 1, 5, 10, 21, 1, 39, 1, 21, 10, 5, 26, 30, 1, 5, 10, 46, 1, 50, 1, 21, 35, 5, 1, 66, 1, 30, 10, 21, 1, 50, 26, 70, 10, 5, 1, 91, 1, 5, 59, 21, 26, 50, 1, 21, 10, 79, 1, 130, 1, 5, 35, 21, 50, 50, 1, 110
OFFSET
1,6
COMMENTS
Sum of squares of divisors of n that are smaller than sqrt(n).
MATHEMATICA
nmax = 80; CoefficientList[Series[Sum[k^2 x^(k (k + 1))/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] // Rest
Table[DivisorSum[n, #^2 &, # < Sqrt[n] &], {n, 80}]
PROG
(PARI) a(n) = sumdiv(n, d, if (d^2 < n, d^2)); \\ Michel Marcus, Dec 02 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 01 2020
STATUS
approved