数学板块学习之莫比乌斯

莫比乌斯函数

对于函数$F(n)={\Sigma }_{d|n}f(d)$ （其中$d|n$表示$d$能够被$n$整除），根据定义有
$F(1)=f(1)$
$F(2)=f(1)+f(2)$
$F(3)=f(1)+f(3)$
$F(4)=f(1)+f(2)+f(4)$
$F(5)=f(1)+f(5)$
$F(6)=f(1)+f(2)+f(3)+f(6)$
$F(7)=f(1)+f(7)$
$F(8)=f(1)+f(2)+f(4)+f(8)$

于是
$f(1)=F(1)$
$f(2)=F(2)-F(1)$
$f(3)=F(3)-F(1)$
$f(4)=F(4)-F(2)$
$f(5)=F(5)-F(1)$
$f(6)=F(6)-F(3)-F(2)+F(1)$
$f(7)=F(7)-F(1)$
$f(8)=F(8)-F(4)$

我们可以由$F(n)={\Sigma }_{d|n}f(d)$ ，如果已知$F(n)$我们就可以反演推出$f(n)$
公式：
$$F(n)={\Sigma }_{d|n}f(d)\Rightarrow f(n)={\Sigma }_{d|n}\mu (d)F(\frac{n}{d})$$
其中的$\mu (d)$即为莫比乌斯函数，定义为

1. 若$d=1$则$\mu (d)=1$
2. 若$d=p_1{p}_2{p}_3\cdots p_k$，其中$p_i$为互异素数，则$\mu (d)=(-1{)}^k$
3. 其他情况$\mu (d)=0$

莫比乌斯函数性质

1.对于任意正整数$n$有
$${\Sigma }_{d|n}\mu (d)=\{\begin{array}{ll}1& \text{n=1}\\ 0& \text{n>1}\end{array}$$
2.对于任意正整数$n$有
$${\Sigma }_{d|n}\frac{\mu (d)}{d}=\frac{\phi (n)}{n}$$
其中的$\phi (n)$是欧拉函数
3.$\mu (d)为积性函数$，也因为它为积性函数所以可以用线性筛求出

求莫比乌斯函数代码

int prime[maxn],mu[maxn],tot = 0;
bool vis[maxn];
void Mobius(int n){
    me(vis,false);
    mu[1] = 1;
    for(int i = 2; i <= n; ++i){
        if(!vis[i]) prime[++tot] = i,mu[i] = -1;
        for(int j = 1; j <= tot & i*prime[j] <= n; ++j){
            vis[i*prime[j]] = true;
            if(i%prime[j] == 0){
                mu[i*prime[j]] = 0;
                break;
            }
            mu[i*prime[j]] = -mu[i];
        }
    }
}

莫比乌斯反演

公式：
$$F(n)={\Sigma }_{d|n}f(d)\Rightarrow f(n)={\Sigma }_{d|n}\mu (d)F(\frac{n}{d})$$
$$F(n)={\Sigma }_{n|d}f(d)\Rightarrow f(n)={\Sigma }_{n|d}\mu (\frac{d}{n})F(d)$$
证明：
$$f(n)={\Sigma }_{d|n}\mu (d)F(\frac{n}{d})={\Sigma }_{d|n}\mu (d){\Sigma }_{k|\frac{n}{d}}f(k)={\Sigma }_{k|n}f(k){\Sigma }_{d|\frac{n}{k}}\mu (d)$$
因为：
$${\Sigma }_{d|n}\mu (d)=\{\begin{array}{ll}1& \text{n=1}\\ 0& \text{n>1}\end{array}$$
所以只有$n=k$：
$${\Sigma }_{k|n}f(k){\Sigma }_{d|\frac{n}{k}}\mu (d)=f(k)=f(n)$$

做题方法(套路)：
一般题目都是将其转换为$gcd(x,y)=1$的形式求解
令$$f(x)=\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==x]$$
$$F(x)=\sum _{x|d}f(d)$$
所以
$$F(x)=\sum _{x|d}\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==d]=\sum _{i=1}^n\sum _{j=1}^m[x|gcd(i,j)]=\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{x}\rfloor$$
从而
$$f(x)=\sum _{x|d}\mu (\frac{d}{x})F(d)=\sum _{x|d}\mu (\frac{d}{x})\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$$
一般求$f(1)$
$$f(1)=\sum _{d=1}^{min(n,m)}\mu (d)\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$$

套路2：
由莫比乌斯定理得
$$[gcd(i,j)==1]=\sum _{d|gcd(i,j)}\mu (d)=\sum _{d|i,d|j}\mu (d)$$
所以
$$ans=\sum _{i=1}^n\sum _{j=1}^m\sum _{d|i,d|j}\mu (d)$$
提贡献：
$$ans=\sum _{d=1}^{min(n,m)}\mu (d)\sum _{d|i}^n\sum _{d|j}^{m}1$$
即
$$ans=\sum _{d=1}^{min(n,m)}\mu (d)\lfloor \frac{n}{d}\rfloor \lfloor \frac{n}{d}\rfloor$$
小性质：
$$d(nm)=\sum _{i|n}\sum _{j|m}[gcd(i,j)==1]$$
其中$d(x)表示x的约数的个数$