莫比乌斯反演总结

参考文档：http://download.csdn.net/detail/williamsun0122/9845697

我是看上面的文档学会的，大家可以参考参考，讲的还比较详细。
没看懂莫比乌斯反演证明，所以没写

莫比乌斯反演

公式

有两种形式

$$(约数和)\\ F(n)=\sum_{d|n} f(d) \iff f(n)=\sum_{d|n} \mu(d)F(\frac{n}{d})\\ (倍数和)\\ F(n)=\sum_{n|d} f(d) \iff f(n)=\sum_{n|d} \mu(\frac{d}{n})F(d)\\ (其中d|n表示n\%d=0即d为n的因子，\mu()是莫比乌斯函数)$$

应用

对于一些函数f(n)，如果我们很难直接求出它的值，而容易求出倍数和或约数和F(n)，那么我们可以通过莫比乌斯反演来求得f(n)的值。
之后会写一些例题的博客，看了就差不多懂了。

莫比乌斯函数

定义

$$设n=p_1^{k_1}*p_2^{k_2}*...*p_m^{k_m},其中p_i为互异素数\\ \mu(n)= \begin{cases} 1 & \text{n=1} \\ (-1)^m & \text{$\prod_{i=1}^m k_i=1即n的所有质因子次数均为1$}\\ 0 & \text{others} \end{cases}$$

性质

1. 对于任意正整数n
   

   $$\sum_{d|n} \mu(d)= \begin{cases} 1 & \text{n=1} \\ 0 & \text{n>1} \end{cases}$$
   证明：
   （1）n=1时显然
   （2）n>1时，n可以分解为$n=p_1^{a_1}*p_2^{a_2}*...*p_m^{a_m}$。根据定义知在n的所有因子d中，$\mu(d)$的值不为0时，只有d的所有质因子次数都为1（即$d=p_1,p_1*p_2,...,或p_1*p_2*...*p_m$）。那么质因子个数为r的d值有$C_m^r$个。所以：
   $$\sum_{d|n} \mu(d)=C_k^0-C_k^1+C_k^2+...+(-1)^k C_k^k=\sum_{i=0}^k (-1)^i C_k^i$$
   现在只需证明$\sum_{i=0}^n (-1)^i C_n^i=0$即可
   由二项式定理：$(x+y)^n=\sum_{i=0}^n C_n^i x^i y^{n-i}$，令x=-1，y=1得证。

2. 莫比乌斯函数是积性函数：$\mu(a*b)=\mu(a)*\mu(b),gcd(a,b)=1$

---

以上是莫比乌斯反演的一些知识，具体怎么应用我之后会在一些例题中说明。

下面我们运用莫比乌斯函数的性质写一个求莫比乌斯函数的模板：

const int maxn = 1e4+5    //所求莫比乌斯函数的范围，我以10000为例
int mu[maxn];        //莫比乌斯函数值
int vis[maxn];      //vis[i]=1表示i不是素数，为0表示i是素数
int prime[maxn];    //存放素数

void Moblus()
{
    memset(vis,0,sizeof(vis);
    int num=0;  //素数的个数
    mu[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!vis[i])   //i是素数
        {
            mu[i]=-1;  //根据定义知单个素数其莫比乌斯函数值为(-1)^1=-1
            prime[num++]=i;  //存进prime[]中
        }
        for(int j=0;j<num && prime[j]*i<maxn;j++)
        {
            vis[prime[j]*i]=1;  //素数的倍数不为素数
            if(i%prime[j]==0)
            {
                mu[prime[j]*i]=0;   //质因子的次数不为1
                break;  //这里为什么跳出我只有自己的理解，不确定对不对，大家就自己在纸上推一推。知道的请给我留言，谢谢。
            }
            else mu[prime[j]*i] = -mu[i];   //根据莫比乌斯函数的性质mu[prime[j]*i]=mu[prime[j]]*mu[i]=-mu[i]，因为prime[j]为单素数
        }
    }
}