首页 > 其他分享 >莫反

莫反

时间:2024-12-24 17:19:26浏览次数:2  
标签:lfloor frac sum rfloor mu 莫反 prod

整除分块和狄利克雷卷积没啥说的。

规定莫比乌斯函数 \(\mu(i)\) 满足 \(i\) 被表示为 \(x\) 个单个质因子的积时返回值为 \((-1)^x\)。其余时为 \(0\),\(\mu(1)=1\)。

有重要性质

\[\sum_{d|n}\mu(d)=[n=1] \]

证明:不妨设

\[n=\prod^m p_i^{a_i},n'=\prod^m p_i \]

则应满足

\[\sum_{d|n}\mu(d)=\sum_{d|n'}\mu(d) \]

后面一部分的每个 \(\mu(d)\) 都是有返回值的,\(n'\) 有保证了质因子不同所以对 \(x=i\) 个质因子的情况有 \(\binom{m}{i}\) 种取法,返回值全为 \((-1)^i\)

二项式定理:

\[\sum_{d|n}\mu(d)=\sum_{d|n'}\mu(d)=\sum_{i=1}^m\binom{m}{i}(-1)^i=(1+(-1))^m=1 \]

然后好了。

有一个很典的应用,即莫反:

\[[gcd(i,j)=1]=\sum_{d|gcd(i,j)}\mu(d)=\sum_{d=1}^n[d|i][d|j]\mu(d) \]

YY的GCD

\[Ans=\sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)\in prime]\\ \]

\[=\sum_{d=1}^n[d\in prime]\sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)=d]\\ \]

\[=\sum_{d\in prime}\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{m}{d}}\sum_{x=1}^{\frac{n}{d}}[x|i][x|j]\mu(x)\\ \]

\[=\sum_{d\in prime}\sum_{x=1}^{\frac{n}{d}}\mu(x)\lfloor\frac{n}{dx}\rfloor\lfloor\frac{m}{dx}\rfloor \]

这种多变量难以整除分块就设 \(T=dx,x=\frac{T}{d}\)
。原式变为

\[=\sum_{d\in prime}\sum_{d|T}^n\mu(\frac{T}{d})\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor\\ \]

\[=\sum_{T=1}^n\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor\sum_{d|T,d\in prime}\mu(\frac{T}{d}) \]

然后前半整除分块后半调和级数预处理,有的时候后半满足积性函数,可以直接在线筛的时候处理出来。复杂度 \(O(\sqrt n)\)

其他题都大同小异,不多说了。

P3704

为啥学校oj写东西不显示latex??

就是说,把求和改成指数形式。

\[Ans=\prod_{i=1}^n\prod_{j=1}^mf(gcd(i,j))\\ \]

\[=\prod_{d=1}^nf^{\sum_{i=1}^{\lfloor \frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor \frac{j}{m}\rfloor}[gcd(i,j)=d]}(d)\\ \]

\[=\prod_{d=1}^nf^{\sum_{x=1}^{\lfloor \frac{n}{d}\rfloor}\mu(x)\lfloor\frac{n}{dx}\rfloor\lfloor\frac{m}{dx}\rfloor}(d) \]

设 \(T=dx,x=\frac{T}{d}\)。

\[Ans=\prod_{d=1}^nf^{\sum_{d|T}^n\mu(\frac{T}{d})\lfloor\frac{n}{T}\rfloor\lfloor \frac{m}{T}\rfloor}(d)\\ \]

\[=\prod_{T=1}^n\prod_{d|T}f^{\mu(\frac{T}{d}\lfloor\frac{n}{T}\rfloor\lfloor \frac{m}{T}\rfloor)}(d)\\ \]

\[=\prod_{T=1}^n(\prod_{d|T}f^{\mu(\frac{T}{d})}(d))^{\lfloor\frac{n}{T}\rfloor\lfloor \frac{m}{T}\rfloor} \]

设 \(F(x)=\prod_{d|T}f^{\mu(\frac{T}{d})}(d)\),可以调和级数复杂度算出来。

\[Ans=\prod_{T=1}^nF^{\lfloor\frac{n}{T}\rfloor\lfloor \frac{m}{T}\rfloor}(T) \]

相当于把数论分块挪到指数上,想一下就是维护前缀积算的时候拿逆元和快速幂算。

标签:lfloor,frac,sum,rfloor,mu,莫反,prod
From: https://www.cnblogs.com/Claimo0611/p/18628264

相关文章

  • 筛子与莫反
    1.Miller-RabinMiller-Rabin是一种接受随机性的正确性较高的素数检验方法,它有一定概率将合数判断为素数,但不会将素数判断为合数。其基本判定思路是,检测素数都具有但合数不具有的特殊性质,如众所周知的费马小定理\(x^{p-1}\equiv1\pmodp\)。1.1费马素性检验费马小定理的逆......
  • 莫反小练
    P1829[国家集训队]Crash的数字表格/JZPTAB不妨假设\(n<m\)。\[\begin{aligned}\sum_{i=1}^n\sum_{j=1}^m\dfrac{ij}{\gcd(i,j)}&=\sum_{d=1}^n\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=d]\dfrac{ij}{d}\\\\&=\sum_{d=1}^n\sum_{i......
  • 萌新的莫反练习笔记
    萌新的莫反练习笔记简单的数论函数恒等函数\(I(n)=1\)。元函数\(e(n)=[n=1]\)。单位函数\(id(n)=n\)。狄利克雷卷积我们设\(f\)和\(g\)的卷积\(f\astg=F\)。卷积还是一个函数。那么,\(F(n)=\sum_{d|n}f(d)g(\frac{n}{d})\)。这就是卷积。显然,\(e\astf=f\)。所以......
  • (未完工)莫反与杜教筛
    \(p\)是质数。\(p^k\)是质数的幂。广告https://github.com/August-Light/NTFuncs这是一个关于各种数论函数的Python库!前置知识数论分块模板题1:[UVA11526]H(n)模板题2:[LuoguP2261][CQOI2007]余数求和线性筛&线性筛数论函数模板题1:[LuoguP3383]【模板......
  • 莫反乱做
    莫比乌斯函数重要性质\([n=1]=\sum\limits_{d|n}\mu(d)\)应用:\(\gcd(a,b)=1\Leftrightarrow\sum\limits_{d|\gcd(a,b)}\mu(d)\)P3327[SDOI2015]约数个数和首先有:\[d(ij)=\sum_{a|i}\sum_{b|j}[\gcd(a,b)=1]\]证明参考题解1,很神秘的构......
  • 重学莫反
    cc0000想学会莫(魔)反(法)!关于莫反:最没用的东西:若\(f(x)=\sum\limits_{d|x}g(d)\),则\(g(x)=\sum\limits_{d|x}\mu(d)f(d)\)。这是莫比乌斯反演的定义式。关于一车有......
  • 莫反的一些练习题(2)
    7.P2231[HNOI2002]跳蚤原题链接题目大意求\(\suma_i\cdotx_i=1\)\((1\leN\leM\le10^8,M^N\le10^{16})\)解的个数分析一道比较有意思的题根据裴蜀定理,\(\s......
  • 莫反的一些练习题(1)
    1.树林里的树.TreesinaWood原题链接题目大意求\(\frac{\sum_{i=-a}^{a}\sum_{j=-b}^{b}\lbrack\lvertgcd(i,j)\rvert=1\rbrack}{(2a+1)(2b+1)-1}(a\le2000,b\le2......
  • 容斥与简单莫反
    容斥与莫比乌斯函数容斥原理:介绍:设集合\(S_1\simS_n\),记\(|S_i|\)表示集合\(S_i\)的大小,设\(\cup\)表示集合的并集运算,\(\cap\)表示集合的交集运算,则\[\left|\bigcup_......