莫比乌斯反演、杜教筛

莫比乌斯反演

莫比乌斯反演基本形式：

对于一个函数$f(x)$

设$g(x)={\sum }_{x\mid d}f(d)$，那么

$f(x)={\sum }_{x\mid d}\mu (\frac{d}{x})\cdot g(d)$

---

$1$：$f(x)=1$

$e$：$f(x)=[x=1]$

$id$：$f(x)=x$

---

推式子：

1. 1到$n$与$n$互质的数的个数：${\sum }_{i=1}^n[gcd(i,n)=1]$

   令$f(x)={\sum }_{i=1}^n[gcd(i,n)=x]$

   $$\begin{gathered} g(x)={\sum }_{x\mid d}f(d)={\sum }_{x\mid d}{\sum }_{i=1}^n[gcd(i,n)=d]={\sum }_{i=1}^n[x\mid gcd(i,n)] \\ =\{\begin{array}{ll}\lfloor \frac{n}{x}\rfloor ,& x\mid n\\ 0,& x̸\mid n\end{array} \end{gathered}$$

   由莫比乌斯反演得：

   $$\begin{gathered} f(x)={\sum }_{x\mid d}\mu (\frac{d}{x})\cdot g(d)=\{\begin{array}{ll}{\sum }_{x\mid d}\mu (\frac{d}{x})\cdot \lfloor \frac{n}{d}\rfloor ,& d\mid n\\ 0,& d̸\mid n\end{array} \\ ={\sum }_{d\mid n}{\sum }_{x\mid d}\mu (\frac{d}{x})\cdot \lfloor \frac{n}{d}\rfloor \end{gathered}$$

   所求式子即为$f(1)={\sum }_{d\mid n}\mu (d)\cdot \lfloor \frac{n}{d}\rfloor$

   根据欧拉函数定义，所求即为$\phi (n)$。

   因此$\phi (n)={\sum }_{d\mid n}\mu (d)\cdot \frac{n}{d}={\sum }_{d\mid n}\mu (\frac{n}{d})\cdot d$

   ---

2. 1D gcd sum：${\sum }_{i=1}^{n}gcd(i,n)={\sum }_{d\mid n}\phi (d)\lfloor \frac{n}{d}\rfloor$

   ${\sum }_{i=1}^{n}gcd(i,n)={\sum }_{d=1}^{n}d{\sum }_{i=1}^n[gcd(i,n)=d]$

   设$f(x)={\sum }_{i=1}^n[gcd(i,n)=x]$

   $$\begin{gathered} g(x)={\sum }_{x\mid d}f(d)={\sum }_{x\mid d}{\sum }_{i=1}^n[gcd(i,n)=d] \\ ={\sum }_{i=1}^n[x\mid gcd(i,n)]=\{\begin{array}{ll}\lfloor \frac{n}{x}\rfloor ,& x\mid n\\ 0,& x̸\mid n\end{array} \end{gathered}$$

   由莫比乌斯反演得：

   $$\begin{gathered} f(x)={\sum }_{x\mid d}\mu (\frac{d}{x})\cdot g(d)=\{\begin{array}{ll}{\sum }_{x\mid d}\mu (\frac{d}{x})\cdot \lfloor \frac{n}{d}\rfloor ,& d\mid n\\ 0,& d̸\mid n\end{array} \\ ={\sum }_{d\mid n}{\sum }_{x\mid d}\mu (\frac{d}{x})\cdot \lfloor \frac{n}{d}\rfloor \end{gathered}$$

   代回得：
   $$\begin{gathered} {\sum }_{i=1}^{n}gcd(i,n)=\{\begin{array}{ll}{\sum }_{d=1}^{n}d{\sum }_{d\mid x}\mu (\frac{x}{d})\cdot \lfloor \frac{n}{x}\rfloor ,& x\mid n\\ 0,& x̸\mid n\end{array} \\ ={\sum }_{x\mid n}{\sum }_{d\mid x}d\cdot \mu (\frac{x}{d})\cdot \lfloor \frac{n}{x}\rfloor ={\sum }_{x\mid n}\phi (x)\cdot \lfloor \frac{n}{x}\rfloor \end{gathered}$$

   (参考$\phi (n)={\sum }_{d\mid n}\mu (\frac{n}{d})\cdot d$)

   即${\sum }_{i=1}^{n}gcd(i,n)={\sum }_{d\mid n}\phi (d)\lfloor \frac{n}{d}\rfloor$

   ---

3. 2D gcd sum：${\sum }_{i=1}^n{\sum }_{j=1}^{m}gcd(i,j)$

   ${\sum }_{i=1}^n{\sum }_{j=1}^{m}gcd(i,j)={\sum }_{d=1}^{n}d{\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=d]$

   设$f(x)={\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=x]$

   $$\begin{gathered} g(x)={\sum }_{x\mid d}f(d)={\sum }_{x\mid d}{\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=d] \\ ={\sum }_{i=1}^n{\sum }_{j=1}^m[x\mid gcd(i,j)]=\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{x}\rfloor \end{gathered}$$

   由莫比乌斯反演得：

   $f(x)={\sum }_{x\mid d}\mu (\frac{d}{x})\cdot g(d)={\sum }_{x\mid d}\mu (\frac{d}{x})\cdot \lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$

   代回得：
   $$\begin{gathered} {\sum }_{d=1}^{n}d\cdot f(d)={\sum }_{d=1}^{n}d{\sum }_{d\mid x}\mu (\frac{x}{d})\cdot \lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{x}\rfloor \\ ={\sum }_{x=1}^n{\sum }_{d\mid x}d\cdot \mu (\frac{x}{d})\cdot \lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{x}\rfloor ={\sum }_{x=1}^n\phi (x)\cdot \lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{x}\rfloor \end{gathered}$$

   ---

4. 1…n和1…m中互质的数的对数

   ${\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=1]$

   设$f(x)={\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=x]$

   由莫比乌斯反演得：

   $f(x)={\sum }_{x\mid d}\mu (\frac{d}{x})\cdot g(d)={\sum }_{x\mid d}\mu (\frac{d}{x})\cdot \lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$

   所求即${\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=1]=f(1)={\sum }_{d=1}^n\mu (d)\cdot \lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$

---

杜教筛

求积性函数$f(i)$的前缀和$S(i)$

$S(n)={\sum }_{i=1}^{n}f(i)$

杜教筛：

找一个积性函数$g(i)$，做$g$和$f$的狄利克雷卷积

$(g\ast f)(i)={\sum }_{d\mid i}g(d)f(\frac{i}{d})$

再做一下卷积的前缀和

${\sum }_{i=1}^n(g\ast f)(i)={\sum }_{i=1}^n{\sum }_{d\mid i}g(d)f(\frac{i}{d})$

把$d$提出：${\sum }_{d=1}^{n}g(d){\sum }_{d\mid i}f(\frac{i}{d})={\sum }_{d=1}^{n}g(d){\sum }_{i=1}^{\frac{n}{d}}f(i)={\sum }_{d=1}^{n}g(d)S(\frac{n}{d})$

$g(1)S(n)={\sum }_{i=1}^{n}g(i)S(\frac{n}{i})-{\sum }_{i=2}^{n}g(i)S(\frac{n}{i})={\sum }_{i=1}^n(g\ast f)(i)-{\sum }_{i=2}^{n}g(i)S(\frac{n}{i})$

前半部分是狄利克雷卷积的前缀和。

如果狄利克雷卷积的前缀和非常好算，那么可以对后半部分进行数论分块，递归计算(记忆化)。

---

栗子：${\sum }_{i=1}^n\mu (i)$

$g(1)S(n)={\sum }_{i=1}^n(g\ast f)(i)-{\sum }_{i=2}^{n}g(i)S(\frac{n}{i})$

看到$\mu$想到$(\mu \ast 1)=e$即${\sum }_{d\mid n}\mu (d)=[n=1]$

$e$的前缀和就是1，十分好计算

因此取$g(x)=1$

$S(n)=1-{\sum }_{i=2}^{n}S(\frac{n}{i})$

先筛出一部分$\mu$的前缀和，然后记忆化递归计算。

---

栗子2：${\sum }_{i=1}^n\phi (i)$

看到$\phi$想到$(\phi \ast 1)=id$即${\sum }_{d\mid n}\phi (d)=n$

取$g(x)=1$

$S(n)=\frac{n(n+1)}{2}-{\sum }_{i=2}^{n}S(\frac{n}{i})$