首页 > 其他分享 >莫比乌斯反演

莫比乌斯反演

时间:2023-03-24 20:22:31浏览次数:44  
标签:frac 乌斯 sum mu 反演 莫比

莫比乌斯反演

引入

莫比乌斯反演用处:对于一些函数 \(f(n)\),如果比较难以求出它的值,但容易求出其倍数和或约束和 \(g(n)\),则可以通过莫比乌斯反演简化运算。

莫尼乌斯函数

定义

定义 \(\mu\) 为莫比乌斯函数,

\[\mu(n) = \begin{cases} 1 & n=1\\ 0& n含平方因子\\ (-1)^k & k为n的本质不同质因子个数\\ \end{cases} \]

详细解释一下后两条,我们令 \(n=\prod_{i=1}^{k}p_i^{c_i}\),也就是将 \(n\) 分解质因数,\(p_i\) 为质因子,\(c_i\ge 1\),则凡是有 \(c_i>1\),\(\mu(n) = 0\),而当任意的 \(c_i\) 都等于 \(1\),\(\mu(n) = (-1) ^ k\)。

性质

\[\sum_{d|n}\mu(d) = \begin{cases} 1 &n=1\\ 0 & n \ne1 \end{cases} \]

证明

既然 \(n=\prod_{i=1}^{k}p_i^{c_i}\) ,那么我们令 \(d = \prod_{i=1}^k p_i^{\beta_i}\),其中 \(0 \le \beta_i \le c_i\)

对于存在 \(\beta_i>=2\) 的情况,\(\mu(d)=0\),我们可以不管。

对于所有 \(\beta_i\) 都小于等于 1 的情况,很显然我们按照\(\beta_i=1\)有几个给它分一下组,也就是按照选了几个质因数给它分组,这样,很显然 \(\mu(d)\) 的值是 \(\C_{k}^0\times(-1)^0 + \C_{k}^1\times(-1)^1+\cdots+\C_{k}^k\times(-1)^k\) 也就是 \(\sum_{i=0}^k\C_k^i\times(-1)^i\)

看到这个形式,我们可以想起来二项式定理:

\[(a+b)^k=\C_k^0a^kb^0+\C_k^1a^{k-1}b^1+\cdots+\C_k^ka^0b^k \]

我们可以令 \(a=1,b=-1\),代入二项式定理正好就是我们刚才推出的式子。

所以

\[\sum_{d|n} \mu(d) = \sum_{i=0}^k\C_k^i\times(-1)^i=(1-1)^k=0 \]

莫比乌斯反演

形式一

定义在正整数域上的两个函数,若

\[F(n) = \sum_{d|n}f(d) \]

\[f(n) = \sum_{d|n}\mu(d)F(\frac nd) \]

莫比乌斯反演相关的题都是去套用这个定理来简化 \(f(n)\) 的计算。

证明

首先将 第一个式子代入第二个式子,消去 \(F(n)\):

\[\sum_{d|n}\mu(d)F(\frac nd)=\sum_{d|n}\mu(d)\sum_{i|\frac nd}f(i) \\ \]

这其实就相当于一个二重循环:

for d|n
    for i|(n/d) 
        sum += mu(d) * f(i)

考虑将循环顺序颠倒,没有影响,\(i\) 可遍历到 \(n\) 的任何因数。

再考虑对于 \(\mu(n)\) 的遍历,既然 \(i|\frac nd\),那么 \(di|n\),那么 \(d|\frac ni\)

则:

\[\sum_{d|n}\mu(d)F(\frac nd)=\sum_{d|n}\mu(d)\sum_{i|\frac nd}f(i) \\=\sum_{i|n}f(i)\sum_{d|\frac ni}\mu(d) \]

前面已经证明过 \(\sum_{d|n}\mu(d) = \begin{cases} 1 &n=1\\ 0 & n \ne1 \end{cases}\),则当只有 \(\frac ni = 1\) 时才会对答案有贡献,其他时候都是 0,则 \(i=n\) 时,答案为 \(f(n)\),证毕。

形式二

我们一般会用到莫比乌斯反演的另外一种形式:

\[若 F(n) = \sum_{n|d}f(d)\\ 则 f(n) = \sum_{n|d}\mu(\frac dn) F(d) \]

证明

形式二的证明与形式一略有不同,但是大致一样。

首先依旧是将第一个式子代入第二个式子。

\[\sum_{n|d}\mu(\frac dn) F(d) = \sum_{n|d}\mu(\frac dn)\sum_{d|i}f(i) \]

设 \(d'=\frac dn\),\(d = d'n\),因为 \(d|i\),所以 \(d'n|i\),所以 \(d'|\frac in\)

\[\sum_{n|d}\mu(\frac dn) F(d) = \sum_{n|d}\mu(\frac dn)\sum_{d|i}f(i) \\ = \sum_{n|i}f(i)\sum_{d'|\frac in} \mu(d') \]

与前面同理,\(\mu(n)\) 只有在 \(n=1\) 时才是 1,所以最终答案为 \(f(n)\).

例题

[HAOI2011]Problem B

[SDOI2015]约数个数和

VLATTICE - Visible Lattice Points

(太懒了懒得写了回头补

标签:frac,乌斯,sum,mu,反演,莫比
From: https://www.cnblogs.com/djc01/p/17253225.html

相关文章

  • 二项式反演
    学习参照cmd的博客,知乎,oi-wiki,某神仙的博客组合恒等式\[\binom{n}{k}=\binom{n-1}{k}+\binom{n-1}{k-1}\\\binom{n}{n_1,n_2,\ldots,n_k}=\frac{(n_1+n_2+\ldots+n_k......
  • 莫比乌斯
         ......
  • 【应用】Lagrange 反演应用
    证明鸽了,所以先开始应用篇。对于一元多项式\(F,G\)我们有Lagrange反演公式:\[n[x^n]F^k=k[w^{-k}]G^{-n}\]绝大多数情况我们都取\(k=1\)。其中多项式\(G\)为\(F......
  • 【洛谷】P2257 YY的GCD(莫比乌斯反演)
    原题链接题意\(T\)组询问,每次询问求:\[\sum_{i=1}^{n}\sum_{i=1}^{m}[\gcd(i,j)\inprime]\]\(T=10^4,n,m\leq10^7\)。思路不难想到枚举质数,将原式化简为:\[\sum......
  • 莫比乌斯反演 & 狄利克雷卷积
    大家好,我不会数学实锤了。文章内容较杂,分章节叙述了的大部分有关内容。为什么把这俩放一起?我不知道。积性函数积性函数:\(\foralla,b\),\(a\perpb\),如果一个函数\(f\)......
  • 快速莫比乌斯/沃尔什变换 (FMT/FWT) 学习笔记
    引入考虑一个基本问题:给定序列\(a_n,b_n\),求出序列\(c_n\),满足\(c_i=\sum_{j\oplusk=i}a_jb_k\),其中\(\oplus\)是一种二元运算符,形如上式的问题一般被称为卷积。......
  • 浅析排列组合、斯特林数、贝尔数、二项式定理与推论及其反演、子集反演、广义容斥
    浅析排列组合、斯特林数、贝尔数、二项式定理与推论及其反演、子集反演、广义容斥目录浅析排列组合、斯特林数、贝尔数、二项式定理与推论及其反演、子集反演、广义容斥更......
  • <学习笔记> 关于二项式反演
    1容斥原理的式子:\[|A1∪A2∪...∪An|=\sum_{1≤i≤n}|Ai|−\sum_{1≤i<j≤n}|Ai∩Aj|+...+(−1)^{n−1}×|A1∩A2∩...∩An|\]一般来说不会直接用容斥原理这个式子,而是......
  • 反演随记
    零、反演的本质$$A\vec{x}=\vec{y}\iffB\vec{y}=\vec{x}$$其中\(\vec{x},\vec{y}\)为列向量,\(A,B\)为任意矩阵。所以反演的证明即证明\(A,B\)互逆,可以通过......
  • 各种反演技巧
    二项式反演\[\sum_{i=0}^n\binom{n}{i}(-1)^i=[n=0]\]所以得到\[f_n=\sum_{i=0}^n\binom{n}{i}g_i\\g_n=\sum_{i=0}^n\binom{n}{i}(-1)^{n-i}f_i\]考虑每一项的贡......