吐槽

额其实这个东西的话。。好像缠着机房里面的dalao们给我讲过好多遍了然后。。

拖到现在才搞懂也是服了qwq(可能有个猪脑子)

感觉就是主要几条式子然后疯狂换元换着换着就化简运算了?

草稿纸杀手qwq


莫比乌斯反演公式

$F(n)$和f(n)是定义在非负整数集合上面的两个函数,并且满足条件$F(n) = \sum\limits_{d\mid n}f(d)$,那么
$$
f(n) = \sum\limits_{d\mid n}\mu(d)F(\frac{n}{d})
$$
这条式子还有另一种描述

$F(n)$和f(n)$足条件F(n) = \sum\limits_{n\mid d}f(d)$,那么
$$
f(n)=\sum\limits_{n\mid d}\mu(\frac{d}{n})F(d)
$$

上面的公式有个$\mu$函数,定义如下:

1. 若 $d=1$,那么$\mu(d) =1$
2. 若$d=\prod\limits_{i=1}^{k}p_i$,且$p_i$均为互异素数,那么$\mu(d) =(-1)^k$
3. 其他情况($d$有平方因子)$\mu(d)=0$

$\mu$的常见性质

对于任意正整数$n$满足
$$
\sum\limits_{d\mid n}\mu(d) = [n=1]
$$

$$
\sum\limits_{d\mid n}\frac{\mu(d)}{d}=\frac{\phi(n)}{n}
$$

证明

啊。。是证明莫比乌斯反演公式啦。。不是上面两条qwq
$$
\sum\limits_{d\mid n}\mu(d)F(\frac{n}{d})=\sum\limits_{d\mid n}\mu(d)\sum\limits_{d'\mid \frac{n}{d}}f(d')=\sum\limits_{d'\mid n}f(d')\sum\limits_{d\mid \frac{n}{d'}}\mu(d)=f(n)
$$

应用

目前做到的几题都是。。各种换元然后优化式子?

总之大概就是上面四条式子(性质两条+反演两条)+各种玄学换来换去,最后好像。。目前做的几题都是化成了一个带有$g(T) = \sum\limits_{d\mid T}f(d)\mu(\frac{T}{d})$的式子,然后就想办法把$g(x)$筛出来

最后的求解基本上是要用到一个(类似)分块的方法用前缀和在根号的时间内把式子里面的其他一些奇奇怪怪的部分求出来

主要题做的也不多qwq大概就先这样吧qwq

(所以说过了这么久才更博肯定不是因为懒嗯)

【learning】莫比乌斯反演的更多相关文章

  1. 【Learning】 莫比乌斯反演

    莫比乌斯反演 ​ 对于两个定义域为非负整数的函数\(F(n)\)和\(f(n)\) ​ 若满足:\(F(n)=\sum\limits_{d|n}f(d)\),则反演得到\(f(n)=\sum\limi ...

  2. hdu1695 GCD(莫比乌斯反演)

    题意:求(1,b)区间和(1,d)区间里面gcd(x, y) = k的数的对数(1<=x<=b , 1<= y <= d). 知识点: 莫比乌斯反演/*12*/ 线性筛求莫比乌 ...

  3. BZOJ 2154: Crash的数字表格 [莫比乌斯反演]

    2154: Crash的数字表格 Time Limit: 20 Sec  Memory Limit: 259 MBSubmit: 2924  Solved: 1091[Submit][Status][ ...

  4. BZOJ2301: [HAOI2011]Problem b[莫比乌斯反演 容斥原理]【学习笔记】

    2301: [HAOI2011]Problem b Time Limit: 50 Sec  Memory Limit: 256 MBSubmit: 4032  Solved: 1817[Submit] ...

  5. Bzoj2154 Crash的数字表格 乘法逆元+莫比乌斯反演(TLE)

    题意:求sigma{lcm(i,j)},1<=i<=n,1<=j<=m 不妨令n<=m 首先把lcm(i,j)转成i*j/gcd(i,j) 正解不会...总之最后化出来的 ...

  6. 莫比乌斯函数筛法 & 莫比乌斯反演

    模板: int p[MAXN],pcnt=0,mu[MAXN]; bool notp[MAXN]; void shai(int n){ mu[1]=1; for(int i=2;i<=n;++i ...

  7. 【BZOJ-2440】完全平方数 容斥原理 + 线性筛莫比乌斯反演函数 + 二分判定

    2440: [中山市选2011]完全平方数 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2371  Solved: 1143[Submit][Sta ...

  8. POI2007_zap 莫比乌斯反演

    题意:http://hzwer.com/4205.html 同hdu1695 #include <iostream> #include <cstring> #include & ...

  9. hdu.5212.Code(莫比乌斯反演 && 埃氏筛)

    Code Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submi ...

  10. CSU 1325 莫比乌斯反演

    题目大意: 一.有多少个有序数对(x,y)满足1<=x<=A,1<=y<=B,并且gcd(x,y)为p的一个约数: 二.有多少个有序数对(x,y)满足1<=x<=A ...

随机推荐

  1. 如何为MySQL服务器和客户机启用SSL

    本文来自我的github pages博客http://galengao.github.io/ 即www.gaohuirong.cn 摘要: mysql5.7后有ssl新特性 自己搭建mysql ent ...

  2. HTTP协议学习笔记

    一.什么是HTTP协议 HTTP协议是指计算机通信网络中两台计算机之间进行通信所必须共同遵守的规定或规则,超文本传输协议(HTTP)是一种通信协议,它允许将超文本标记语言(HTML)文档从Web服务器 ...

  3. Java经典编程题50道之十

    一球从100米高度自由落下,每次落地后反跳回原高度的一半:再落下……求它在第10次落地时,共经过多少米?第10次反弹多高? public class Example10 {    public sta ...

  4. 通过核心概念了解webpack工作机制

    webpack 是一个现代 JavaScript 应用程序的静态模块打包器(module bundler).当 webpack 处理应用程序时,它会递归地构建一个依赖关系图(dependency gr ...

  5. Servlet中文乱码问题解决办法

    首先对于源jsp网站和servlet里面的字符集要一样,一般支持中文的字符集为UTF-8最好采用这个字符集(除此之外还有gb2312); 对于源jsp文件的代码中需要设置 设置你的page里面的字符集 ...

  6. C语言_了解一下C语言中的四种存储类别

    C语言是一门通用计算机编程语言,应用广泛.C语言的设计目标是提供一种能以简易的方式编译.处理低级存储器.产生少量的机器码以及不需要任何运行环境支持便能运行的编程语言. C语言中的四种存储类别:auto ...

  7. (MonoGame从入门到放弃-3)-放弃MonoGame

    又一段时间过去了,这一章没内容了.我真的已经放弃MonoGame的学习了,MonoGame用起来感觉就是在自己实现2d游戏引擎一样,好多现代游戏引擎有的内容都没有...,我只是想做游戏,而不是给引擎添 ...

  8. UVA - 1371 Period 二分+dp

    思路:设字符串x的长度为n,y的长度为m,那么答案一定在[0, m]之间,那么可以二分求答案. d(i, j)表示第一个串前i个字符至少需要经过多少次才能的到第二个串的前j个字符,转移方程d(i, j ...

  9. Enable multi-tenancy on ironic

    Multi-tenancy 是openstack ironic从Ocata版本开始支持的新特性,通过network-generic-switch插件控制交换机,Ironic可以实现在不同租户间机网络隔 ...

  10. 源码编译安装Apache-附一键部署脚本

    1.进入apache官网https://www.apache.org/,点击Download 2.如图选择 3.选择httpd 4.下载两个包,2.2为CentOS6使用,2.4为CentOS7使用 ...