【Mobius绮丽的辗转】莫比乌斯反演

Problem

对于给出的n个询问，每次求有多少个数对(x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y) = k，gcd(x,y)函数为x和y的最大公约数。
$1≤n≤50000，1≤a≤b≤50000，1≤c≤d≤50000，1≤k≤50000$

Sub problem

设$Ans(i,j)$表示有多少个数对(x,y)，满足x≤i，c≤y≤j，且gcd(x,y) = k。
我们可以先求出$Ans(b,d),Ans(b,c-1),Ans(a-1,d),Ans(a-1,c-1)$，
然后$ans=Ans(b,d)-Ans(b,c-1)-Ans(a-1,d)+Ans(a-1,c-1)$。

那么问题就变成了如何求$Ans(n,m)$。

Discuss

讨论一下$Ans(n,m)$如何求，其中$n<m$。
先设$f(k)$，表示有多少个数对(x,y)，满足x≤n，c≤y≤m，且gcd(x,y) = k。
显然$Ans(n,m)=f(k)$。

再设$g(k)$，表示多少个数对(x,y)，满足x≤n，c≤y≤m，且k|gcd(x,y)
因为$k|gcd(x,y)$，所以设$x=k*x',y=k*y'$；
由于$x'$只能取$1...\lfloor\frac{n}{k}\rfloor$，$y'$只能取$1...\lfloor\frac{m}{k}\rfloor$

所以

$$g(k)=\lfloor\frac{n}{k}\rfloor*\lfloor\frac{m}{k}\rfloor$$

同时，我们会有

$$g(k)=\sum^{\lfloor\frac{n}{k}\rfloor}_{i=1}f(i*k)\tag1$$

此时，我们将$g(k)$用$f(k)$表示，并且$g(k)$是容易求出结果的。

---

Mobius

正片开始

我们非常功利地得出结论：
正当我们遇到这种式子时，

$$g(i)=\sum_{d|i}f(d)\tag2$$ 或
$$g(i)=\sum_{j=1}^{\lfloor\frac{n}{k}\rfloor}f(i*j)\tag3$$
当 $g[d]$是 积性函数，我们可以将上式转化为，
$$f(i)=\sum_{d|i}g(d)*\mu(\frac{i}{d})$$
$$f(i)=\sum_{j=1}^{\lfloor\frac{n}{k}\rfloor}g(i*j)*\mu(j)$$
其中
$$\mu(x)=\begin{cases} 1,&x=1\\ (-1)^k,&x=\prod_{i=1}^{k}p_i,p_i\in \mathbb P\\ 0,&otherwise \end{cases}$$

---

Discuss：$\mu$的性质

（1）$\mu$是积性函数
可以证明，$\mu$函数也是积性函数，所以$\mu$可以通过线性筛法预处理，如下代码。

miu[1]=1;
for (i=2;i<maxn;i++){
    if (!bz[i]){
        p[++p[0]]=i;
        miu[i]=-1;
    }
    for (j=1;j<=p[0];j++){          
        k=i*p[j];
        if (k>=maxn) break;
        bz[k]=true;
        if (i%p[j]==0){
            miu[k]=0;
            break;
        }else miu[k]=-miu[i];
    }
}

（2）$\mu$的“和性质”

$$\sum_{d|n}\mu(d)= \begin{cases} 0,&n=1\\ 1,&n>1 \end{cases}$$

Back to the Problem

题目的式子是

$$g(k)=\sum^{\lfloor\frac{n}{k}\rfloor}_{i=1}f(i*k)\tag1$$
跟 $(2)$有异曲同工之妙，
所以
$$\begin{align} f(k)&=\sum_{i=1}^{\lfloor\frac{n}{k}\rfloor}g(i*k)*\mu(i)\\ &=\sum_{i=1}^{\lfloor\frac{n}{k}\rfloor}\lfloor\frac{n}{ik}\rfloor*\lfloor\frac{m}{ik}\rfloor*\mu(i) \end{align}$$
然而，这并没有什么卵用，我们仍然过不了。
还能优化？？

Deeplier discuss

我们发现，
其实$\lfloor\frac{n}{ik}\rfloor*\lfloor\frac{m}{ik}\rfloor$很多时候是相同的取值。
所以我们可以将相同值的$\lfloor\frac{n}{ik}\rfloor*\lfloor\frac{m}{ik}\rfloor$合并一起来计算，来优化时间复杂度。
显然$\lfloor\frac{n}{ik}\rfloor$的取值最多有$2*\sqrt n$种，
所以可以把时间复杂度优化到$O(2*\sqrt n+2*\sqrt m)$一次询问。

Ending

至此，我们已解决了这道题。
原题Code。

Proving

$\mu$的“和性质”

求证：

$$\sum_{d|n}\mu(d)= \begin{cases} 0,&n=1\\ 1,&n>1 \end{cases}$$
证明：
$n=1$时显然；
讨论 $n>1$的情况，
因为 $\mu$的定义，
$$\mu(x)=\begin{cases} 1,&x=1\\ (-1)^k,&x=\prod_{i=1}^{k}p_i,p_i\in \mathbb P\\ 0,&otherwise \end{cases}$$
所以 $\sum_{d|n}\mu(d)$中，只有当 $d$的任意质因子的指数不能超过$1$时， $\mu(d)$才会对 和产生贡献。
我们设 $n$的质因子个数为$q$个。
那么，
$$\sum_{d|n}\mu(d)=\sum_{i=0}^qC^i_q*(-1)^i$$
我们观察一下杨辉三角：
$$1\\ 1\ 1\\ 1\ 2\ 1\\ 1\ 3\ 3\ 1\\ 1\ 4\ 6\ 4\ 1\\ 1\ 5\ 10\ 10\ 5\ 1\\ 1\ 6\ 15\ 20\ 15\ 6\ 1\\ ...$$
显然的是，当 $q$是偶数时，由杨辉三角的对称性， $$\sum_{d|n}\mu(d)=\sum_{i=0}^qC^i_q*(-1)^i=0$$
现在考虑 $q(q>1)$是奇数的情况，
$$\begin{align} \sum_{i=0}^qC^i_q*(-1)^i &=C^0_q-C^q_q+\sum_{i=1}^{q-1}C^i_q*(-1)^i\\ &=C^0_q-C^q_q+\sum_{i=1}^{q-1}(C^{i-1}_{q-1}+C^{i}_{q-1})*(-1)^i\\ &=C^0_q-C^q_q+\sum_{i=1}^{q-1}C^{i-1}_{q-1}*(-1)^i+\sum_{i=1}^{q-1}C^{i}_{q-1}*(-1)^i\\ &=\sum_{i=0}^{q-1}C^i_{q-1}*(-1)^{i+1}+\sum_{i=0}^{q-1}C^{i}_{q-1}*(-1)^i\\ \end{align}$$
由 $q-1$是偶数，综上， $$\sum_{d|n}\mu(d)=\sum_{i=0}^qC^i_q*(-1)^i=0$$
得证。

证明反演

求证：

$$g(i)=\sum_{d|i}f(d)\Rightarrow f(i)=\sum_{d|i}g(d)*\mu(\frac{i}{d})$$
证明：
$$\begin{align} \sum_{d|i}g(d)*\mu(\frac{i}{d}) &=\sum_{d|i}\mu(\frac{i}{d})\sum_{d'|d}f(d')\\ \end{align}$$
这里经历一个重要的过程：转换主体，
感性地想，所有的 $\mu(\frac{i}{d})$都与 $f(d')$相乘过，其中 $d'|d$；
反过来，那么所有的 $f(d')$都与 $\mu(\frac{i}{d})$相乘过，其中 $d'|d$。
所以，
$$\begin{align}\sum_{d|i}\mu(\frac{i}{d})\sum_{d'|d}f(d') &=\sum_{d'|i}f(d')\sum_{d'|d,d|i}\mu(\frac{i}{d})\\ \end{align}$$
令 $x=\frac{i}{d}$，则 $d=\frac{i}{x}$，那么
$$\begin{align} \sum_{d'|i}f(d')\sum_{d'|d,d|i}\mu(\frac{i}{d}) &=\sum_{d'|i}f(d')\sum_{d'|\frac{i}{x},\frac{i}{x}|i}\mu(x)\\ &=\sum_{d'|i}f(d')\sum_{x|\frac{i}{d'}}\mu(x)\\ \end{align}$$
由 $\mu$的“和性质”，
当 $d'!=i$时，则 $\frac{i}{d'}>1$，所以 $\sum_{x|\frac{i}{d'}}\mu(x)=0$；
当 $d'=i$时，则 $\frac{i}{d'}=1$，所以 $\sum_{x|\frac{i}{d'}}\mu(x)=1$。
所以 $$\sum_{d|i}f(d)*\mu(\frac{i}{d})=\sum_{d'|i}f(d')\sum_{x|\frac{i}{d'}}\mu(x)=f(i)$$
综上， $$f(i)=\sum_{d|i}g(d)*\mu(\frac{i}{d})$$
得证。

---

另一个变式$(3)$类似。

True Ending

至此，Mobius反演已证明完毕。