狄利克雷卷积及莫比乌斯反演

1.数论函数

定义两个数论函数\(f(n)\)与\(g(n)\)

则\((f+g)(n)=f(n)+g(n)\)

2.狄利克雷卷积

定义两个数论函数的狄利克雷卷积\(*\)

定义数论函数\(t=f*g\)

则\(\mathbf t(n)=\sum_{ij=n}\mathbf f(i)\mathbf g(j)\)

显然,\(f*g=g*f,f*(g*k)=(f*g)*k,f*(g+k)=f*g+f*k\)

定义数论函数\(f\)的单位元\(?\),使得\(\epsilon\ast\mathbf f=\mathbf f*?\)

不难看出,当\(i\)等于\(1\)时,\(?(i)=1\),否则,\(?(i)=0\)

定义数论函数\(g\)为函数\(f\)的逆元,即\(f*g=?\)

则\(g(n)=\dfrac{1}{f(1)}\left([n==1]?\sum\limits_{i|n,i\ne1}f(i)g(\dfrac{n}{i})\right)\)

3.积性函数

如果一个数论函数\(f\)当\(n⊥m\)的时候满足\(f(nm)=f(n)f(m)\),则称此函数为积性函数

常见的积性函数:\(*σ*_0\)(\(n\)的因数个数),\(\varphi(n)\)(\([1,n]\)中与\(n\)互质的数的个数)

结论:

1.两个积性函数的狄利克雷卷积是积性函数

证明如下:

设\(t=f*g\),则
\[
\begin{aligned}\mathbf t(nm)&=\sum_{d\mid nm}\mathbf f(d)\mathbf g\left(\frac{nm}d\right)\\&=\sum_{a\mid n,b\mid m}\mathbf f(ab)\mathbf g\left(\frac{nm}{ab}\right)\\&=\sum_{a\mid n,b\mid m}\mathbf f(a)\mathbf f(b)\mathbf g\left(\frac na\right)\mathbf g\left(\frac mb\right)\\&=\left(\sum_{a\mid n}\mathbf f(a)\mathbf g\left(\frac na\right)\right)\left(\sum_{b\mid m}\mathbf f(b)\mathbf g\left(\frac mb\right)\right)\\&=\mathbf t(n)\mathbf t(m)\end{aligned}
\]
2.一个积性函数的逆也是积性函数

证明如下:

设\(g\)为\(f\)的逆元

1.当\(nm=1\)时,则\(g(1)=1\),结论成立

2.当\(nm>1\)时,则有
\[
\begin{aligned}\mathbf g(nm)&=-\sum_{d\mid nm,d\neq1}\mathbf f(d)\mathbf g\left(\frac{nm}d\right)\\&=-\sum_{a\mid n,b\mid m,ab\neq1}\mathbf f(ab)\mathbf g\left(\frac{nm}{ab}\right)\\&=-\sum_{a\mid n,b\mid m,ab\neq1}\mathbf f(a)\mathbf f(b)\mathbf g\left(\frac na\right)\mathbf g\left(\frac mb\right)\\&=\mathbf f(1)\mathbf f(1)\mathbf g(n)\mathbf g(m)-\sum_{a\mid n,b\mid m}\mathbf f(a)\mathbf f(b)\mathbf g\left(\frac na\right)\mathbf g\left(\frac mb\right)\\&=\mathbf g(n)\mathbf g(m)-\left(\sum_{a\mid n}\mathbf f(a)\mathbf g\left(\frac na\right)\right)\left(\sum_{b\mid m}\mathbf f(b)\mathbf g\left(\frac mb\right)\right)\\&=\mathbf g(n)\mathbf g(m)-\epsilon(n)\epsilon(m)\\&=\mathbf g(n)\mathbf g(m)\end{aligned}
\]
Tips:对于一个积性函数,显然\(f(1)=1\),\(f(1)\ne1\)的情况暂时不加讨论

运用:

对于一个积性函数,可以用线性筛的方法来快速递推得到

我们可以把一个积性函数\(f\),可以把\(f(n)\)分解为\(\mathbf f(n)=\prod_{i=1}^t\mathbf f(p_i^{k_i})\),即是把它质因数分解

由于在线性筛素数的时候我们可以顺便求出每个数的最小质因数\(p_1\),最小质因数的次数\(k_1\)以及\(n/p_1^{k_1}\),就可以利用递推式\(\mathbf f(n)=\mathbf f(p_1^{k_1})\mathbf f(n/p_1^{k_1})\)直接计算 \(\mathbf f\)了

由于$σ_0 $和 \(\varphi\)在素数幂处的值很容易得到,\(k>0\)时,\(\sigma_0(p^k)=k+1,\varphi(p^k)=p^{k-1}(p-1)\)

从而推出公式:
\[
\sigma_0(n)=\prod_{i=1}^t(k_i+1), \varphi(n)=\prod_{i=1}^np_i^{k_i-1}(p_i-1)=n\prod_{i=1}^t\left(1-\frac1{p_i}\right)
\]

4.莫比乌斯反演

定义\(\mathbf1\)的逆是\(\mu\)

这样对于\(\mathbf g=\mathbf f\ast\mathbf1\),就有\(\mathbf f=\mathbf f\ast\mathbf1\ast\mu=\mathbf g\ast\mu\)

即如果\(\mathbf g(n)=\sum_{d\mid n}\mathbf f(d)\),则\(\mathbf f(n)=\sum_{d\mid n}\mu\left(\frac nd\right)\mathbf g(d)\)

这就证明了莫比乌斯反演

对于\(\mu\)的求法

显然,由于\(\mathbf1\)是积性的,而且\(\mu\)是\(\mathbf 1\)的逆,所以\(\mu\)也是积性的

可得:
\[
\mu(p^k)=\begin{cases}1&k=0\\-1&k=1\\0&k>1\end{cases}
\]
从而得出:
\[
\mu(n)=\begin{cases}(-1)^t&n=p_1p_2\dots p_t\text{且$ p_i $互不相同}\\0&n\text{不满足上述条件}\end{cases}
\]
Tips:对于\(n=1\),则\(t=0\)

反演推论:

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

从\(\mu\)是1的逆显然很好理解

原文地址:https://www.cnblogs.com/magicduck/p/11423193.html

时间: 2024-10-09 20:47:47

狄利克雷卷积及莫比乌斯反演的相关文章

狄利克雷卷积和莫比乌斯反演

狄利克雷卷积是函数与函数之间的运算 运算定义为:$(f*g)(x) = \sum_{d|x}{f(d) * g(\frac{x}{d})}$ 该运算有交换律:$f*g=g*f$ 有结合律:$f*(g*h)=(f*g)*h$ 有分配率:$f*(g+h)=f*g+f*h$ 在这种新定义的运算中有许多特殊的元素 $e$ 满足 $e(n) = [n==1]$ 即当只有$n=1$时$e(n)=1$,否则为$0$ (对命题$P$,引入符号$[P]$,如果命题成立则值为$1$,否则为$0$.) $d$ 满足$

Dirichlet卷积和莫比乌斯反演

半夜不睡写博客 1.Dirichlet卷积 定义2个数论函数f,g的Dirichlet卷积$(f*g)n=\sum_{d|n}f(d)g(\frac{n}{d})$ Dirichlet卷积满足交换律,结合律,加分分配律,若f,g为积性函数,则f*g也为积性函数 2.莫比乌斯反演 如果有2个函数f,g满足$f(n)=\sum_{d|n}g(d)$ 则有$g(n)=\sum_{d|n}\mu(d)f(\frac{n}{d})$ 两者可以互相推得,即$f=g*1 \Leftrightarrow g=\

数论入门——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛

一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数\(\mu(n)\)为: 当n有平方因子时,\(\mu(n)=0\). 当n没有平方因子时,\(\mu(n)=(-1)^{\omega(n)}\),\(\omega(n)\)表示n不同质因子的个数. 性质1: \(\sum_{d|n}\mu(d)=[n=1]\) 证明:我们把n分解质因数,则原式\(=(-1+1)^{\omega(n)}=0\). 因为对于不同的质因子,只有选和不选两种方案,这是一个组合数相加的形式,偶数加奇数减,根据二项式

狄利克雷卷积 && 莫比乌斯反演

狄利克雷卷积 && 莫比乌斯反演 狄利克雷卷积 数论函数及其运算 数论函数是指定义域是正整数,值域是一个数集的函数. 加法,逐项相加,即\((f+h)(n)=f(n)+h(n)?\): 数乘,这个数和每一项都相乘,即 \((xf)(n)=x·f(n)?\) 狄利克雷卷积 定义两个数论函数的狄利克雷卷积 \(*:?\) 若\(t=f*g?\),则\(t(n)=\sum_{d|n}^{}f(d)·g(\frac{n}{d})?\),又或者写成\(t(n)=\sum_{ij=n}f(i)\cdo

[莫比乌斯反演]【学习笔记】[更新中]

参考资料: [大部分还没看完,目前主要看了popoqqq那篇 orz] http://wenku.baidu.com/link?url=Kzzxkk64CFU7sfDeJbGKNpZpFJzJY1ZwNoaPgGo7tPSpv4KJvGAkStkpzytG46gjQuqNX7NB0merxfS4knD2H5fw7s4oHu1o1-6p16_VbEm http://wenku.baidu.com/view/77396ebb27d3240c8547ef2e.html?re=view 浅谈一类积性函数

【bzoj4407】于神之怒加强版 莫比乌斯反演+线性筛

题目描述 给下N,M,K.求 输入 输入有多组数据,输入数据的第一行两个正整数T,K,代表有T组数据,K的意义如上所示,下面第二行到第T+1行,每行为两个正整数N,M,其意义如上式所示. 输出 如题 样例输入 1 2 3 3 样例输出 20 题解 莫比乌斯反演+线性筛 $\sum\limits_{i=1}^n\sum\limits_{j=1}^m\gcd(i,j)^k\\=\sum\limits_{d=1}^{\min(n,m)}d^k\sum\limits_{i=1}^n\sum\limits

【Luogu3768】简单的数学题(莫比乌斯反演,杜教筛)

[Luogu3768]简单的数学题(莫比乌斯反演,杜教筛) 题面 洛谷 \[求\sum_{i=1}^n\sum_{j=1}^nijgcd(i,j)\] $ n<=10^9$ 题解 很明显的把\(gcd\)提出来 \[\sum_{d=1}^nd\sum_{i=1}^n\sum_{j=1}^nij[gcd(i,j)==d]\] 习惯性的提出来 \[\sum_{d=1}^nd^3\sum_{i=1}^{n/d}\sum_{j=1}^{n/d}ij[gcd(i,j)==1]\] 后面这玩意很明显的来一发

【Project Euler】530 GCD of Divisors 莫比乌斯反演

[题目]GCD of Divisors [题意]给定f(n)=Σd|n gcd(d,n/d)的前缀和F(n),n=10^15. [算法]莫比乌斯反演 [题解]参考:任之洲数论函数.pdf 这个范围显然杜教筛也是做不了的,而且考虑直接化简f(n)也遇到了困难,所以考虑将前缀和的Σ一起化简. $$F(n)=\sum_{i=1}^{n}\sum_{d|i}(d,\frac{i}{d})$$ 这一步很常见的是第一重改为枚举倍数,但这样化简后面就推不下去了. 这道题必须最后转成$\sigma_0(n)$才

莫比乌斯反演&#183;学习记录

莫比乌斯反演·学习记录 cyw在6.8左右学的莫比乌斯反演,记录一下 这个东西感觉不大好描述,我一开始也不知道这玩意能干嘛(其实现在也不知道) CYW认为,对关于一些因数/倍数关系进行操作的行为,可以用莫比乌斯反演来解决 莫比乌斯函数 这并不是什么高大上的东西,但是很有用 对于 莫比乌斯函数 的定义是 $d=1,\mu(d)=1 $ \(d=\prod_{k=1}^n p_k\;(k\in prime),\mu(d)=(-1)^k\) 即数\(d\)可以被表示为若干互异素数相乘的形式(指数不超过