[OI] 莫比乌斯函数与莫比乌斯反演
9. 莫比乌斯函数与莫比乌斯反演
9.1 莫比乌斯函数
9.1.1 定义
设 \(\mu\) 为莫比乌斯函数,则有:
\]
直观地说,只要 \(x\) 的某个质因子出现的次数超过一次,则有 \(\mu(x)=0\),否则,\(mu(x)\) 即为 \((-1)^{t}\),其中 \(t\) 为 \(x\) 的质因数总数.
9.1.2 性质
\(No.9121\) \(\forall (i,j)=1\rightarrow\mu(i\times j)=\mu(i)\times \mu(j)\)
\]
若 \(\mu(A)=1\),则有 \(\mu(A)\times \mu(B)=1\times \mu(B)=\mu(1\times B)\),结果成立
若 \(\mu(A)=0\),则有 \(\mu(A)\times\mu(B)=0\),结果成立
否则设 \(A=\prod_{i\in prime}i,B=\prod_{j\in prime}j\ (\forall i\neq j)\),则有 \(\mu(A)\times\mu(B)=\mu(A\times B)\)
\(No.9122\)
\]
\]
\(n=0\) 时,原式等于 \(0\)
\(n=1\) 时,原式等于 \(\mu(1)=1\)
否则,设 \(n\) 的唯一分解式为 \(n=\prod_{1\le i\le k}p_{i}^{a_{i}}\),由 \(\mu(n)\neq 0\) 知 \(\forall a_{i}=1\)
显然,\(\forall d\mid n\),有 \(i\) 个质因子的数 \(d\) 有 \(C^{i}_{k}\) 种组合情况,每种组合的值为 \((-1)^{i}\),即:
\]
使用二项式定理转化:
\]
发现此时在定义域 \((k\neq 0)\) 内恒为 \(0\),易知当 \(n\gt 1\) 时,\(k_{n}\ge 1\),因此 \(\mu(n)=0\),证毕.
9.2 莫比乌斯反演
9.2.1 不完全结论
设
\]
对于
\]
考虑到
\]
根据 \(No.9122\) 得出推导式:
\]
代入有:
\]
即:
\]
将原式等价转换,可以得到
\]
其中
\]
(可以发现,在这里仅仅是将判断条件与枚举条件互相变换了,而答案是不变的)
在区间 \([1,\lfloor{\frac{n}{k}}\rfloor]\) 中,满足条件 \(d\mid i\) 的 \(i\) 有 \(\lfloor{\frac{n}{kd}}\rfloor\) 个,因此原式可化为
\]
可以证明,当 \(d\gt \min(\lfloor{\frac{n}{k}}\rfloor,\lfloor{\frac{m}{k}}\rfloor)\) 时不存在解.
9.2.2 完全结论
设两个函数 \(f,g\) 满足如下等式:
\]
则有:
\]
特别地:对于莫比乌斯反演还有如下等式:
设两个函数 \(f,g\) 满足如下等式:
\]
则有:
\]
[OI] 莫比乌斯函数与莫比乌斯反演的更多相关文章
- 莫比乌斯函数筛法 & 莫比乌斯反演
模板: int p[MAXN],pcnt=0,mu[MAXN]; bool notp[MAXN]; void shai(int n){ mu[1]=1; for(int i=2;i<=n;++i ...
- bzoj 2440 完全平方数 【莫比乌斯函数】
题目 题意:第Ki 个不是完全平方数的正整数倍的数. 对于一个数t,t以内的数里的非完全平方数倍数的个数:num=1的倍数的数量−一个质数平方数(9,25,49...)的倍数的数量+两个质数的积平方数 ...
- 51nod 1244 莫比乌斯函数之和 【杜教筛】
51nod 1244 莫比乌斯函数之和 莫比乌斯函数,由德国数学家和天文学家莫比乌斯提出.梅滕斯(Mertens)首先使用μ(n)(miu(n))作为莫比乌斯函数的记号.具体定义如下: 如果一个数包含 ...
- BZOJ 1101 莫比乌斯函数+分块
思路: 题目中的gcd(x,y)=d (x<=a,y<=b)可以转化成 求:gcd(x,y)=1 (1<=x<=a/d 1<=y<=b/d) 设 G(x,y)表示x ...
- 莫比乌斯函数&莫比乌斯反演
莫比乌斯函数:http://wenku.baidu.com/view/fbec9c63ba1aa8114431d9ac.html Orz PoPoQQQ
- hdu 1965 (莫比乌斯函数 莫比乌斯反演)
GCD Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submiss ...
- 【bzoj4804】欧拉心算 莫比乌斯反演+莫比乌斯函数性质+线性筛
Description 给出一个数字N 求\(\sum_{i=1}^{n}\sum_{j=1}^{n}\varphi(gcd(i,j))\) Input 第一行为一个正整数T,表示数据组数. 接下来T ...
- 【bzoj3601】一个人的数论 莫比乌斯反演+莫比乌斯函数性质+高斯消元
Description Sol 这题好难啊QAQ 反正不看题解我对自然数幂求和那里是一点思路都没有qwq 先推出一个可做一点的式子: \(f(n)=\sum_{k=1}^{n}[(n,k)=1]k^d ...
- 莫比乌斯函数 && HDU-1695
莫比乌斯函数定义: $$\mu(d)=\begin{cases}1 &\text{d = 1}\\(-1)^r &\text{$d=p_1p_2...p_r,其中p_i为不同的素数$} ...
- BZOJ 2440 莫比乌斯函数+容斥+二分
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 5473 Solved: 2679[Submit][Sta ...
随机推荐
- 靶机: EvilBox---One
靶机: EvilBox---One 准备工作 靶机地址: https://download.vulnhub.com/evilbox/EvilBox---One.ova MD5 校验:c3a65197b ...
- 初读nginx.conf
学成在线里教的部署门户,记录一下NGINX的部分配置是如何使用的 location块用于定义如何处理不同URL请求 1.代理: server { listen 80; server_name file ...
- Jmeter函数助手36-P
P函数用于获取jmeter属性值.类似property函数 属性名称:填入jmeter的属性名称 默认值:缺省值,当获取属性值为空时则返回该值 1.填入属性名称获取属性值${__P(language, ...
- tar命令备份压缩7天生产日志
[root@localhost logs]# cat tar_7day.sh #!/bin/bash #压缩日期[当天的前一天] todayStamp_1=`date -d "-1 day& ...
- 【SpringMVC】10 对Ajax的应用
编写一个AjaxController package cn.dai.controller; import org.springframework.web.bind.annotation.GetMapp ...
- Streamlit运行出现ModuleNotFoundError: No module named ‘altair.vegalite.v4‘ —— ModuleNotFoundError: No module named 'altair.vegalite.v4'
参考: https://blog.csdn.net/ikun_King/article/details/131852167 解决方法: pip install altair=4.2.2
- xshell打开vim后颜色异常——xshell连接ubuntu打开vim后界面覆盖一层绿色
参考原文: https://blog.csdn.net/Blank_Shen/article/details/106527312 =================================== ...
- 花店橱窗(线性DP)
线性DP--花店橱窗 谨以此题解献给线性dp最后一道题 题目大致 Description xq和他的老婆xz最近开了一家花店,他们准备把店里最好看的花都摆在橱窗里.但是他们有很多花瓶,每个花瓶都具有各 ...
- 利用XtraBackup对MGR集群进行扩容
运行了一段时间以后MGR集群,需要扩容节点,这是一个常见的需求.很多时候我们都喜欢用mysqldump工具来进行,因为这个工具有一个很好用的参数叫做master-data以及single-transa ...
- zabbix4.0配置短信报警
1.准备工作 #访问短信网址:172.16.98.1,网线插LAN口 #账号&密码:admin 安装ubuntu系统模拟http请求工具(命令行模式) # apt-get install ht ...