OFFSET
1,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
Vaclav Kotesovec, Plot of Sum_{k=1..n} a(k) / n^(3/2) for n = 1..10^7
FORMULA
G.f.: Sum_{k>=2} (k^2-1) * x^(k^2-1)/(1 - x^(k^2-1)).
Sum_{k=1..n} a(k) ~ zeta(3/2)*n^(3/2)/3. - Vaclav Kotesovec, Jan 21 2023
MATHEMATICA
Table[Sum[If[IntegerQ[Sqrt[d+1]], d, 0], {d, Divisors[n]}], {n, 1, 100}] (* Vaclav Kotesovec, Jan 21 2023 *)
PROG
(PARI) a(n) = sumdiv(n, d, issquare(d+1)*d);
(PARI) my(N=100, x='x+O('x^N)); concat([0, 0], Vec(sum(k=2, sqrtint(N+1), (k^2-1)*x^(k^2-1)/(1-x^(k^2-1)))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 20 2023
STATUS
approved