莫比乌斯反演

反演初步和其他

↘锁芯ラ 提交于 2019-11-28 17:42:48
参考了刘汝佳老师《算法艺术与信息学竞赛》。 离散变换与反演 有些时候,我们所求解问题的答案可以表示成一个类似前缀和的形式: \(f(x)=\sum a_ig(i)\) 。 更多时候, \(g(x)\) 是我们所需的答案,不易求出;而 \(f(x)\) 则可以通过另一种更好的方式求出来。如何根据 \(f(x)=\sum a_ig(i)\) 来推导出 \(g(x)\) 的解析式呢?这个反推的过程,我们可以叫做 反演 。我们用一个例子来具体说明。 我们已经用容斥原理得出了错排问题的答案,即在 \(n\) 个元素中的所有排列中,第 \(i(1 \leq i \leq n)\) 个元素不是 \(i\) 的排列数目 \(D_n\) 。现在我们来尝试用反演的方法。 我们反过来求:有且仅有 \(k\) 个元素处在正确位置的方案数为 \(\dbinom{n}{k}D_{n-k}\) ,即剩余的 \(n-k\) 个元素不允许再排在原来的位置。 易知,所有排列里,一定存在 \(0,1,\cdots,n\) 个元素恰好排在其位置的方案。我们把他们加起来得到: \[ P_n=n!=\sum_{k=0}^n \dbinom{n}{k}D_{n-k} \tag{*} \] 反演(解出 \(D_n\) )得到: \[ D_n=n!\sum_{r=0}^n{(-1)^r \frac{1}{r!}} \tag{**}

莫比乌斯反演入门

梦想与她 提交于 2019-11-28 17:42:28
      转载自----- http://blog.csdn.net/qw4990/article/details/14055183 这个文章主要讲一下ACM中1个常用的莫比乌斯反演公式,看到很多博客上面公式是有,但是都没证明,《组合数学》上的证明又没看懂, 就自己想了种证明方法,觉得比《组合数学》的证明简单些,就写一下,希望对初学莫比乌斯反演的同学有帮助。 PS:下面公式出现的sigma是累加,另外建议大家看的时候 把公式在纸上写出来! 一:什么是莫比乌斯反演 简单点的说,就是先给出一个函数 F(n) ,然后再由 F(n)定义一个新函数 G(n) 其中 G(n) = sigma(F(d)) (其中d被“包含”于n) 然后 现在我们不知道 F(n)的值 , 却知道 G(n), 接着我们就可以通过 反演由G(n)反向得到F(n) 什么叫 (其中d被“包含”于n) ?以及怎么理解反演? 通过下面的几个例子说明 例1: 我们直接定义 G(n)=sigma(F(i)) (1<=i<=n) {这里的每个F(i),相对于G(n)实际上就是一种包含关系了!!} 然后我们现在已经知道 G(n)=n*(n+1)/2; 接下来 我们要通过 G(n)反向得到F(n) 的过程,就是反演 当然,这个问题很简单,很容易都可以看出来 F(n)=n ~~ ----------------------------

莫比乌斯反演

久未见 提交于 2019-11-28 17:42:16
莫比乌斯反演 (PS:在评论区中众多dalao的催促下,我认真的写了三天三夜写完了这篇 杜教筛 ,保证是精品!) 前言 (这大概是我第一次写学习笔记吧OvO) 可能每一个刚开始接触莫比乌斯反演的OIer,起初都会厌恶这个神奇的东西。 (我也一样233) 每一个人厌恶的原因有许多,可能是这个烦人的式子,也可能仅仅只是因为不理解 \(\mu\) 函数而感到不爽。当然,莫比乌斯反演有一个小小的预备知识: 整除分块 那么我们先从莫比乌斯反演中最基础的莫比乌斯函数 \(\mu\) 开始说起: 莫比乌斯函数 首先,我们可以先明确一点,莫比乌斯函数并不是什么很高大上的东西,它其实只是一个由容斥系数所构成的函数。 \(\mu(d)\) 的定义是: 当 \(d=1\) 时, \(\mu(d)=1\) ; 当 \(d=\Pi_{i=1}^{k}p_i\) 且 \(p_i\) 为互异素数时, \(\mu(d)=(-1)^k\) 。(说直白点,就是 \(d\) 分解质因数后,没有幂次大于平方的质因子,此时函数值根据分解的个数决定); 只要当 \(d\) 含有任何质因子的幂次大于等于2,则函数值为0. 当然,莫比乌斯函数也有很多有趣的性质: 对于任意正整数 \(n\) , \(\sum_{d|n}\mu(d)=[n=1]\) 。( \([n=1]\) 表示只有当 \(n=1\) 成立时,返回值为 \(1\)

莫比乌斯反演入门

拥有回忆 提交于 2019-11-27 13:59:00
。。。。您点进来我的博客,想必是信任我这个新初一大菜鸡能讲清楚吧!可是我来写这篇博客,完全是因为考试的时候看到数学题推到类似∑∑【gcd(i,j)==1】的式子就推不下去了。其实莫比乌斯反演就是求这样的东西的。 对于这样的式子,我经常说一句话:当年bmh201708普及组的时候提前30分钟出考场吊打提高组,就是冲着数学题去AK的 对于这样的人,我们用3个字形容他:金艺轲! 首先声明一点:我不会markdown 我zbs用莫比乌斯定理不是问题:【n==1】=∑mu(d)(d|n) 如果这个式子带进去,将绝杀。而且可以带 单走一个交换枚举项,不用说,他死定了:∑∑[gcd(i,j)==1] =∑∑∑mu(d)(d|gcd(i,j))=∑mu(d)∑∑(d|i&&d|j) (到这里,有没有发现gcd(i,j)的所有因数都是i和j的因数?) 给莫比乌斯倒一杯卡布奇诺! 原式=∑mu(d)[n/d]*[m/d] 后面整数分块处理O(sqrt(n)),前面莫比乌斯函数因为是积性函数,所以直接上欧拉筛。 你卡掉我试试!你能把我zbs卡掉?卡掉我当场吃屎! https://www.luogu.org/record/22884149 来源: https://www.cnblogs.com/zbsakioi/p/11366682.html

莫比乌斯反演学习笔记

邮差的信 提交于 2019-11-26 17:42:00
昨天知道了莫比乌斯函数,今天就来学学莫比乌斯反演: 莫比乌斯反演: 若 \[F(n)=\sum_{x\mid n}f(x)\] 可以得到一个结论: \[f(n)=\sum_{x\mid n}\mu (x)F(\left \lfloor \frac{n}{d} \right \rfloor)\] 若 \[F(n)=\sum_{n\mid x}f(x)\] 又可以得到一个结论: \[f(n)=\sum_{n\mid x}\mu (\frac{x}{n})F(x)\] 证明我们可以用 狄利克雷乘积 来证明。 算法应用: 此处感谢我们学校的教练——symbol,谢谢您的课件,这给了我很大帮助。 终于可以学习杜教筛了,Yes! 来源: https://www.cnblogs.com/zhouyifei/p/11329483.html

狄利克雷卷积和莫比乌斯反演

[亡魂溺海] 提交于 2019-11-26 14:26:29
积性函数(前置知识) 积性函数定义: 函数 \(f(x)\) 满足 \(gcd(a, b) = 1\) 时, \(f(ab) = f(a)f(b)\) 则 \(f(x)\) 为积性函数 常见积性函数 (具体证明可以百度=w=) 欧拉函数 \(ϕ(n)\) \(ϕ(n) = n ∗∏(pi − 1)/pi\) 莫比乌斯函数 \(µ(n)\) 当 \(n\) 有平方因子 如 \(n=∏_{i=1}^{t}pi^{ci}(表示n有t个互不相同质因子pi,每个pi的次数是ci)\) 当某 \(ci>=2(即有平方因子)\) ,则 \(µ(n) = 0\) 否则,若 \(n\) 为 \(k\) 个不同质数的乘积, \(µ(n) = (−1)^k\) 除数函数 \(σ_k(n)\) 表示所有正因子的 \(k\) 次幂和 \(σ_0(n) = d(n)\) 表示正因子的个数 \(σ_1(n) = σ(n)\) 表示正因子的和 完全积性函数 幂函数 \(id_k(n) = n^k\) \(id_0(n) = 1(n) =1\) \(id_1(n) = id(n) = n\) 单位函数 \(ϵ(n) = [n = 1]\) : \(即ϵ(n)仅当n=1时值为1,其它都为0\) 狄利克雷卷积 Dirichlet 卷积 对两个数论函数 \(f, g\) ,定义其 Dirichlet 卷积为新函数 \(f