n∑k=0(nk)kd
where d is some fixed positive integer. Is this a well known sum that has a faster-than-O(n) evaluation? It looks similar to Faulhaber's formula, except with binomial coefficients.
Answer
Suppose we seek to evaluate
Sd(n)=n∑k=0(nk)kd.
using
kd=d!2πi∫|z|=ϵ1zd+1exp(kz)dz.
This yields for the sum
d!2πi∫|z|=ϵ1zd+1n∑k=0(nk)exp(kz)dz=d!2πi∫|z|=ϵ1zd+1(1+exp(z))ndz=d!2πi∫|z|=ϵ1zd+1(2+exp(z)−1)ndz=d!2πi∫|z|=ϵ1zd+1n∑q=0(nq)2n−q(exp(z)−1)qdz=n∑q=0(nq)2n−q×q!×{dq}
for a final answer of
d∑q=0(nq)2n−q×q!×{dq}.
This is a polynomial of degree d in n.
The change of the upper limit is justified as follows: if n>d
then the Stirling number is zero for n≥q>d and we may replace
n by d. On the other hand if n<d the binomial coefficient is
zero for n<q≤d and we may again replace n by d.
Here we have used the species equation for labeled set partitions
P(UP≥1(Z))
which gives the bivariate generating function of the Stirling numbers
of the second kind exp(u(exp(z)−1)).
No comments:
Post a Comment