数论函数——莫比乌斯反演

一些函数的一些性质

取整函数 \(\lfloor x \rfloor\)

(一)\(\lfloor x \rfloor <= x < \lfloor x \rfloor +1\)

(二)对任意x与正整数a,b\(\lfloor \lfloor \frac{x}{a} \rfloor /b\rfloor=\lfloor \frac{x}{ab}\rfloor\)

(三)对于正整数n,1 -- n中d的倍数个数为 \(\lfloor \frac{n}{d} \rfloor\)

(四)若n为正整数,\(\lfloor \frac{n}{d}\rfloor\)不同取值个数不超过\(2\times\sqrt{n}种\)

证明:
\((1)若d \leq{\sqrt{n}},\lfloor \frac{n}{d}\rfloor只有不超过\sqrt{n}种\)

\((2)若d>\sqrt{n},\lfloor \frac{n}{d} \rfloor \leq \frac{n}{d} \leq \sqrt{n},\lfloor \frac{n}{d}\rfloor 不超过\sqrt{n}种\)

\(综上,\lfloor \frac{n}{d}\rfloor 不超过2\times{\sqrt{n}}种\)

调和数

定义 \[Hn=\sum\limits_{k=1}^{n}\frac{1}{k}\]运算得\[
Hn=ln(n)+r+o(1) \]其中r为欧拉马歇罗尼常数,r约为0.577 ,因此\[
\sum\limits_{d=1}^{n} \lfloor \frac{n}{d}\rfloor = o(n\times{logn)} \]

素数计数函数

\[ \pi(n)\sim \frac{n}{ln (n)}\]\[
因此n附近素数密度近似为\frac{1}{ln(n)}\]\[
第n个素数pn\sim n\times{ln(n)}\]

数论函数

积性函数

\[ f为数论函数,对互质的正整数a,b,f(a\times{b})=f(a)\times{f(b})\]

完全积性函数

\[ f为数论函数,对任意的正整数a,b,f(a\times{b})=f(a)\times{f(b})\]
\(若f为积性函数,\) \[
n=p1^{a1}\times{p_2^{a_2}}\times{p_3^{a_3}}......\times{p_s^{a_s}}\]\[
f(n)=f(p_1^{a_1})\times{f(p_2^{a_2})}\times{f(p_3^{a_3})}......\times{f(p_s^{a_s})}\]

单位函数

\[\epsilon(n)=[n==1]= \left\{
\begin{aligned}
1&,n=1\ 0&,n!=1\\end{aligned}
\right. \]

除数函数

\(\delta_k(n)表示n因子的k次方之和\)
\[
\delta_k(n)=\sum\limits_{d|n}d^k\]

Euler函数:\(\phi(n)\)

\(Euler函数表示不超过n且与n互质的正整数个数\)
性质:
\[ n=\sum\limits_{d|n}\phi(d) \]
证明:

\(若gcd(n,i)=d,gcd(\frac{n}{d},\frac{i}{d})=1\)

\(其中\frac{i}{d}是不超过\frac{n}{d}的整数,由欧拉函数的定义,i的个数为\phi(\frac{n}{d})个\)

\(对于所有的d|n,n=\sum\limits_{d|n}\phi(\frac{n}{d})=\sum\limits_{d|n}\phi(d)\)

Dirichlet 卷积

\(f,g为数论函数,数论函数h满足\) \[
h(n)=\sum\limits_{d|n}f(d)g(\frac{n}{d})\]
\(则h为f与g的Dirichlet卷积,记为\)
\[h=f \ast g\]

性质

\((1)单位函数\epsilon为Dirichlet卷积的单位元,即对于任意函数f,有\)
\[\epsilon \ast f=f \ast \epsilon =f\]

\((2)满足交换律和结合律\)

\((3) 若f,g为积性函数,f \ast g也为积性函数\)

\((4) 逆函数:f \ast f_逆=\epsilon\)

\(定义幂函数:Id_k(n)=n^k,Id=Id_1\)

\(所以除数函数: \delta_k=1 \ast Id_k\)

\(Euler函数: \phi(n) = 1 \ast Id\)

计算Dirichlet卷积:

\(f,g为数论函数,则 f \ast g在n处的值需要枚举n的所有约数\)

\(计算f \ast g的前n项,枚举1到n中每个数的倍数\)

Mobius 函数

\[\mu(n)= \left\{
\begin{aligned}
&1&n=1 \&(-1)^s&n=p_1\times{p_2}\times{p_3}......\times{p_s}\&0&otherwise \\end{aligned}
\right. \]
\(其中p_1,p_2,p_3为素数\)

性质

\[\sum\limits_{d|n}\mu(n)= \epsilon(n) \Rightarrow \mu \ast 1 = \epsilon\]

证明:

\(n=1,显然成立\)

\(n>1\)

\(设n有s个不同的素因子,由Mobius函数定义,\)

\(当且仅当d无平方因子的时候,\mu(d)!=0\)

\(于是d中的每一个素因子的指数只能是0或1\)

\(故由二项式定理\)

\[\sum\limits_{d|n}\mu(d)=\sum\limits_{k=0}^s (-1)^k(s,k)=(1-1)^s=0\]

\(得证\)

Mobius变换

\(设f为数论函数,定义函数g满足:\)
\[g(n)=\sum\limits_{d|n}f(d)\]
\(则称g为f的Mobius变换,f是g的Mobius逆变换\)

\(Dirichlet卷积\) \[g=f \ast 1\]

Mobius反演

\(g(n)=\sum\limits_{d|n}f(d)的充要条件为f(n)=\sum\limits_{d|n}g(d)\mu(\frac{n}{d})\)

证明:

\[g=f \ast 1\]
\[f=f \ast \epsilon =f \ast 1 \ast \mu =g \ast \mu\]

得证。

应用

利用Dirichlet卷积可以解决一系列求和问题。
常用一个Dirichlet卷积替换求和式中的一部分,然后调换求和顺序,最终降低时间复杂度。

常用:
\(\mu \ast 1= \epsilon\)
\(Id = 1 \ast \phi\)

原文地址:https://www.cnblogs.com/liuquanxu/p/11741198.html

时间: 2024-10-10 20:27:23

数论函数——莫比乌斯反演的相关文章

数论函数与莫比乌斯反演

数论函数 取整函数 定义 对于实数 \(x\),记 \(?x?\) 为不超过 \(x\) 的最大整数. \(\lfloor x \rfloor\) 也是满足如下关系的唯一整数: \(\lfloor x \rfloor ≤x<\lfloor x \rfloor+1\) 对于正整数 \(n\),\(1\) 到 \(n\) 中 \(d\) 的倍数有 \(?\frac{n}{d}?\) 个 性质1 对于任意的 \(x\) 与正整数 \(a\),\(b\),我们均有: \[ ??\frac{x}{a}?/

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=\

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

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

莫比乌斯反演(入门)

细细算来,看反演已经有一两个星期了 刚开始的时候也是走了不少的弯路 和其他的算法一样,只要你懂了,就会有一种不过如此的感觉(误 感觉反演还是刚入门,不过还是先写一篇不完全的总结吧,不然过段时间就要忘记了 虽说看反演看了好久才懂,但是现在回头看看,其实很多时间还是花费在弯路上,真正的输出时间不过是最近的一两天 建议的前置技能: 容斥的简单应用 唯一分解定理 欧拉函数的定义 积性函数的定义 然后就可以看反演啦 反演推荐的资料还是贾志鹏线性筛 不要在电脑上直接过一遍就算看过了,这个资料上面有一些不加证

【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)$才

初探莫比乌斯反演

我们首先从这个函数出发:\(\mu(n)\),它被称作莫比乌斯函数 它的定义式是这样的 \[ \mu(n)= \begin{cases} 1& n=1\(-1)^r& n=p_1p_2\ldots p_r,p_1,p_2,\ldots p_r\in prime\0& p^k|n,p\in prime,k>1 \end{cases} \] 即:在对\(n\)进行质因子分解之后,所有质因子的幂次都不大于1的话,\(\mu(n)=(-1)^r\),r为n互不相同的质因子的个数.若n

莫比乌斯函数与莫比乌斯反演

莫比乌斯函数 定义 莫比乌斯函数\(\mu(n)\),当\(n=1\)时,\(\mu(n)=1\):当\(n>1\)时,设\(n\)的唯一分解式为\(n=p_1^{c_1}\cdots p_k^{c_k}\),则\(\mu(n)\)定义为 \(\mu(n)= \begin{cases} (-1)^k,c_1=c_2=\cdots=c_k=1 \\ 0, \exists\, c_i>1(1\leq i\leq k)\\ \end{cases}\) 性质 \(\sum\limits_{d|n}\m

莫比乌斯反演详解

话说这是我打的第一篇算法博客2333 话不多说直接进入正题 一.莫比乌斯函数μ 什么是μ? μ(n)={ 1,n=1: (-1)k,n=p1*p2*p3*……*pk ,pn为互不相等的质数: 0,else: } μ的性质: 1.积性函数:线性筛 在线性筛质数的基础上加几句处理就好了 void get() { mu[1]=1; for(int i=2;i<=10000000;i++) { //cout<<i<<endl; if(!p[i]) { mu[i]=-1; prime[

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

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