login
A327238
Expansion of Sum_{k>=1} ((1 + k * x^k)^k - 1).
7
1, 4, 9, 20, 25, 63, 49, 160, 108, 350, 121, 940, 169, 1225, 1475, 2304, 289, 7560, 361, 8025, 12446, 7139, 529, 58192, 3750, 13858, 61965, 102655, 841, 191181, 961, 318464, 220704, 40460, 354172, 1304370, 1369, 63175, 629863, 4012608, 1681, 1916733, 1849
OFFSET
1,2
LINKS
FORMULA
a(n) = Sum_{d|n} (n/d)^d * binomial(n/d,d).
a(p) = p^2, where p is prime.
MATHEMATICA
nmax = 43; CoefficientList[Series[Sum[((1 + k x^k)^k - 1), {k, 1, nmax}], {x, 0, nmax}], x] // Rest
Table[DivisorSum[n, (n/#)^# Binomial[n/#, #] &], {n, 1, 43}]
PROG
(PARI) a(n)={sumdiv(n, d, (n/d)^d * binomial(n/d, d))} \\ Andrew Howroyd, Sep 14 2019
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Ilya Gutkovskiy, Sep 14 2019
STATUS
approved