首页 > 其他分享 >数论 Part : Dirichlet 卷积 & 莫比乌斯反演 & 杜教筛

数论 Part : Dirichlet 卷积 & 莫比乌斯反演 & 杜教筛

时间:2024-08-25 12:37:51浏览次数:5  
标签:bf Dirichlet 函数 text sum 杜教 mu 反演 alpha

\(\text{-1 前言}\)

\(\text{-1.0 日志}\)

24.08.24:启动本文企划,正式着笔。

\(\text{-1.1 本文记号说明}\)

本文使用 \(\cdot\) 表示乘号,\(*\) 表示卷积,\(\mathbb{P}\) 表示质数集。

\(\text{0 基础函数科技}\)

  1. 单位函数 \({\bf 1}(x)=1\)。
  2. 幂函数 \(id^k(x)=x^k\)。
  3. 恒等函数(幂函数的一种) \(id(x)=x\)。
  4. 单位元函数 \(\epsilon(x)=[x=1]\)
  5. 欧拉函数 \(\varphi(x)=\sum_{i=1}^x [\gcd(i,x)=1]\)
  6. 莫比乌斯函数 \(\mu(x)=\left\{ \begin{aligned}1, & & x=1,\\0, & & x\text{ 具有平方质因子,} \\ (-1)^k,& &k\text{为不同质因子个数.} \end{aligned}\right.\)
  7. 约数幂函数 \(\sigma^k(x)=\sum_{d|x}d^k\)
  8. 约数个数函数 \(d(x)=\sigma^0(x)=\sum_{d|x}1\)。
  9. 约数和函数 \(\sigma(x)=\sum_{d|x}d\)

\(\text{1 Dirichlet 卷积}\)

\(\text{1.0 定义}\)

两个数论函数 \(f(x)\) 和 \(g(x)\) 的 \(\text{Dirichlet}\) 卷积定义为:

\[h(x)=(f*g)(x)=\sum_{d|n} f(d)g\left(\frac{n}{d}\right)=\sum_{i\cdot j=n} f(i)g(j). \]

后者形式通常用于证明某些命题。

\(\text{1.1 性质}\)

  1. 交换律:\(f*g=g*f\),显而易见的,读者自证不难。
  2. 结合律:\(f*(g*h)=(f*g)*h\)。下文给出证明。
  3. 分配律:\((f+g)*h=f*h+g*h\),读者自证不难。
  4. 两个积性函数的卷积仍是积性函数。

性质二的证明:

\[(f*(g*h))(n)=\sum_{i|n} \left(f(i) \sum_{j\left|\frac{n}{i}\right.}g(j)h\left(\frac{n}{ij}\right) \right) = \sum_{i\cdot j\cdot k=n} f(i)g(j)h(k). \]

可发现任意的顺序运算最终结果均为此。

\(\text{1.2 代数结构}\)

若用 \(\mathbb{F}\) 表示数论积性函数集,则由她和 \(\text{Dirichlet}\) 卷积构成的代数结构 \((\mathbb{F},*)\) 是一个阿贝尔群。

  1. 封闭性:两个数论函数的 \(\text{Dirichlet}\) 卷积显然仍是一个数论函数。
  2. 单位元:该群的单位元是 \(\epsilon\),即任意数论函数 \(f(x)\) 都有 \(f*\epsilon=f\)。
  3. 结合律:于 \(1.1\) 节已证。
  4. 逆元:对于一个数论积性函数 \(f\),满足 \(f*g=\epsilon\) 的函数 \(g\) 称之为 \(f\) 的 \(Dirichlet\) 逆,记作 \(f^{-1}\)。可以证明,任意积性数论函数均有其逆元。
  5. 交换律:于 \(1.1\) 节已证。

莫比乌斯函数除了 \(0\) 节的第二定义以外,它的第一定义为 \(\mu=e^{-1}\),即莫比乌斯函数是单位函数的 \(\text{Dirichlet}\) 逆。可以得到 \(\sum_{d|n}\mu(d)=[n=1]\)。

\(\text{1.3 基础函数的 Dirichlet 卷积}\)

  1. \(\epsilon * f=f\)。由 \(1.2\) 节已证。
  2. \(id * {\bf 1}=\sigma\)。证明见下文。
  3. \({\bf 1}*{\bf 1}=d\)。证明见下文。
  4. \({\bf 1} * \varphi = id\)。证明见下文。
  5. \(\mu*id=\varphi\)。证明见下文。

对上述第二点的证明:

\[(id*{\bf 1})(n) =\sum_{d |n}id(d){\bf 1}\left(\frac{n}{d}\right) = \sum_{d|n}d=\sigma(n). \]

对上述第三点的证明:

\[({\bf 1}*{\bf 1})(n)=\sum_{d|n}{\bf 1}(d){\bf 1}\left(\frac{n}{d}\right)=\sum_{d|n}1=d(n).\]

对上述第四点的证明:要证它,即证 \(\phi(n)=\sum_{d|n}\varphi(d)=n\)。

  1. 当 \(n=1\) 时,等式显然成立。
  2. 当 \(n\in \mathbb{P}\) 时,等式左边为 \(\varphi(1)+\varphi(n)=1+n-1=n\),等式成立。
  3. 当 \(n=p^\alpha,p\in\mathbb{P},\alpha\in \mathbb{N^*}\) 时,等式左边等于 \(\sum_{i=0}^\alpha \varphi(p^i)=1+\sum_{i=1}^\alpha p^i-p^{i-1}=p^\alpha-p^{\alpha-1}+p^{\alpha-1}-p^{\alpha-2}+\cdots+p^1-p^0+1=p^\alpha=n\),等式成立。
  4. 当 \(n=\prod_{i=1}^m p_i^{\alpha_i},p_i\in\mathbb{P},\alpha_i\in\mathbb{N^*}\) 时,由 \(\text{Dirichlet}\) 卷积的性质三,\(\phi(n)\) 也是一个积性函数,则 \(\phi(n)=\prod_{i=1}^m \phi(p_i^{\alpha_i})=\prod_{i=1}^n p_i^{\alpha_i}=n\)。
  5. 证毕。

对上述第五点的证明:根据第四点,等式两边同卷积 \(\mu\) 得 \(\mu * {\bf 1} * \varphi = id * \mu\),因为 \(\mu\) 和 \(\bf 1\) 互为逆元,所以 \(\varphi = id * \mu\),得证。

\(\text{2 莫比乌斯反演}\)

定理:对于两个数论函数 \(f(x)\),\(g(x)\),若 \(f(n)=\sum_{d|n} g(n)\),则 \(g(n)=\sum_{d|n}\mu(d)f\left(\frac{n}{d}\right)\),此时称 \(f\) 为 \(g\) 的莫比乌斯变换,\(g\) 为 \(f\) 的莫比乌斯反演。

证明:使用 \(\text{Dirichlet}\) 卷积。容易发现,\(f=g*{\bf 1}\),所以将两边同时卷积一个莫比乌斯函数 \(\mu\),可以得到 \(f*\mu=g*{\bf 1}*\mu\),由于 \(\mu\) 是单位函数 \(\bf 1\) 的 \(\text{Dirichlet}\) 逆,则 \(f*\mu=g\),得证。

\(\text{3 杜教筛}\)

杜教筛是一种用于求解积性函数前缀和的思想。

\(\text{3.0 欧拉函数前缀和}\)

定义:

\[\Phi(n)=\sum_{i=1}^n \varphi(i) \]

\(\text{3.0.0 引理}\)

\(\text{3.0.1 证明}\)

\(\text{3.1 莫比乌斯函数前缀和}\)

\(\text{3.0.0 引理}\)

\(\text{3.0.1 证明}\)

\(\text{4 后记}\)

标签:bf,Dirichlet,函数,text,sum,杜教,mu,反演,alpha
From: https://www.cnblogs.com/LaDeX-Blog/p/18378825

相关文章

  • 【2】容斥与二项式反演
    【2】容斥与二项式反演1.1容斥原理容斥原理基于的是下面的恒等式:\[\sum\limits_{i=0}^n\dbinom{n}{i}(-1)^i=0^n=[n=0]\]这个式子有什么意义呢?我们考虑一个长度为\(N\)的序列,并且要求其中每个元素都满足某个限制,计算满足这个条件的序列数量。每个元素都满足限制\(\Leftri......
  • 9. 容斥与反演
    9.容斥与反演容斥原理:\[|\bigcup_{i=1}^nP_i|=\sum_{S\subseteqU}(-1)^{|S|-1}|\bigcap_{s\inS}P_s|\]感性理解:\(P_i\):”满足某种性质的元素的集合“;左边:具有任意一种性质的元素的并,右边:至少具备多个性质的元素。证明:考虑一个元素\(x\),设其包含在\(k\)个集合内,那么当......
  • 2024牛客多校第九场 C.Change Matrix 欧拉反演
    这题是欧拉反演的应用,之前没学过欧拉函数和欧拉反演,傻傻对着\(gcd(i,j)\)不知道怎么化简。首先对原来的矩阵进行转化,拆成\(n\)个小矩阵因为\(gcd(i,j)=\sum_{x|i,x|j}\phi(x)\)这是因为对于任意的正整数\(n\)都有\(n=\sum_{d|n}\phi(d)\),证明见oiwiki:https://oi-wi......
  • 反演(2)
    CP4反演与共轴圆系还是有很大关联的。我们说,共轴圆系反演后还是共轴圆系,理由如下:对于有两个交点的共轴圆系,反演后的所有圆还是过这两个点(的对应点),所以还是共轴圆系对于切于某点的共轴圆系,由反演的保相切,它们依旧相切与一点对于无交点的共轴圆系,我们找到与它共轭的共轴圆......
  • 反演(1)
    反演是一种几何变换。在给出它的具体变换前,需要明确几个概念:直线是一种退化的圆,我们将直线与圆统称为广义圆所有直线交于一个点,即无穷远点\(P_\infty\)需要指出的是,反演中所述的无穷远点只有一个,这与射影几何中无穷个的无穷远点有一定区别上述的定义可以给出广义圆的相......
  • 【无人机】基于动态反演和扩展状态观测器的无人机鲁棒姿态控制研究(Matlab代码实现)
     ......
  • lg容斥与反演
    容斥与反演容斥之前从没有搞清楚的:容斥是一种方法,为了做到不重复计数,先算总和再去除重复的方法。所以我们可以计算任意具备一种性质的元素个数(并),通过计算“至少具备了某些元素的个数”(交)。另一种形式:总数-不满足所有性质的元素=任意满足一种性质的元素此时,不满足所有性质即......
  • 二项式定理+二项式反演
    序(感谢9G对本博客证明的大力支持)前置知识1:排列组合2:多步容斥\[\dbinom{n}{m}=\binom{n}{n-m}=\mathrm{C}_n^m=\mathrm{C}_n^{n-m}\]\[\dbinom{n}{m}*\binom{m}{s}=\dbinom{n}{s}*\binom{n-s}{n-m}\]\[\dbinom{n}{x-1}+\binom{n}{x}=\dbinom{n+1}{x}\]证明:\(\dbinom{n}{x......
  • BZOJ2839/LG10596 集合计数 题解(二项式反演+扩展欧拉定理)
    题目大意:一个有\(N\)个元素的集合有\(2^N\)个不同子集(包含空集),现在要在这\(2^N\)个集合中取出若干集合(至少一个),使得它们的交集的元素个数为\(K\),求取法的方案数,答案模\(10^9+7\)。为表述方便,不妨设这\(i\)个元素分别为\(1\simn\)。前置知识:二项式反演。考虑设\(g(......
  • 容斥定理及二项式反演
    二项式定理:\[(a+b)^n=\sum_{i=0}^{n}\binom{n}{i}a^ib^{n-i}\]很好理解。我们经常会使用的式子:\[(1+x)^n=\sum_{i=0}^{n}x^i\binom{n}{i}\]容斥定理:\[\begin{split}\left|\bigcup_{i=1}^{n}S_i\right|=&\sum_{i}|S_i|-\sum_{i<j}|S_i\capS_j|+\sum_{i<j&......