Friday, 15 March 2019

combinatorics - Can this binomial polynomial sum be simplified?



nk=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)=nk=0(nk)kd.


using

kd=d!2πi|z|=ϵ1zd+1exp(kz)dz.



This yields for the sum
d!2πi|z|=ϵ1zd+1nk=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+1nq=0(nq)2nq(exp(z)1)qdz=nq=0(nq)2nq×q!×{dq}



for a final answer of

dq=0(nq)2nq×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 nq>d and we may replace
n by d. On the other hand if n<d the binomial coefficient is
zero for n<qd and we may again replace n by d.




Here we have used the species equation for labeled set partitions
P(UP1(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

real analysis - How to find limhrightarrow0fracsin(ha)h

How to find limh0sin(ha)h without lhopital rule? I know when I use lhopital I easy get $$ \lim_{h\rightarrow 0}...