莫比乌斯反演入门
莫比乌斯反演
前置知识:容斥 整除分块
莫比乌斯函数
设正整数 \(N\) 按照算术基本定理分解质因数为 \(N=\Pi_{i=1}^m p_i^{c_i}\) ,定义函数:
我们称 \(\mu(N)\) 为莫比乌斯函数
通俗地讲,当 \(N\) 包含相等的质因子时, \(\mu(N)=0\)
当 \(N\) 的所有质因子各不相同时,若 \(N\) 有偶数个质因子,\(\mu(N)=1\)
若 \(N\) 有奇数个质因子,\(\mu(N)=-1\)
它有一些性质:
-
对于任意正整数 \(n\) ,\(\sum_{d|n}\mu(d)=[n=1]\)
( \([n=1]\) 表示只有当 \(n=1\) 成立时,返回值为 \(1\);否则,值为 \(0\))
简单地证明一下:
\(n=1\) 时显然成立
对于 \(n\neq1\) 的情况:
证明一:对于一个 \(p_i\) 强制选或不选,其他随便,会产生 奇数/偶数 个质因子两种情况,抵消为 \(0\)
证明二:若 \(d\) 中含 \(p_i^{a_i},a_i>=2\) ,即 \(\mu(d)=0\) 忽略
那么, \(\mu\) 的取值仅仅与质因数个数有关
\((i)\) \(n\) 有奇数个(\(x\)个)质因数
\(\sum_{d|n}\mu(d)=-1\times(C_x^1+C_x^3+C_x^5...+C_x^x)+(C_x^2+C_x^4+...+C_x^{x-1})\)
又有:\(C_x^y=C_x^{x-y}(y<x)\) 全部抵消
\((ii)\) \(n\) 有偶数个质因数:同理。
- 对于任意正整数 \(n\),\(\sum_{d|n}\frac{μ(d)}d=\frac{ϕ(n)}n\)
从略。(不会)
程序实现并不难,我们可以在线性筛素数的程序上略作修改,便可以筛出 \(μ\) 函数
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];
}
}
}
莫比乌斯反演
定理:\(F(n)\) 和 \(f(n)\) 是定义在非负整数集合上的两个函数,并且满足条件:\(F(n)=\sum_{d|n}f(d)\)
那么存在一个结论:\(f(n)=\sum_{d|n}\mu(d)F(\lfloor \frac nd \rfloor)\)
即莫比乌斯反演定理
我们可以通过定义证明它:
它的另外一种形式是:当 \(F(n)\) 和 \(f(n)\) 满足:\(F(n)=\sum_{n|d}f(d)\)
可以推出:\(f(n)=\sum_{n|d}\mu(\frac dn)F(d)\)
感觉这个式子,可能在莫比乌斯反演中更加好用。
题目
YY的gcd
[POI2007]ZAP-Queries
[SDOI2015]约数个数和
[HAOI2011]Problem b
[洛谷P1829 国家集训队]Crash的数字表格
原文地址:https://www.cnblogs.com/wyb-sen/p/15027818.html
- 各大公司移动端页面 - 导航的实现
- JavaScript 运行机制之执行顺序详解
- Math对象面试题目
- Highcharts AJAX JSON JQuery 实现动态数据交互显示图表 柱形图
- 企业支付宝账号开发接口教程--JAVA-UTF-8(实际操作------SpringMVC+JSP)
- 用pandas 进行投资分析
- 【专业技术】android 应用程序如何获取root权限
- Nginx+Keepalived(双机热备)搭建高可用负载均衡环境(HA)
- SpringMVC+MongoDB+Maven整合(微信回调Oauth授权)
- ZeroClipboard实现多个浏览器兼容的复制文本到剪贴板的功能
- Shiro 权限框架使用总结
- Apriori算法介绍(Python实现)
- linux学习第六十二篇:添加自定义监控项目,配置邮件告警,测试告警,不发邮件的问题处理
- Entity Framework Core 2.0 入门
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- Reactive-MongoDB异步Java Driver解读
- 解密Go协程的栈内存管理
- 深入浅出mongodb之实战
- 想成为可视化高手?这篇合集就够了 | Vue
- 谈谈Vue开发过程中用到的插件
- 快速入门使用tikz绘制深度学习网络图
- why哥被一道基础面试题给干懵了,一气之下写出万字长文。
- 结构与算法(03):单向链表和双向链表
- 有赞crash平台符号化实践
- 什么是时间分片(Time Slicing)?
- 逐行分析鸿蒙系统的 JavaScript 框架
- 48张小图带你领略Flex 布局之美
- 怎样设计一个 JavaScript 插件系统
- 一道 React 面试题:在浏览器、组件和元素中都渲染了些什么?
- 「新手入门福利」一张脑图带你掌握Git命令