一类恒等式的应用(范德蒙德卷积与超几何函数)

你可以在这里找到一个PDF版本


翻到了2年前的一篇日报,一类恒等式的应用 -- foreverlastnig 的博客 ,里面指出了一个恒等式的应用:

k=0n(rk)(snk)=(r+sn)

这个恒等式就是大名鼎鼎的范德蒙德卷积,它最早是由中国人朱世杰于1303年发现的,法国人范德蒙德在18世纪重新发现了它。本文尝试从超几何函数的角度更进一步探究其在组合数的恒等变形中的重要意义。本文内容大部分来源于对《具体数学》第五章其中有关超几何函数的整理。

本文均只在实数范围内讨论,不考虑复数。

约定与前置知识

连续阶乘(广义阶乘函数):

最开始,大家熟知的阶乘定义是:

n!=k=1nk

欧拉将其拓展到了全体实数上(也适用于部分复数,超出了本文的范围):

x!=0txetdt

同样满足阶乘的基本性质,即:

x!=x(x1)!

欧拉还定义了一个相仿的函数,称为 Γ 函数,满足:

Γ(x+1)=x!Γ(x+1)=xΓ(x)(x)!Γ(x)=πsin(πx)

通过广义阶乘,我们可以定义广义上的二项式系数(取恰当极限):

(ab)=limralimsbr!s!(rs)!

通过广义阶乘,定义下降幂和上升幂:

xn=Γ(x+n)/Γ(x)xn=x!/(xn)!

对于广义阶乘的乘法逆元,欧拉还证明了,

1x!=limn(n+xn)nx

由此可知,对于任何负整数 x,都有

1x!=0

二项式系数:

基本的二项式恒等式应当被熟练使用,并了解多项式推理法,可以参考二项式系数 -- zhiyangfan 的小窝。下面列出一些常用的供读者快速索引,

(nk)=(nnk),nN,kZ(a+b)r=k(rk)akbrk(rk)=rk(r1k1),kZ,k0(rs)=(r1k)+(r1k1),kZ(rs)=(1)s(s1rs),sZ

生成函数:

基本的生成函数知识会被使用,可以参考铃悬的数学小讲堂——生成函数初步 -- 铃悬 的博客

超几何函数

化简求和式子的时候,我们通常可以将原式转化为下面形式

k0tk

考虑这个求和相邻两项的比值 tk+1/tk,若为一个常数 c,则有

tk=t0ck

运用等比数列的求和公式,

k0tk=k0t0ck=t0k0ck=t01c

上面运用了几何级数的求和公式。然而大部分情况相邻两项的比值并非一个常数,举个例子,求下面式子的封闭形式,其中 n0

A=k=0n(nk)2

翻转求和顺序,并重写求和范围,

A=k0(nnk)2

tk=(nnk)2,其相邻两项的比值

tk+1tk=(nnk1)2(nnk)2=(n!(nk1)!(k+1)!/n!(nk)!k!)2=(kn)2(k+1)2

是关于 k 的一个有理函数(两个关于 k 的多项式的商),由此考虑构造级数

S=k0((n)kk!)2

其相邻两项的比值也为 (kn)2(k+1)2。这个级数 S 的常数项的值为 1,巧的是 A 的常数项也为 1,这使得我们断定这两个级数相等,即

A=t0S=S

通过后面要讲述的方法,我们将可以轻松的求出 S=(2nn),从而求出了 A 的封闭形式。

为了更加通用和方便地通过 tk+1/tk 构造出形如 S 的级数,我们使用一般的超几何函数,它是关于 x 的带有 n+m 个参数的幂级数,定义为

F(a1,,amb1,,bn|x)=k0a1kamkxkb1kbnkk!

(在行内,我们则使用记号 F(a1,,am;b1,,bn;x) 表示)不难发现,超几何级数的常数项总为 1,而其相邻两项的比值为

tk+1tk=a1k+1amk+1b1k+1bnk+1b1kbnka1kamkk!xk+1(k+1)!xk=(k+a1)(k+am)x(k+b1)(k+bn)(k+1)

这是关于 k 的一个有理函数。因此,相邻两项的比值是有理函数(且可以因式分解成上面的形式)的级数总可以写成超几何函数的形式。使得许多组合恒等式成为超几何函数的特例而已。在上面的例子里,S 其实是超几何函数的一个特殊例子,即

S=F(n,n1|1)

高斯超几何函数与范德蒙德卷积

一个拥有两个上参数和一个下参数的超几何函数被称为高斯超几何函数(又称普通超几何函数),它的形式是

F(a,bc|x)=k0akbkzkckk!

高斯超几何函数取 x=1 时与范德蒙德卷积有着密切的联系。考虑范德蒙德卷积,我们知道

k=0n(rk)(snk)=(r+sn)

另一方面,

L.H.S.=k0(rk)(snk)

其第 k 项为 tk=(rk)(snk),相邻两项的比值

tk+1tk=(rk+1)(snk1)(rk)(snk)=r!s!k!(rk)!(nk)!(sn+k)!r!s!(k+1)!(rk1)!(nk1)!(sn+k+1)!=(rk)(nk)(k+1)(sn+k+1)=(kr)(kn)(k+sn+1)(k+1)

由此可知,

L.H.S.=t0F(r,nsn+1|1)=(sn)F(r,nsn+1|1)=R.H.S

代入范德蒙德卷积右边,并将 t0 移项,得到

F(r,nsn+1|1)=(r+sn)(sn)F(r,nsn+1|1)=(r+s)!(sn)!(r+sn)!s!

a=r,b=n,c=sn+1,整理得到

F(a,bc|1)=Γ(cab)Γ(c)Γ(ca)Γ(cb)

高斯证明了,这个等式成立需满足 bZ,b0 或者 c>a+b 时成立,否则等式左边的级数不收敛。这个等式从某种意义上说是范德蒙德卷积的另一种表现形式,这种形式能够更好的展现各种组合和式与范德蒙德卷积之间的关系。下面举几个应用,

例题1 (平行求和法)求下面式子的封闭形式,其中整数 r,n0

k=0n(r+kk)

令原式为 S,翻转求和顺序

S=k=0n(r+nknk)=k0(r+nknk)

级数的第 ktk=(r+nknk),其相邻两项的比值

tk+1tk=(r+nk1nk1)/(r+nknk)=knknr=(k+1)(kn)(knr)(k+1)

由此可知,

S=t0F(1,nnr|1)=(r+nn)F(1,nnr|1)

其中 n0,应用范德蒙德卷积,

S=(r+nn)Γ(nr(n)1)Γ(nr)Γ(nr(n))Γ(nr1)=(r+nn)Γ(r1)Γ(rn)Γ(r)Γ(rn1)=(r+nn)r+n+1r+1=(r+n+1n)

例题2 求下面式子的封闭形式

k0(n+k2k)(2kk)(1)kk+1

令原式为 S,级数的第 ktk=(n+k2k)(2kk)(1)kk+1,其相邻两项的比值

tk+1tk=(n+k+12k+2)(2k+2k+1)(1)(k+1)k+2(n+k2k)(2kk)(1)kk+1=(k+n+1)(kn)(k+2)(k+1)

由此可知,

S=t0F(n+1,n2|1)=F(n+1,n2|1)

其中 2>n+1+(n),应用范德蒙德卷积,

S=Γ(2n1+n)Γ(2)Γ(2n1)Γ(2+n)=1Γ(1n)Γ(2+n)=sin(πn)n(n+1)π

特别的,当 n 为自然数时,S=[n=0]

例题3 求下面式子的封闭形式,其中整数 nm0

k=0m(mk)/(nk)

令原式为 S,取消求和上界以把 S 写成级数的形式,

S=k0(mk)/(nk)

级数的第 ktk=(mk)/(nk),其相邻两项的比值

tk+1tk=(mk+1)(nk)(mk)(nk+1)=kmkn=(k+1)(km)(kn)(k+1)

由此可知,

S=t0F(1,mn|1)=F(1,mn|1)

其中整数 m0,应用范德蒙德卷积,

S=Γ(n1+m)Γ(n)Γ(n1)Γ(n+m)=n+1n+1m

例题4 求下面式子的封闭形式,其中整数 nm0

k=0n1(km)1nk

令原式为 S,翻转求和顺序,

S=k=0n1(n1km)1k+1=k0(n1km)k+1

到目前为止一切正常,级数的第 ktk=(n1km)k+1,我们似乎只用算出它的相邻两项的比值即可

tk+1tk=(k+1)(k+n2m)(k+2)(k+n1m)=(k+1)(k+mn+1)(k+2)(kn+1)

这显然不是一个高斯超几何函数相邻两项的比值,它对应的是,

S=t0F(mn+1,1,1n+1,2|1)

我们尚且还不知道这样的超几何函数如何求解。退一步,我们注意到似乎只要在求和中移动 k,就可以使比值的 (k+1)(k+2) 变为 k(k+1)。具体而言,考虑

S=k0(n1km)k+1=k1(nkm)k=limϵ0(nm)ϵ+k0(nkm)k+ϵ

后半部分级数的第 ktk=(nkm)k+ϵ,其相邻两项的比值

tk+1tk=limϵ0(nk1m)k+ϵ+1(nkm)k+ϵ=limϵ0(nk1m)k+1(nkm)k+ϵ=limϵ0(k+ϵ)(nk1m)(k+1)(nkm)=limϵ0(k+ϵ)(k+mn)(kn)(k+1)

由此而言,

S=limϵ0(nm)ϵ+t0F(ϵ,mnn|1)=limϵ0(nm)ϵ(F(ϵ,mnn|1)1)

其中整数 mn0,应用范德蒙德卷积,

S=limϵ0(nm)ϵ(Γ(mϵ)Γ(n)Γ(nϵ)Γ(m)1)=(nm)limϵ01ϵ(m!(n+ϵ)!sin(mπ)sin(nπ+ϵπ)n!(m+ϵ)!sin(mπ)sin(mπ+ϵπ)1)=(nm)limϵ01ϵ(m!(n+ϵ)!sin(mπ)(sin(nπ)cos(ϵπ)+sin(ϵπ)cos(nπ))n!(m+ϵ)!sin(mπ)(sin(mπ)cos(ϵπ)+sin(ϵπ)cos(mπ))1)=(nm)limϵ01ϵ(m!(n+ϵ)!sin(mπ)sin(nπ)n!(m+ϵ)!sin(mπ)sin(mπ)1)=(nm)limϵ01ϵ(m!(n+ϵ)!n!(m+ϵ)!1)

对于后半部分的极限,可以直接运用洛必达法则求值,得到

limϵ01ϵ(m!(n+ϵ)!n!(m+ϵ)!1)=HnHm

综上所述,

S=(nm)(HnHm)

库默尔公式与迪克逊公式

有几个公式与范德蒙德卷积十分相似的,属于同一类型,其中一个是

k=02n(1)k(rk)(r2nk)=(1)n(rn)

一方面,考虑生成函数 (1x)r(1+x)r=(1x2)r 两侧 xn 的系数相等,原等式显然成立。
另一方面,从超几何函数的角度考虑,

L.H.S.=k0(1)k(rk)(r2nk)=(r2n)F(2n,r2n+r+1|1)=R.H.S.

代入原等式的右侧,并移项,

F(2n,r2n+r+1|1)=(1)n(rn)(r2n)F(2n,r2n+r+1|1)=(1)n(2n)!(2nr)!n!(rn)!

a=2n,b=r,整理得到

F(a,b1+ba|1)=(b/2)!(ba)b/2b!

这就是库默尔公式,由库默尔证明于1836年。

例题5 求下面式子的封闭形式,其中整数 a,b0

k(1)k(a+ba+k)(b+ab+k)

令原式为 S,确定 S 的求和范围,写出它的超几何函数形式,并对于这个式子使用库默尔公式,

S=k=aa(1)k(a+ba+k)(b+ab+k)=k0(1)ka(a+bk)(a+bba+k)=(1)a(a+bba)F(ab,2a1a+b|1)=π(4)a(a+bba)Γ(a+b+1)b!Γ(12a)=(a+b)!a!b!

考虑例题5中的恒等式的三项版本,我们可以用其导出一个与库默尔公式相似的超几何函数恒等式

kZ(1)k(a+ba+k)(b+cb+k)(c+ac+k)=(a+b+c)!a!b!c!

这个恒等式对于整数 a,b,c0 成立,证明参见 A Proof of Dixon's Identity -- Journal of Integer Sequences, Vol. 19 (2016)。我们从超几何函数的角度讨论,

确定求和范围,不妨设 0abc,则

L.H.S.=k=aa(1)k(a+ba+k)(b+cb+k)(c+ac+k)=k0(1)ka(a+bk)(b+cb+ka)(a+cc+ka)=(1)a(b+cba)(a+cca)F(2a,ab,aca+b+1,a+c+1|1)=R.H.S.

代入原等式的右侧,并移项,

F(2a,ab,aca+b+1,a+c+1|1)=(1)a(a+b+c)!a!b!c!(b+cba)(a+cca)

a=2a,b=ab,c=ac,整理得到

F(a,b,c1+ca,1+cb|1)=(c/2)!(ca)c/2(cb)c/2c!(cab)c/2

这就是迪克逊公式,迪克逊最早证明了上面等式在 a+b<1+c/2 时成立。

拓展范德蒙德卷积

关于三个二项式系数的卷积十分罕见而重要,考虑整数 m,n,有

k(mr+sk)(n+rsnk)(r+km+n)=(rm)(sn)

对此证明的关键点是注意到 (r+km+n)=j(rm+nj)(kj),然后使用范德蒙德卷积,即可堆出右边的表达式。

另一方面,

L.H.S.=k0(mr+sk)(n+rsnk)(r+km+n)=(n+rsn)(rm+n)F(1+r,rms,nrs+1,rmn|1)=R.H.S.

代入原等式的右侧,并移项,

F(1+r,rms,nrs+1,rmn|1)=(rm)(sn)(n+rsn)(rm+n)

a=1+r,b=rms,c=rs+1,整理得到

F(a,b,nc,a+bcn+1|1)=(ca)n(cb)ncn(cab)n

其中等式成立的条件是整数 n0。它通常被称为Saalschütz恒等式,最早由普法夫发表于1797年。从某种意义上说,这个恒等式是拓展范德蒙德卷积,因为当 n 时的极限就是范德蒙德卷积。

例题6 求下面式子的封闭形式,其中整数 m,n0

k0(n+k2k)(2kk)(1)kk+m+1

求出原级数的超几何表达,并应用拓展范德蒙德卷积,发现这是 c=1 时的特殊情况。

超几何变换

我们已经拥有了若干超几何函数与组合数之间的恒等式,但这个话题还有更多讨论内容,因为超几何函数也有它们自己的恒等式,称为超几何变换。这样一来,即使一个超几何函数的值尚不知道,我们也可以尝试先利用超几何变换将其变形为一个已知的形式。

我们在研究例题4时注意到对超几何级数移项有时是一种有效的处理方法,我们现在考虑一般情况下,

F(a1,,amb1,,bn|x)=k0a1kamkxkb1kbnkk!=1+k1a1kamkxkb1kbnkk!=1+a1amxb1bnk0(a1+1)k(am+1)kxk(b1+1)k(bn+1)k(k+1)!=1+a1amxb1bnF(a1+1,,am+1,1b1+1,,bn+1,2|x)

由此我们得到右移项变换

F(a1,,amb1,,bn|x)=1+a1amxb1bmF(a1+1,,am+1,1b1+1,,bn+1,2|x)

这是最基本的一个超几何变换。如果我们在一点知道它,那么就不用在解决例题4时绕弯路了。

1797年,德国数学家J.F.普法夫(Johann Friedrich Pfaff)发表了一个惊人的反射定律

1(1x)aF(a,bc|x1x)=F(a,cbc|x)

例题7 求下面式子的封闭形式,其中整数 n0

k(1)k2k(nk)(n+kk)/(n+b+kk)

连接反射定律和库默尔公式,我们可以得到

F(a,1a1+ba|12)=2aF(a,b1+ba|1)=2a(b/2)!(ba)b/2b!

将这个等式转换为关于二项式系数的恒等式,有

k(1)k2k(nk)(n+kk)/(n+b+kk)=2n(b/2)!(b+n)!b!(b/2+n)!

例题8 求下面式子的封闭形式,

k=0m2m+11+k(2)k(mk)

令原式为 S,翻转求和顺序,

S=k02m+1m+1k(2)mk(mmk)=2m+1m+1(2)mF(m1,11|12)

应用反射定律,

F(m1,11|12)=12m+1F(m1,01|1)

注意到对于 k0,有

0k=Γ(k)Γ(0)=[k=0]

因此,

F(m1,01|1)=1

代回原式,得到

S=(1)m(2m+1)2(m+1)

例题9 求下面式子的封闭形式,其中整数 m0

k=0m2m+12m+1k(2)k(mk)

令原式为 S,翻转求和顺序,依次应用反射定律和库默尔公式

S=k02m+1m+1+k(2)mk(mmk)=2m+1m+1(2)mF(m,m+1m+2|12)=2m+1m+1(1)mF(m,1m+2|1)=(1)m2m+1m+1(1/2)!(m+1)1/2=(1/2m)1

这真是个不错的结果!不过我们不应该止步于此,是时候该找出更多超几何变换了。通过二项式定理,我们很容易注意到,对于整数 m,有

k=0m(m+rk)xkymk=k=0m(rk)(x)k(x+y)mk

为了证明这个恒等式,我们首先假定 mr0,随后使用多项式推理法。显然有

(x+y)m+ryr=(x+y)m+r(x+(x+y))r

运用二项式定理,

L.H.S.=yrk=0m+r(m+rk)xkym+rk=k=0m+r(m+rk)xkymk

其中 m+rm,而当 k>m+r 时始终有 (m+rk)=0,故可以重写求和上界为 m

L.H.S.=k=0m(m+rk)xkymk

同理,

R.H.S.=(x+y)m+rk=0r(rk)(x)k(x+y)rk=k=0r(rk)(x)k(x+y)mk=k=0m(rk)(x)k(x+y)mk

代入 L.H.S=R.H.S,对于 mr0m+1r 的取值,有

k=0m(m+rk)xkymk=k=0m(rk)(x)k(x+y)mk

注意到等式两边均为关于 rm 次或更低次多项式,在上述 m+1 个不同取值处相等说明等式在一般情况下成立。

另一方面,我们尝试将其转化为超几何函数的一个恒等式。具体而言,考虑将等式两边关于 y 偏微分 n 次,同时除以 n!,并用 mnk 代替 k

dndynn!L.H.S.=dndynn!k=0m(m+rk)xkymk=1n!k=0mdndyn(m+rk)xkymk=1n!k=0m(m+rk)xk(mk)nymnk=k0(m+rmkn)(n+k)nn!xmnkyk=k0(m+rmnk)(n+kn)xmnkyk

dndynn!R.H.S.=dndynn!k=0m(rk)(x)k(x+y)mk=1n!k=0mdndyn(rk)(x)k(x+y)mk=1n!k=0m(rk)(x)k(mk)n(x+y)mnk=k0(rmnk)(n+k)nn!(x)mnk(x+y)k=k0(rmnk)(n+kn)(x)mnk(x+y)k

代入 dndynn!L.H.S.=dndynn!R.H.S.

k0(m+rmnk)(n+kn)xmnkyk=k0(rmnk)(n+kn)(x)mnk(x+y)k

x=1,将等式两边写成超几何函数的形式并整理,最终得到一条超几何变换式,对于整数 n0

F(a,nc|x)=(ac)n(c)nF(a,n1+a+nc|1x)

这就是范德蒙德变换,当 x=1 时它退化为范德蒙德卷积。

微分方程与高斯恒等式

为了简便起见,我们简记关于 x 的求导算子为 D,即

D=ddx

上面推理范德蒙德变换所使用的微分方法看起来十分有用,事实上,这个例子是具有指导意义的。我们来看看一个一般的超几何函数关于 x 求导时会发生什么:

DF(a1,,amb1,,bn|x)=k1a1kamkxk1b1kbnk(k1)!=k+11a1k+1amk+1xkb1k+1bnk+1k!=k0a1(a1+1)kam(am+1)kxkb1(b1+1)kbn(bn+1)kk!=a1amb1bnF(a1+1,,am+1b1+1,,bn+1|x)

这样参数被移出来了一项并且原参数都增加了 1。当然,我们尝试是否可以用微分法调整超几何函数中的一个参数而保持其他参数不变。为此,我们考虑算子

ϑ=xD

ϑ 算子先对函数求关于 x 的导数,然后再把结果乘 x。这个算子给出

ϑF(a1,,amb1,,bn|x)=k0ka1kamkzkb1kbnkk!

这个式子乍一看不太有用,但是我们尝试用 ϑ 和某个上参数(如 a1)的和乘以 F(a1,am;b1,bn;x),得到

(ϑ+a1)F(a1,,amb1,,bn|x)=k0(a1+k)a1kamkzkb1kbnkk!=k0a1(a1+1)ka2kamkzkb1kbnkk!=a1F(a1+1,a2,amb1,,bn|x)

这样,只有选定的一个上参数改变了。同时,一个类似的技巧也可以运用于下参数,不过这将把下参数减少而不是增加:

(ϑ+b11)F(a1,,amb1,,bn|x)=k0(k+b11)a1kamkzkb1kbnkk!=k0(b11)(a1+1)kamkzk(b11)kbnkk!=(b11)F(a1,,amb11,b2,bn|x)

现在我们尝试把这两个式子结合起来看(简记 F(a1,am;b1,bn;x)=F),我们有

(ϑ+a1)(ϑ+am)F=a1amF(a1+1,,am+1b1,,bn|x)(ϑ+b11)(ϑ+bn1)F=(b11)(bn1)F(a1,,amb11,,bn1|x)

我们容易发现,上面一行右边的超几何函数是下面一行右边的超几何函数的导数,由此我们得到一个一般的超几何函数都满足的微分方程

D(ϑ+b11)(ϑ+bn1)F=(ϑ+a1)(ϑ+am)F

这个似乎让人眼花缭乱式子的式子就是超几何微分方程。我们不妨从一个基本的例子出发,只考虑普通超几何函数 F(x)=F(a,b;c;x)。代入上面的微分方程,我们有

D(ϑ+c1)F=(ϑ+a)(ϑ+b)F

代入算子 Dϑ 的含义,并整理成通常的微分方程的记号,

x(1x)F(x)=(c(a+b+1)x)F(x)abF(x)=0

尝试求解这个微分方程会发现,F(a,b;c;x) 是唯一满足这个微分方程且常数项为 1 的形式幂级数。事实上,这不是一个巧合。我们尝试从任意一个超几何微分方程回到幂级数:假设 F(x)=k0tkxk 是某个满足超几何微分方程的幂级数。从超几何微分方程的形式出发,必定有

tk+1tk=(k+a1)(k+am)(k+b1)(k+bn)(k+1)

从而解得,F(x)=cF(a1,,am;b1,,bn;x)。因此,我们证明了超几何级数 F(a1,,am;b1,,bn;x) 是唯一满足超几何微分方程且常数项为 1 的形式幂级数。

我们还可以用 x 乘超几何微分方程的等式两边,这将给出一个相当有启发性的只包含算子 ϑ 式子:

ϑ(ϑ+c1)F=x(ϑ+a)(ϑ+b)F

这个式子的启发性意义在于其算子与超几何级数的相邻两项的比值一一对应:左边第一个因子 ϑ=(ϑ+11) 对应相邻两项的比值中的 k+1,与一般的超几何级数中的分母的 k! 对应;其他因子 (ϑ+bj1)bjk 相对应,而右边的 zzk 相对应, (ϑ+aj) 则与 ajk 相对应。

例题10 (二项式定理)求证:

(a+b)r=k=0r(rk)akbrk

这不是什么新鲜的式子,不过从超几何的角度来研究一下总是有益的。首先将等式右边重写为超几何函数的形式,

R.H.S.=brF(r|ab)

x=a/b,得到

R.H.S.=brF(r|x)

其中 F(r;;x) 满足超几何微分方程

ϑF=x(ϑr)F

另一方面 f(x)=(1x)r 也满足上面的超几何微分方程,由此可知

(1x)r=F(r|x)

代回原式,

L.H.S.=(bx+b)r=br(1x)r=brF(r|x)=R.H.S.

由此,原恒等式成立。

言归正传,这个微分理论的一个用途是证明新的超几何变换,因为每个确定的超几何微分方程只有一个(常数项为 1 的)解,从而如果两个形式上不同的超几何函数满足同一个超几何微分方程,那么它们一定总是相等。考虑下面两个超几何级数

F(2a,2ba+b+12|x),F(a,ba+b+12|4x(1x))

它们都满足微分方程

x(1x)F(x)+(a+b+12)(12x)F(x)4abF(x)=0

由此可知,

F(2a,2ba+b+12|x)=F(a,ba+b+12|4x(1x))

这就是高斯恒等式,由德国数学家高斯在1813年发表。特别地,当 x=1/2 时,可以连接高斯恒等式和范德蒙德卷积,

F(2a,2ba+b+12|12)=F(a,ba+b+12|1)=Γ(1/2)Γ(a+b+1/2)Γ(a+1/2)Γ(b+1/2)=πab+1/2(a1/2)!

运用这个公式时必须注意原超几何级数只在 a+b+1/2 是正整数时收敛。

例题11 求下面式子的封闭形式,其中整数 mn0

k=0m(1)k2k(mkn)(m+n+1k)

令原式为 S。注意到只有满足 0kmn 的项是非零的,重写求和上界,

S=k0(1)k2k(mkn)(m+n+1k)

计算其相邻两项的比值

tk+1tk=(k+m1n)(m+n+1k+1)2(mkn)(m+n+1k)=(kmn1)(km+n)2(km)(k+1)

由此可知,

S=(mn)F(mn1,nmm|12)=[m+n是偶数](mn)π(mn12)n+1m2(mn1)!=[m+n是偶数]2nm((m+n)/2n)

例题12 求下面式子的封闭形式,其中整数 n0n2(mod3)

k8k(nk)(13n16k)/(43n23k)

令原式为 S。将 S 写成超几何函数的形式,并应用高斯恒等式和库默尔公式,

S=F(1613n,n2343n|8)=F(1323n,2n2343n|1)=(2nn)/(43n23n)

这实在是一个令人感到不可思议的结果,这个恒等式即使是在小数值时也不容易手工计算。然而这个恒等式完全没有用(除了用作一道例题之外),它几乎肯定不会在实际问题中出现。

小结与闲话

那本文到此所有的知识点就都讲解完毕了~,写了超过2.6w字,总时长四天,还是很希望能让大家(当然还有我自己)学到新东西的。在《具体数学》上还有部分超几何求和Gosper算法)和机械求和法的讲解作为超几何方法的更进一步拓展。另一方面,根据Wolfram Mathematica的提示(DifferenceRootReduce[]),超几何函数与生成函数、整式递推和矩阵还有更深刻的联系,我以后也许会单独谈谈(其实现在我还没想清楚这个联系是什么,但是可以确定它是存在的,参考一个三重求和下的组合数问题 -- 方而静的笔记)。

写这篇文章时我每个式子都在纸上推导后又用数学软件验证,并与书上结果核对了,但难免仍有错漏,欢迎读者指出。祝您早日成为数数大师!

posted @   方而静  阅读(1169)  评论(1编辑  收藏  举报
相关博文:
阅读排行:
· DeepSeek 开源周回顾「GitHub 热点速览」
· 物流快递公司核心技术能力-地址解析分单基础技术分享
· .NET 10首个预览版发布:重大改进与新特性概览!
· AI与.NET技术实操系列(二):开始使用ML.NET
· .NET10 - 预览版1新功能体验(一)
点击右上角即可分享
微信分享提示