标签:lfloor frac 乌斯 sum mid rfloor mu 反演 莫比
积性函数
积性函数:对于任意互质的整数 \(a,b\) 有 \(f(ab)=f(a)f(b)\) 则称 \(f(x)\) 的数论函数。
完全积性函数:对于任意整数 \(a,b\) 有 \(f(ab)=f(a)f(b)\) 的数论函数。
-
常见的积性函数:\(\varphi,\mu,\sigma,d\)
-
常见的完全积性函数:\(\epsilon,I,id\)
\(\epsilon(n) = [n=1], I(n) = 1, id(n) = n\)
狄利克雷卷积
设 \(f, g\) 是两个数论函数,它们的狄利克雷卷积卷积是:\((f*g)(n) = \sum \limits _{d | n} f(d) g(\frac{n}{d})\)
性质:满足交换律,结合律
单位元:\(\epsilon\) (即 \(f*\epsilon=f\))
结合狄利克雷卷积得到的几个性质:
- \(\mu * I = \epsilon\)
- \(\varphi * I = id\)
- \(\mu * id = \varphi\)
莫比乌斯函数
定义
\[\mu(n)=\left\{ \begin{aligned} 1\ &(n=1)\\ (-1)^m\ &(m 为质数个数)\\ 0\ &(含有平方因子) \end{aligned} \right. \]性质
- \[\lfloor\frac{n}{xy}\rfloor=\lfloor\frac{\lfloor\frac{n}{x}\rfloor}{y}\rfloor \]
- \[\sum_{d|n}\mu(d)=[n=1] \]
- \[\sum_{d\mid n}\frac{\mu(d)}{d}=\frac{\varphi(n)}{n} \]
-
\[\sigma_k(nm)=\sum_{x\mid n}\sum_{y\mid n}[\gcd(x,y)=1](\frac{n}{x}\cdot y)^k
\]
\(\sum_{i=1}^n\sum_{j=1}^m\sigma_k(ij)\) 的推式子过程
公式
\(f(i), g(i)\) 是两个数论函数:
或者
\[\begin{aligned} &f(i) = \sum_{d = 1}^{\left\lfloor\frac n i\right\rfloor}g(d\cdot i)\\ \Rightarrow &g(i) = \sum_{d = 1}^{\left\lfloor\frac n i\right\rfloor}f(d\cdot i)\mu(d) \end{aligned} \]例题
求 \(\gcd(i,j)=k\) 的个数
\[\sum_{i=1}^n\sum_{j = 1}^m[\gcd(i,j)=k] \]设 \(f(k)\) 为 \(\gcd(i,j)=k\) 的个数,\(g(k)\) 为 \(\gcd(i,j)\) 是 \(k\) 的倍数的个数,则:
\[g(k) = \sum_{d = 1}^{\left\lfloor\frac n k\right\rfloor}f(d\cdot x) \]而 \(g(k)=\lfloor\frac{n}{k}\rfloor\cdot\lfloor\frac{m}{k}\rfloor\),故直接莫反+数论分块求 \(f(k)\) 即可。
求 \(\gcd(i,j)\) 的 \(k\) 次幂
\[\sum_{i=1}^n\sum_{j=1}^m\gcd(i,j)^k\bmod{10^9+7} \]用之前的做法将得到:
\[\sum_{d=1}^nd^k\sum_{x=1}^{\lfloor\frac{n}{d}\rfloor}\mu(x)\lfloor\frac{n}{dx}\rfloor\lfloor\frac{m}{dx}\rfloor \]这样可以做到 \(O(n^{\frac{3}{4}})\),但是还有更快的做法,枚举 \(dx\),记为 \(T\):
\[\sum_{T=1}^n\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor\sum_{d\mid T} d^k\cdot\mu(\frac{T}{d}) \]设 \(f(n)=\sum_{d\mid n}d^k\mu(\frac{n}{d})\),则只要求出 \(f\) 的前缀和,就可以 \(O(\sqrt n)\) 解决单次询问。
Sengxian 说
这个交换枚举顺序,然后变成积性函数预处理+数论分块的技巧很常见。
标签:lfloor,frac,乌斯,sum,mid,rfloor,mu,反演,莫比 来源: https://www.cnblogs.com/Callis/p/16147563.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。