|
|
A320568
|
|
a(n) = Sum_{k=1..n} (-1)^(n-k)*binomial(n,k)*sigma(k).
|
|
3
|
|
|
1, 1, -2, 5, -14, 40, -111, 293, -731, 1726, -3882, 8408, -17796, 37423, -79337, 170917, -373812, 823585, -1809844, 3934750, -8424747, 17749392, -36874749, 75862016, -155359339, 318331170, -655146929, 1356952623, -2828151136, 5920984735, -12420140296, 26036563525
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
COMMENTS
|
Inverse binomial transform of A000203.
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1/(1 + x))*Sum_{k>=1} k*x^k/((1 + x)^k - x^k).
L.g.f.: Sum_{k>=1} sigma(k)*x^k/(k*(1 + x)^k) = Sum_{n>=1} a(n)*x^n/n.
|
|
MAPLE
|
seq(add((-1)^(n-k)*binomial(n, k)*sigma(k), k=1..n), n=1..32); # Paolo P. Lava, Oct 25 2018
|
|
MATHEMATICA
|
Table[Sum[(-1)^(n - k) Binomial[n, k] DivisorSigma[1, k], {k, n}], {n, 32}]
nmax = 32; Rest[CoefficientList[Series[1/(1 + x) Sum[k x^k/((1 + x)^k - x^k), {k, 1, nmax}], {x, 0, nmax}], x]]
|
|
PROG
|
(GAP) Flat(List([1..10], n->Sum([1..n], k->(-1)^(n-k)*Binomial(n, k)*Sigma(k)))); # Muniru A Asiru, Oct 15 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|