首页 > 其他分享 >【Coel.学习笔记】莫比乌斯反演

【Coel.学习笔记】莫比乌斯反演

时间:2022-10-18 19:55:15浏览次数:69  
标签:乌斯 sum Coel mu 反演 text 莫比

闲话

记得在差不多一年前写扩展欧拉定理的时候我提了一句

这周终于把古代猪文搞定了,数论这块的内容就只剩个博弈论了
别提莫比乌斯反演之类的东西,我不想搞

甚至刚开始写的时候还笔误把莫反写成了莫队……

转眼一年过去了,来填上这个大坑吧!

前言

莫比乌斯反演是一个基于莫比乌斯函数 \(\mu(n)\) 的反演算法。对于某些求和函数,如果直接求解的时间复杂度很高,而倍数或因数之和求解较为容易,则可以通过莫比乌斯反演简化运算,提高效率。

前置知识

这些前置知识可以算到莫比乌斯反演的一部分,不过也可以作为一个独立的内容(在进阶指南里其实就是这样划分的)。

莫比乌斯函数的定义

记正整数 \(x\) 的质因数分解结果为 \(x=p_1^{\alpha_1}p_2^{\alpha_2}...p_k^{\alpha_n}\)(\(p_i\) 为质数且 \(\alpha_i\geq 1\)),则

\[\mu(x)=\begin{cases} 0 & \exists i,\alpha_i>1\text{(即含有平方因子)} \\ (-1)^k & \forall i,\alpha_i=1\text{ (即不含有平方因子) } \end{cases} \]

特殊地,\(\mu(1)=1\)。

莫比乌斯函数的性质

记 \(n\) 的约数的莫比乌斯函数之和 \(S_n = \sum_{d | n}\mu(d)\),则

\[S_n=\begin{cases}1 & n=1 \\ 0 & \text{otherwise.}\end{cases} \]

\(n=1\) 显然正确,其他情况证明利用二项式定理容易得到,这里略去。

顺带一提,这个性质可以说明 \(\mu\) 在狄利克雷生成函数中为黎曼函数 \(\zeta\) 的逆元,不过对于反演没啥用(

反演定理

若对于函数 \(F(n),f(n)\),\(F(n)=\sum_{d|n}f(d)\),则 $$f(n)=\sum_{d|n}\mu(d)F(\dfrac{n}{d})$$

下面是证明。我们证明得到的式子左边等于右边即可:

\[\begin{aligned}\sum_{d|n}\mu(d)F(\dfrac{n}{d}) &= \sum_{d|n}\mu(d)\sum_{i|\frac{n}{d}}f(i) & \text{【代入函数 F】}\\ &=\sum_{i|n}f(i)\sum_{d|\frac{n}{i}}\mu(d)& \text{【进行数论变换】}\\ &=f(n)& \text{【利用前面提到的性质】} \end{aligned} \]

\(\text{Q.E.D.}\)


另一种反演定理的形式为:若 \(F(n)=\sum_{n|d}f(d)\),则 \(f(n)=\sum_{n|d}\mu(\dfrac{d}{n})F(d)\)。

类比上一个的证明过程也不难得证:

\[\begin{aligned}\sum_{n|d}\mu(\dfrac{d}{n})F(d) &=\sum_{n|d}\mu(\dfrac{d}{n})\sum_{d|i}f(i) \\ &=\sum_{n|i}f(i)\sum_{d^\prime|\frac{i}{n}}\mu(d^\prime)\\ &= \sum_{n|i}f(i)S(\dfrac{i}{n})\\ &=f(n)\end{aligned} \]

\(\text{Q.E.D.}\)

利用这两个定理,我们可以把从 \(f(n)\) 到 \(F(n)\) 翻转过来求解,这个过程就叫做莫比乌斯反演

例题

待补充……

标签:乌斯,sum,Coel,mu,反演,text,莫比
From: https://www.cnblogs.com/Coel-Flannette/p/16803866.html

相关文章

  • 莫比乌斯
    莫比乌斯函数定义\[\mu(n)=\begin{cases}1&n=1\\0&n含有平方因子\\(-1)^k&k为本质不同的质因子个数\end{cases}\]性质莫比乌斯函数是积性函数,且......
  • 算法学习笔记(数学):数论分块 + 容斥原理 + 莫比乌斯函数
    算法学习笔记(数学):数论分块+容斥原理+莫比乌斯函数这篇文章主要是要讲一道题目(链接在这里)以及梳理一下数论分块,莫比乌斯函数,容斥原理这些知识。先介绍下知识点吧qwq......
  • 【Coel.学习笔记】基环树动态规划
    引入基环树(又称环套树)是一种特殊的图,在原有的树形结构上添加一条边,就会形成一个环,看起来就像从环延伸出树。特别地,对于有向图而言,环上点所连接的边指向环外为外向树,反之为......
  • luogu P5518 幽灵乐团 / 莫比乌斯反演基础练习题
    重新学习了一下莫比乌斯反演,再来写一次这道题。Part1\[\prod_{i=1}^A\prod_{j=1}^B\prod_{k=1}^C\dfrac{\text{lcm}(i,j)}{\gcd(i,k)}\]\[=\prod_{i=1}^A\prod_{j=1}^B......
  • AGC038C LCMs 详解(莫比乌斯反演好题)
    ProblemAGC038C给定一个长为\(n\)的序列\(A_1,A_2,\cdots,A_n\),求\(\sum_{i=1}^{n}{\sum_{j=i+1}^{n}{lcm(A_i,A_j)}}\bmod998244353\)\(n\leq2\times10^5,A_i......
  • 【Coel.学习笔记】特殊的图 - 仙人掌与圆方树
    你是什么仙人?引入仙人掌是一种特殊的无向图,它的任意一条边至多只出现在一条简单回路(每个点只出现一次的回路是简单回路,特殊地,自环不算简单回路)。这里借用一下[SHOI2006......
  • 【Coel.学习笔记】随机化算法:模拟退火与爬山法
    简介模拟退火(\(\text{SimulateAnneal}\))和爬山法是随机化算法,二者的原理都在于通过随机生成答案并检查,把答案逐步缩小在一个可行的区间,尽可能地靠近正确答案。在考场......
  • 容斥与反演技巧(一)
    二项式反演我们直接上式子好了一般有两种形式,第一种是\[f(n)=\sum_{i=0}^n\binom{n}{i}g(i)\iffg(n)=\sum_{i=0}^n(-1)^{n-i}\binom{n}{i}f(i)\]第二种是\[f(n)=\sum......
  • 【Coel.解题报告】【没事找事】CSP-S2 真题解析
    昨天刚考完CSP-S1,反正没什么想做的(最近好颓废…),来复盘一下。本次比赛评价(转载):CSP-S1是由CCF自主研发的一款全新开放世界冒险游戏。游戏发生在一个被称作「基数排序......
  • 【瞎口胡】单位根反演
    单位根反演是用单位根来表示整除关系的东西。定义式\[\left[k\midn\right]=\dfrac{1}{k}\sum\limits_{i=0}^{k-1}\omega_k^{in}\]如果\(k\midn\),那么\(\omeg......