吐槽

额其实这个东西的话。。好像缠着机房里面的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. PHP读取大文本文件并处理数据的思路

    //处理文件 $file = fopen($filename, "r") or exit("Unable to open file!"); $total_lin ...

  2. memcached 的实践操作

    memcached安装和使用   yum install -y libevent  memcached  libmemcached   启动命令:   /etc/init.d/memcached st ...

  3. 4、flask之分页插件的使用、添加后保留原url搜索条件、单例模式

    本篇导航: flask实现分页 添加后保留原url搜索条件 单例模式 一.flask实现分页 1.django项目中写过的分页组件 from urllib.parse import urlencode ...

  4. CentOS6.9安装mysql5.7

    1.查看系统版本 uname -a 2.卸载old mysql rpm -qa|grep -i mysql rpm -ev MySQL-client-5.5.25a-1.rhel5 如果提示依赖错误, ...

  5. 报错信息 The jsp:param action must not be used outside the jsp:include, jsp:forward, or jsp:params elements 的原因及解决办法

    如果你的代码是这样的话就会报错 <jsp:forward page="02.jsp"></jsp:forward> <jsp:param value= ...

  6. asp.net core 一 Centos 环境部署

             .netcore的运行环境,创建asp.net core 项目    CentOS 7 ,dotnet-sdk-2.0.0-2.0.0-1.x86_64 直接在liunx创建项目并运 ...

  7. hihoCoder 1523 数组重排2 贪心

    题意:给定一个1-N的排列A1, A2, - AN,每次操作小Hi可以选择一个数,把它放到数组的最左边. 请计算小Hi最少进行几次操作就能使得新数组是递增排列的. 思路:最后的序列是递增的,那么必定满 ...

  8. hihoCoder 403 Forbidden 字典树

    题意:给定个规则,个ip,问这些ip是否能和某个规则匹配,如果有多个规则,则匹配第一个.如果没能匹配成功,则认为是"allow",否则根据规则决定是"allow" ...

  9. react——一个todolist的demo

    代码如下: function ToDoListHeader(props) { return <h1 className={props.className}>ToDoList</h1& ...

  10. [php] in_array 判断问题(坑)

    <?php $arr = array("Linux"); if (in_array(0, $arr)) { echo "match"; } ?> 执 ...