解析数论之数论函数【其一】
@Coding: Typora+LaTeX
@Author : DorinXL(博客)
@Time : 2023/7/24 16:54
目录
- Chapter1 数论函数介绍
- Chapter2 数论函数的狄利克雷乘积
- Chapter3 莫比乌斯反转
Chapter1 数论函数介绍
在正整数上定义的实数或复数的函数称为数论函数。
Section1 莫比乌斯函数μ(n)μ(n)
-
Definition:
根据算术基本定理,将大于1的自然数n分解为若干个质数乘积形式n=pα11pα22⋯pαkkμ(n)={1,if n=1(−1)k,if α1=α2=⋯=αk=10,otherwise根据算术基本定理,将大于1的自然数n分解为若干个质数乘积形式n=pα11pα22⋯pαkkμ(n)=⎧⎨⎩1,if n=1(−1)k,if α1=α2=⋯=αk=10,otherwise对于出现莫比乌斯函数的计算来说,我们更加关注于组成n的若干素数乘积,他们的幂是否都是1。
-
Examples:
n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
---|---|---|---|---|---|---|---|---|---|---|
μ(n)μ(n) | 1 | -1 | -1 | 0 | -1 | 1 | -1 | 0 | 0 | 1 |
-
Theorem:
- Th1:如果n≥1n≥1,我们有:∑d|nμ(d)=[1n]={1,if n=10,if n>1∑d|nμ(d)=[1n]={1,if n=10,if n>1
Proof:
对于n=1n=1,显然成立。
对于n>1n>1,将nn分解为n=pα11⋯pαkkn=pα11⋯pαkk,我们知道pαpα的因子只能是1,p,p2,⋯,pα1,p,p2,⋯,pα,而nn由很多素数幂组成,在此基础上,我们只考虑那些会让μμ不为0的展开项,所以我们展开求和公式:
∑d|nμ(d)=μ(1)+μ(p1)+⋯+μ(pk)+μ(p1p2)+⋯+μ(p1p2p3)+⋯+μ(p1⋯pk)=1+C1k(−1)+C2k(−1)2+⋯+Ckk(−1)k=(1−1)k(二项式定理)=0∑d|nμ(d)=μ(1)+μ(p1)+⋯+μ(pk)+μ(p1p2)+⋯+μ(p1p2p3)+⋯+μ(p1⋯pk)=1+C1k(−1)+C2k(−1)2+⋯+Ckk(−1)k=(1−1)k(二项式定理)=0 - Th1:如果n≥1n≥1,我们有:
-
Note:
二项式定理:
(x+y)n=(n0)xny0+(n1)xn−1y1+⋯+(nn)x0yn(x+y)n=(n0)xny0+(n1)xn−1y1+⋯+(nn)x0yn
Section2 欧拉函数φ(n)φ(n)
-
Definition:
欧拉函数φ(n)φ(n)被定义为不大于n并且与n互素的数的个数。
-
Examples:
n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
---|---|---|---|---|---|---|---|---|---|---|
φ(n)φ(n) | 1 | 1 | 2 | 2 | 4 | 2 | 6 | 4 | 6 | 4 |
-
Theorem:
- Th1:如果n≥1n≥1,我们有:∑d|nφ(d)=n∑d|nφ(d)=n
Proof:
用SS表示不大于nn的集合{1,2,⋯,n}{1,2,⋯,n},接下来将SS分解为若干个不相交的集合AA,A(d)={k|(k,n)=d,1≤k≤n}A(d)={k|(k,n)=d,1≤k≤n}。集合A里面装的是对于nn的所有因子,在SS的范围内与nn的最大公因数为相同nn的因子的数放在一起。例如对于n=10n=10:
S={1,2,3,4,5,6,7,8,9,10}A(1)={1,3,7,9}A(2)={2,4,6,8}A(5)={5}A(10)={10}S={1,2,3,4,5,6,7,8,9,10}A(1)={1,3,7,9}A(2)={2,4,6,8}A(5)={5}A(10)={10}现在我们想用一个符号表示AA中的符号,我们选择f(d)f(d)来表示A(d)A(d)中的个数,那么就有:
∑d|nf(d)=n∑d|nf(d)=n现在我们把(k,n)=d,0<k≤n(k,n)=d,0<k≤n转化为(kd,nd)=1,0<kd≤nd(kd,nd)=1,0<kd≤nd,如此一来我们就可以找到A(d)A(d)的数量与φ(kd)φ(kd)之间的关系:f(d)=φ(nd)f(d)=φ(nd),于是
∑d|nφ(nd)=n=∑d|nφ(d)∑d|nφ(nd)=n=∑d|nφ(d)- Th2:如果n≥1n≥1,我们有:φ(n)=∑d|nμ(d)ndφ(n)=∑d|nμ(d)nd
Proof:
基于欧拉函数的定义我们可以将其改写为:
φ(n)=n∑k=1[1(n,k)]φ(n)=n∑k=1[1(n,k)]我们使用莫比乌斯函数的定理来改写这个式子:
φ(n)=n∑k=1[1(n,k)]=n∑k=1∑d|(n,k)μ(d)=n∑k=1∑d|nd|kμ(d)φ(n)=n∑k=1[1(n,k)]=n∑k=1∑d|(n,k)μ(d)=n∑k=1∑d|nd|kμ(d)如何解读上面的求和条件呢?对于n的一个固定的因子d,我们需要满足k,1≤k≤nk,1≤k≤n是d的倍数求和。所以我们用k=qdk=qd代替:
φ(n)=∑d|nnd∑q=1μ(d)=∑d|nμ(d)nd∑q=11=∑d|nμ(d)ndφ(n)=∑d|nnd∑q=1μ(d)=∑d|nμ(d)nd∑q=11=∑d|nμ(d)nd- Th3:如果n>1n>1,我们有:φ(n)=n∏p|n(1−1p)φ(n)=n∏p|n(1−1p)
Proof:
对于n=1n=1,没有素数整除11,这个式子没有意义。
对于n>1n>1,令p1,p2,⋯,prp1,p2,⋯,pr为n的不同素因子,那么我们有:
∏p|n(1−1p)=n∏i=1(1−1pi)=(1−1p1)(1−1p2)⋯(1−1pr)=1−∑1pi+∑1pipj−∑1pipjpk+⋯+(−1)kpip2⋯pk(这一行的分母可以看作n的因子,分子可以看作μ(d),因为如果出现素数平方幂会等于0)=∑d|nμ(d)d∏p|n(1−1p)=n∏i=1(1−1pi)=(1−1p1)(1−1p2)⋯(1−1pr)=1−∑1pi+∑1pipj−∑1pipjpk+⋯+(−1)kpip2⋯pk(这一行的分母可以看作n的因子,分子可以看作μ(d),因为如果出现素数平方幂会等于0)=∑d|nμ(d)d所以有
φ(n)=∑d|nμ(d)nd=n∑d|nμ(d)d=n∏p|n(1−1p)φ(n)=∑d|nμ(d)nd=n∑d|nμ(d)d=n∏p|n(1−1p)- Th4:欧拉函数有如下性质:
- 对于素数PP与α≥1α≥1,有φ(Pα)=pα−pα−1φ(Pα)=pα−pα−1。
- φ(mn)=φ(m)φ(n)(dφ(d))φ(mn)=φ(m)φ(n)(dφ(d)),这里d=(m,n)d=(m,n)。
- φ(mn)=φ(m)φ(n)φ(mn)=φ(m)φ(n),如果(m,n)=1(m,n)=1。
- a|ba|b得出φ(a)|φ(b)φ(a)|φ(b)
- 当n≥3n≥3时,φ(n)φ(n)是偶数,而且,如果nn有rr个不同的奇素因子时,2r|φ(n)2r|φ(n)。
Proof:
4.1:在φ(n)=n∏p|n(1−1p)φ(n)=n∏p|n(1−1p)中取n=Pαn=Pα得证。
4.2:假设有m,nm,n两个整数,mnmn积的每一个素因数也是mm或者nn的素因数,我们将p|m,p|np|m,p|n的素因子pp遍历出来,会出现重复的因子,为了防止多余计算,我们将多出来的这些出现过的因子除去,于是结合Th3就有了下面的式子:
φ(mn)mn=∏p|mn(1−1p)=∏p|m(1−1p)∏p|n(1−1p)∏p|(m,n)(1−1p)=φ(m)mφ(n)nφ(d)dφ(mn)mn=∏p|mn(1−1p)=∏p|m(1−1p)∏p|n(1−1p)∏p|(m,n)(1−1p)=φ(m)mφ(n)nφ(d)d4.3:是4.2的特殊情况
4.4:由a|ba|b我们得出b=ac,1≤c≤bb=ac,1≤c≤b,如果c=bc=b,那么a=1a=1,φ(a)|φ(b)φ(a)|φ(b)成立。
如果c<bc<b:
φ(b)=φ(ac)=φ(a)φ(c)dφ(d)=dφ(a)φ(c)φ(d),d=(a,c)φ(b)=φ(ac)=φ(a)φ(c)dφ(d)=dφ(a)φ(c)φ(d),d=(a,c)接下来使用归纳法,假设对所有小于bb的整数,φ(a)|φ(b)φ(a)|φ(b)成立,那么作为小于bb的cc自然满足这一式子,那么既然d=(a,c),d|cd=(a,c),d|c,于是就有φ(d)|φ(c)φ(d)|φ(c),于是上面式子的右侧就变成了φ(a)φ(a)的倍数。于是φ(a)|φ(b)φ(a)|φ(b)成立。
4.5:我们假设n是偶数,n=2α,α≥2n是偶数,n=2α,α≥2,那么由4.1我们知道φ(n)φ(n)肯定是偶数。如果n至少有一个奇数素因子,我们写:
φ(n)=n∏p|np−1p=n∏p|np∏p|n(p−1)=C(n)∏p|n(p−1)φ(n)=n∏p|np−1p=n∏p|np∏p|n(p−1)=C(n)∏p|n(p−1)对于上面的式子,C(n)C(n)是一个整数,而∏p|n(p−1)∏p|n(p−1)是一个偶数(因为有至少一个素因子的贡献),所以φ(n)φ(n)是偶数。
如果nn有rr个不同的奇素因子时,每个素因子都能在上面的式子∏p|n(p−1)∏p|n(p−1)中提供一个因子2,于是就有2r|φ(n)2r|φ(n)
- Th1:如果n≥1n≥1,我们有:
Section3 恒等函数I(n)I(n)、单位函数u(n)u(n)
-
Definition:
I(n)=[1n]={1,n=10,n>1u(n)=1I(n)=[1n]={1,n=10,n>1u(n)=1 -
Theorem:
- Th1:(该定理需要前置知识狄利克雷卷积)对所有的ff,我们有I∗f=f∗I=fI∗f=f∗I=f
Proof:
(f∗I)(n)=∑d|nf(d)I(nd)=∑d|nf(d)[dn]=f(n)(f∗I)(n)=∑d|nf(d)I(nd)=∑d|nf(d)[dn]=f(n)
Section4 Mangoldt(曼戈尔特函数)Λ(n)Λ(n)
-
Definition:对每一个整数n≥1n≥1,我们定义
Λ(n)={logp如果n=pm,p为素数,m≥10其他Λ(n)={logp如果n=pm,p为素数,m≥10其他 -
Examples:
n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
---|---|---|---|---|---|---|---|---|---|---|
Λ(n)Λ(n) | 0 | log2log2 | log3log3 | log2log2 | log5log5 | 0 | log7log7 | log2log2 | log3log3 | 0 |
-
Theorem:
- Th1:若n≥1n≥1,我们有logn=∑d|nΛ(d)logn=∑d|nΛ(d)
Proof:
n=1n=1,两边都是0,成立
n>1n>1,算术基本定理:n=r∏k=1pαkkn=r∏k=1pαkk,两边取对数:
logn=r∑k=1αklogpklogn=r∑k=1αklogpk现在我们关注要证明的式子的右端logn=∑d|nΛ(d)logn=∑d|nΛ(d),对于Λ()Λ()来说,非零的项来自pmk(m=1,2,⋯,αk;k=1,2,⋯,r)pmk(m=1,2,⋯,αk;k=1,2,⋯,r),于是
∑d|nΛ(d)=r∑k=1αk∑m=1Λ(pmk)=r∑k=1αk∑m=1logpk=r∑k=1αklogpk=logn∑d|nΛ(d)=r∑k=1αk∑m=1Λ(pmk)=r∑k=1αk∑m=1logpk=r∑k=1αklogpk=logn- Th2:若n≥1n≥1,我们有Λ(n)=∑d|nμ(d)lognd=−∑d|nμ(d)logdΛ(n)=∑d|nμ(d)lognd=−∑d|nμ(d)logd
Proof:
对上面的定理使用莫比乌斯反转:
logn=∑d|nΛ(d)⇕Λ(d)=∑d|nμ(d)log(nd)=logn∑d|nμ(d)−∑d|nμ(d)logd=logn⋅I(n)−∑d|nμ(d)logdlogn=∑d|nΛ(d)⇕Λ(d)=∑d|nμ(d)log(nd)=logn∑d|nμ(d)−∑d|nμ(d)logd=logn⋅I(n)−∑d|nμ(d)logd对于所有的n,logn⋅I(n)=0logn⋅I(n)=0,所以证明完毕。
- Th1:若n≥1n≥1,我们有
Section5 Liouville(刘维尔函数)λ(n)λ(n)
-
Definition:我们规定λ(1)=1λ(1)=1,如果n=pα11⋯pαkkn=pα11⋯pαkk,我们规定
λ(n)=(−1)α1+⋯+αkλ(n)=(−1)α1+⋯+αk -
Theorem:
- Th1:对每一个n>1n>1,我们有∑d|nλ(d)={1n是平方数0其他∑d|nλ(d)={1n是平方数0其他
Proof:
令g(n)=∑d|nλ(d)g(n)=∑d|nλ(d),gg是积性的。运用算术基本定理我们只需要确定g(pα)g(pα):
g(pα)=∑d|pαλ(d)=1+λ(p)+λ(p2)+⋯+λ(pα)=1−1+1−⋯+(−1)α={1α是奇数0α是偶数g(pα)=∑d|pαλ(d)=1+λ(p)+λ(p2)+⋯+λ(pα)=1−1+1−⋯+(−1)α={1α是奇数0α是偶数在这种情况下,n=r∏k=1pαkk,g(n)=r∏k=1g(pαkk)n=r∏k=1pαkk,g(n)=r∏k=1g(pαkk),如果有指数αα是奇数,那么g(n)=0g(n)=0,如果所有的αα都是偶数,那么g(n)=1g(n)=1。
- Th1:对每一个n>1n>1,我们有
Section6 除数函数σα(n)σα(n)
-
Definition:对于实数或复数αα以及任意整数n>1n>1,我们规定:
σα(n)=∑d|ndασα(n)=∑d|ndα当α=0α=0时,σ0(n)σ0(n)是n的因子个数,常用d(n)d(n)表示。
当α=1α=1时,σ1(n)σ1(n)是n的因子之和,常用σ(n)σ(n)表示。
-
Theorem:
- Th1:对n≥1n≥1,我们有σ−1α(n)=∑d|ndαμ(d)μ(nd)σ−1α(n)=∑d|ndαμ(d)μ(nd)
Proof:
证明需要狄利克雷卷积相关知识:
σα(n)=Nα∗uσα(n)=Nα∗u
σ−1α(n)=(Nα∗u)−1=(μNα)∗u−1=(μNα)∗μσ−1α(n)=(Nα∗u)−1=(μNα)∗u−1=(μNα)∗μ- Th2:σα(pa)={a+1if α=1,pα(a+1)−1pα−1if α≠1σα(pa)=⎧⎪⎨⎪⎩a+1if α=1,pα(a+1)−1pα−1if α≠1
Proof:
σα(pa)=1α+pα+p2α+⋯+paα={a+1if α=1,pα(a+1)−1pα−1if α≠1σα(pa)=1α+pα+p2α+⋯+paα=⎧⎪⎨⎪⎩a+1if α=1,pα(a+1)−1pα−1if α≠1 - Th1:对n≥1n≥1,我们有
Chapter2 数论函数的狄利克雷乘积
Section1 狄利克雷卷积
-
Definition:如果ff和gg是两个数论函数,我们规定他们的狄利克雷卷积由下面的等式确定:
h(n)=(f∗g)(n)=∑d|nf(d)g(nd)h(n)=(f∗g)(n)=∑d|nf(d)g(nd)上面的式子还可以改写成其他形式:
h(n)=∑d|nf(d)g(nd)=∑d|nf(nd)g(d)=∑xy=nf(x)g(y)h(n)=∑d|nf(d)g(nd)=∑d|nf(nd)g(d)=∑xy=nf(x)g(y) -
Examples:
-
Ex1:μ∗u=Iμ∗u=I
Proof:
μ∗u=∑d|nμ(d)u(nd)=∑d|nμ(d)=Iμ∗u=∑d|nμ(d)u(nd)=∑d|nμ(d)=I -
Ex2:μ∗N=φμ∗N=φ
Proof:
μ∗N=∑d|nμ(d)N(nd)=∑d|nμ(d)nd=φμ∗N=∑d|nμ(d)N(nd)=∑d|nμ(d)nd=φ -
Ex3:u∗N=σ1u∗N=σ1
Proof:
u∗N=∑d|nu(d)N(nd)=∑d|n(d)=σ1u∗N=∑d|nu(d)N(nd)=∑d|n(d)=σ1 -
Ex4:u∗u=σ0u∗u=σ0
Proof:
u∗u=∑d|nu(d)u(nd)=∑d|n1=σ0u∗u=∑d|nu(d)u(nd)=∑d|n1=σ0 -
Ex5:u∗φ=Nu∗φ=N
Proof:
u∗φ=∑d|nu(d)φ(nd)=∑d|nφ(d)=Nu∗φ=∑d|nu(d)φ(nd)=∑d|nφ(d)=N
-
-
Theorem:
- Th1:狄利克雷卷积满足交换律,即f∗g=g∗ff∗g=g∗f
Proof:
f∗g=∑xy=nf(x)g(y)=∑xy=ng(x)f(y)=f∗gf∗g=∑xy=nf(x)g(y)=∑xy=ng(x)f(y)=f∗g- Th2:狄利克雷卷积满足结合律,即(f∗g)∗h=f∗(g∗h)(f∗g)∗h=f∗(g∗h)
Proof:
(f∗g)∗h=∑xy=n(f∗g)(x)h(y)=∑xy=n∑ab=xf(a)g(b)h(y)=∑abc=nf(a)g(b)h(c)=f∗(g∗h)(展开过程同上)(f∗g)∗h=∑xy=n(f∗g)(x)h(y)=∑xy=n∑ab=xf(a)g(b)h(y)=∑abc=nf(a)g(b)h(c)=f∗(g∗h)(展开过程同上)- Th3:狄利克雷卷积满足分配律,即f∗(g+h)=f∗g+f∗hf∗(g+h)=f∗g+f∗h
Proof:
f∗(g+h)=∑xy=nf(x)(g+h)(y)=∑xy=nf(x)(g(y)+h(y))=∑xy=nf(x)g(y)+∑xy=nf(x)h(y)=f∗g+f∗hf∗(g+h)=∑xy=nf(x)(g+h)(y)=∑xy=nf(x)(g(y)+h(y))=∑xy=nf(x)g(y)+∑xy=nf(x)h(y)=f∗g+f∗h
Section2 狄利克雷逆
-
Definition:如果ff是一个数论函数且f(1)≠0f(1)≠0,则存在唯一的一个称为ff的狄利克雷逆函数的数论函数,使
f∗−1∗f=f−1∗f=If∗−1∗f=f−1∗f=I -
Theorem:
- Th1:我们计算f−1f−1的递推公式:f−1={1f(1)n=1−1f(1)∑d|nd<nf(nd)f−1(d)otherwisef−1=⎧⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪⎩1f(1)n=1−1f(1)∑d|nd<nf(nd)f−1(d)otherwise
Proof:
对n=1n=1,(f∗f−1)(1)=I(1)=1⇒f−1=1f(1)(f∗f−1)(1)=I(1)=1⇒f−1=1f(1)
对n≥1n≥1,
(f∗f−1)(n)=∑d|nf(d)f−1(nd)=0f(1)f−1(n)+∑d|nd<nf(nd)f−1(d)=0f−1(n)=−1f(1)∑d|nd<nf(nd)f−1(d)(f∗f−1)(n)=∑d|nf(d)f−1(nd)=0f(1)f−1(n)+∑d|nd<nf(nd)f−1(d)=0f−1(n)=−1f(1)∑d|nd<nf(nd)f−1(d)- Th2:(f∗g)−1=f−1∗g−1(f∗g)−1=f−1∗g−1
Proof:
(f∗g)−1=f−1∗g−1(f∗g)∗(f∗g)−1=f−1∗g−1∗(f∗g)I=f−1∗f∗g∗g−1I=I∗I=I(f∗g)−1=f−1∗g−1(f∗g)∗(f∗g)−1=f−1∗g−1∗(f∗g)I=f−1∗f∗g∗g−1I=I∗I=I - Th1:我们计算f−1f−1的递推公式:
Chapter3 莫比乌斯反转
-
Definition:ff,gg是两个数论函数,如果有:
f(n)=∑d|ng(d)⇕g(n)=∑d|nf(d)μ(nd)f(n)=∑d|ng(d)⇕g(n)=∑d|nf(d)μ(nd)Proof:
证明⇓⇓:
f(n)=∑d|ng(d)即f=g∗uf∗μ=g∗u∗μ=g∗(u∗μ)=g∗I=gf(n)=∑d|ng(d)即f=g∗uf∗μ=g∗u∗μ=g∗(u∗μ)=g∗I=g证明⇑⇑:
g(n)=∑d|nf(d)μ(nd)即g=f∗μu∗g=u∗f∗μ=(u∗μ)∗f=I∗f=fg(n)=∑d|nf(d)μ(nd)即g=f∗μu∗g=u∗f∗μ=(u∗μ)∗f=I∗f=f
本文作者:DorinXL
本文链接:https://www.cnblogs.com/DorinXL/p/17577874.html
版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步