A326121 Expansion of Sum_{k>=1} k^2 * x^(2*k) / (1 - k * x^k).

0, 1, 1, 5, 1, 18, 1, 33, 28, 58, 1, 246, 1, 178, 369, 577, 1, 1539, 1, 2774, 2531, 2170, 1, 16706, 3126, 8362, 20413, 35366, 1, 116444, 1, 135425, 178479, 131362, 94933, 1110999, 1, 524650, 1596521, 2530946, 1, 7280892, 1, 8403734, 16364457, 8389138, 1, 78568322, 823544, 43420683

Offset

1,4

Links

Table of n, a(n) for n=1..50.

Formula

a(n) = Sum_{d|n, d>1} (n/d)^d = Sum_{d|n, d<n} d^(n/d).

a(p) = 1, where p is prime.

a(n) = A055225(n) - n.

Mathematica

nmax = 50; CoefficientList[Series[Sum[k^2 x^(2 k)/(1 - k x^k), {k, 1, nmax}], {x, 0, nmax}], x] // Rest
Table[DivisorSum[n, (n/#)^# &, # > 1 &], {n, 1, 50}]

Prog

(PARI) a(n)={sumdiv(n, d, if(d > 1, (n/d)^d))} \\ Andrew Howroyd, Sep 10 2019

Crossrefs

Cf. A055225, A078308, A087909.

Sequence in context: A121755 A104174 A050400 * A008971 A151335 A297174

Adjacent sequences: A326118 A326119 A326120 * A326122 A326123 A326124

Keyword

nonn

Author

Ilya Gutkovskiy, Sep 10 2019

Status

approved