莫比乌斯反演

前置技能

数论基础

约定

1.　$gcd(a,b)=(a,b)$
2.　$[a]=\{\begin{array}{l}1,a=true\\ 0,a=false\end{array}$

数论分块

已知 $n,k$，则暴力地求 $\begin{array}{r}S=\sum _{i=1}^{{}^{{}^{{}^{{}^{}}}}{n}^{{}^{{}^{{}^{}}}}}\lfloor \frac{n}{i}\rfloor \end{array}$ 需要 $O(n)$ 的时间。我们把表 $(k=10)$ 列出来看看。

$i$	$1$	$2$	$3$	$4$	$5$	$6$	$7$	$8$	$9$	$10$
$\frac{n}{i}$	$10$	$5$	$3$	$2$	$2$	$1$	$1$	$1$	$1$	$1$

其实最多只有 $2\sqrt{n}$ 种取值。
设当前枚举的下标为 $i$，当前位置的取值为 $\lfloor \frac{n}{i}\rfloor$，则相同取值的最右端点下标为 $\lfloor \frac{n}{\lfloor \frac{n}{i}\rfloor }\rfloor$。
因而可以写出如下代码：

for(int l=1,r;l<=n;l=r+1){
    r=n/(n/l);
    ans+=(r-l+1)*(n/l);
}

时间复杂度$O(\sqrt{n})$

---

变形1　已知 $n,k$，求 $\begin{array}{r}S=\sum _{i=1}^{{}^{{}^{{}^{}}}{n}^{{}^{{}^{}}}}\lfloor \frac{k}{i}\rfloor \end{array}$
解： 计算右端点时特判一下是否越界即可。

int len=min(n,k);
for(int l=1,r;l<=len;l=r+1){
    r=min(n,k/(k/l));
    if(r>n) r=n;//防止越过n
    ans+=(long long)(k/l)*(r-l+1);
}

---

变形2　已知 $n,a,b$，求 $\begin{array}{r}S=\sum _{i=1}^{{}^{{}^{{}^{}}}{n}^{{}^{{}^{}}}}\lfloor \frac{a}{i}\rfloor \lfloor \frac{b}{i}\rfloor \end{array}$
解： 分段的时候分得细一些即可，这对时间复杂度的影响是常数级别的。

int len=min(n,min(a,b));
for(int l=1,r;l<=len;l=r+1){
    r=min(a/(a/l),b/(b/l));
    if(r>n) r=n;
    ans+=(r-l+1)*(a/l)*(b/l);
}

---

变形3　已知 $n,a,b,k$，求 $\begin{array}{r}S=\sum _{i=1}^{{}^{{}^{{}^{}}}{n}^{{}^{{}^{}}}}\lfloor \frac{a}{ik}\rfloor \lfloor \frac{b}{ik}\rfloor \end{array}$
解： 首先，我们发现 $\sum$ 的上界可以调紧一些。
$$S=\sum _{i=1}^{min(n,a,b)}\lfloor \frac{a}{ik}\rfloor \lfloor \frac{b}{ik}\rfloor$$

由于 $\lfloor \frac{a}{ik}\rfloor \lfloor \frac{b}{ik}\rfloor =\lfloor \frac{\lfloor \frac{a}{k}\rfloor }{i}\rfloor \lfloor \frac{\lfloor \frac{b}{k}\rfloor }{i}\rfloor$，因此把 $k$ 提出来：
$$S=\sum _{i=1}^{min(n,\frac{a}{k},\frac{b}{k})}\lfloor \frac{\lfloor \frac{a}{k}\rfloor }{i}\rfloor \lfloor \frac{\lfloor \frac{b}{k}\rfloor }{i}\rfloor$$

这样就可以了。

int len=min(n,min(a/=k,b/=k));
for(int l=1,r;l<=len;l=r+1){
    r=min(a/(a/l),b/(b/l));
    if(r>n) r=n;
    ans+=(r-l+1)*(a/l)*(b/l);
}

变形4　已知 $n$，求 $\begin{array}{r}S=\sum _{i=1}^{{}^{{}^{{}^{}}}{n}^{{}^{{}^{}}}}i\times \lfloor \frac{n}{i}\rfloor \end{array}$
解： 乘法分配律即可。

for(int l=1,r;l<=n;l=r+1){
	r=n/(n/l);
	ans+=(l+r)*(r-l+1)*(n/l)/2;
}

---

积性函数

例子

1. 常函数 $1(n)=1$
2. 单位元 $ϵ(n)=\{\begin{array}{l}1,n=1\\ 0,n\ne 1\end{array}$
3. 约数个数函数 $\begin{array}{r}d(n)=\sum _{x|n}^{{}^{{}^{{}^{{}^{{}^{}}}}}}1\end{array}$
4. 设 $p$ 的唯一分解为 $\begin{array}{r}p=\prod p_i^{k_i}(\forall i\ne j,p_i\ne p_j)\end{array}$，则 $\begin{array}{r}d(p)=\prod (k_i+1)\end{array}$。（积性证明）
5. 约数和函数 $\begin{array}{r}\sigma (n)=\sum _{x|n}^{{}^{{}^{{}^{{}^{}}}}}x\end{array}$
6. 约数k次幂函数 $\begin{array}{r}{\sigma }_k(n)=\sum _{x|n}^{{}^{{}^{{}^{{}^{}}}}}{x}^k\end{array}$
7. 欧拉函数 $\begin{array}{r}\phi (n)=\sum _1^{{}^{{}^{}}{n}^{{}^{}}}[(i,n)=1]\end{array}$
8. 莫比乌斯函数 $\mu (x)=\{\begin{array}{ll}1,& x=1\\ (-1{)}^k,& x=p_1{p}_2...p_k,\forall p_i\ne p_j(p_i为x的唯一分解)\\ 0,& others\end{array}$
9. $id(n)=n$

性质

如果 $f$ 和 $g$ 是积性函数，则下列函数也是积性函数。
1.　$h(x)=f(x^p)$
2.　$h(x)=f^p(x)$
3.　$h(x)=f(x)g(x)$
4.　$h(x)=\sum _{d|x}f(d)g(\frac{x}{d})$

部分积性证明

欧拉函数

先说一下，网上那些根据通式证明的大多都是错的，因为通式需要积性证明。
结论　若$(n,m)=1$，则$\phi (nm)=\phi (n)\phi (m)$
证明　构造矩阵
　　　$\left[\begin{array}{ccccccc}& 1& 2& \cdots & r& \cdots & m\\ & m+1& m+2& \cdots & m+r& \cdots & 2m\\ & 2m+1& 2m+2& \cdots & 2m+r& \cdots & 3m\\ & ⋮& ⋮& \ddots & ⋮& \ddots & ⋮\\ & (n-1)m+1& (n-1)m+2& \cdots & (n-1)m+r& \cdots & nm\end{array}\right]$
　　　可以发现，$\phi (nm)$ 指的是矩阵中与 $nm$ 互质的数的个数。
　　　由于 $(n,m)=1$，因此与 $nm$ 互质的数相当于与 $n$ 互质且和 $m$ 互质的数的个数。
　　　$\because (km+r,m)=(r,m)$
　　　$\therefore$ 每一列的 $n$ 个元素同时与 $m$ 互素当且仅当 $(r,m)=1$。
　　　$\therefore$ 与 m互素的列共有 $\phi (m)$ 列
　　　设 $r$ 满足 $(r,m)=1$，则第 $r$ 列的元素为 $r,m+r,\cdots ,(n-1)m+r$
　　　用反证法可以得出：它们模 $n$ 的余数互不相同

---

证明　设存在 $i$ 和 $j$ 满足 $i\not\equiv j(\mathrm{mod}n)$ 且 $im+r\equiv jm+r(\mathrm{mod}n)$
　　　则根据同余的性质得 $im\equiv jm(\mathrm{mod}n)$
　　　又 $\because (n,m)=1$
　　　$\therefore$ 存在 $m^{-1}$，满足 $m\times m^{-1}\equiv 1(\mathrm{mod}n)$
　　　$\therefore i\equiv j(\mathrm{mod}n)$
　　　与假设矛盾。得证。

---

    \qquad\,\,\, 因此 $r,m+r,\cdots ,(n-1)m+r$ 这 $n$ 个数取遍了 $0$ 到 $n-1$，其中和 $n$ 互质的数的个数为 $\phi (n)$。综上所述，共有 $\phi (m)$ 列，每列有 $\phi (n)$ 个数，因此总共有 $\phi (n)\times \phi (m)$ 个。
　　　故$\phi (nm)=\phi (n)\phi (m)$

莫比乌斯函数

结论　若 $(n,m)=1$ 则 $\mu (nm)=\mu (n)\mu (m)$。
证明　其实根据定义就可以很简单地得到。
　　　1. 当 $n$ 存在某个大于 $1$ 的因子为完全平方数时，$\mu (n)=0$
　　　　$nm$ 的因子必然也有这个完全平方数，故 $\mu (nm)=\mu (n)\mu (m)=0$
　　　2. 当 $n$ 不存在大于 $1$ 的因子为完全平方数时
　　　　$\because (n,m)=1$，即 $n$ 和 $m$ 不含有相同的因子
　　　　$\therefore$ $nm$ 的质因子个数等于 $n$ 的质因子个数加 $m$ 的质因子个数
　　　　又 奇数$+$奇数$=$偶数，奇数$+$偶数$=$奇数，偶数$+$偶数$=$偶数
　　　　对应 $(-1)\times (-1)=1,1\times (-1)=-1,1\times 1=1$
　　　　$\therefore \mu (nm)=\mu (n)\mu (m)$，得证。

筛积性函数

欧拉函数

欧拉函数有下面的性质。
性质１　若 $m\in prime,n\equiv 0(\mathrm{mod}m)$ 则 $\phi (n\times m)=\phi (n)\times m$
证明　　构造矩阵
$$\left[\begin{array}{ccccccc}& 1& 2& \cdots & r& \cdots & n\\ & n+1& n+2& \cdots & n+r& \cdots & 2n\\ & 2n+1& 2n+2& \cdots & 2n+r& \cdots & 3n\\ & ⋮& ⋮& \ddots & ⋮& \ddots & ⋮\\ & (m-1)n+1& (m-1)n+2& \cdots & (m-1)n+r& \cdots & mn\end{array}\right]$$　　由于 $(n,m)=1$，因此与 $nm$ 互质的数相当于与 $n$ 互质且与 $m$ 互质的数的个数。又由于 $$(kn+r,n)=(r,n)$$　　因此每一列的 $m$ 个元素同时与 $n$ 互素当且仅当 $(r,n)=1$。因此共有 $\phi (n)$ 列满足，其余列均不与 $n$ 互质。   \, 其中满足 $(r,n)=1$ 的第 $r$ 列第 $k$ 行的数为 $kn+r$ 又$$n\equiv 0(\mathrm{mod}m)$$　　因此 $$kn+r\equiv 0\times k+r\equiv r(\mathrm{mod}m)$$因此 $$(kn+r,m)=(r,m)$$又 $$m|n,(r,n)=1$$因此 $(r,m)=1$ 故 $(kn+r,m)=1$
$\therefore$ 满足 $(r,n)=1$ 的第 $r$ 列的所有 $m$ 个数均与 $n$ 互质。
$\therefore \phi (n\times m)=\phi (n)\times m$

性质２　若 $p\in prime,i\not\equiv 0(\mathrm{mod}p)$ 则 $\phi (i\times p)=\phi (i)\times (p-1)$。
证明　　$\because$ 欧拉函数是积性函数，$(i,p)=1$
　　　　$\therefore \phi (i\times p)=\phi (i)\times \phi (p)$
　　　　又 $p\in prime$
　　　　$\therefore \phi (p)=p-1$
　　　　$\therefore \phi (i\times p)=\phi (i)\times (p-1)$
代码：

void euler_phi(int n){
    memset(phi,0,sizeof phi);//欧拉函数表
    memset(vis,0,sizeof vis);//是否是质数
	phi[1]=1;
    for(int i=2;i<=n;i++){
	    if(!vis[i]){//找到质数
		    p[++tot]=i;//保存质数 
		    phi[i]=i-1;//质数的欧拉函数为它自身减1
		}
	    for(int j=1;j<=tot&&i*p[j]<=n;j++){
		    vis[i*p[j]]=1;//筛去合数
		    if(i%p[j]==0){
			    phi[i*p[j]]=phi[i]*p[j];//性质1
			    break;//继续下去没有意义，见欧拉筛法 
			}else
			    phi[i*p[j]]=phi[i]*(p[j]-1);//性质2
		}
	}
}

莫比乌斯函数

void euler_mu() {
	memset(mu,0,sizeof mu);//莫比乌斯函数表
	memset(vis,0,sizeof vis);//是否是质数
    mu[1]=1;//特判
    for(int i=2;i<=n;++i) {
	    if(!vis[i]) p[++tot]=i,mu[i]=-1;//是质数
	    for(int j=1;j<=tot&&i*p[j]<=n;++j) {
		    vis[i*p[j]]=1;
		    if(i%p[j]==0) {
			    mu[i*p[j]]=0;//有完全平方因子
			    break;
			}
		    mu[i*p[j]]=-mu[i];//积性函数(p是质数，mu[p[j]]=-1)
		}
	}
}

一般积性函数

%% 　　如果想要筛出函数 $f$ ，需要能快速求出下面三个位置的函数取值。
$$f(1),f(p),f(p^k)$$　　其实就是含有的质因子数小于等于1的所有数对应的函数值。其中 $p$ 是质数，$k$ 是正整数。
　　对于一个常见积性函数，以上所有函数值都会给出，但如果是一个一般积性函数，就需要自行计算，上面的取值，现在假设我们已经求出了上面三种位置的函数值。
考虑在筛的时候，外层循环枚举到 $i$，内层循环枚举到 $pri[j]$，也就是说，当前正在筛 $i\times pri[j]$ 并计算 $f(i\times pri[j])$。
考虑对 $i$ 唯一分解： 　$$\begin{array}{r}i=\prod _{i=1}^k{p}_i^{a_i}=p_1^{a_1}{p}_2^{a_2}\dots p_k^{a_k}\end{array}$$　　由欧拉筛法的精华所在（当 $i\equiv 0(\mathrm{mod}p[j])$ 时可以跳出内层循环），我们可以发现，内层循环满足 $pri[j]⩽p_1$（当且仅当跳出内层循环时取 $“=”$）。下面开始分类讨论：
Part 1: 当 $pri[j]<p_1$ 时，必然有 $(i,pri[j])=1$，根据积性函数的定义有 $f(i\times pri[j])=f(i)\times f(pri[j])$。
Part 2:当 $pri[j]=p_1$ 时，我们还需要记录下每个数的最小值因子的指数次幂 $low[i]=p_1^{a_1}$，这个数值可以递推。有了这个数值以后，可以发现， $\begin{array}{r}\frac{i}{low[i]}\end{array}$ 不含有 $p_1$ 因子，而 $low[i]\times pri[j]$ 只含有 $p_1$ 因子（$pri[j]=p_1$），因而有 $$\begin{array}{r}\left(\frac{i}{low[i]},low[i]\times pri[j]\right)=1\end{array}$$

故$$f\left(\frac{i}{low[i]}\right)\times f(low[i]\times pri[j])=f\left(i\times pri[j]\right)$$

至此，我们求出了所有位置的函数值。

void sieve(long long n){
    vis[1]=low[1]=1;
    f[1]=对1直接定义;
    for (long long i=2;i<=n;i++){
        if(!vis[i]) low[i]=pri[++tot]=i,f[i]=对质数直接定义;
        for(long long j=1;j<=tot&&i*pri[j]<=n;j++){
            vis[i*pri[j]]=1;
            if(i%pri[j]==0){
                low[i*pri[j]]=low[i]*pri[j];
                if(low[i]==i) f[i*pri[j]]=对质数的若干次幂进行定义(一般由f[i]递推);
                else f[i*pri[j]]=f[i/low[i]]*f[low[i]*pri[j]];
                break;
            }
			low[i*pri[j]]=pri[j];
            f[i*pri[j]]=f[i]*f[pri[j]];
        }
    }
}

Dirichlet 卷积

定义

数论函数　定义域为正整数，陪域¹（不懂）为复数的函数称为数论函数。
$\text{Dirichlet}$ 卷积　定义两个数论函数 $f,g$ 的 $\text{Dirichlet}$ 卷积为 $$(f\ast g)(n)=\sum _{d|n}f(d)g(\frac{n}{d})$$　　也叫做狄利克雷卷积。可以验证，它的单位元为 $ϵ=e$。

性质

⒈　$f\ast g=g\ast f$
⒉　$(f\ast g)\ast h=f\ast (g\ast h)$
⒊　$(f+g)\ast h=f\ast h+g\ast h$
⒋　$f\ast ϵ=f$
⒌　$d=1\ast 1$ （代入定义式即得）
⒍　$\sigma =d\ast 1$ （代入定义式即得）
⒎　$ϵ=\mu \ast 1$
⒏　$id=\phi \ast 1$
⒐　$\phi =\mu \ast id$
⒑　若 $f,g$ 均为积性函数，则 $f\ast g$ 为积性函数
⒒　若 $c\ast a=c\ast b$ 且 $c(1)\ne 0$，则 $a=b$

证明

性质3证明
$\begin{array}{rl}(f+g)\ast h& =\sum _{d|n}(f+g)(d)h(\frac{n}{d})\\ & =\sum _{d|n}(f(d)+g(d))h(\frac{n}{d})\\ & =\sum _{d|n}f(d)h(\frac{n}{d})+g(d)h(\frac{n}{d})\\ & =\sum _{d|n}f(d)h(\frac{n}{d})+\sum _{d|n}g(d)h(\frac{n}{d})\\ & =f\ast h+g\ast h\end{array}$

性质7证明
性质7等价于 $\begin{array}{r}\sum _{d|n}^{{}^{{}^{{}^{{}^{}}}}}\mu (d)=[n=1]\end{array}$
设 $\begin{array}{r}n=\prod _{i=1}^k{p}_i^{c_i},n^{\prime }=\prod _{i=1}^k{p}_i\end{array}$
则有 $\begin{array}{r}\sum _{d|n}\mu (d)=\sum _{d|n^{\prime }}\mu (d)=\sum _{i=0}^k{C}_k^i\times (-1{)}^i\end{array}$
当 $k=0$ 时显然成立，否则有二项式定理

二项式定理：$\begin{array}{r}(a+b{)}^n=\sum _{i=0}^{{}^{{}^{{}^{}}}{n}^{{}^{{}^{}}}}{C}_n^i{a}^{n-i}{b}^i\end{array}$

当 $a=1,b=-1$ 时，有 $\begin{array}{r}(1-1{)}^n=\sum _{i=0}^{{}^{{}^{{}^{}}}{n}^{{}^{{}^{}}}}(-1{)}^i\times C_n^i=0\end{array}$，得证。

性质8证明
性质8等价于 $\begin{array}{r}\sum _{d|n}^{{}^{{}^{{}^{}}}}\phi (d)=n\end{array}$
化简分数 $\frac{1}{n},\frac{2}{n},\frac{3}{n},\cdots ,\frac{n}{n}$
  \\\,
最简分数 $\frac{a}{b}$ 在上面出现的话，当且仅当 $b|n$ 且 $(a,b)=1$
那么不同的分母共有 $\begin{array}{r}\sum _{x|n}^{{}^{{}^{{}^{{}^{}}}}}1=d(n)\end{array}$ 个
其中，以 $b$ 为分母的分数共 $\phi (b)$ 个
故 $\begin{array}{r}\sum _{d|n}^{{}^{{}^{{}^{}}}}\phi (d)=n\end{array}$，得证。

性质9证明
　　易得 $id=id\ast ϵ$
　　$\because \phi \ast 1=id,\mu \ast 1=ϵ$
　　$\therefore \phi \ast 1=id\ast \mu \ast 1$
　　两边同时约去 $1$，得 $\phi =\mu \ast id$，得证。

性质10证明
前置结论：1. 若 $(n,m)=1,d|nm$，则必定存在 $a|n,b|m$ 且满足 $ab=d$。
　　　　　2. 若 $(n,m)=1,a|n,b|m$，则有 $(a,b)=1$。
　　　　　感性理解一下即可。
令 $(n,m)=1$，则有
$$(f\ast g)(nm)=\sum _{d|nm}f(d)g(\frac{nm}{d})$$根据前置结论1，有：
$$(f\ast g)(nm)=\sum _{a|n,b|m}f(ab)g(\frac{nm}{ab})$$根据前置结论2，可得：$(a,b)=1$。又因为 $f,g$ 均为积性函数，因此有：
$$\begin{array}{rl}(f\ast g)(nm)& =\sum _{a|n,b|m}f(a)f(b)g(\frac{n}{a})g(\frac{m}{b})\\ & =\sum _{a|n}f(a)g(\frac{n}{a})\sum _{b|m}f(b)g(\frac{m}{b})\\ & =(f\ast g)(n)\times (f\ast g)(m)\end{array}$$

因而 $f\ast g$ 为积性质函数。
性质11证明
　　很容易发现，本条性质说明狄利克雷卷积满足消去律。
　　要证明性质11，只需要证明在狄利克雷卷积意义下，所有函数都有逆元即可。即对于所有函数 $f$，存在函数 $f^{-1}$，使得 $f\ast f^{-1}=ϵ$。
　　注意到 $ϵ(1)=1$ 因而若函数 $f$ 满足 $f(1)=0$，则这个函数不可能有逆元，因为其与任何函数的狄利克雷卷积在 $1$ 处的取值必然为 $0$。
用数学表示，就是：设有数论函数 $f$，若 $f(1)=0$，则不可能存在函数 $f^{-1}$ 使得 $f\ast f^{-1}=ϵ$。因为当 $n=1$ 时， $(f\ast f^{-1})(n)=f(1)f^{-1}(1)=0\ne ϵ(1)=1$。
　　事实上，对于函数 $f$，当 $f(1)\ne 0$ 时，其狄利克雷卷积意义下的逆元必然存在，为：
$$f^{-1}(n)=\frac{1}{f(1)}(ϵ(n)-\sum _{d|n,d\ne 1}f(d)f^{-1}(\frac{n}{d}\left)\right)$$　　因此，我们只需要证明 $f\ast f^{-1}=ϵ$ 即可。容易得到：
$$\begin{array}{rl}(f\ast f^{-1})(n)& =\sum _{d|n}f(d)f^{-1}\left(\frac{n}{d}\right)\\ & =f(1)f^{-1}(n)+\sum _{d|n,d\ne 1}f(d)f^{-1}\left(\frac{n}{d}\right)\\ & =f(1)\times \frac{1}{f(1)}(ϵ(n)-\sum _{d|n,d\ne 1}f(d)f^{-1}(\frac{n}{d}\left)\right)+\sum _{d|n,d\ne 1}f(d)f^{-1}\left(\frac{n}{d}\right)\\ & =ϵ(n)-\sum _{d|n,d\ne 1}f(d)f^{-1}\left(\frac{n}{d}\right)+\sum _{d|n,d\ne 1}f(d)f^{-1}\left(\frac{n}{d}\right)\\ & =ϵ(n)\end{array}$$　　得证。

中心结论

结论　$\begin{array}{r}\sum _{d|(i,j)}^{{}^{{}^{}}}\mu (d)=[(i,j)=1]\end{array}$
证明　$\because \mu \ast 1=ϵ$
　　　$\begin{array}{r}\therefore \sum _{d|n}\mu (d)=[n=1]\end{array}$
　　　将 $n=(i,j)$ 代入，得证。

莫比乌斯反演

结论　设 $f,g$ 均为数论函数，若满足 $\begin{array}{r}f(n)=\sum _{d|n}^{{}^{{}^{{}^{{}^{}}}}}g(d)\end{array}$，则有 $\begin{array}{r}g(n)=\sum _{d|n}^{{}^{{}^{{}^{{}^{}}}}}\mu (d)f(\frac{n}{d})\end{array}$
证明　原命题等价于：若 $f=g\ast 1$，则 $g=f\ast \mu$
　　　两边同时卷 $\mu$ 得 $f\ast \mu =g\ast 1\ast \mu$
　　　$\because 1\ast \mu =ϵ$
　　　$\therefore f\ast \mu =g\ast 1\ast \mu \Rightarrow f\ast \mu =g$，得证
没错，就这么一点。因为有了上面的证明，这里的证明十分简单。

变形1　若 $\begin{array}{r}F(n)=\sum _{n|d}^{{}^{{}^{{}^{{}^{}}}}}f(d)\end{array}$，则 $\begin{array}{r}f(n)=\sum _{n|d}^{{}^{{}^{{}^{{}^{}}}}}\mu (\frac{d}{n})F(d)\end{array}$

变形2　若 $\begin{array}{r}f(k)=\sum _{{}_{{}_{{}_{{}_{}}}}i=1_{{}_{{}_{{}_{}}}}}^{\lfloor \frac{n}{k}\rfloor }g(ik)\end{array}$，则 $\begin{array}{r}g(k)=\sum _{{}_{{}_{{}_{{}_{}}}}i=1_{{}_{{}_{{}_{}}}}}^{\lfloor \frac{n}{k}\rfloor }f(ik)\mu (i)\end{array}$

题表

题目	题解
[POI2007]ZAP-Queries	莫比乌斯反演
[HAOI2011]Problem b	莫比乌斯反演
[NOI2010]能量采集	欧拉反演

参考资料

[1] 租酥雨,积性函数与线性筛

回收区

无视就好。

**性质１**若 $m\in prime,n≡0 \pmod m$ 则 $\varphi(n\times m)=\varphi(n)\times m$
**证明**构造矩阵
$\begin{bmatrix}
&1&2&\cdots&r&\cdots&n\\
&n+1&n+2&\cdots&n+r&\cdots&2n\\
&2n+1&2n+2&\cdots&2n+r&\cdots&3n\\
&\vdots&\vdots&\ddots&\vdots&\ddots&\vdots\\
&(m-1)n+1&(m-1)n+2&\cdots&(m-1)n+r&\cdots&mn\\
\end{bmatrix}$
由于 $(n,m)=1$，因此与 $nm$ 互质的数相当于与 $n$ 互质且与 $m$ 互质的数的个数。
$\because (kn+r,n)=(r, n)$
$\therefore$ 每一列的 $m$ 个元素同时与 $n$ 互素当且仅当 $(r,n)=1$。
$\therefore$ 共有 $\varphi(n)$ 列，其余列均不与 $n$ 互质。
$\,$其中满足 $(r,n)=1$ 的第 $r$ 列第 $k$ 行的数为 $kn+r$
又 $n≡0\pmod m$
$\therefore kn+r≡0\times k+r≡r\pmod m$
$\therefore (kn+r,m)=(r,m)$
又 $m\in prime,r<m$
$\therefore (kn+r,m)=1$
$\therefore$ 满足 $(r,n)=1$ 的第 $r$ 列的所有 $m$ 个数均与 $n$ 互质。
$\therefore \varphi(n\times m)=\varphi(n)\times m$
　

**性质１**　对于质数 $p$ 和正整数 $k$，有 $\varphi(p^k)=\varphi(p^{k-1})\times p$
**证明**　　由于 $p\in prime$，因此除了 $p$ 的倍数之外，$1\sim p^k$ 之间所有数都与 $p^k$ 互质。也就是 ${\forall (1\leqslant t\leqslant p^k )∧(p\not|\,\;t),\text{gcd}(t,p^k)=1}$。而 $1\sim p^k$之间有 $\dfrac{p^k}{p}$ 个数是 $p$ 的倍数。
　　　　因此有 $\varphi(p^k)=p^k-\dfrac {p^{k}}{p}=p^k\left(1-{\dfrac 1p}\right)$。
　　　　因而 $\varphi(p^k)=p^{k-1}\left(1-{\dfrac 1p}\right)\times p=\varphi(p^{k-1})\times p$
　　　　更进一步，我们还可以得到 $\varphi$ 的通项公式：
　　　　$\begin{aligned}\varphi(p)&=\varphi(p_1^{a_1})\times \varphi(p_2^{a_2})\times \dots\times \varphi(p_k^{a_k})\\
&=n\times (1-\frac {1}{p_1})\times (1-\frac {1}{p_2})\times ……\times (1-\frac {1}{p_k})
\end{aligned}$

---

1. 陪域是映射中的概念，其包含了值域在内，一般就是指值域，但又略有差别。其中值域是上下界绝对严格的，即不会出现边界虚大或虚小的情况，而陪域包括了不包含值域的部分，因而陪域仅仅是绝对包含值域的。 ↩︎