[数论数学]莫比乌斯反演定理

  • Post author:
  • Post category:其他

近期学习了莫比乌斯反演,算是一个学习笔记吧…

在这里首先要说明:

1:本文讨论的所有函数为数论函数,即定义域为

D

=

N

D=N^*

D=N的函数;

2:

d

n

f

(

d

)

\sum \limits_{d|n}f(d)

dnf(d)表示

d

d

d取遍

n

n

n的所有正因子,再将所有的

f

(

d

)

f(d)

f(d)相加,例如当

n

=

6

n=6

n=6时,

d

n

f

(

d

)

=

f

(

1

)

+

f

(

2

)

+

f

(

3

)

+

f

(

6

)

\sum \limits_{d|n}f(d)=f(1)+f(2)+f(3)+f(6)

dnf(d)=f(1)+f(2)+f(3)+f(6)

3:

i

=

1

k

a

i

=

a

1

×

a

2

×

.

.

.

×

a

k

\prod \limits_{i=1}^ka_i=a_1 \times a_2 \times … \times a_k

i=1kai=a1×a2×...×ak,即所有

a

i

a_i

ai的乘积;

4:

[

p

]

[p]

[p](其中p是一个表达式)表示当p成立时,值为1;否则值为0;

5:

(

a

,

b

)

(a,b)

(a,b)表示

g

c

d

(

a

,

b

)

gcd(a,b)

gcd(a,b),即

a

,

b

a,b

a,b的最大公因数;

6:几个在狄利克雷卷积中会出现的基本函数:

  • ε

    (

    n

    )

    =

    [

    n

    =

    1

    ]

    \varepsilon(n)=[n=1]

    ε(n)=[n=1],也就是

    ε

    (

    n

    )

    =

    {

    1

    n

    =

    1

    0

    n

    1

    \varepsilon(n)=\begin{cases}1\qquad n=1\\0 \qquad n \neq 1 \end{cases}

    ε(n)={1n=10n̸=1

  • I

    (

    n

    )

    =

    1

    I(n)=1

    I(n)=1

  • i

    d

    (

    n

    )

    =

    n

    id(n)=n

    id(n)=n

一、积性函数

定义:如果一个函数

f

(

n

)

f(n)

f(n)满足:当

(

a

,

b

)

=

1

(a,b)=1

(a,b)=1时,有

f

(

a

b

)

=

f

(

a

)

×

f

(

b

)

f(ab)=f(a) \times f(b)

f(ab)=f(a)×f(b),则称

f

(

n

)

f(n)

f(n)为积性函数。

进一步,如果对于任意数

a

,

b

a,b

a,b,都有

f

(

a

b

)

=

f

(

a

)

×

f

(

b

)

f(ab)=f(a) \times f(b)

f(ab)=f(a)×f(b),则称

f

(

n

)

f(n)

f(n)为完全积性函数。

积性函数的基本性质:对于任意正整数

n

1

n \neq 1

n̸=1,根据唯一分解定理,可以进行分解素因数:即

n

=

p

1

c

1

.

.

.

p

k

c

k

=

i

=

1

k

p

i

c

i

n=p_1^{c_1}…p_k^{c_k}=\prod \limits_{i=1}^kp_i^{c_i}

n=p1c1...pkck=i=1kpici,显然所有的

p

i

c

i

p_i^{c_i}

pici两两互素,因此有:

f

(

n

)

=

f

(

p

1

c

1

)

×

.

.

.

×

f

(

p

k

c

k

)

=

i

=

1

k

f

(

p

i

c

i

)

f(n)=f(p_1^{c_1}) \times … \times f(p_k^{c_k})=\prod \limits_{i=1}^kf(p_i^{c_i})

f(n)=f(p1c1)×...×f(pkck)=i=1kf(pici)

因此,我们可以用类似线性筛素数的方法线性筛出某一积性函数的值,只需要知道对于每一个素数以及素数的方幂结果为多少即可。

由积性函数以及完全积性函数的定义,

ε

(

n

)

,

I

(

n

)

,

i

d

(

n

)

\varepsilon(n),I(n),id(n)

ε(n),I(n),id(n)显然皆为完全积性函数。

证明积性函数的常见方法是利用上面的唯一分解式,将

a

,

b

a,b

a,b唯一分解,再用

(

a

,

b

)

=

1

(a,b)=1

(a,b)=1的条件可以得到分解出来的所有素数两两不等,从而继续证明。

例如:证明

d

(

n

)

d(n)

d(n)

n

n

n的正因子个数为积性函数。

证明:取

(

a

,

b

)

=

1

(a,b)=1

(a,b)=1,将

a

,

b

a,b

a,b分解得:

a

=

i

=

1

k

p

i

c

i

a=\prod \limits_{i=1}^kp_i^{c_i}

a=i=1kpici,且

b

=

i

=

1

l

q

i

d

i

b=\prod \limits_{i=1}^lq_i^{d_i}

b=i=1lqidi,则

p

i

,

q

i

p_i,q_i

pi,qi两两不等(由

(

a

,

b

)

=

1

(a,b)=1

(a,b)=1可得),于是:

d

(

a

)

=

i

=

1

k

(

c

i

+

1

)

d(a)=\prod \limits_{i=1}^k(c_i+1)

d(a)=i=1k(ci+1),且

d

(

b

)

=

i

=

1

l

(

d

i

+

1

)

d(b)=\prod \limits_{i=1}^l(d_i+1)

d(b)=i=1l(di+1)
由于

p

i

,

q

i

p_i,q_i

pi,qi两两不等,所以

d

(

a

b

)

=

i

=

1

k

(

c

i

+

1

)

i

=

1

l

(

d

i

+

1

)

=

d

(

a

)

×

d

(

b

)

d(ab)=\prod \limits_{i=1}^k(c_i+1)\prod \limits_{i=1}^l(d_i+1)=d(a) \times d(b)

d(ab)=i=1k(ci+1)i=1l(di+1)=d(a)×d(b),即

d

(

n

)

d(n)

d(n)是积性函数。

*积性函数另一个不太常用的性质:若

f

(

n

)

f(n)

f(n)是积性函数,则或者

f

(

n

)

f(n)

f(n)是值为

0

0

0的常值函数,或者

f

(

1

)

=

1

f(1)=1

f(1)=1

证明十分容易,因为

(

n

,

1

)

=

1

(n,1)=1

(n,1)=1,所以

f

(

n

)

=

f

(

1

)

×

f

(

n

)

f(n)=f(1) \times f(n)

f(n)=f(1)×f(n),则

f

(

n

)

=

0

f(n)=0

f(n)=0或者

f

(

1

)

=

1

f(1)=1

f(1)=1

二、莫比乌斯函数

定义:若将

n

n

n分解素因数为

n

=

p

1

c

1

.

.

.

p

k

c

k

=

i

=

1

k

p

i

c

i

n=p_1^{c_1}…p_k^{c_k}=\prod \limits_{i=1}^kp_i^{c_i}

n=p1c1...pkck=i=1kpici,则莫比乌斯函数定义为:

μ

(

n

)

=

{

1

n

=

1

(

1

)

k

c

1

=

.

.

.

=

c

k

=

1

0

m

a

x

(

c

1

,

.

.

.

,

c

k

)

>

1

\mu(n)=\begin{cases}1\qquad\qquad n=1\\ (-1)^k\qquad c_1=…=c_k=1\\ 0 \qquad\qquad max(c_1,…,c_k)>1\end{cases}

μ(n)=1n=1(1)kc1=...=ck=10max(c1,...,ck)>1

(最后一行的表示其实等价于:若

n

n

n有一个因子为完全平方数,则

μ

(

n

)

=

0

\mu (n)=0

μ(n)=0

莫比乌斯函数是积性函数。

证明:取

(

a

,

b

)

=

1

(a,b)=1

(a,b)=1,对

a

,

b

a,b

a,b进行讨论:

  • a

    =

    1

    a=1

    a=1

    b

    =

    1

    b=1

    b=1,不妨设

    a

    =

    1

    a=1

    a=1,则

    μ

    (

    a

    )

    =

    1

    \mu(a)=1

    μ(a)=1,因为

    μ

    (

    b

    )

    =

    μ

    (

    b

    )

    \mu(b)=\mu(b)

    μ(b)=μ(b),所以

    μ

    (

    a

    b

    )

    =

    μ

    (

    a

    )

    ×

    μ

    (

    b

    )

    \mu(ab)=\mu(a) \times \mu(b)

    μ(ab)=μ(a)×μ(b)

  • μ

    (

    a

    )

    =

    0

    \mu(a)=0

    μ(a)=0

    μ

    (

    b

    )

    =

    0

    \mu(b)=0

    μ(b)=0,不妨设

    μ

    (

    a

    )

    =

    0

    \mu(a)=0

    μ(a)=0,则存在

    p

    2

    a

    p^2|a

    p2a,所以

    p

    2

    a

    b

    p^2|ab

    p2ab,所以有

    μ

    (

    a

    b

    )

    =

    0

    \mu(ab)=0

    μ(ab)=0,所以

    μ

    (

    a

    b

    )

    =

    μ

    (

    a

    )

    ×

    μ

    (

    b

    )

    \mu(ab)=\mu(a) \times \mu(b)

    μ(ab)=μ(a)×μ(b)

  • 否则,设

    a

    =

    p

    1

    p

    2

    .

    .

    .

    p

    k

    b

    =

    q

    1

    q

    2

    .

    .

    .

    q

    l

    a=p_1p_2…p_k,b=q_1q_2…q_l

    a=p1p2...pkb=q1q2...ql,则

    μ

    (

    a

    )

    =

    (

    1

    )

    k

    μ

    (

    b

    )

    =

    (

    1

    )

    l

    \mu(a)=(-1)^k,\mu(b)=(-1)^l

    μ(a)=(1)kμ(b)=(1)l,因为

    (

    a

    ,

    b

    )

    =

    1

    (a,b)=1

    (a,b)=1,所以

    p

    i

    ,

    q

    i

    p_i,q_i

    pi,qi两两不等,因此

    μ

    (

    a

    b

    )

    =

    (

    1

    )

    k

    +

    l

    \mu(ab)=(-1)^{k+l}

    μ(ab)=(1)k+l,所以:

    μ

    (

    a

    b

    )

    =

    μ

    (

    a

    )

    ×

    μ

    (

    b

    )

    \mu(ab)=\mu(a) \times \mu(b)

    μ(ab)=μ(a)×μ(b)

综上所述,

μ

(

a

b

)

=

μ

(

a

)

×

μ

(

b

)

\mu(ab)=\mu(a) \times \mu(b)

μ(ab)=μ(a)×μ(b),所以

μ

(

n

)

\mu(n)

μ(n)是积性函数。

前面已经提到,积性函数可以用线性筛的方法筛出,下面就是使用线性筛求

μ

(

1

)

\mu(1)

μ(1)

μ

(

x

)

\mu(x)

μ(x)的代码,关键部分有注释:

void init (int x) {
	mob[1]=1;    //由定义可得mob[1]=1
	for (int i=2;i<=x;i++) {
		if (!p[i]) {
			n[++cnt]=i;
			mob[i]=-1;    //若i是素数,则mob[i]=(-1)^1=-1
		}
		for (int j=1;j<=cnt;j++) {
			p[i*n[j]]=1;
			if (i%n[j]==0) {    //若i是n[j]的倍数
				mob[i*n[j]]=0;    //则i*n[j]有n[j]^2的平方因数,所以函数值为0
				break;
			} else {
				mob[i*n[j]]=-mob[i];    //i*n[j]多了素因子n[j],-1的次数加1
			}
		}
	}
}

下面是一个莫比乌斯函数的重要性质,也是莫比乌斯反演定理的核心内容:

d

n

μ

(

n

)

=

[

n

=

1

]

\sum \limits_{d|n}\mu(n)=[n=1]

dnμ(n)=[n=1]

证明:

  • n

    =

    1

    n=1

    n=1,则上式的值为

    μ

    (

    n

    )

    =

    1

    \mu(n)=1

    μ(n)=1,成立;

  • 否则,设

    n

    n

    n

    k

    k

    k个不同的素因子,则只需考虑这些素因子的组合即可(因为一旦某个素因子的个数大于1,则结果为0,不计入求和答案);那么取奇数个素因数时,莫比乌斯函数的值为-1,取偶数个素因数时,莫比乌斯函数的值为1,因此原式可转化为:

    C

    n

    0

    C

    n

    1

    +

    .

    .

    .

    +

    (

    1

    )

    n

    C

    n

    n

    C_n^0-C_n^1+…+(-1)^nC_n^n

    Cn0Cn1+...+(1)nCnn

  • 再由于

    C

    n

    m

    =

    C

    n

    1

    m

    1

    +

    C

    n

    1

    m

    (

    0

    &lt;

    m

    &lt;

    n

    )

    C_n^m=C_{n-1}^{m-1}+C_{n-1}^m\quad(0&lt;m&lt;n)

    Cnm=Cn1m1+Cn1m(0<m<n),因此原式可以化为:

C

n

0

C

n

1

0

C

n

1

1

+

.

.

.

+

(

1

)

n

1

C

n

1

n

2

+

(

1

)

n

1

C

n

1

n

1

+

(

1

)

n

C

n

n

C_n^0-C_{n-1}^0-C_{n-1}^1+…+(-1)^{n-1}C_{n-1}^{n-2}+(-1)^{n-1}C_{n-1}^{n-1}+(-1)^nC_n^n

Cn0Cn10Cn11+...+(1)n1Cn1n2+(1)n1Cn1n1+(1)nCnn

  • 此时发现中间项全部抵消,剩余部分为:

C

n

0

C

n

1

0

+

(

1

)

n

1

C

n

1

n

1

+

(

1

)

n

C

n

n

=

1

1

+

1

1

=

0

C_n^0-C_{n-1}^0+(-1)^{n-1}C_{n-1}^{n-1}+(-1)^nC_n^n=1-1+1-1=0

Cn0Cn10+(1)n1Cn1n1+(1)nCnn=11+11=0

  • d

    n

    μ

    (

    n

    )

    =

    0

    (

    n

    1

    )

    \sum \limits_{d|n}\mu(n)=0\quad(n \neq 1)

    dnμ(n)=0(n̸=1),得证。

三、狄利克雷卷积

定义:已知函数

f

(

n

)

,

g

(

n

)

f(n),g(n)

f(n),g(n),若函数

h

(

n

)

=

d

n

f

(

d

)

g

(

n

d

)

h(n)=\sum \limits_{d|n}f(d)g(\frac{n}{d})

h(n)=dnf(d)g(dn),记为

h

=

f

g

h=f*g

h=fg,称

h

h

h

f

,

g

f,g

f,g的狄利克雷卷积。

一般,我们还可以将狄利克雷卷积写成这样的形式:

h

(

n

)

=

i

j

=

n

f

(

i

)

g

(

j

)

h(n)=\sum \limits_{ij=n}f(i)g(j)

h(n)=ij=nf(i)g(j) 这更利于我们解决一些问题。

狄利克雷卷积的运算律:

(一)交换律:

f

g

=

g

f

f*g=g*f

fg=gf

证明:因为

(

f

g

)

(

n

)

=

i

j

=

n

f

(

i

)

g

(

j

)

(f*g)(n)=\sum \limits_{ij=n}f(i)g(j)

(fg)(n)=ij=nf(i)g(j)

i

,

j

i,j

i,j地位相等,因此

i

,

j

i,j

i,j位置可交换,满足交换律;

(二)结合律:

(

f

g

)

h

=

f

(

g

h

)

(f*g)*h=f*(g*h)

(fg)h=f(gh)

证明:根据狄利克雷卷积定义,有:

(

(

f

g

)

h

)

(

n

)

=

i

j

=

n

(

k

l

=

i

f

(

k

)

g

(

l

)

)

h

(

j

)

=

k

l

j

=

n

f

(

k

)

g

(

l

)

h

(

j

)

((f*g)*h)(n)=\sum \limits_{ij=n}(\sum \limits_{kl=i}f(k)g(l))h(j)=\sum \limits_{klj=n}f(k)g(l)h(j)

((fg)h)(n)=ij=n(kl=if(k)g(l))h(j)=klj=nf(k)g(l)h(j)同时有:

(

f

(

g

h

)

)

(

n

)

=

i

j

=

n

f

(

i

)

(

k

l

=

j

g

(

k

)

h

(

l

)

)

=

i

k

l

=

n

f

(

i

)

g

(

k

)

h

(

l

)

(f*(g*h))(n)=\sum \limits_{ij=n}f(i)(\sum \limits_{kl=j}g(k)h(l))=\sum \limits_{ikl=n}f(i)g(k)h(l)

(f(gh))(n)=ij=nf(i)(kl=jg(k)h(l))=ikl=nf(i)g(k)h(l)两式相等,代表

(

f

g

)

h

=

f

(

g

h

)

(f*g)*h=f*(g*h)

(fg)h=f(gh)

(三)分配律:

(

f

+

g

)

h

=

f

h

+

g

h

(f+g)*h=f*h+g*h

(f+g)h=fh+gh

将外部求和进行结合,有

f

(

i

)

+

g

(

i

)

=

(

f

(

i

)

+

g

(

i

)

)

\sum f(i)+\sum g(i)=\sum (f(i)+g(i))

f(i)+g(i)=(f(i)+g(i)),由此易证分配率成立(与莫比乌斯反演关系不大,不多做证明了)

狄利克雷卷积有一个结论,虽然这里不会用到,但也非常有用:

两个积性函数的狄利克雷卷积还是积性函数,即若

f

,

g

f,g

f,g是积性函数,则

f

g

f*g

fg也是积性函数

证明:记

h

=

f

g

h=f*g

h=fg,取

(

a

,

b

)

=

1

(a,b)=1

(a,b)=1,则根据定义有:

h

(

a

b

)

=

d

a

b

f

(

d

)

g

(

a

b

d

)

h(ab)=\sum \limits_{d|ab}f(d)g(\frac{ab}{d})

h(ab)=dabf(d)g(dab)由于

(

a

,

b

)

=

1

(a,b)=1

(a,b)=1,因此可以将

d

d

d表示成:

d

=

d

1

×

d

2

(

d

1

a

,

d

2

b

)

d=d_1 \times d_2\quad(d1|a,d2|b)

d=d1×d2(d1a,d2b),而且对于任意一组不同的

(

d

1

,

d

2

)

(d_1,d_2)

(d1,d2),所对应的

d

d

d都不同。因此将原式改写成:

h

(

a

b

)

=

d

1

a

d

2

b

f

(

d

1

d

2

)

g

(

a

b

d

1

d

2

)

h(ab)=\sum \limits_{d_1|a}\sum \limits_{d_2|b}f(d_1d_2)g(\frac{ab}{d_1d_2})

h(ab)=d1ad2bf(d1d2)g(d1d2ab)因为

(

a

,

b

)

=

1

(a,b)=1

(a,b)=1,所以

(

d

1

,

d

2

)

=

1

,

(

a

d

1

,

b

d

2

)

=

1

(d_1,d_2)=1,(\frac{a}{d_1},\frac{b}{d_2})=1

(d1,d2)=1,(d1a,d2b)=1,根据

f

,

g

f,g

f,g是积性函数可以进行拆分:

h

(

a

b

)

=

d

1

a

d

2

b

f

(

d

1

)

f

(

d

2

)

g

(

a

d

1

)

g

(

b

d

2

)

h(ab)=\sum \limits_{d_1|a}\sum \limits_{d_2|b}f(d_1)f(d_2)g(\frac{a}{d_1})g(\frac{b}{d_2})

h(ab)=d1ad2bf(d1)f(d2)g(d1a)g(d2b)将与

d

1

d_1

d1有关的式子提到外层求和:

h

(

a

b

)

=

d

1

a

(

f

(

d

1

)

g

(

a

d

1

)

d

2

b

f

(

d

2

)

g

(

b

d

2

)

)

h(ab)=\sum \limits_{d_1|a}(f(d_1)g(\frac{a}{d_1})\sum \limits_{d_2|b}f(d2)g(\frac{b}{d_2}))

h(ab)=d1a(f(d1)g(d1a)d2bf(d2)g(d2b))内外层求和分别独立,可以拆去括号,此时再由定义即可得到:

h

(

a

b

)

=

d

1

a

(

f

(

d

1

)

g

(

a

d

1

)

)

d

2

b

(

f

(

d

2

)

g

(

b

d

2

)

)

=

h

(

a

)

×

h

(

b

)

h(ab)=\sum \limits_{d_1|a}(f(d_1)g(\frac{a}{d_1}))\sum \limits_{d_2|b}(f(d2)g(\frac{b}{d_2}))=h(a) \times h(b)

h(ab)=d1a(f(d1)g(d1a))d2b(f(d2)g(d2b))=h(a)×h(b)得证。

以函数的狄利克雷卷积类比数的乘法,我们已经得到了一系列运算律,下面做进一步的分析。
在乘法中,如果两个数相乘结果为1,则称它们互为倒数,也称为逆元,我们可以类似定义狄利克雷卷积中的逆元。但在此之前,需要先搞清楚卷积中的”1″是多少,即找到一个函数,使得任意一个函数与它求卷积后仍然等于原函数。

定理:对于任意函数

f

(

n

)

f(n)

f(n)

f

ε

=

f

(

ε

(

n

)

=

[

n

=

1

]

)

f*\varepsilon=f \quad (\varepsilon(n)=[n=1])

fε=f(ε(n)=[n=1])

证明:由狄利克雷卷积定义有:

(

f

ε

)

(

n

)

=

d

n

f

(

d

)

ε

(

n

d

)

(f*\varepsilon)(n)=\sum \limits_{d|n}f(d)\varepsilon(\frac{n}{d})

(fε)(n)=dnf(d)ε(dn),又因为当且仅当

d

=

n

d=n

d=n时,

ε

(

n

d

)

=

1

0

\varepsilon(\frac{n}{d})=1 \neq 0

ε(dn)=1̸=0,因此对答案有贡献的部分只是

d

=

n

d=n

d=n,因此有:

(

f

ε

)

(

n

)

=

f

(

n

)

×

1

=

f

(

n

)

(f*\varepsilon)(n)=f(n) \times 1=f(n)

(fε)(n)=f(n)×1=f(n),得证。

因此我们定义卷积意义下的逆元为:若

f

g

=

ε

f*g=\varepsilon

fg=ε,则称

f

,

g

f,g

f,g互为逆元,记作

f

=

g

1

f=g^{-1}

f=g1

下面我们就莫比乌斯函数的逆元进行讨论。

我们已经证明了

d

n

μ

(

d

)

=

[

n

=

1

]

\sum \limits_{d|n}\mu(d)=[n=1]

dnμ(d)=[n=1],我们可以在每个求和式后乘上一个

1

1

1,将

[

n

=

1

]

[n=1]

[n=1]

ε

(

n

)

\varepsilon(n)

ε(n)代换,因此就得到了:

μ

I

=

ε

\mu*I=\varepsilon

μI=ε,也就是

I

=

μ

1

I=\mu^{-1}

I=μ1,其中

I

(

n

)

=

1

I(n)=1

I(n)=1

下面是一个联系了欧拉函数和莫比乌斯函数的定理:

φ

=

μ

i

d

(

i

d

(

n

)

=

n

)

\varphi=\mu*id\quad(id(n)=n)

φ=μid(id(n)=n)

证明:设将n进行素因数分解后为:

n

=

p

1

c

1

.

.

.

p

k

c

k

n=p_1^{c_1}…p_k^{c_k}

n=p1c1...pkck,由定义的角度出发,欧拉函数指的是小于等于n的与n互素的数的个数,那么我们只需要用n减去与n不互素的数的个数即可。
其中,

n

=

1

×

n

=

μ

(

1

)

×

i

d

(

n

)

n=1 \times n=\mu(1) \times id(n)

n=1×n=μ(1)×id(n),而与n不互素的数,其实就是所有

p

i

p_i

pi的倍数,但是

p

i

p_i

pi倍数组成的数集中有交集,因此我们需要用容斥原理,设

S

i

S_i

Si表示

p

i

p_i

pi倍数组成的集合,

S

|S|

S表示S的元素数量,则:

i

=

1

k

S

i

=

i

=

1

k

S

i

1

i

&lt;

j

n

S

i

S

j

+

.

.

.

+

(

1

)

k

1

i

=

1

k

S

i

|\bigcup\limits_{i=1}^kS_i|=\sum\limits_{i=1}^k|S_i|-\sum\limits_{1 \leq i&lt;j \leq n}|S_i\cap S_j|+…+(-1)^{k-1}|\bigcap\limits_{i=1}^kS_i|

i=1kSi=i=1kSi1i<jnSiSj+...+(1)k1i=1kSi而事实上,

(

1

)

k

(-1)^k

(1)k可以直接由莫比乌斯函数提供,而

S

i

=

n

p

i

|S_i|=\frac{n}{p_i}

Si=pin,类似的,用

n

n

n去除若干

p

i

p_i

pi的积就可以得到交集的元素个数了,而这就是卷积式:

d

n

,

d

1

(

(

1

)

×

μ

(

d

)

×

n

d

)

\sum \limits_{d|n,d \neq 1}((-1) \times\mu(d) \times\frac{n}{d})

dn,d̸=1((1)×μ(d)×dn)注意莫比乌斯函数中-1的次数与容斥原理中的次数相差1,而且若d有平方因子,则不会计入答案。因此我们就可以得到:

φ

(

n

)

=

n

d

n

,

d

1

(

(

1

)

×

μ

(

d

)

×

n

d

)

=

d

n

μ

(

d

)

i

d

(

n

d

)

=

(

μ

i

d

)

(

n

)

\varphi(n)=n-\sum \limits_{d|n,d \neq 1}((-1) \times\mu(d) \times\frac{n}{d})=\sum\limits_{d|n}\mu(d)id(\frac{n}{d})=(\mu*id)(n)

φ(n)=ndn,d̸=1((1)×μ(d)×dn)=dnμ(d)id(dn)=(μid)(n)

φ

=

μ

i

d

\varphi=\mu*id

φ=μid,得证。

由上面这个定理,我们可以得到,莫比乌斯函数与容斥原理息息相关,事实上我们在证明

μ

I

=

ε

\mu*I=\varepsilon

μI=ε时也运用到了这一性质。

在此介绍这个定理,是因为使用这个定理再加上下面即将介绍的莫比乌斯反演定理就可以得到:

i

d

=

φ

I

id=\varphi*I

id=φI,也就是:

n

=

d

n

φ

(

d

)

n=\sum \limits_{d|n}\varphi(d)

n=dnφ(d)

四、莫比乌斯反演定理

下面进入最后一个版块…真正的莫比乌斯反演定理:

F

(

n

)

=

d

n

f

(

n

)

 

 

f

(

n

)

=

d

n

F

(

d

)

μ

(

n

d

)

F(n)=\sum\limits_{d|n}f(n)\ \Leftrightarrow\ f(n)=\sum\limits_{d|n}F(d)\mu(\frac{n}{d})

F(n)=dnf(n)  f(n)=dnF(d)μ(dn)

其实只要有了前面的基础,这个结论就非常浅显简单,需要三个之前提到和一个之前没提到的基础知识作为引理:

  1. (

    f

    g

    )

    h

    =

    f

    (

    g

    h

    )

    (f*g)*h=f*(g*h)

    (fg)h=f(gh),记为

    f

    g

    h

    f*g*h

    fgh

  2. f

    ε

    =

    f

    f*\varepsilon=f

    fε=f

  3. μ

    I

    =

    ε

    \mu*I=\varepsilon

    μI=ε

  4. f

    =

    g

    f=g

    f=g,则

    f

    h

    =

    g

    h

    f*h=g*h

    fh=gh

最后一条结论非常显然,因此之前没有说明。

下面就开始我们的证明:

先证必要性(

\Rightarrow

):

F

(

n

)

=

d

n

f

(

n

)

F(n)=\sum\limits_{d|n}f(n)

F(n)=dnf(n)可以写作

F

=

f

I

F=f*I

F=fI,则根据引理4,将两边同时卷上

μ

\mu

μ

F

μ

=

(

f

I

)

μ

F*\mu=(f*I)*\mu

Fμ=(fI)μ右侧可以通过引理1将括号加在右边:

(

f

I

)

μ

=

f

(

I

μ

)

(f*I)*\mu=f*(I*\mu)

(fI)μ=f(Iμ)而通过引理3,我们知道:

f

(

I

μ

)

=

f

ε

f*(I*\mu)=f*\varepsilon

f(Iμ)=fε再使用引理2可得:

f

ε

=

f

f*\varepsilon=f

fε=f再看最初的等式,其实已经得到了结论:

f

=

F

μ

f

(

n

)

=

d

n

F

(

d

)

μ

(

n

d

)

f=F*\mu\qquad f(n)=\sum\limits_{d|n}F(d)\mu(\frac{n}{d})

f=Fμf(n)=dnF(d)μ(dn)再证充分性(

\Leftarrow

):

f

(

n

)

=

d

n

F

(

d

)

μ

(

n

d

)

f(n)=\sum\limits_{d|n}F(d)\mu(\frac{n}{d})

f(n)=dnF(d)μ(dn)可以写作

f

=

F

μ

f=F*\mu

f=Fμ,类似刚才,我们两边同时卷上

I

I

I:

f

I

=

(

F

I

)

μ

=

F

(

I

μ

)

f*I=(F*I)*\mu=F*(I*\mu)

fI=(FI)μ=F(Iμ)通过引理3和引理2可得:

F

(

I

μ

)

=

F

ε

=

F

F*(I*\mu)=F*\varepsilon=F

F(Iμ)=Fε=F于是我们就得到了结论:

F

=

f

I

F

(

n

)

=

d

n

f

(

d

)

F=f*I\qquad F(n)=\sum\limits_{d|n}f(d)

F=fIF(n)=dnf(d)

对于莫比乌斯反演定理的另一种理解是:将

I

I

I视为

μ

1

\mu^{-1}

μ1,或者将

μ

\mu

μ视为

I

1

I^{-1}

I1,这时结论更加显然,只需在两边同时卷上

I

I

I

μ

\mu

μ,同侧逆元抵消即可。
这个过程可以看做移项后乘法变除法的过程。

另一个莫比乌斯反演定理的形式:

F

(

n

)

=

n

d

f

(

d

)

 

 

f

(

n

)

=

n

d

F

(

d

)

μ

(

d

n

)

F(n)=\sum\limits_{n|d}f(d)\ \Rightarrow\ f(n)=\sum\limits_{n|d}F(d)\mu(\frac{d}{n})

F(n)=ndf(d)  f(n)=ndF(d)μ(nd)(说明一下,我并没有在网上找到太多关于这个公式的描述,特别是有关逆定理的证明并没有找到,而且我也不会证 ,所以这里并没有加上反方向的箭头)

这个式子难以用卷积来直接证明,这里直接用了求和的恒等变换的方法,可供参考:

k

=

d

n

k=\frac{d}{n}

k=nd,则要证的等式右侧为:

k

N

F

(

k

n

)

μ

(

k

)

=

k

N

μ

(

k

)

k

n

c

f

(

c

)

\sum\limits_{k \in N^*} F(kn)\mu(k)=\sum\limits_{k \in N^*}\mu(k)\sum\limits_{kn|c}f(c)

kNF(kn)μ(k)=kNμ(k)kncf(c)这里可以用类似于狄利克雷卷积结合律的性质交换求和号位置(有点玄学,需要理解一下):

k

N

μ

(

k

)

k

n

c

f

(

c

)

=

n

c

f

(

c

)

k

c

n

μ

(

k

)

\sum\limits_{k \in N^*}\mu(k)\sum\limits_{kn|c}f(c)=\sum\limits_{n|c}f(c)\sum\limits_{k|\frac{c}{n}}\mu(k)

kNμ(k)kncf(c)=ncf(c)kncμ(k)此时后面一个是我们熟悉的卷积式,即

k

c

n

μ

(

k

)

=

ε

(

c

n

)

\sum\limits_{k|\frac{c}{n}}\mu(k)=\varepsilon(\frac{c}{n})

kncμ(k)=ε(nc),当且仅当

c

=

n

c=n

c=n时计入答案:

n

c

f

(

c

)

k

c

n

μ

(

k

)

=

n

c

f

(

c

)

ε

(

c

n

)

=

f

(

n

)

\sum\limits_{n|c}f(c)\sum\limits_{k|\frac{c}{n}}\mu(k)=\sum\limits_{n|c}f(c)\varepsilon(\frac{c}{n})=f(n)

ncf(c)kncμ(k)=ncf(c)ε(nc)=f(n)至此,我们已经证明了上述定理2。

以上就是莫比乌斯反演的理论部分了…这还有很多应用,包括与狄利克雷卷积的组合应用,各种求和问题以及杜教筛都是基于这一板块的内容的。


版权声明:本文为ix35_原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。