🆙 你好,我是DorinXL,|

DorinXL

园龄:5年4个月粉丝:22关注:7

几个函数的平均阶

几个函数的平均阶

目录

  • Chapter1 d(n)的平均阶
  • Chapter2 除数函数σα(n)的平均阶
  • Chapter3 μ(n)Λ(n)的平均阶

Chapter1 d(n)的平均阶

  • Theorem:

    • Th1:对所有x1,我们有

    nxd(n)=xlogx+(2C1)x+O(x)

    其中C是欧拉常数

    Proof:

    我们知道,d(n)=d|n1,所以我们有

    nxd(n)=nxd|n1

    我们把n=qd,那么对于q,d展开就变成了

    nxd(n)=q,dqdx1

    结合图像我们可以形象的明白,我们取的点即是在qd平面上的一些格点。

    image-20230817171835151

    也就是当n=1,2,...,[x]时,曲线qd=n上的格点的个数。

    我们固定每一个小于x的d,计算水平上q可以取的格点的数目,然后对于所有的dx求和,最终得到:

    nxd(n)=dxqxd1

    根据欧拉求和公式,我们可以得到

    qxd1=xd+O(1)

    继续用欧拉求和公式,得到、

    nxd(n)=dxqxd1=dx{xd+O(1)}=xdx1d+O(x)=x{logx+C+O(1x)}+O(x)=xlogx+O(x)

    由此得出nxd(n)xlogx,x

    我们可以加强这个式子,还是从图片入手,我们换个思路:利用对称性,格点的总数等于在直线q=d下面的格点数的2倍加上这条线上的格点数:

    image-20230817173942375

    既然qd=n,我们更关注q=d这条线,那么遍历d从1到[x],每固定d,那么q就有[xd]个格点可以选择,那么阴影区域的格点数为[xd]d个。

    于是我们得到新的式子

    nxd(n)=2dx{[xd]d}+[x]

    利用[y]=y+O(1)和求和公式的推论得到

    nxd(n)=2dx{[xd]d}+[x]=2dx{xdd+O(1)}+O(x)=2xdx1d2dxd+O(x)=2x{logx+C+O(1x)}2{x2+O(x)}+O(x)=xlogx+(2C1)x+O(x)

    这就是除数函数d(n)的部分和的迪利克雷渐进公式。

    Dirichlet's asymptotic formula for the partial sums of the divisor function d(n).

Chapter2 除数函数σα(n)的平均阶

  • 我们在上面讨论了α=0的情况,接下来我们讨论α=1α>1的情况。

  • Theorem:

    • Th1:对所有的x1我们有

    nxσ1(n)=12ζ(2)x2+O(xlogx)

    注:能够证明ζ(2)=π26,所以σ1(n)的平均阶是π2n12

    Proof:

    与上一个chapter相似,我们写成下面的这种形式:

    nxσ1(n)=nxq|nq=q,dqdxq=dxqxdq=dx{12(xd)2+O(xd)}=x22dx1d2+O(xdx1d)=x22dx{1x+ζ(2)+O(1x2)}+O(xlogx)=12ζ(2)x2+O(xlogx)

    配合欧拉求和完成求证。

    • Th2:如果x1α>0,α1,我们有

    nxσα(n)=ζ(α+1)α+1xα+1+O(xβ)

    其中,β=max{1,α}

    Proof:

    配合欧拉求和公式的推论,我们开始证明:

    nxσα(n)=nxq|nqα=dxqxdqα=dx(1α+1(dx)α+1+O(xαdα))=xα+1α+1dx1dα+1+O(xαdx1dα)=xα+1α+1{xαα+ζ(α+1)+O(xα1)}+O(xα{x1α1α+ζ(α)+O(xα)})=ζ(α+1)α+1xα+1+O(x)+O(1)+O(xα)=ζ(α+1)α+1xα+1+O(xβ)

    其中,β=max{1,α}

    • Th3:对于负数ασα(n)的平均阶,我们让α=β,β>0

      如果β>0,令δ=max{0,1β},当x>1我们有

    nxσβ(n)={ζ(β+1)x+O(xδ),β1ζ(2)x+O(logx),β=1

    其中,β=max{1,α}

    Proof:

    我们有

    nxσβ(n)=nxd|n1dβ=dx1dβqxd1=dx1dβ{xd+O(1)}=xdx1dβ+1+O(dx1dβ)

    最后一项当β=1时是O(logx),当β1时是O(xδ)

    xdx1dβ+1=x1ββ+ζ(β+1)x+O(xβ)=ζ(β+1)x+O(x1β)

Chapter3 μ(n)Λ(n)的平均阶

  • Theorem:

    • Th1:如果h=fg,令

    H(x)=nxh(n),F(x)=nxf(n),G(x)=nxg(n)

    则有

    H(x)=nxf(n)G(xn)=nxg(n)F(xn)

    Proof:

    U={0,0,x<11,x1

    于是根据广义卷积得到:F=fU,G=gU

    使用结合律得到

    fG=f(gU)=(fg)U=HgF=g(fU)=(gf)U=H

    • Th2:如果F(x)=nxf(n),则有

    nxd|nf(d)=nxf(n)[xn]=nxF(xn)

    Proof:

    我们把上面证明出来的定理中,对所有的n,g(n)=1,那么G(x)=[x]

    于是这个定理对应上一个定理的结论。

    • Th3:对所有的x1,我们有

    nxμ(n)[xn]=1nxΛ(n)[xn]=log[x]!

    ​ 当且仅当x<2时,等号成立。

    Proof:

    根据定理2:

    nxμ(n)[xn]=nxd|nμ(d)=nx[1n]=1nxΛ(n)[xn]=nxd|nΛ(d)=nxlogn=log[x]!

    • Th4:对所有的x1,我们有

    |nxμ(n)n|1

    ​ 当且仅当x<2时,等号成立。

    Proof:

    如果x<2,那么这个求和里面只有一个μ(1)=1

    假设x2,对每个实数y,我们有{y}=y[y],于是有

    1=nxμ(n)[xn]=nxμ(n)(xn{xn})=xnxμ(n)nnxμ(n)xn

    因为0{y}<1,所以

    x|nxμ(n)n|=|1+nxμ(n)xn|1+nxxn=1+x+2nxxn<1+x+[x]1=x

    两边除x得到不等号

    • Th5:(勒让德等式 Legendre's identity)

      对每一个x1,我们有

    [x]!=pxpα(p)

    ​ 其中,

    α(p)=m=1[xpm]

    ​ 注:α(p)的和式是有限的,因为对于p>x,[x]pm=0

    Proof:

    如果n不是素数幂,那么有Λ(n)=0

    对于素数幂的n,有Λ(pm)=logp

    所以

    log[x]!=nxΛ(n)[xn]=pxm=1[xpm]logp=pxα(p)logp

    下面利用欧拉求和公式去确定log[x]!的渐进公式。

    • Th6:如果x2,我们有

    log[x]!=xlogxx+O(logx)

    同时就有

    nxΛ(n)[xn]=xlogxx+O(logx)

    Proof:

    在欧拉求和公式里取f(t)=logt可以得到:

    nxlogn=1xlogtdx+1xt[t]tdt(x[x])logx=xlogxx+1+1xt[t]tdt+O(logx)

    1xt[t]tdt=O(1x1tdt)=O(logx)

    接下来是上一个定理的推论:

    • Th7:如果x2,我们有

    px[xp]logp=xlogx+O(x)

    Proof:

    如果n不是素数幂,那么有Λ(n)=0

    nx[xn]Λ(n)=pm=1[xpm]Λ(pm)=pm=1[xpm]logp=px[xp]logp+pxm=2[xpm]logp

    接下来证明等式最后一项为O(x)

    pxm=2[xpm]logp=pxlogpm=2[xpm]pxlogpm=2xpm=xpxlogpm=2(1p)m=xpxlogpp(p1)xn=2lognn(n1)=O(x)

    于是证明了

    nx[xn]Λ(n)=px[xp]logp+O(x)

    结合上一个定理

    nxΛ(n)[xn]=xlogxx+O(logx)

    可以得到

    px[xp]logp+O(x)=xlogxx+O(logx)px[xp]logp=xlogx+O(x)

本文作者:DorinXL

本文链接:https://www.cnblogs.com/DorinXL/p/17724738.html

版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。

posted @   DorinXL  阅读(32)  评论(0编辑  收藏  举报
历史上的今天:
2022-09-23 使用mixamo为obj文件绑骨并导入unity
👨‍💼
点击右上角即可分享
微信分享提示
评论
收藏
关注
推荐
深色
回顶
收起
  1. 1 404 not found REOL
404 not found - REOL
00:00 / 00:00
An audio error has occurred.

作曲 : Reol

作词 : Reol

fade away...do over again...

fade away...do over again...

歌い始めの一文字目 いつも迷ってる

歌い始めの一文字目 いつも迷ってる

どうせとりとめのないことだけど

伝わらなきゃもっと意味がない

どうしたってこんなに複雑なのに

どうしたってこんなに複雑なのに

噛み砕いてやらなきゃ伝わらない

ほら結局歌詞なんかどうだっていい

僕の音楽なんかこの世になくたっていいんだよ

Everybody don't know why.

Everybody don't know why.

Everybody don't know much.

僕は気にしない 君は気付かない

何処にももういないいない

Everybody don't know why.

Everybody don't know why.

Everybody don't know much.

忘れていく 忘れられていく

We don't know,We don't know.

目の前 広がる現実世界がまた歪んだ

目の前 広がる現実世界がまた歪んだ

何度リセットしても

僕は僕以外の誰かには生まれ変われない

「そんなの知ってるよ」

気になるあの子の噂話も

シニカル標的は次の速報

麻痺しちゃってるこっからエスケープ

麻痺しちゃってるこっからエスケープ

遠く遠くまで行けるよ

安定なんてない 不安定な世界

安定なんてない 不安定な世界

安定なんてない きっと明日には忘れるよ

fade away...do over again...

fade away...do over again...

そうだ世界はどこかがいつも嘘くさい

そうだ世界はどこかがいつも嘘くさい

綺麗事だけじゃ大事な人たちすら守れない

くだらない 僕らみんなどこか狂ってるみたい

本当のことなんか全部神様も知らない

Everybody don't know why.

Everybody don't know why.

Everybody don't know much.

僕は気にしない 君は気付かない

何処にももういないいない

Everybody don't know why.

Everybody don't know why.

Everybody don't know much.

忘れていく 忘れられていく

We don't know,We don't know.