Learning：数论（四）莫比乌斯反演（+题集）

莫比乌斯函数的定义：

$\mu (x)=\{\begin{array}{l}1,n=1\\ (-1{)}^k,n=p_1{p}_2...p_k，其中p_1,p_2,...,p_k为互不相同的素数\\ 0,其他情况\end{array}$
自己稍微推一推就可以知道这是一个积性函数。
所以很显然我们可以线性筛$\mu$
由定义得：
$\mu (1)=1$
$\mu (p)=-1$
$\mu (p^k)=0(k>1)$
所以就可以直接筛了，代码如下：

void init(){
    mu[1]=1;
    for(int i=2;i<N;i++){
        if(!vis[i]){
            prm[cnt++]=i;
            mu[i]=-1;
        }
        for(int j=1;j<cnt&&i*prm[j]<N;j++){
            vis[i*prm[j]]=1;
            if(!(i%prm[j])){
                mu[i*prm[j]]=0;
                break;
            }
            mu[i*prm[j]]=-mu[i];
        }
    }
    return;
}

莫比乌斯反演

莫比乌斯反演有一个非常非常非常重要的性质：
${\sum }_{d|n}\mu (d)=[n=1]$
怎么证？
我们设$n={\prod }_{i=1}^k{p}_i^{{\alpha }_i}$
考虑$n$的所有因数，他们的素因子都是$p_0,p_1,...,p_{k-1}$中的，而我们会发现只有当因数的所有的素因子次数为$1$时对答案的贡献才不为$0$，所以我们可以认为$n$的每个对答案的贡献不为$0$因数都是在$p_0,p_1,...,p_{k-1}$中取若干个素数相乘得到的。特殊的，因数1我们可以当做一个素数都不取。
当$n=1$时，显然$原式=1$。
当$n>1$时，
若$k$为奇数，
首先易证$C_n^m=C_n^{n-m}$
选奇数个素数相乘得到的因数个数为$C_k^1+C_k^3+...+C_k^k$，这种因数的函数值为$-1$
选偶数个素数相乘得到的因数个数为$C_k^0+C_k^2+...+C_k^{k-1}$，这种因数的函数值为$1$
首先显然$C_n^m=C_n^{n-m}$
所以易证$C_k^1+C_k^3+...+C_k^k=C_k^0+C_k^2+...+C_k^{k-1}$，所以因数的函数值加起来为$0$。
若$k$为偶数，
选奇数个素数相乘得到的因数个数为$C_k^1+C_k^3+...+C_k^{k-1}$，这种因数的函数值为$-1$
选偶数个素数相乘得到的因数个数为$C_k^0+C_k^2+...+C_k^k$，这种因数的函数值为$1$
首先显然$C_n^m=C_{n-1}^{m-1}+C_{n-1}^m(m̸=0),C_n^0=C_{n-1}^0$
所以$C_k^0+C_k^2+...+C_k^k=C_{k-1}^0+C_{k-1}^1+C_{k-1}^2+...+C_{k-1}^{k-2}+C_{k-1}^{k-1}=2^{k-1}$
$C_k^1+C_k^3+...+C_k^{k-1}=C_k^1+C_k^2+C_k^3+...+C_k^k-(C_k^0+C_k^2+...+C_k^k)$
$=2^k-(C_k^0+C_k^2+...C_k^k)=2^k-2^{k-1}=2^{k-1}=C_k^0+C_k^2+...C_k^k$
所以因数的函数值加起来为$0$。
综上，${\sum }_{d|n}\mu (d)=[n=1]$

这个性质有多重要，我们接下来看一道题（前置知识：数论分块）：
有$t(1\le t\le 10000)$组数据，每组数据给定$n,m,(1\le n\le m\le 10^7)$求$$\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=1]$$
O(nm)秒了
我们会发现暴力的时间复杂度是$O(nm)$，显然会TLE，所以我们考虑优化。
$[gcd(i,j)=1]$这个式子很神奇？所以就从这里入手。
仔细分析就会发现$$[gcd(i,j)=1]=\sum _{d|gcd(i,j)}\mu (d)$$，因为由前面的莫比乌斯函数的性质可得当且仅当$gcd(i,j)=1$时，${\sum }_{d|gcd(i,j)}\mu (d)$才会等于$1$，否则等于$0$
所以原式就可以变成：$$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|gcd(i,j)}\mu (d)$$
然后枚举$d$，此时$i,j$就是$d$的倍数，所以小于$n$的$i$就有$\lfloor \frac{n}{d}\rfloor$个，$j$同理。所以我们就可以轻松得到：$$\sum _{d=1}^n\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor \mu (d)$$
线性筛出$\mu$的前缀和，然后数论分块$O(\sqrt{n}+\sqrt{m})$求解。
这就是莫比乌斯反演（其实有个公式的，但是用这个公式大部分题目都很难推，所以不如直接用莫比乌斯函数的性质来搞）
莫比乌斯反演只有多推才能熟练，所以
我们继续！

有$t(1\le t\le 10000)$组数据，每组数据给定$n,m,k(1\le n\le m\le 10^7)$.求:$$\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=k]$$（HDU1695数据范围稍有增加）
咋一看和刚刚的怎么长得那么像？
我们会发现只有当$i,j$都是$k$的倍数时，且$gcd(\frac{i}{k},\frac{j}{k})=1$时才能成立
所以原式转化为$$\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{i=1}^{\lfloor \frac{m}{k}\rfloor }[gcd(i,j)=1]$$
这不就转化为原来的那道题了吗？得到原式等于$$\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\lfloor \frac{n}{dk}\rfloor \lfloor \frac{m}{dk}\rfloor \mu (d)$$
直接数论分块

接着下一题：
有$t(1\le t\le 10000)$组数据，每组数据给定$n,m(1\le n\le m\le 10^7)$求$$\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)为素数]$$
（bzoj2820）
咋一看怎么和刚刚的还是那么像？
记素数集为$prime$
所以我们直接枚举素数，得到$$\sum _{p\in prime}\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=p]$$
由上题可得：原式=$$\sum _{p\in prime}\sum _{d=1}^{\frac{n}{p}}\lfloor \frac{n}{pd}\rfloor \lfloor \frac{m}{pd}\rfloor \mu (d)$$
然后我们记$k=pd$，枚举$k$，得到：
$$\sum _{k=1}^n\lfloor \frac{n}{k}\rfloor \lfloor \frac{m}{k}\rfloor \sum _{p\in prime,p|k}\mu (\frac{k}{p})$$
设函数$f(n)={\sum }_{p\in prime,p|n}\mu (\frac{n}{p})$
所以原式=$$\sum _{k=1}^n\lfloor \frac{n}{k}\rfloor \lfloor \frac{m}{k}\rfloor f(k)$$
我们可以在线性筛的时候处理$f(n)$。
当$n$为素数时，显然$f(n)=1$
当$n$为合数时，
记$prm[j]$为$n$的最小素因子（线性筛中当前枚举的素数），$i=\frac{n}{p}$
当$prm[j]|i$（$n$的最小素因子$prm[j]$的次数大于$1$）时，
当$i$无多个次数大于1的素因子时，当且仅当$p=prm[j]$时$\mu (\frac{n}{p})̸=0$，此时$f(n)=\mu (\frac{n}{prm[j]})=\mu (i)$
当$i$有多个次数大于1的素因子时，无论$prm[j]$等于多少，$\mu (\frac{n}{p})=0$，又因为易知$\mu (i)=0$，所以$f(n)=\mu (\frac{n}{p})=0=\mu (i)$，
所以当$prm[j]|i$时$f(n)=\mu (i)$
当$prm[j]̸|i$（$n$的最小素因子$prm[j]$的次数等于$1$）时，
此时，$n$就比$i$多了一个素因数$prm[j]$，因为
$$f(i)=\sum _{p|i,p\in prime}\mu (\frac{i}{p})$$
$$f(n)=\sum _{p|n,p\in prime}\mu (\frac{n}{p})=\sum _{p|n,p\in prime}\mu (\frac{i\ast prm[j]}{p})$$
因为$i$不含素因子$prm[i]$，所以
$$\mu (\frac{i\ast prm[j]}{p})=-\mu (\frac{i}{p})$$
又因为$n$比$i$多一个素因子$prm[j]$，所以
$$f(n)=-\sum _{p|n,p\in prime,p̸=prm[j]}\mu (\frac{i}{p})+\mu (\frac{n}{prm[j]})$$
$$=\mu (i)-\sum _{p|i,p\in prime}\mu (\frac{i}{p})$$
$$=\mu (i)-f(i)$$
这样我们就可以在线性筛的时候预处理$f$，然后直接数论分块解决。

继续
有$t(1\le t\le 1000000)$组数据，每组数据给定$n,m,k(1\le n\le 10^7)$求$$\sum _{i=1}^{n}gcd(i,n)$$
前置知识：狄利克雷卷积
首先枚举$gcd(i,n)$：
$$\sum _{d|n}d\sum _{i=1}^n[gcd(i,n)=d]$$
除掉$d$得到：
$$\sum _{d|n}d\sum _{i=1}^{\frac{n}{d}}[gcd(i,\frac{n}{d})=1]$$
反演：
$$\sum _{d|n}d\sum _{i=1}^{\frac{n}{d}}\sum _{k|i,k|\frac{n}{d}}\mu (k)$$
枚举$k$：
$$\sum _{d|n}\sum _{k|\frac{n}{d}}d\mu (k)\frac{n}{dk}$$
我们记$T=dk$，然后枚举T：
$$\sum _{T|n}\frac{n}{T}\sum _{d|T}d\mu (\frac{T}{d})$$
记$g(n)=n,f(n)={\sum }_{d|n}d\mu (\frac{n}{d})$
显然$g$是积性函数
然后我们发现$f$是$g$和$\mu$的狄利克雷卷积，所以也是积性函数
然后我们发现原式其实也就是$f$和$g$的狄利克雷卷积，所以也是一个积性函数
直接线性筛，然后$O(1)$询问。

其实还有一种更简单的做法，需要用到欧拉函数的性质${\sum }_{d|n}\varphi (d)=n$
利用这个性质我们就可以吧原式转化为：
$$\sum _{i=1}^n\sum _{d|i,d|n}\varphi (d)$$
直接枚举$d$，就可以直接得到
$$\sum _{d|n}\varphi (d)\frac{n}{d}$$
这显然是狄利克雷卷积的形式，是一个积性函数，直接线性筛就OK了。

再看一道题
有$t(1\le t\le 10000)$组数据，每组数据给定$n,m,k(1\le n\le m\le 10^7)$.求:$$\sum _{i=1}^n\sum _{j=1}^{m}gcd(i,j)$$
这道题和刚刚有点像？
还是先枚举$gcd(i,j)$：
$$\sum _{d=1}^{n}d\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=d]$$
通过之前第二题的结论可以得到：
$$\sum _{d=1}^{n}d\sum _{k=1}^n\lfloor \frac{n}{dk}\rfloor \lfloor \frac{m}{dk}\rfloor \mu (k)$$
我们令$T=dk$，然后直接枚举$T$：
$$\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor \sum _{d|T}d\mu (\frac{T}{d})$$
记$f(n)={\sum }_{d|T}d\mu (\frac{T}{d})$，显然是狄利克雷卷积的形式，是一个积性函数，可以直接线性筛。
所以原式可以转化为：
$$\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor f(T)$$
直接数论分块$O(\sqrt{n}+\sqrt{m})$
这道题同样也可以用欧拉函数，而且同样简单很多
原式可以直接转化为
$$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|i,d|j}\varphi (d)$$
直接枚举$d$：
$$\sum _{d=1}^n\varphi (d)\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$$
然后数论分块就OK了。
就是这么简单。
那干嘛要讲莫比乌斯反演的做法？
练手感…因为很多题都是只能用莫比乌斯反演做的，比如说下一题：
有$t(1\le t\le 1000000)$组数据，每组数据给定$n,m,k(1\le n\le 10^7)$求$$\sum _{i=1}^{n}lcm(i,n)$$
先讲一种用欧拉函数的$O(nlogn)$做法（当然在这道题会TLE，但是可以了解一下，这种做法非常神奇）
$$\sum _{i=1}^{n}lcm(i,n)$$
$$=n\sum _{i=1}^n\frac{i}{gcd(i,n)}$$
$$=n\sum _{d|n}\frac{1}{d}\sum _{i=1}^{n}i[gcd(i,n)=d]$$
$$=n\sum _{d|n}\frac{1}{d}\sum _{i=1}^{\frac{n}{d}}id[gcd(i,\frac{n}{d})=1]$$
$$=n\sum _{d|n}\sum _{i=1}^{\frac{n}{d}}i[gcd(i,\frac{n}{d})=1]$$
$$=n\sum _{d|n}\frac{d\varphi (d)+[d=1]}{2}$$
前面几行应该都懂了吧，都做了那么多题了
问题是最后一行…
我们只考虑后面的式子：${\sum }_{i=1}^{\frac{n}{d}}i[gcd(i,\frac{n}{d})=1]$
因为$d$是$n$的因数，所以我们可以用$d$代替$\frac{n}{d}$：
$$\sum _{i=1}^{d}i[gcd(i,d)=1]$$
其实就是求小于等于$d$与$d$互质的数的和
然后我们会发现如果$x\perp d(x<d)$，那么$(d-x)\perp d$
所以我们可以把与$d$互质的数首尾两两分组，他们的和就是$d$，然后他们所有的和就是$\frac{d\varphi (d)}{2}$，如果$x=d-x$也不影响，特判一下d=1的情况就OK了。
然后线性筛出$\varphi$，调和级数时间预处理答案。
但是这道题要$O(n)$才行啊，有更优秀的做法吗？
莫比乌斯反演
$$\sum _{i=1}^{n}lcm(i,n)$$
$$=n\sum _{i=1}^n\frac{i}{gcd(i,n)}$$
$$=n\sum _{d|n}\frac{1}{d}\sum _{i=1}^{n}i[gcd(i,n)=d]$$
$$=n\sum _{d|n}\frac{1}{d}\sum _{i=1}^{\frac{n}{d}}id[gcd(i,\frac{n}{d})=1]$$
$$=n\sum _{d|n}\sum _{i=1}^{\frac{n}{d}}i[gcd(i,\frac{n}{d})=1]$$
$$=n\sum _{d|n}\sum _{i=1}^{\frac{n}{d}}i\sum _{k|i,k|\frac{n}{d}}\mu (k)$$
$$=n\sum _{d|n}\sum _{k|\frac{n}{d}}\mu (k)k\sum _{i=1}^{\frac{n}{dk}}i$$
$$=n\sum _{k|n}\mu (k)k\sum _{d|\frac{n}{k}}\sum _{i=1}^{d}i$$
后面那个式子用等差数列求和：
$$=n\sum _{k|n}\mu (k)k\sum _{d|\frac{n}{k}}\frac{d(d+1)}{2}$$
$$=\frac{n}{2}\sum _{k|n}\mu (k)k\sum _{d|\frac{n}{k}}(d^2+d)$$
我们记$f(n)={\sum }_{k|n}\mu (k)k{\sum }_{d|\frac{n}{k}}{d}^2,g(n)={\sum }_{k|n}\mu (k)k{\sum }_{d|\frac{n}{k}}d$
可以轻松证明这两个都是积性函数，可以线性筛。
然后就OK了。
留一道题给读者思考：
【BZOJ2693】jzptab