莫比乌斯反演,狄利克雷卷积,杜教筛

本文深入探讨了数论中的莫比乌斯反演及其应用,通过实例解释了莫比乌斯函数的性质,并证明了莫比乌斯反演的两个重要式子。接着介绍了狄利克雷卷积,包括其定义、性质以及与莫比乌斯反演的关系。最后,简要提及了杜教筛在高效计算数列求和中的作用。
摘要由CSDN通过智能技术生成

Part 1 莫比乌斯反演

首先定义莫比乌斯函数\(\mu(n)\)
\(n=\prod_{i=1}^mp_i^{k_i}\),其中\(p_{1-n}\)为互不相等的质数,则

\[\mu(n)=\begin{cases}0 & \exists k_i>1(i\in[1,m]) \\(-1)^m & else \end{cases} \]

通俗地讲就是:如果\(n\)能被一个质数的平方整除,则\(\mu(n)=0\),否则若\(n\)\(m\)个质因数,则\(\mu(n)=(-1)^m\)

由定义可以得到几个结论:
\(1.\mu(1)=(-1)^0=1\)
\(2.\mu(p)=-1\)\(p\)为质数)
\(3.\sum_{d|n}\mu(d)=[n=1]\)。其中\([]\)是条件运算符,当里面的条件成立时值为\(1\),否则为\(0\)
下面来证明第三个结论。
首先,若\(n=1\),结论显然成立。
否则,我们同样把\(n\)写成\(\prod_{i=1}^mp_i^{k_i}\)的形式,那么结论的式子可以展开为

\[\sum_{d_1=0}^{k_1}\sum_{d_2=0}^{k_2}···\sum_{d_m=0}^{k_m} \mu(\prod_{i=1}^mp_i^{d_i}) \]

可以发现当存在一个\(d_i>1\)时,后面的\(\mu\)的值为0。所以我们只需要考虑\(d_i\)\(0\)\(1\)的情况。这时可以把我们式子看成:在一个大小为\(m\)的集合中选任意个数,这些数的乘积的\(\mu\)值之和。又由于这些数都是互不相等的质数,所以乘积的\(\mu\)值只与选了多少数有关。那么枚举选了多少个数,式子就可以写成:

\[\sum_{i=0}^m{m\choose i}(-1)^i \\=\sum_{i=0}^m{m\choose i}(-1)^i\times1^{m-i} \\=0^m \]

所以只有当\(m=0\),即\(n=1\)时,\(\sum_{d|n}\mu(d)=1\),其余都为0。

接下来就是莫比乌斯函数的一种应用,被称为莫比乌斯反演。
通常有两种式子:

\[若F(n)=\sum_{d|n}f(d) \\则f(n)=\sum_{d|n}\mu(\dfrac{n}{d})F(d) \]

证明很简单,把\(F\)展开即可:

\[f(n)=\sum_{d|n}\mu(\dfrac{n}{d})\sum_{r|d}f(r) \\=\sum_{r|n}f(r)\sum_{d|\frac{n}{r}}\mu(\dfrac{n}{dr}) \\=\sum_{r|n}f(r)\sum_{d|\frac{n}{r}}\mu(d) \]

根据上面的性质可以知道,只有当\(r=n\)时,第二个\(\sum\)才有值,为1,此时等号右边等于\(f(n)\)
第二种式子为:

\[若F(n)=\sum_{n|d}f(d) \\则f(n)=\sum_{n|d}\mu(\dfrac{d}{n})F(d) \]

证明和上面类似:

\[f(n)=\sum_{n|d}\mu(\dfrac{d}{n})\sum_{d|r}f(r) \\=\sum_{n|r}f(r)\sum_{d|\frac{r}{n}}\mu(d) \]

同理,只有当\(r=n\)时,第二个\(\sum\)有值,为1。
下面总结的狄利克雷卷积有更好的方式证明第一种式子。

Part 2 狄利克雷卷积

首先总结一下积性函数和几个常见的积性函数。
积性函数的定义为:若\(n,m\)互质,且\(f(n)\times f(m)=f(nm)\),则\(f\)为积性函数。
另外有个完全积性函数的定义:若\(f(n)\times f(m)=f(nm)\),则\(f\)为完全积性函数。
下面总结一下常见的积性函数:
\(1.e(n)=[n=1]\)元函数。
\(2.I(n)=1\)常值函数。
\(3.id(n)=n\)不知道叫什么函数。
\(4.p_k(n)=n^k\)幂函数
\(4.\mu(n)\)莫比乌斯函数,定义在上面。
\(5.\varphi(n)=\sum_{i=1}^n[gcd(i,n)=1]\),即欧拉函数,意义是\(n\)以内与\(n\)互质的数的个数。
\(6.d(n)=\sum_{d|n}1\)约数个数函数,值为\(n\)的约数个数。
可以注意到前四个实际上还是完全积性函数。
积性函数可以通过欧拉筛(线筛)来求。

狄利克雷卷积是一种运算,符号为\(*\),定义两个函数的狄利克雷卷积\((f*g)(n)=\sum_{d|n}f(d)g(\dfrac{n}{d})\)
狄利克雷卷积有几个常用的式子:
对于任意函数\(f\),都有\((f*e)(n)=f(n)\)。证明就是把式子展开。
\((\mu*I)(n)=e\)。展开之后就是莫比乌斯函数的结论。
\((id*id)(n)=n\times d(n)\)还是展开再证明
\((\varphi*I)(n)=id(n)\)
最后一个式子可以展开为\(\sum_{d|n}\varphi(d)=n\)
这个式子有一个巧妙的证明和一个中规中矩的证明。

巧妙的证明:
尝试列出\(\dfrac{1}{n},\dfrac{2}{n},···,\dfrac{n}{n}\),并把他们化为最简分数。设化简后的分母集合为\(S\)
引理\(1\)\(S\)恰好为\(n\)的约数集。
证明:若存在一个\(n\)的约数\(d\)不在\(S\)中,则\(\dfrac{n/d}{n}\)不在最初的数列中。不成立。
\(x\)不是\(n\)的约数。由于最初的数列中没有一个数能化简得到分母为\(x\)的分数,所以\(x\)不可能在\(S\)中。得证。
引理\(2\):对于\(S\)中的任意元素\(x\),化简后的数列中所有以\(x\)为分母的分数的分子都与\(x\)互质,且分子集合恰好是所有\(x\)以内与\(x\)互质的数。
证明:分子分母显然,因为是最简分数。若存在一个数\(y\leq x\)且与\(x\)互质,且不在分子集中出现过,则分数\(\dfrac{y\times\frac{n}{x}}{n}\)不在最初的数列中出现,不成立。得证。
因此,化简后以约数\(d\)为分母的分数一共有\(\varphi(d)\)个。而分数一共有\(n\)个,因此\(\sum_{d|n}\varphi(d)=n\)

中规中矩的证明:
引理:\(f(n)=\sum_{d|n}\varphi(d)\)为积性函数。
证明:

\[f(n)\times f(m)=(\sum_{d|n}\varphi(d))(\sum_{r|m}\varphi(r)) \\=\sum_{d|n}\varphi(d)\sum_{r|m}\varphi(r) \\=\sum_{d|n}\sum_{r|m}\varphi(d)\varphi(r) \\=\sum_{dr|nm}\varphi(dr) \\=f(nm) \]

这里涉及到\(\varphi\)为积性函数。不懂就搜。

\[f(n)=\prod_{i=1}^mf(p_i^{k_i}) \\=\prod_{i=1}^m(1+\sum_{j=1}^{k_i}(p_i^j-p_i^{j-1})) \\=\prod_{i=1}^mp_i^{k_i}=n \]

狄利克雷卷积也有运算律:
·交换律:\(f*g=g*f\)。根据狄利克雷卷积定义式可知其显然成立。
·结合律:\((f*g)*h=f*(g*h)\)
证明:

\[\sum_{d|n}(\sum_{r|d}f(r)g(\frac{d}{r}))h(\frac{n}{d}) \\=\sum_{d|n}\sum_{r|d}f(r)g(\frac{d}{r})h(\frac{n}{d}) \\=\sum_{r|n}f(r)\sum_{d|\frac{n}{r}}g(d)h(\frac{n}{dr}) \]

·分配律:\((f+g)*h=f*h+g*h\)。把左边展开即得证。

上面说狄利克雷卷积可以证明莫比乌斯反演的第一个式子,现在填坑:
\(F(n)=\sum_{d|n}f(d)\),则\(f(n)=\sum_{d|n}\mu(d)F(\dfrac{n}{d})\)
证明:\(F(n)=\sum_{d|n}f(d)\)可以写成:\(F=f*I\)
两边同时卷一个\(\mu\)可以得到\(F*\mu=f*I*\mu=f*(I*\mu)=f*e=f(n)\)。即得证。

同时还有一个常用的式子:\(\sum_{d|n}d\times\mu(\dfrac{n}{d})=\varphi(n)\)
证明:把左边写成\(id*\mu\)。由于\(id=\varphi*I\),所以\(id*\mu=\varphi*I*\mu=\varphi*e=\varphi(n)\)

Part 3 杜教筛

可以在\(O(n^{\frac{2}{3}})\)的复杂度内求出\(\sum_{i=1}^nf(i)\)
设要求的为\(S(n)\),并设\(h=f*g\)\(g\)是我们选的一个与\(f\)有联系的函数,最后会讲怎么选。则

\[\sum_{i=1}^nh(i)=\sum_{i=1}^n\sum_{d|i}g(d)f(\frac{i}{d}) \\=\sum_{d=1}^ng(d)\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}f(i) \\=\sum_{d=1}^ng(d)S(\lfloor\frac{n}{d}\rfloor) \]

然后把\(d=1\)的一项留着,其余的减过去,得到

\[g(1)S(n)=\sum_{i=1}^nh(i)-\sum_{d=2}^ng(d)S(\lfloor\frac{n}{d}\rfloor) \]

这个就可以整除分块然后递归了。这就是杜教筛的套路式。
观察这个式子,可以发现杜教筛的使用条件:
\(h\)\(g\)的前缀和必须很好求。比如\(f(i)=i\)这样的函数。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值