首页 > 其他分享 >「学习笔记」莫比乌斯反演套路总结

「学习笔记」莫比乌斯反演套路总结

时间:2023-01-18 22:02:30浏览次数:70  
标签:lfloor frac 套路 sum rfloor times mu 反演 莫比

基本套路

Crash的数字表格 / JZPTAB
求:

\[\sum_{i=1}^n\sum_{j=1}^m\mathrm{lcm}(i,j) \]

\[\begin{aligned} &\sum_{i=1}^n\sum_{j=1}^m\mathrm{lcm}(i,j)\\ =&\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{\gcd(i,j)}\\ =&\sum_{d=1}^{\min(n,m)}\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}ijd[\gcd(i,j)=1] &\text{套路1:枚举}\gcd(i,j)\\ =&\sum_{d=1}^{\min(n,m)}d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}ij\sum_{k|\gcd(i,j)}\mu(k) &\text{套路2:拆开}[\gcd(i,j)=1]\\ =&\sum_{d=1}^{\min(n,m)}d\sum_{k=1}^{\frac{\min(n,m)}{d}}k^2\mu(k)\sum_{i=1}^{\lfloor\frac{n}{kd}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{kd}\rfloor}ij\\ =&\sum_{d=1}^{\min(n,m)}d\sum_{k=1}^{\frac{\min(n,m)}{d}}k^2\mu(k)\sum_{i=1}^{\lfloor\frac{n}{kd}\rfloor}i\sum_{j=1}^{\lfloor\frac{m}{kd}\rfloor}j &\text{套路3:枚举}k\\ =&\sum_{d=1}^{\min(n,m)}d\sum_{k=1}^{\frac{\min(n,m)}{d}}k^2\mu(k)\frac{\lfloor\frac{n}{kd}\rfloor(\lfloor\frac{n}{kd}\rfloor+1)}{2}\times\frac{\lfloor\frac{m}{kd}\rfloor(\lfloor\frac{m}{kd}\rfloor+1)}{2}\\ =&\sum_{T=1}^{\min(n,m)}\sum_{k|T}\frac{T}{k}\times k^2\mu(k)\frac{\lfloor\frac{n}{T}\rfloor(\lfloor\frac{n}{T}\rfloor+1)}{2}\times\frac{\lfloor\frac{m}{T}\rfloor(\lfloor\frac{m}{T}\rfloor+1)}{2} &\text{套路4:设}T=kd\\ =&\sum_{T=1}^{\min(n,m)}T\times \frac{\lfloor\frac{n}{T}\rfloor(\lfloor\frac{n}{T}\rfloor+1)}{2}\times\frac{\lfloor\frac{m}{T}\rfloor(\lfloor\frac{m}{T}\rfloor+1)}{2}\sum_{k|T}k\mu(k)\\ \end{aligned}\]

\(O(n\log n)\) 暴力计算出 \(g(T)=\sum_{k\mid T}k\mu(k)\) 的前缀和, 再进行数论分块即可。

拆函数

\[\varphi(ij)=\frac{\varphi(i)\varphi(j)}{\varphi(\gcd(i,j))}\times \gcd(i,j) \]

\[\mu(ij)=\mu(i)\mu(j)[\gcd(i,j)=1] \]

\[d(ij)=\sum_{x|i}\sum_{y|j}[\gcd(x,y)=1] \]

线性筛求函数前缀和

于神之怒加强版
前部分省略。
最后需要求一个函数的前缀和:

\[g(T)=\sum_{d|T}d^k\mu(\frac{T}{d}) \]

两个积性函数的卷积也是积性函数,可以将 \(g(T)\) 拆为 \(\prod{g(p_i^{c_i})}\)。

\[\begin{aligned} g(T)&=\prod g(p_i^{c_i})\\ &=\prod (p_i^{k\times (c_i-1)}\times \mu(p_i)+p_i^{k\times c_i}\times \mu(i))\\ &=\prod (p_i^{k\times (c_i-1)}\times (p_i^k-1)) \end{aligned}\]

因为线性筛是不断乘质因子的过程,所以这样就可以线性筛求出这个函数。

标签:lfloor,frac,套路,sum,rfloor,times,mu,反演,莫比
From: https://www.cnblogs.com/apjifengc/p/17060676.html

相关文章

  • 反演
    本人刚学OI一秒,对了解众多反演之后对反演的本质有一些看法本人现阶段学过的反演和类似反演的包括:容斥,二项式反演,莫比乌斯反演,子集和反演。Min-Max容斥单位根反演斯特林反......
  • 莫比乌斯反演学习笔记
    莫比乌斯函数定义\[\mu(n)=\begin{cases}1&n=1\\0&n\text{含有平方因子}\\(-1)^k&\text{其中}k\text{为}n\text{本质不同的质因子个数}\end{cases}......
  • 莫比乌斯反演
    莫比乌斯反演考虑狄利克雷前缀和的式子,\(g(n)=\sum\limits_{d\midn}f(d)\)。知\(f\)求\(g\)是naive的;知\(g\)求\(f\)呢?首先上式等价于\(g=f\astI\)......
  • 欧拉函数与莫比乌斯函数的一些性质
    前置知识翡蜀定理与算数基本定理的证明积性函数若有一个数论函数\(f\)满足以下性质:\(\left(1\right)f\left(1\right)=1\)\(\left(2\right)\)若\(a,b\)互质,那么\(f\le......
  • 学习开源项目的几个实用套路
    记得我的leader之前说过,很多人工作之后就丧失了钻研技术的热情,这个确实,我发现自己多少也有这个问题。转眼已经毕业一年多了,回想这一年,有些惭愧,感觉不仅技术能力上并没有......
  • P3704 [SDOI2017]数字表格——莫比乌斯反演
    莫比乌斯反演\(\color{red}{f(n)=\sum\limits_{d|n}}g(d)\Leftrightarrowg(n)=\sum\limits_{d|n}\mu(d)f(\dfrac{n}{d})\)例题:P3704[SDOI2017]数字表格题意:给出\(n,......
  • 狄利克雷卷积和莫比乌斯反演初探(施工中)
    0.前置知识瞅这1.狄利克雷卷积定义定义域为\(\mathbb{N_+}\)的函数称为数论函数。对于两个数论函数\(f,g\),其狄利克雷卷积为\(h(n)=\sum\limits_{d\midn}f(d)......
  • 莫比乌斯反演草稿纸
    这只小蒟蒻做莫比乌斯反演的练习题时用$\LaTeX$打了一些草稿,又不舍得扔掉,于是放在这个博客里供大家欣(tu)赏(cao)。P2257YY的GCD$$\text{求}\sum_{x=1}^n\sum_{y=1}^m[g......
  • HDU 6439 2018CCPC网络赛 Congruence equationI(杜教筛 + 莫比乌斯反演 + 伯努利数)
      大致题意:给你一个长度为k的序列a。对于序列c,当  时,;当时,取[0,m)中任意一个数字。令  表示满足  的序列c的方案数。现在让你求 。          ......
  • 莫比乌斯反演
    莫比乌斯函数 设正整数$x=p_1^{c_1}p_2^{c_2}...p_m^{c_m}$,定义函数 $\left(\begin{array}{ccc}1&2&3\\4&5&6\\7&8&9\\0&0&0\end{array}\right)$......