首先给大家介绍一下莫比乌斯函数吧,其实这个函数挺好理解的,只是一个容斥系数
μ(d)的定义是:
当d=1时,μ(d)=1;
当d=Πki=1pi且pi为互异素数时,μ(d)=(−1)k。(说直白点,就是d分解质因数后,没有幂次大于平方的质因子,此时函数值根据分解的个数决定);
只要当d含有任何质因子的幂次大于等于2,则函数值为0.

莫比乌斯函数的性质

1、对于任意正整数n,∑d|nμ(d)=[n=1]。([n=1]表示只有当n=1成立时,返回值为1;否则,值为0;(这个就是用μ是容斥系数的性质可以证明)(PS:这一条性质是莫比乌斯反演中最常用的)
2、对于任意正整数n,∑d|nμ(d)d=ϕ(n)n。(这个性质很奇妙,它把欧拉函数和莫比乌斯函数结合起来)

介绍完性质现在再来介绍一下那个莫比乌斯函数的公式
定理:F(n)和f(n)是定义在非负整数集合上的两个函数,并且满足条件:
F(n)=∑d|nf(d)
那么存在一个结论:

f(n)=∑d|nμ(d)F(⌊nd⌋)
这个定理就称作莫比乌斯反演定理。

这个公式的证明:

对于莫比乌斯函数的线性筛法(其实就是在筛质数的前提下进行筛,如果这个数是质数,那么他的莫比乌斯函数就为-1,如果是能被i*(primes[j])消掉的话,他的莫比乌斯函数就为0,否则就为-mobius[i]
下面请看代码:

void get_mu(int n)
{mu[1]=1;for(int i=2;i<=n;i++){if(!vis[i]){prim[++cnt]=i;mu[i]=-1;}for(int j=1;j<=cnt&&prim[j]*i<=n;j++){vis[prim[j]*i]=1;if(i%prim[j]==0)break;else mu[i*prim[j]]=-mu[i];}}}

至于狄利克雷卷积法,我学了之后再补充吧哈哈哈

注意

当d是n的倍数的时候也满足莫比乌斯反演定理,推导过程如下


这个μ是莫比乌斯函数的意思

莫比乌斯反演的证明(非狄利克雷卷积法)相关推荐

  1. matlab狄利克雷函数,数论入门1——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛...

    数论入门1 一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数$\mu(n)$为: 当n有平方因子时,$\mu(n)=0$. 当n没有平方因子时,$\mu(n)=(-1)^{\omega(n)} ...

  2. 狄利克雷卷积与莫比乌斯函数

    狄利克雷卷积与莫比乌斯函数 数论函数与积性函数 狄利克雷卷积 莫比乌斯函数 相关性质 基础问题与公式 若想使用莫比乌斯反演,熟练掌握狄利克雷卷积包括定义.记号以及相关的性质.证明等是非常有好处的. 数 ...

  3. 数学--数论--莫比乌斯反演

    一.莫比乌斯反演涉及知识 1.莫比乌斯函数 2.莫比乌斯的线性筛法 3.狄利克雷卷积 4.莫比乌斯反演详解 5.整除法分块 6.杜教筛 二.μ 莫比乌斯函数定义 μ(n)={1n=1(−1)kn= P ...

  4. 从积性函数到莫比乌斯反演

    积性函数 积性函数:对于数论函数 fff ,若任意互质的 p,qp,\ qp, q 都有 f(pq)=f(p)f(q)f(pq)=f(p)f(q)f(pq)=f(p)f(q) ,则称 fff 是积性函 ...

  5. [数论]莫比乌斯反演1

    索引 莫比乌斯反演1 定理 莫比乌斯反演2 证明 莫比乌斯反演3 技巧 前言 本篇内容全部为定理,无证明 定义 莫比乌斯函数的符号为\(\mu\),通俗的来讲 \[ \mu(n) = \left\{ ...

  6. 算法学习-莫比乌斯反演

    写在前面 必须把更多的精力放在文化课上了, 所以这段时间的学习和数学相关的比较多, 希望可以对文化课有帮助. 莫比乌斯反演公式 g(n)=∑d|nf(d)⇒f(n)=∑d|nμ(d)g(nd) g(n ...

  7. 莫比乌斯反演入门讲解

    莫比乌斯反演实际上是一两个公式定理的运用,自认为想要掌握它的话,其中的证明还是有必要了解的.看过网上一些博客,感觉都只证明了一半,没看到有人将这个定理完全证明出来.然而我最近在正好在学习初等数论,发现 ...

  8. 狄利克雷卷积莫比乌斯反演证明

    狄利克雷卷积简介 卷积这名字听起来挺学究的,今天学了之后发现其实挺朴实hhh. 卷积: "(n)"表示到n的一个范围. 设\(f,g\)是两个数论函数(也就是说,以自然数集为定义域 ...

  9. 狄利克雷卷积_莫比乌斯反演及狄利克雷卷积

    参考文档: https://wenku.baidu.com/view/fbec9c63ba1aa8114431d9ac.html 假设$F(n)=sum_{d|n}f(d)$,那么$f(n)=sum_ ...

最新文章

  1. [Python_7] Python Socket 编程
  2. 2021辽宁省大学生程序设计竞赛(热身赛)【完结】
  3. 【干货】产品经理处理工作中突发事件的新技能Get!
  4. linux命令:until循环
  5. nssl1167-桐人的约会【最短路】
  6. 操作系统【八】文件管理
  7. iis5.0+php5.0+mysql5.0配置完全手册_IIS5.0+PHP5.0+MySQL5.0配置完全手册
  8. python天天向上的力量三天打鱼两天晒网_017 示例3-天天向上的力量-Go语言中文社区...
  9. 95-235-040-源码-task-Flink 对用户代码异常处理
  10. ICE通信之IceGrid服务(一)
  11. 电子商务概论(农)之章节课后题
  12. 结构相似度索引(SSIM)全攻略:理论+代码(PyTorch)
  13. JS实现四舍五入保留两位小数
  14. [求助]python模拟布丰投针
  15. Admob反馈利诱性质流量(无效流量)解决方案
  16. 强化学习--蒙特卡洛法
  17. java程序设计实用教程高飞pdf_普通高等教育“计算机类专业”规划教材:Java程序设计实用教程习题集 pdf epub mobi txt 下载...
  18. 第四章 06 打印星号
  19. 【Rust日报】 2019-02-23
  20. 交换机的工作原理是什么?

热门文章

  1. [imx8mpevk] pcie endpoint test
  2. 初学者制作VMOS场效应管小功放
  3. Dolphinscheduler 学习笔记
  4. java枚举类构造方法默认是,枚举类型的构造方法_Java语言程
  5. 【前端】性能优化 - WebP
  6. 极限交付:软件项目外包成功的保障
  7. 微信签到抽奖程序java源码_某宝买的微信抽奖签到墙源码,年会,学习源码--已经配置完成了,源码没问题,可直接用...
  8. 3D空间堆叠PCB设计 ——PCB线路优化项目总结
  9. 【CV/Matlab系列】基于matlab GUI的视频监控界面
  10. kumo词云使用io_Microsoft Kumo Search可以在Google上使用。 再次。