莫比乌斯反演

最新推荐文章于 2022-05-03 02:59:51 发布

ezoiHY

最新推荐文章于 2022-05-03 02:59:51 发布

阅读量2.7k

点赞数

分类专栏：学习笔记数学莫比乌斯反演

本文链接：https://blog.csdn.net/EZharry/article/details/81734896

版权

学习笔记同时被 3 个专栏收录

15 篇文章 0 订阅

订阅专栏

数学

7 篇文章 0 订阅

订阅专栏

莫比乌斯反演

1 篇文章 0 订阅

订阅专栏

莫比乌斯反演是一个十分玄幻的东西，它可以把 $o(n^2)$ 的时间复杂度降到 $o(n\sqrt{n})$ 甚至更低

1.公式

这是莫比乌斯反演最基本的东西，两个定义在正整数集上的函数 $F(n)$ 和 $f(n)$
若满足这个式子

F (n) = \sum d | n f (d)

$F(n)=\sum_{d|n}f(d)$
则会有

f (n) = \sum d | n μ (d) F (n d)

$f(n)=\sum_{d|n}\mu(d)F(\frac{n}{d})$

2.莫比乌斯函数 $\mu$

$\mu$ 是个神奇的东西，它是一个定义在正整数集上的函数

μ (n) = ⎧ ⎩ ⎨ ⎪ ⎪ 1 (- 1) k 0 n = 1 n = p 1 p 2 p 3 . . . p k (其 中 p i 为 互 异 质 数) 其 他 情 况

$\mu(n)=\left\{ \begin{array}{rcl} 1 & & n=1\\ (-1)^k & & n=p_1p_2p_3...p_k (其中p_i为互异质数)\\ 0 & & 其他情况\\ \end{array} \right.$
这是前50个莫比乌斯函数的取值
这里写图片描述

它有两个性质
性质一：

\sum d | n μ (d) = {10 n = 1 n! = 1

$\sum_{d|n}{\mu(d)}=\left\{ \begin{array}{rcl} 1 & & n=1\\ 0 & & n!=1 \end{array} \right.$
性质二：

\sum d | n μ ( d ) d = ϕ ( n ) n

$\sum_{d|n}\frac{\mu(d)}{d}=\frac{\phi(n)}{n}$

3.莫比乌斯反演的实现

其实线性筛就能求出 $\mu$
代码：

void get(int n){
    m[1]=1;
    for(int i=2;i<=n;i++){
        if(!vis[i]){
            pri[++cnt]=i;
            m[i]=-1;
        }
        for(int j=1;j<=cnt and pri[j]*i<=n;j++){
            int k=pri[j]*i;
            vis[k]=1;
            if(i%pri[j]==0){
                m[k]=0;
                break;
            }else m[k]=-m[i];
        }
    }
}