莫比乌斯反演+常见数论函数的性质+狄利克雷卷积+数论分块+杜教筛学习笔记
一些约定
- 本文中所有未知数如没有特别说明，均为整数
- \(\gcd(a,b)\)表示a,b的最大公约数
- \(a|b\)表示a能整除b
- \([表达式]\)表示表达式成立的时候为1，不成立的时候为0.如\([n=1]\)在\(n=1\)的时候为1，否则为0
- \(\omega(n)\)表示n本质不同的质因子个数
- 在没有说明的情况下,所有除法都看做下取整
预备知识
积性函数
首先，我们定义:定义域和值域都在整数集合中的函数称为数论函数(这个定义不严谨，但不影响等下的推倒)
定义:
若数论函数\(f\)是积性函数，那么对于\(\forall \gcd(a,b)=1,f(ab)=f(a)f(b)\)
若数论函数\(f\)是完全积性函数，那么对于\(\forall a,b \in N,f(ab)=f(a)f(b)\)
常见的积性函数
1.莫比乌斯函数
\[\mu(n)=\begin{cases} 1,n=1 \\ (-1)^k,n=\prod_{i=1}^k p_k^{a_k},a_k=1 (即n的所有k个质因子指数均为1的时候) \\ 0，otherwise\end{cases}\]
其中第二种情况严谨的说，应该写成\((-1)^{\omega(n)}\)
证明：
我们分类讨论\(\mu\)的符号
(1)当\(\mu(a)=1\)时，\(a=1\),\(\mu(ab)=\mu(b)=\mu(a)\mu(b)\)
(2)当\(\mu(a)=0\)时，a有指数>1的质因子，乘上b后，ab一定有指数>1的质因子,因此\(\mu(ab)=0=\mu(a)\mu(b)\),
(3)当\(\mu(a)=-1\)时，\(\gcd(a,b)=1\),则说明a,b没有公共质因子
(3.1)当\(\mu(b)=1\)时，同(1)得证
(3.2)当\(\mu(b)=0\)时，同(3)得证
(3.3)当\(\mu(b)=-1\)时，a,b的质因子指数都为1，因为a,b没有公共质因子，因此ab的质因子指数也都为1,则\(\mu(ab)=(-1)^{\omega(a)+\omega(b)}=(-1)^{\omega(a)}(-1)^{\omega(b)}=\mu(a)\mu(b)\)
因此，\(\mu\)是积性函数
性质1：若p是质数，\(\mu(p)=-1\)
证明：满足第二种情况，\(\mu(p)=(-1)^1=-1\)
性质2：若p是质数,\(i \in N\)且\(i\)不整除p，\(\mu(pi)=-\mu(i)\)
证明：
(1)当\(\mu(i)=0或1\)的时候，显然得证
(2)当\(\mu(i)=-1\)时i的质因数分解里显然没有p，那么pi的质因子个数比i多1个，会使得\((-1)^k\)的符号改变
根据性质1,2,我们可以类似线性筛素数的方法，写出线性筛莫比乌斯函数的代码
void sieve(int n){
mu[1]=1;
for(int i=2;i<=n;i++){
if(!vis[i]){
prime[++cnt]=i;
mu[i]=-1;//性质1
}
for(int j=1;j<=cnt&&i*prime[j]<=n;j++){
vis[i*prime[j]]=1;
if(i%prime[j]==0){//prime[j]|i,i*prime[j]中，质因子prime[j]的指数>1,mu值为0
mu[i*prime[j]]=0;
break;
}else mu[i*prime[j]]=-mu[i];//性质2
}
}
}
2.欧拉函数
\[\varphi(n)=\sum_{i=1}^n [gcd(i,n)=1]\]
表示1~n中与n互质的数
性质1：若p是质数，则\(\varphi(p)=p-1\)
证明：显然1~p-1的数都与p互质
性质2：若m是正整数，则\(a^{\varphi(m)} \equiv 1 (\mathrm{mod}\ m)\)
证明：这就是著名的欧拉定理，它的证明需要用到剩余系相关知识，不属于(wo)本文(ye)的讨论(bu)范围(hui)，请读者自行使用搜索引擎
性质3：若n的质因数分解式为\(\prod_{i=1}^k p_i^{a_i}\),则\(\varphi(n)=n\prod_{i=1}^k(1-\frac{1}{p_i})\)
证明：
设p,q均为n的质因子，1~n中p的倍数有\(\frac{n}{p}\)个，q的倍数有\(\frac{n}{q}\)个，既是p的倍数也是q的倍数的数有\(\frac{n}{pq}\)个。根据容斥原理，既不是p的倍数也不是q的倍数的数有\(n-\frac{n}{p}-\frac{n}{q}+\frac{n}{pq}=n(1-\frac{1}{p})(1-\frac{1}{q})\)个
两个数的结论显然可以推广到k个数，得证
由性质3，我们显然可以证明\(\varphi\)是积性函数
3.约数个数函数\(d(n)\)
若n的质因数分解式为\(\prod_{i=1}^k p_i^{a_i}\),则\(d(n)=\prod_{i=1}^k(a_i+1)\)
4.约数和函数\(\sigma\)
若n的质因数分解式为\(\prod_{i=1}^k p_i^{a_i}\),则\(\sigma(n)=\prod_{i=1}^k \sum_{j=0}^{a_i}p_i^j)\)
常见的完全积性函数
1.常函数\(I(n)=1\)
2.恒等函数\(id(n)=n\)
3.单位函数\(\varepsilon(n)=[n=1]\)
狄利克雷卷积
定义：如果函数\(F,f,g\)满足\[F(n)=\sum_{d|n} f(d)g(\frac{n}{d})\]，则F是f和g的狄利克雷卷积，记作\(F=(f*g)\),或\(F(n)=(f*g)(n)\)
由定义知，狄利克雷卷积满足交换律和结合律
性质1：\(f*\varepsilon=f\)(任意函数卷积上单位函数都等于本身)
证明：
\((f*\varepsilon)(n)=\sum_{d|n}f(d)\varepsilon(\frac{n}{d})\)
显然只有\(\frac{n}{d}=1\)的时候乘积才不为0，\((f*\varepsilon)(n)=f(n)\varepsilon(1)=f(n)\)
性质2:\(\mu *I=\varepsilon\)
证明：
当\(n=1\)时，\(\mu(n)=\varepsilon(n)=1\)
当\(n \neq 1\)时，我们要求证\(\varepsilon(n)=\sum_{d|n} \mu(d)\)
设n的质因数分解式为\(\prod_{i=1}^k p_i^{a_i}\),则d表示在这里面选择一些质因子乘起来的情况。只有选的质因子指数为1的时候才会对答案有1或-1的贡献，取决于。那么问题就转换成在k个数中选0,1,2,3...个数的方案数。
\[\sum_{d|n} \mu(d)=\sum_{i=0}^k C_k^i (-1)^i=\sum_{i=0}^k C_k^i (-1)^i 1^{k-i}=(-1+1)^k=0\]
倒数第二步用到了二项式定理.此时\(\varepsilon(n)=0\)，得证
性质3：\(\varphi * I=id\)
证明:
\[id(n)=n=\sum_{i|n}\sum_{j=1}^n[gcd(j,n)=i]\](把1~n的数按gcd分类)
\[=\sum_{i|n}\sum_{j=1}^{\frac{n}{i}}[gcd(j,\frac{n}{i})=1]\](\(gcd(ka,kb)=kgcd(a,b)\))
\[=\sum_{i|n} \varphi(\frac{n}{i})\](欧拉函数的定义)
\[=\sum_{i|n} \varphi(i)\]
\[=(\varphi*I)(n)\]
性质4：\(\mu*id=\varphi\)
证明：
由性质2知\(\mu * I=\varepsilon\)
由性质3知\(id=\varphi *I\)
那么\(id * \varepsilon=id=\varphi * I\)
代入\(\mu * I=\varepsilon\)，\(id * \mu * I=\varphi *I\)
约掉\(I\)，得\(\mu*id=\varphi\)
性质5：\(d=I*I\)
证明:\(d(n)=\sum_{d|n}1 \times 1\),n的约数个数就是整除n的数的个数，显然得证
性质6: \(\sigma=I*id\)
证明:\(\sigma(n)=\sum_{d|n}d=\sum_{d|n}1 \times id(d)\),显然得证
我们再来总结一下比较重要的几条性质，因为在莫比乌斯反演里会用到:
\(\mu *I=\varepsilon\)
\(\varphi * I=id\)
\(f * \varepsilon =f\)
莫比乌斯反演
莫比乌斯反演是一种简化计算的方法
比如，我们可以用函数f表示F，比如\(F(n)=\sum_{d|n}f(d)\)那么能不能用F来表示f呢？
莫比乌斯反演的约数形式
定理:若\(F(n)=\sum_{d|n}f(d)\)，则\(f(n)=\sum_{d|n}F(d) \mu(\frac{n}{d})\)
这就是莫比乌斯反演的约数形式(第一种形式)，观察式子得到，我们用\(F,\mu\)来表示了\(f\)
如果我们要求f，但f比较难求解。而f求和可以得到一个函数F，F很好求，那么就可以用莫比乌斯反演求解出f.因此用莫比乌斯反演做题的时候，重点是找出好求解的函数F，来优化f的求解过程．
关于这个定理的证明，可以不依靠狄利克雷卷积直接证明，但式子较繁琐。这里用卷积证明：
把已知条件写成卷积形式，\(F=f*I\)
两边卷积上\(\mu\),得\(F*\mu=f*I*\mu\)
由\(\mu *I=\varepsilon\)，得\(F*\mu=f *\varepsilon=f\)
再写成求和形式，\(f(n)=\sum_{d|n}F(d) \mu(\frac{n}{d})\)
莫比乌斯反演的倍数形式
定理，若\(F(n)=\sum_{n|d}f(d)\),则\(f(n)=\sum_{d|n}F(d) \mu(\frac{d}{n})\)
这就是莫比乌斯反演的倍数形式(第二种形式)。观察式子得到，这是一个无限求和，但往往题目中其实只需要求到一个上界。
这种形式的证明见:https://oi-wiki.org/math/mobius/
为了加深理解，可以看一道例题:[BZOJ 2301] [HAOI 2011] Problem b (莫比乌斯反演)(有证明)
数论分块
数论分块其实和数据结构中的分块关系不大
在莫比乌斯反演中，我们经常需要求解这种形式的式子\(\sum_{i=1}^n f(i) \lfloor \frac{n}{i} \rfloor\)
暴力枚举求解i的复杂度仍然较高
观察到d比较大的时候n/d的值变化不大，如n=6时n/4=n/5=n/6.
因此我们可以找出所有下取整值相同的区间[l,r],再用\(\sum \lfloor\frac{n}{l} \rfloor \sum_{i=l}^r f(i)\)就可以求出答案，求和部分可以前缀和预处理。对于一个块，假设它的起始位置的下标为l，那么可以得到的是，它的结束位置的下标为n/(n/l)
如果我们预处理出f的前缀和，那么数论分块的复杂度是\(O(\sqrt{n})\)的。
下面的代码求解了\(\sum_{d=1}^n \mu(d) \lfloor \frac{n}{d} \rfloor \lfloor \frac{m}{d} \rfloor\)
long long ans=0;
int l,r;
for(l=1;l<=m;l=r+1){
r=min(n/(n/l),m/(m/l));
ans+=(sum_mu[r]-sum_mu[l-1])*(n/l)*(m/l);//sum_mu为mu的前缀和
}习题
[BZOJ 2301] [HAOI 2011] Problem b (莫比乌斯反演)(有证明)