|
|
A363661
|
|
a(n) = Sum_{d|n} (n/d)^d * binomial(d+n,n).
|
|
4
|
|
|
2, 12, 32, 150, 282, 1890, 3488, 21582, 54650, 282612, 705564, 4072224, 10400782, 55006530, 158987232, 790611350, 2333606526, 11573213196, 35345264180, 168673694070, 540848064614, 2500462200182, 8233430728152, 37445946291600, 126411051769652
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
a(n) = [x^n] Sum_{k>0} (1/(1 - k*x^k)^(n+1) - 1).
|
|
MAPLE
|
f:= proc(n) local d;
add((n/d)^d * binomial(n+d, n), d = numtheory:-divisors(n))
end proc:
|
|
MATHEMATICA
|
a[n_] := DivisorSum[n, (n/#)^# * Binomial[# + n, n] &]; Array[a, 30] (* Amiram Eldar, Jul 05 2023 *)
|
|
PROG
|
(PARI) a(n) = sumdiv(n, d, (n/d)^d*binomial(d+n, n));
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|