具体数学组合数习题选做(Genshining)

8.计算

k(nk)(1)k(1kn)nk(nk)(1)k(1kn)n

解:

考虑

nk=0(nk)(1)kf(k)=(1)nn![kn]f(k)nk=0(nk)(1)kf(k)=(1)nn![kn]f(k)

其中 f(k)f(k)nn 次多项式。

显然可以设 f(k)=(1kn)nf(k)=(1kn)n

容易发现,[xn]f(k)=(1)nnn[xn]f(k)=(1)nnn

故原式 =n!nn=n!nn

根据斯特林近似,n!nn2πnenn!nn2πnen

15.求

k(nk)3(1)kk(nk)3(1)k

可以证明,该式不可使用超几何项部分求和。

考虑机械求和法,考虑:

k(a+ba+k)(b+cb+k)(c+ac+k)(1)k=(a+b+ca,b,c),a,b,cZk(a+ba+k)(b+cb+k)(c+ac+k)(1)k=(a+b+ca,b,c),a,b,cZ

不妨设 2|n2|n。此时取 a=b=c=n2a=b=c=n2,得原式 ==

(3n/2)!(n/2)!3(3n/2)!(n/2)!3

(2n1/2n)(2n1/2n)(2n1/22n)(2n1/22n) 的简单关系。

考虑

(rk)(r1/2k)=(2r2k)(2kk)22k(rk)(r1/2k)=(2r2k)(2kk)22k

同时化开两者,得:

LHS=(4n2n)22n,RHS=(4n2n)22nLHS=(4n2n)22n,RHS=(4n2n)22n

显然。

证明:

1z!=limn(n+mn)nm,zN1z!=limn(n+mn)nm,zN

这是阶乘的典型定义。事实上可以证明 zZzZ 的情况。

RHS=nnm+1inz!RHS=nnm+1inz!

有限项相乘,所以答案是 1z!1z!

22.证明:

x!(x1/2)!=(2x)!(1/2)!22xx!(x1/2)!=(2x)!(1/2)!22x

化为 22x(2x)!=(1/2)!x!(x1/2)!22x(2x)!=(1/2)!x!(x1/2)!

考虑阶乘定义。

RHS=limn(n+xn)(n+x1/2n)n2x(n1/2n)=limn(2n+2x2n)n2xRHS=limn(n+xn)(n+x1/2n)n2x(n1/2n)=limn(2n+2x2n)n2x

LHS=limn(2n)2x22x=RHSLHS=limn(2n)2x22x=RHS

23.求 F(n,1;;1)F(n,1;;1)

这里我们认为错排数 DnDn 是封闭形式。考虑其式子:

Dn=n!nk=0(1)kk!Dn=n!nk=0(1)kk!

这个超几何函数是

k0(n)¯kk0(n)¯¯¯k

因此答案是 (1)nDn(1)nDn

24.求:

k(nm+k)(m+k2k)4kk(nm+k)(m+k2k)4k

机械求和,启动!

超几何形式是

(nm)F(mn,m;1/2;1)(nm)F(mn,m;1/2;1)

考虑高斯超几何函数的公式:

F(a,b;c;1)=Γ(cab)Γ(c)Γ(ca)Γ(cb)F(a,b;c;1)=Γ(cab)Γ(c)Γ(ca)Γ(cb)

答案等于

(2n2m)(2n2m)

25.证明

(a1b1)F(a1,a2,,am;b1+1,,bn;z)=a1F(a1+1,a2,,am;b1+1,,bn;z)b1F(a1,a2,,am;b1,,bn;z)(a1b1)F(a1,a2,,am;b1+1,,bn;z)=a1F(a1+1,a2,,am;b1+1,,bn;z)b1F(a1,a2,,am;b1,,bn;z)

F=F(a1,a2,,am;b1+1,,bn;z;)F=F(a1,a2,,am;b1+1,,bn;z;)

于是可以得出其算子等价形式:(a1b1)=(ϑ+a1)(ϑ+b1)(a1b1)=(ϑ+a1)(ϑ+b1),这是显然的。

类似的,有

(a1a2)F(a1,a2,,am;b1,,bn;z)=a1F(a1+1,a2,,am;b1,,bn;z)a2F(a1,a2+1,,am;b1,,bn;z)(a1a2)F(a1,a2,,am;b1,,bn;z)=a1F(a1+1,a2,,am;b1,,bn;z)a2F(a1,a2+1,,am;b1,,bn;z)

以及其他线性组合。

28.证明欧拉恒等式:

F(a,b;c;z)=(1z)cabF(ca,cb;c;z)F(a,b;c;z)=(1z)cabF(ca,cb;c;z)

考虑反射定律,第一次应用得到:

F(a,b;c;z)=(1z)aF(a,cb;c;z1z)F(a,b;c;z)=(1z)aF(a,cb;c;z1z)

注意到上参数可以交换,且 z1zz1z 算子的平方是 11

再一次应用得到:

(1z)ab+cF(cb,ca;c;z)(1z)ab+cF(cb,ca;c;z)

29.证明:合流超几何函数满足:

ezF(a;b;z)=F(ba;b;z)ezF(a;b;z)=F(ba;b;z)

展开:

(k0zkk!)(k0a¯k(1)kzkb¯kk!)=k0(ba)¯kzkb¯kk!(k0zkk!)(k0a¯¯¯k(1)kzkb¯¯¯kk!)=k0(ba)¯¯¯kzkb¯¯¯kk!

考虑两边取 [zk][zk]

ki=0(1)ia¯ib¯ii!(ki)!=(ba)¯kb¯kk!ki=0(1)ia¯ib¯ii!(ki)!=(ba)¯¯¯kb¯¯¯kk!

考虑左边化为超几何形式。然后应用高斯超几何函数的公式。

LHS=F(a,k;b;1)k!=(ba)¯kk!b¯kLHS=F(a,k;b;1)k!=(ba)¯¯¯kk!b¯¯¯k

得证。

30.解超几何微分方程

zF(z)+F(z)=11zzF(z)+F(z)=11z

常数项很讨厌,考虑再次求导。

z(1z)F(z)+(23z)F(z)F(z)=0(1)z(1z)F′′(z)+(23z)F(z)F(z)=0(1)

考虑高斯超几何函数的微分方程:

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

考虑展开。

(zF(z)+(c1)F(z))=(ϑ+a)(zF(z)+bF(z))(zF(z)+(c1)F(z))=(ϑ+a)(zF(z)+bF(z))

zF(z)+F(z)+(c1)F(z)=azF(z)+abF(z)+z(zF(z)+F(z)+bF(z))zF′′(z)+F(z)+(c1)F(z)=azF(z)+abF(z)+z(zF′′(z)+F(z)+bF(z))

z(z1)F(z)+(az+bz+zc)F(z)+abF(z)=0(2)z(z1)F′′(z)+(az+bz+zc)F(z)+abF(z)=0(2)

(1)(1) 式与 (2)(2) 式相等,得到:

(a+b+1)z+c=23z,ab=1(a+b+1)z+c=23z,ab=1

于是

a=1,b=1,c=2a=1,b=1,c=2

F=F(1;1;2)F=F(1;1;2)

33.利用 Gosper 算法求 δkk21δkk21

t(k+1)t(k)=(k+1)(k1)(k+2)k=p(k+1)q(k)p(k)r(k+1)t(k+1)t(k)=(k+1)(k1)(k+2)k=p(k+1)q(k)p(k)r(k+1)

初始 p(k)=1,q(k)=(k+1)(k1),r(k)=(k+1)(k1)p(k)=1,q(k)=(k+1)(k1),r(k)=(k+1)(k1)

发现 (k1)|r(k),(k+1)|q(k)(k1)|r(k),(k+1)|q(k)

p(k)(k+11)21_p(k)p(k)(k+11)21–––p(k)p(k)=k,q(k)=(k1),r(k)=(k+1)p(k)=k,q(k)=(k1),r(k)=(k+1)

T(k)=(k+1)s(k)(k+1)(k1)=s(k)k1T(k)=(k+1)s(k)(k+1)(k1)=s(k)k1

p(k)=q(k)s(k+1)r(k)s(k)p(k)=q(k)s(k+1)r(k)s(k)

2p(k)=Q(k)(s(k)+s(k+1))+R(k)(s(k+1)s(k)),Q(k)=q(k)r(k)=2,R(k)=q(k)+r(k)=2k2p(k)=Q(k)(s(k)+s(k+1))+R(k)(s(k+1)s(k)),Q(k)=q(k)r(k)=2,R(k)=q(k)+r(k)=2k

[kd+d1]RHS=2d+4[kd+d1]RHS=2d+4。若其不为 00,则 degs=degpdegQ=1degs=degpdegQ=1

s(k)=ak+bs(k)=ak+b,则 k=(k1)(ak+a+b)(k+1)(ak+b)k=(k1)(ak+a+b)(k+1)(ak+b)

解得 s(k)=k+12s(k)=k+12

那么 T(k)=δkk21=2k+1k(2k2)T(k)=δkk21=2k+1k(2k2)

34.证明:

kcF(a1,,an;b1,,bm;z)k=limζ0F(c,a1,,an;ζc,b1,,bn;z)kcF(a1,,an;b1,,bm;z)k=limζ0F(c,a1,,an;ζc,b1,,bn;z)

并计算 kv(nk)(1)kkv(nk)(1)k

第一个式子显然;在 kckc 的时候其取值不变,在 k>ck>c 时上指标取 00

t(k)=(nk)(1)kt(k)=(nk)(1)k

kc(nk)(1)k=kcF(n;;z)k=limζ0F(ζc,n;c;z)kc(nk)(1)k=kcF(n;;z)k=limζ0F(ζc,n;c;z)

=limζ0Γ(ζm)Γ(ζ+n)Γ(ζ)Γ(ζm+n)=limζ0Γ(ζm)Γ(ζ+n)Γ(ζ)Γ(ζm+n)

=(1)m(n1m)=(1)m(n1m)

36.库默尔定理模板题

38.证明:

nN,a,b,c,s.t.n=(a1)+(b2)+(c3),a<b<c,a,b,cNnN,a,b,c,s.t.n=(a1)+(b2)+(c3),a<b<c,a,b,cN

cc

s.t.,(c3)n,(c+13)>ns.t.,(c3)n,(c+13)>n

此时

n(c3)<(c+13)(c3)=(c2)n(c3)<(c+13)(c3)=(c2)

再同样方法取 bb,已证 b<cb<caa 同理。

39.证明下降幂\上升幂二项式定理(这里就挂上升幂)

(x+y)¯n=k(nk)x¯ky¯nk(x+y)¯¯¯n=k(nk)x¯¯¯ky¯¯¯¯¯¯¯¯¯¯nk

暴力展开,斯特林数启动!

ni=0[ni](x+y)i=ni=0(ni)(ik=0xk[ik])(nik=0yk[ink])ni=0[ni](x+y)i=ni=0(ni)(ik=0xk[ik])(nik=0yk[ink])

考虑提取两边系数。

[xayb]LHS=[na+b](a+ba),[xayb]RHS=k(nk)[ka][nkb][xayb]LHS=[na+b](a+ba),[xayb]RHS=k(nk)[ka][nkb]

下证明其等价形式:斯特林数范德蒙德卷积

[nl+m](l+ml)=k[kl][nkm](nk)[nl+m](l+ml)=k[kl][nkm](nk)

注意到第一类斯特林数列的 OGF 很糟糕,取指数生成函数。

[nl+m](l+ml)xnn!=k[kl][nkm](nk)xnn![nl+m](l+ml)xnn!=k[kl][nkm](nk)xnn!

k[kl][nkm](nk)xnn!=[xn](k[km]xkk!)(k[kl]xkk!)k[kl][nkm](nk)xnn!=[xn](k[km]xkk!)(k[kl]xkk!)

RHS=[xn](ln(1x))mm!×(ln(1x))ll!=(ln(1x))m+lm!l!RHS=[xn](ln(1x))mm!×(ln(1x))ll!=(ln(1x))m+lm!l!

RHS(l+ml)=[xn](ln(1x))m+l(m+l)!=[nl+m]=LHS(l+ml)RHS(l+ml)=[xn](ln(1x))m+l(m+l)!=[nl+m]=LHS(l+ml)

LHS=RHS,QED.LHS=RHS,QED.

于是我们证到了原结论。

具体数学上的证明是啥神秘方法,没看懂。

45.求

kn(2kk)4kkn(2kk)4k

封闭形式。

考虑

(rk)(r1/2k)=(2kk)(2r2k)22k(rk)(r1/2k)=(2kk)(2r2k)22k

r=kr=k,则:

(k1/2k)=(2kk)22k(k1/2k)=(2kk)22k

平行求和原式

kn(k1/2k)=(n+1/2n)kn(k1/2k)=(n+1/2n)

46.求

k(2k1k)(4n2k12nk)(1)k1(2k1)(4n2k1),n>0k(2k1k)(4n2k12nk)(1)k1(2k1)(4n2k1),n>0

封闭形式。

考虑

B2(z)1=k0(2k+1k)zk2k+1B2(z)1=k0(2k+1k)zk2k+1

B2(z)1=k0(2k1k)(1)k+1zk2k1B2(z)1=k0(2k1k)(1)k+1zk2k1

乘起来二者

[z2n]B2(z)1B2(z)1=k(2k1k)(4n2k12nk)(1)k+1(2k1)(4n2k1)[z2n]B2(z)1B2(z)1=k(2k1k)(4n2k12nk)(1)k+1(2k1)(4n2k1)

考虑 LHSLHS 的封闭形式。

B2(z)=zB2(z)2+1B2(z)=zB2(z)2+1

B2(z)=114z2zB2(z)=114z2z

B12(z)=1+1+4z2B12(z)=1+1+4z2

LHS=1+14z2×1+1+4z2LHS=1+14z2×1+1+4z2

=116z24+B12(z)2+B2(z)12z+14=116z24+B12(z)2+B2(z)12z+14

116z24=k0(1/2k)(16)kz2k116z24=k0(1/2k)(16)kz2k

考虑加倍公式,是

k0(1/2k)(16)kz2k=n(2nn)4nz2n2n1k0(1/2k)(16)kz2k=n(2nn)4nz2n2n1

1/41/4 对答案没有影响。加上中间两项,取个相反数,答案就是:

(2nn)4n12n1+(4n12n)4n1(2nn)4n12n1+(4n12n)4n1

  1. 给定 n,zn,z,求常数 θθ 使得

(nk)zk(k+θ)δk(nk)zk(k+θ)δk

可超几何求和。

t(k+1)t(k)=(kn)(k+θ+1)(z)(k+1)(k+θ)t(k+1)t(k)=(kn)(k+θ+1)(z)(k+1)(k+θ)

这里不必按传统方法;

q(k)=k+θ,q(k)=z(kn),r(k)=k2z(kn)=(zkk+nz)(s+Es)+(zk+k+nz)Δsdegs=0let s=αα=1z+1,θ=nzz+1(nk)zk(knz1+z)δk=nz+1(n1k1)zk+Cq(k)=k+θ,q(k)=z(kn),r(k)=k2z(kn)=(zkk+nz)(s+Es)+(zk+k+nz)Δsdegs=0let s=αα=1z+1,θ=nzz+1(nk)zk(knz1+z)δk=nz+1(n1k1)zk+C

59.对于 m,nZ+m,nZ+,求以下式子的封闭形式。

k1(nlogmk)k1(nlogmk)

k1(nlogmk)=p0(np)k1[logmk=p]=p0(np)(mp+11mp+1)=(m1)p0(np)mp=(m1)(m+1)pk1(nlogmk)=p0(np)k1[logmk=p]=p0(np)(mp+11mp+1)=(m1)p0(np)mp=(m1)(m+1)p

60.估计 (n+mn)(n+mn),并特殊地,求出 n=mn=m 时的形式。

n!nn2πnen.(n+mn)(n+m)n+mn+m2πen+mnnenmmemnm(2π)2=(n+m)n+mnnmm2π1n+1mwhen n=m(2nn)=22n4πnn!nn2πnen.(n+mn)(n+m)n+mn+m2πen+mnnenmmemnm(2π)2=(n+m)n+mnnmm2π1n+1mwhen n=m(2nn)=22n4πn

61.证明卢卡斯定理。

my blog

62.求 (npmp)modp2,pprime,m,nZ+(npmp)modp2,pprime,m,nZ+ 的简单形式。

容易证明推广范德蒙德卷积:

k=n(riki)=(rn)k=n(riki)=(rn)

反向应用,

(npmp)=k=mp(pk1)(pk2)(pkn)(pi)[i=0]+[i=p](modp)k=mp(pk1)(pk2)(pkn)modp2=(nm)(npmp)=k=mp(pk1)(pk2)(pkn)(pi)[i=0]+[i=p](modp)k=mp(pk1)(pk2)(pkn)modp2=(nm)

为什么?考虑唯一产生贡献的就是 mmkikipp,其他是 00 的项。

63.计算下式子的封闭形式

nk=0(4)k(n+k2k)nk=0(4)k(n+k2k)

首先由

Bt(z)r1t+tBt(z)1=k0(tk+rk)zkBt(z)r1t+tBt(z)1=k0(tk+rk)zk

可以得到:

B2(z)r14z=k2k+rkzkB1(z)r+11+4z=k(rkk)zkB2(z)r14z=k2k+rkzkB1(z)r+11+4z=k(rkk)zk

然后,

kn(nkk)zk=B1(z)n+1(z)n+1B2(z)n+11+4zproof :for k>n[zk](z)n+1B2(z)n+11+4z=(1)n+1(1)kn1[zkn1]B2(z)n+114z=(1)k(2kn1kn1)=(1)k(2kn1k)=(nkk)=[zk]B1(z)n+11+4zkn(nkk)zk=B1(z)n+1(z)n+1B2(z)n+11+4zproof :for k>n[zk](z)n+1B2(z)n+11+4z=(1)n+1(1)kn1[zkn1]B2(z)n+114z=(1)k(2kn1kn1)=(1)k(2kn1k)=(nkk)=[zk]B1(z)n+11+4z

所以可以去掉高次项。现回到原式。

nk=0(4)k(n+k2k)=nk=0(n+knk)(4)k=nk=0(2nkk)(4)nk=(4)nnk=0(2nkk)(14)knk=0(4)k(n+k2k)=nk=0(n+knk)(4)k=nk=0(2nkk)(4)nk=(4)nnk=0(2nkk)(14)k

我们发现这时 1+4z0,考虑🐨使用极限。

那么展开那个公式:

kn(nkk)zk=11+4z((1+1+4z2)n+1(11+4z2)n+1)

考虑取极限 1+4z0

limz14kn(nkk)zk=limx01x((1+x2)n+1(1x2)n+1)=limx0n+1i=0(n+1i)(xi(x)i)x2n+1

发现分子的常数项是 0,那么我们只用关心分子的 [x1] 项。

就是

nk=0(4)k(n+k2k)=(4)n2n+122n=(1)n(2n+1)

感觉推导好恶臭,想扔联考(

事实上,具体数学给出了这样的推导过程:

let nk=0(4)k(n+k2k)=SnSn=2Sn1Sn2 Sn=(1)n(2n+1)

🐨 <= 🤡

64.求以下式子的封闭形式:

nk=0(nk)/k+12

先拆开高斯记号。

=k0(n2k)+(n2k+1)k+1=k0(n+12k+1)k+1=k02n+1(n+22k+2)

下证明:

k(n2k)=2n1

证:令原式等于 Bn

Bn=k(n2k)=k((n12k1)+(n12k))=2n1Bn1+Bn1=2n1

则:

k02n+1(n+22k+2)=2n+1(Bn+21)=2n+22n+1

65.证明:

k(n1k)nk(k+1)!=n

k(n1k)nk(k+1)!=nk(n1k)nnk+1(k+1)!=nnk(n1k)nk(nk)!=nn

下试求 k(n1k)nk(nk)! 的封闭形式。

t(k)=(n1k)nk(nk)!=p(k+1)p(k)q(k)r(k+1)

可得:p(k)=kn,q(k)=n,r(k)=k

2(kn)=(nk)(Es+s)+(n+k)Δs

所以 degs=0s=1

T(k)=r(k)s(k)t(k)p(k)=k(n1k)nk(nk)!kn=k(n1k)nk(nk1)!=(n1)!nk(k1)!

代入原式:

k(n1k)nk(nk)!=n1k=0t(k)=T(n)T(0)=nn(n1)!(n1)!=nn

66.对于整数 m,求以下式子的封闭形式:

0jk(1jkj)(jm)2j

枚举 L=kj

L,j0(1jL)(jm)2j

枚举 l=L

=L0j(1jL)(jm)2j=l0j(1jl)(jm)2j(2l+1)=j(jm)2jl0(1jl)(2l+1)=j(jm)2jjl=0(1)jl(2l+1)=j(jm)2j(1)jjl=0(1)l(2l+1)

考虑后面那个怎么求。考虑有限微积分:

Δ(1)x=(2)(1)x(1)xδx=(1)x+12x(1)xδx=xδ(1)x+12=x(1)x+12(1)x2δx=x(1)x+1(1)x+122=(1)x+1(2x1)4

所以后面那个式子就出来了。

jl=0(1)l(2l+1)=(1)j(j+1)

代入原式得:

j(jm)2j(1)jjl=0(1)l(2l+1)=j(jm)2j(j+1)=(m+1)j(j+1m+1)2j=2(m+1)j(jm)2j

熟知:

k0(kn)zk=zn(1z)n+1

代入 z=12,得:

=2(m+1)j(jm)2j=2(m+1)(1/2)m(1/2)m+1=4m+4

所以:

0jk(1jkj)(jm)2j=4m+4

67.求以下式子的封闭形式

nk=0((k2)2)(2nkn)

考虑化简那个嵌套组合数。

((k2)2)=(k2k)(k2k2)8=(k+1)(k2)(k1)k8=3(k+14)

代入原式

nk=0((k2)2)(2nkn)=3k(k+14)(2nkn)=3(2n+2n+5)

68.求以下式子的封闭形式:

minki=nki0mi=1(ki2)

mi=1(ki2)=k2in2

容易知道,对两个 x,y 执行 x1x,y+1y 更优的充分必要条件是 xy>1

那么最优状态一定满足

i,j[1,m],|kikj|1

可以构造:

k1=k2==knmnm=nm+1,knmnm+1=,km=nm

下证明其唯一性。设:

k1=k2=ap=K+1,ap+1==am=K

那么:

(K+1)p+(mp)K=np+mK=n

0pm,故 p=nmodK

而此时 K=nm

得证。

答案计算是 trival 的。为

m(nm2)+(nmodm)nm

70.求以下式子的封闭形式:

k(nk)(2kk)(1/2)k,nZ+

t(k)=(nk)(2kk)(1/2)k

那么:

t(k+1)t(k)2(kn)(k+1/2)(k+1)2

则原式等于 F(n,1/2;1;2)

发现这个形式不好处理,但是如何做 1/21/2 是我们熟知的。

那么令 nkk,得到另一个超几何形式:

(2)n(2nn)F(n,n;1/2n;1/2)

根据高斯恒等式,原式等于

(2)n(2nn)F(n/2,n/2;1/2n;1)

然后就等于

(2)n(2nn)Γ(1/2n)Γ(1/2)Γ2(1n2)=(2)n(2nn)π(1/2nn+12)

n 为奇数时上式二项式系数下指标为负整数,原式为 0

n 为偶数的时候:

我们必须知道 Γ 函数乘法公式。

Γ(z)Γ(z+1/2)=212zπΓ(2z)

那么

Γ(n+1/2)=(2n)!πn!4n

于是等于:

(2n)!((n/2)!)2((n)!)3=(1)n(2n)!((n/2)!)2((n)!)3=(1)n(nn/2)(2nn)

解释一下第二个等号:

1z!=(n+zn)nz

在同一极限过程中,把阶乘化开,把组合数变成连乘除以下指标阶乘,抵消一下。

考虑对于每一个原来化的阶乘取相反数变成熟知的正的,考察发现只有连乘小于零的部分变了。

然后上下除一下,发现在极限意义下只有符号变了。然后把正的值乘上 (1)n 就行了。

后补:当时我好像不知道余元公式,成为🤡

也可以使用反射定律然后库莫尔公式。

于是答案是

2n(nn/2)

最后翻书《Mathematics for the Analysis of Algorithms》发现直接生成函数就可以了。傻逼具体数学。

71.设

Sn=k0(n+km+2k)ak,A(z)=k0akzk

请用 A 表示 S(z)=k0S(k)zk。(后面还有个简单应用)

S(z)=i0k0(i+km+2k)akzi=k0aki(i+km+2k)zi=k0akzki(im+2k)zi=k0akzkzm+2k(1z)m+2k+1=k0zm+k(1z)m+2k+1ak=zm(1z)m+1A(z(1z)2)

72.求证对于 n,m,kN+

(mnk)n2kv(k)Z

v(k)k 的二进制下 1 的个数。

容易发现,(mnk)n2kv(k)=m(mn)(m2n)(mkn+n)nkv(k)k!

此时考虑对于 p|n,pprime

Vp(k)V2(k)=i1k2i=i1kkmod2i2i=i1k2ii1kmod2i2i=ki1kmod2i2i

考虑 k=(bnbn1b1b0)2

i1kmod2i2i=i1i1p=0bp2p2i=p0bp2pip+112i=p0bp2p12p=p0bp=v(k)

Vp(nkv(k))=V2(k!)Vp(k!)

故对于 p|n,Vp(m(mn)(m2n)(mkn+n)nkv(k)k!)0

对于 pn,显然 (n,p)=1;设 nn 的逆元,

m(mn)(m2n)(mkn+n)nk(mn)k_=k!(modp)

则对于 pn,Vp(m(mn)(m2n)(mkn+n)nkv(k)k!)=00

得证。

77.求对于 m>1,下式的封闭形式。

0kin1j<m(kj+1kj)

拆阶乘然后抵消一下:

1j<m(kj+1kj)=(kmk1,k2k1,k3k2,,kmkm1)

考虑确定 km 时的值。此时考虑生成函数。

多项式系数的下面其实就是 m 个数随便取非负值,取 x,乘上贡献为 1x!

km![xkm](1+x+x22!+)m=km![xkm]emx=mkm

然后对 km0n 求和得到:

0kin1j<m(kj+1kj)=nkm=0mkm=mn+11m1

78.求以下式子封闭形式:

2m2k=0(kmodm(2k+1)mod(2m+1))

试了一下上指标和下指标的关系,没有成功。考虑容斥。扩展上指标:

2m2k=0(kmodm(2k+1)mod(2m+1))=2m2+m1k=0(kmodm(2k+1)mod(2m+1))2m2+m1k=2m2+1(kmodm(2k+1)mod(2m+1))=2m2+m1k=0(kmodm(2k+1)mod(2m+1))m1k=1(k2k+2)=2m2+m1k=0(kmodm(2k+1)mod(2m+1))

注意到 gcd(m,2m+1)=1,那么 (kmodm,(2k+1)mod(2m+1)) 两两不同。

于是可以变成

2m2+m1k=0(kmodm(2k+1)mod(2m+1))=m1i=02mj=0(ij)=m1i=02i=2m1

a.求

gcd((2n1),(2n3),,(2n2n1))

考虑和式。

k(2n2k+1)=22n1

那么

gcd((2n1),(2n3),,(2n2n1))=2min(V2((2n1)),V2((2n3)),,V2((2n2n1)))=2V2(2n)

b.证明

lcm((n0),(n1),,(nn))=lcm(1,2,,n+1)n+1

考虑 Vk,若 krn+1<kr+1,则取 (nkr1),根据库莫尔定理,这是最优的。

然后从 n+1krlowbit 开始有贡献,即要减去 Vk(n+1kr)=Vk(n+1),得证。

80.证明

k,nN,(nk)(enk)k

容易证明:k!(k/e)k

nk_nk,即证。

81.证明

for θ(0,1),x[0,1],l,m,nN,m<n(1)nm1k(lk)(m+θn+k)xk>0

容易发现:

fl,m,n(1)=(1)nm1(l+m+θl+n)

这里面的负的全都抵消了:组合数展开有 nm1 个负数项。

容易发现,fl,m,n(x)=lfl1,m,n+1(x),归纳法知是负的。

得证。

82.证明:

gcd((n1k1),(nk+1),(n+1k))=gcd((n1k),(n+1k+1),(nk1))

m=nk,即:(只写了第一项)

p,min(Vp(m)Vp(m+k),)=min(Vp(k)Vp(m+k),)

我们知道,如果 Vp(a)<Vp(b),那么 Vp(a)=Vp(|a±b|)

所有如果 Vp(m)Vp(k),其更小者必定等于 0

简写为 min(x1,y1,z1)=min(x2,y2,y3),而 x1+y1+z1=x2+y2+z2

我们假设 x1x2,y,z 类似。(否则可直接去除掉)

则两边都必须是 0,此时符合;否则肯定一边全是 0,另一边全部不是,不符合和相等。

83.证明

j,k(1)j+k(j+kk+l)(rj)(nk)(s+njkmj)=(1)(n+rn+l)(srmnl)

借助生成函数,我们可以把 j,k 分开。

LHS =[xlym]j,k(1)j+k(1+x)j+kxk(rj)(nk)(1+y)s+njkyj=[xlym](1)r(xy1)n+r(1+y)srxn

这是容易的。

85.证明:

nm=1(1)m1k1<k2<<kmn(2n+k3in)=(1)n(n!)3(2nn)

即:

δ[1:n](1)δi(2n+δi8in)=(1)n(n!)3,δi{0,1}

可以证明:对于不多于 n 次多项式 f,有:

δ[1:n](1)δif(δixi)=(1)nann!xi

据说可以归纳法,但我不会。

87.证明:

ζ=eπi/mmN+

kn/m(nmkk)zmk=Bm(zm)n+1(m+1)Bm(zm)m0j<m(ζ2j+1zB1+1/m(ζ2j+1z)1/m)n+1(m+1)B1+1/m(ζ2j+1z)11

具体数学上的过程有点小问题😅

“辱骂和恐吓绝不是战斗”——鲁迅

这个式子巨长,但是没有什么厉害的地方。我们熟知:

Bt(z)rtB1t(z)+1t=(tk+rk)zk

然后:

Bm(zm)n+1(m+1)Bm(zm)m=Bm(zm)n(m+1)mBm(zm)1

这不就是 t=m,r=n 的式子吗。

=k0(nmkk)zmk

我们发现后面的式子是一样的:

0j<m(ζ2j+1zB1+1/m(ζ2j+1z)1/m)n+1(m+1)B1+1/m(ζ2j+1z)11=0j<m(ζ2j+1z)n+1m(ζ2j+1zB1+1/m(ζ2j+1z)1/m)n+1(1+1/m)B1+1/m(ζ2j+1z)1(1+1/m)+1=0j<m(ζ2j+1z)n+1mk0((n+1)/m+(1+1/m)kk)(ζ2j+1z)k=0j<m1mk>n((1+1/m)kn1kn1)(ζ2j+1z)k=1mk>n((1+1/m)kn1kn1)zk0j<m(ζ2j+1)k

后面这个式子是可以求的。

0j<m(ζ2j+1)k=ζk0j<mζ2jk=ζk0j<m(ζ2k)j=ζkζ2mk1ζ2k1=(isinkπm+coskπm)isin2kπ+cos2kπ1isin2kπm+cos2kπm1

发现分子是 0,那么原式是 0km,因为此时分母不是 0。而考虑 km 会发生什么:

(isinkπm+coskπm)isin2kπ+cos2kπ1isin2kπm+cos2kπm1=(1)k/misin2kπisin2kπm=(1)k/m2kπi2kπmi=m(1)k/m

所以

0j<m(ζ2j+1)k=(1)k/mm[km]

那么我们可以化简了!

1mk>n((1+1/m)kn1kn1)zk0j<m(ζ2j+1)k=k>n/m((m+1)kn1k)(zm)k=k>n/m(nmkk)zmk

完 全 胜 利

mkn(nmkk)zmk=Bm(zm)n+1(m+1)Bm(zm)m0j<m(ζ2j+1zB1+1/m(ζ2j+1z)1/m)n+1(m+1)B1+1/m(ζ2j+1z)11=k0(nmkk)zmkk>n/m(nmkk)zmk=kn/m(nmkk)zmk

k0(rk)/(sk)=s+1sr+1

成立条件。

显然,根据高斯超几何函数的公式,需要有:

Real(r)>Real(s+1)

91.92. 利用微分方程证明超几何变换,繁而不难,略去。(92. 应用了微分有限函数之积微分有限)

后面可能会补 86,但之后的主要是机械求和,不写了。

研究题很抽象,不知道看不看

posted @   British_Union  阅读(160)  评论(2编辑  收藏  举报
相关博文:
阅读排行:
· 地球OL攻略 —— 某应届生求职总结
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 提示词工程——AI应用必不可少的技术
· .NET周刊【3月第1期 2025-03-02】
点击右上角即可分享
微信分享提示