莫比乌斯反演的证明(非狄利克雷卷积法)
首先给大家介绍一下莫比乌斯函数吧,其实这个函数挺好理解的,只是一个容斥系数
μ(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的倍数的时候也满足莫比乌斯反演定理,推导过程如下
这个μ是莫比乌斯函数的意思
莫比乌斯反演的证明(非狄利克雷卷积法)相关推荐
- matlab狄利克雷函数,数论入门1——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛...
数论入门1 一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数$\mu(n)$为: 当n有平方因子时,$\mu(n)=0$. 当n没有平方因子时,$\mu(n)=(-1)^{\omega(n)} ...
- 狄利克雷卷积与莫比乌斯函数
狄利克雷卷积与莫比乌斯函数 数论函数与积性函数 狄利克雷卷积 莫比乌斯函数 相关性质 基础问题与公式 若想使用莫比乌斯反演,熟练掌握狄利克雷卷积包括定义.记号以及相关的性质.证明等是非常有好处的. 数 ...
- 数学--数论--莫比乌斯反演
一.莫比乌斯反演涉及知识 1.莫比乌斯函数 2.莫比乌斯的线性筛法 3.狄利克雷卷积 4.莫比乌斯反演详解 5.整除法分块 6.杜教筛 二.μ 莫比乌斯函数定义 μ(n)={1n=1(−1)kn= P ...
- 从积性函数到莫比乌斯反演
积性函数 积性函数:对于数论函数 fff ,若任意互质的 p,qp,\ qp, q 都有 f(pq)=f(p)f(q)f(pq)=f(p)f(q)f(pq)=f(p)f(q) ,则称 fff 是积性函 ...
- [数论]莫比乌斯反演1
索引 莫比乌斯反演1 定理 莫比乌斯反演2 证明 莫比乌斯反演3 技巧 前言 本篇内容全部为定理,无证明 定义 莫比乌斯函数的符号为\(\mu\),通俗的来讲 \[ \mu(n) = \left\{ ...
- 算法学习-莫比乌斯反演
写在前面 必须把更多的精力放在文化课上了, 所以这段时间的学习和数学相关的比较多, 希望可以对文化课有帮助. 莫比乌斯反演公式 g(n)=∑d|nf(d)⇒f(n)=∑d|nμ(d)g(nd) g(n ...
- 莫比乌斯反演入门讲解
莫比乌斯反演实际上是一两个公式定理的运用,自认为想要掌握它的话,其中的证明还是有必要了解的.看过网上一些博客,感觉都只证明了一半,没看到有人将这个定理完全证明出来.然而我最近在正好在学习初等数论,发现 ...
- 狄利克雷卷积莫比乌斯反演证明
狄利克雷卷积简介 卷积这名字听起来挺学究的,今天学了之后发现其实挺朴实hhh. 卷积: "(n)"表示到n的一个范围. 设\(f,g\)是两个数论函数(也就是说,以自然数集为定义域 ...
- 狄利克雷卷积_莫比乌斯反演及狄利克雷卷积
参考文档: https://wenku.baidu.com/view/fbec9c63ba1aa8114431d9ac.html 假设$F(n)=sum_{d|n}f(d)$,那么$f(n)=sum_ ...
最新文章
- [Python_7] Python Socket 编程
- 2021辽宁省大学生程序设计竞赛(热身赛)【完结】
- 【干货】产品经理处理工作中突发事件的新技能Get!
- linux命令:until循环
- nssl1167-桐人的约会【最短路】
- 操作系统【八】文件管理
- iis5.0+php5.0+mysql5.0配置完全手册_IIS5.0+PHP5.0+MySQL5.0配置完全手册
- python天天向上的力量三天打鱼两天晒网_017 示例3-天天向上的力量-Go语言中文社区...
- 95-235-040-源码-task-Flink 对用户代码异常处理
- ICE通信之IceGrid服务(一)
- 电子商务概论(农)之章节课后题
- 结构相似度索引(SSIM)全攻略:理论+代码(PyTorch)
- JS实现四舍五入保留两位小数
- [求助]python模拟布丰投针
- Admob反馈利诱性质流量(无效流量)解决方案
- 强化学习--蒙特卡洛法
- java程序设计实用教程高飞pdf_普通高等教育“计算机类专业”规划教材:Java程序设计实用教程习题集 pdf epub mobi txt 下载...
- 第四章 06 打印星号
- 【Rust日报】 2019-02-23
- 交换机的工作原理是什么?
热门文章
- [imx8mpevk] pcie endpoint test
- 初学者制作VMOS场效应管小功放
- Dolphinscheduler 学习笔记
- java枚举类构造方法默认是,枚举类型的构造方法_Java语言程
- 【前端】性能优化 - WebP
- 极限交付:软件项目外包成功的保障
- 微信签到抽奖程序java源码_某宝买的微信抽奖签到墙源码,年会,学习源码--已经配置完成了,源码没问题,可直接用...
- 3D空间堆叠PCB设计 ——PCB线路优化项目总结
- 【CV/Matlab系列】基于matlab GUI的视频监控界面
- kumo词云使用io_Microsoft Kumo Search可以在Google上使用。 再次。