数论函数与莫比乌斯反演

数论函数

取整函数   

定义

对于实数 \(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}?/b? = ?\frac{x}{ab}?
\]

  性质2

\(?\frac{n}{d}?\) 可能的取值不超过 2√n 种。

  证明

对于正整数 \(n\),考虑当 \(1≤d≤n\) 时,\(?\frac{n}{d}?\) 的不同的取值个数。
若 \(d≤\sqrt n\),则能得到的 \(?\frac{n}{d}?\) 只有不超过 \(\sqrt n\) 种。
若 \(d>\sqrt n\),则 \(?\frac{n}{d}?≤\frac{n}{d}<\sqrt n\),又因为 \(?\frac{n}{d}?\)是正整数,故此时
可能的取值也不超过 \(\sqrt n\) 种。
综上,\(?\frac{n}{d}?\) 可能的取值不超过 \(2\sqrt n\) 种。

调和数 

\[
\begin{align}
&H_{n}=\sum_{k=1}^{n}\frac{1}{k}\&=\ln n+\gamma+o(1)
\end{align}
\]

可以推出:
\[
\begin{align}
\sum_{d=1}^{n}\lfloor \frac{n}{d} \rfloor=\Theta(n logn)
\end{align}
\]

素数计数函数

  定义&素数定理

令素数计数函数 \(\pi(n)\) 表示不超过 \(n\) 的素数个数。我们有如下的素数定理:

\[
\pi(n)\sim\frac{n}{\ln n}
\]

  推论:  

\(n\) 附近的素数密度近似是 \(\frac{1}{\ln n}\) 。

第 \(n\) 个素数 \(p_{n}\sim n \ln n\) 。

积性函数

  定义

设 \(f\) 是数论函数,若对任意互质的正整数 \(a, b\) ,都有 \(f(ab) = f(a)f(b)\) ,则称 $f $ 是积性函数。

若对任意的正整数 \(a, b\) ,都有 \(f(ab) = f(a)f(b)\) ,则称 $ f $ 是完全积性的。

单位函数

  定义

单位函数 \(\epsilon (n)\) 定义为:
\[
\begin{align}
\epsilon(n)=[n=1]=\left\{
\begin{matrix}
&1,n=1;\\
&0,n\neq1.
\end{matrix}\right.
\end{align}
\]

除数函数

  定义

除数函数 \(\sigma_{k}\) 用来表示 \(n\) 的因子的 \(k\) 次方之和:  

\[
\begin{align}
\sigma_{k}(n)=\sum_{d|n}d^{k}
\end{align}
\]

约数个数 \(\sigma_{0}(n)\) 常记为 \(d(n)\) ,约数和 \(\sigma_{1}(n)\) 常记为 \(\sigma(n)\) 。
除数函数都是积性函数。

\(Euler\) 函数

  定义:

\(Euler\) 函数 \(φ(n)\) 表示不超过 \(n\) 且与 \(n\) 互质的正整数的个数。

由 \(n\) 的标准分解并结合容斥原理,我们可以得到 \(Euler\) 函数的表达式:\[\varphi(n)=n\cdot\prod_{i=1}^{s}(1-\frac{1}{p_{i}})\]

其中 \(n = p_{1}^{\alpha_1}p_{2}^{\alpha_2} · · · p_{s}^{\alpha_s}\) 是 \(n\) 的标准分解。
由此易见 \(Euler\) 函数是积性函数。

  性质1

对于任意 $ n\(,\)Euler$ 函数有如下性质:

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

  证明1

将 \(1\) 到 \(n\) 中的所有整数按与 \(n\) 的最大公约数分类。

若 \(gcd(n, i) = d\),那么 \(gcd(\frac{n}{d} , \frac{i}{d} ) = 1\) 。而又 \(\frac{i}{d}\) 是不超过 \(\frac{n}{d}\) 的整数,故这样的 \(i\) 有 \(φ(\frac{n}{d})\) 个。
考虑所有 \(d | n\),我们也就考虑到了所有 \(1\) 到 \(n\) 之间的 \(n\) 个整数,因此有\[n=\sum_{d|n}\varphi(\frac{n}{d})=\sum_{d|n}\varphi(d)\]即:
\[
Id=\varphi *1
\]

证明2

可以先证明 \(f(n)=\sum_{d|n}\varphi (d)\) 为积性函数,然后再证对于质数 \(p\),有
\[
f(p^c)=\sum_{d|p^c}\varphi(d)=\varphi(1)+\varphi(p)+\varphi(p^2)+...+\varphi(p^{c-1})=p^c
\]
可以通过定义与等比数列求和得出,然后结论易得

小技巧,研究一个积性函数,先研究其在质数的幂时的表现。

性质2

\(p|n\) ,则 \(\varphi(np)=\varphi(n)p\)。

使用时有:

\(\varphi(p^k)=(p-1)p^{(k-1)}\)

代码中:

phi[t]=phi[i]*(i%p[j]?p[j]-1:p[j]);

\(Mobius\) 函数

莫比乌斯函数的定义:
\[
\mu(n)=
\begin{aligned}
&1 && n=1 \&0 && 有完全平方因子 \&(-1)^p && 是p个不同素因子积\\end{aligned}
\]
第一次看到这个时我很懵逼,感觉这个函数不那么自然,不知道为什么有这个函数……

但实际上,它是常函数 \(1\) 的逆,即 $\mu * 1=\epsilon $。

而具体推导可以看这里。(其实好像算是一种容斥?)

而这满足了进行莫比乌斯反演的需要。

积性函数的逆:

若 \(f*g=\epsilon\),则 \(f\) 与 \(g\) 互逆。

莫比乌斯反演

现有关系:

\[
\begin{align}
F(n)=\sum_{d|n}f(d)
\end{align}
\]
即:
\[
F=f*1
\]
如果我们易求 \(f\),那么就可以轻松求出 \(F\),反之,若 \(F\) 易求,我们如何求出 \(f\)?
\[
F*\mu=f*1*\mu=f*\epsilon=f
\]
倍数的莫比乌斯反演:

若:
\[
F(n)=\sum_{n|d}f(d)
\]
则:
\[
f(n)=\sum_{n|d}\mu(d/n)F(d)
\]
又作:

若:
\[
F(n)=\sum_{k=1}^{\infty}f(kn)
\]
则:
\[
f(n)=\sum_{k=1}^{\infty}\mu(k)f(kn)
\]

技巧:
\[
[gcd(i,j)=1]=\sum_{d|gcd(i,j)}\mu(d)
\]
证明:

由 \(\mu * 1=\epsilon\),

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

将 \(n\) 替换成 \(gcd(i,j)\) 就是上式了。

然后 \(d|gcd(x,y)\) 可以转化为 \(d|x,d|y\),

然后经常可以枚举 \(d\) 来根据 \(i\),\(j\) 的贡献来分块求答案。

主要思想:

使用“交换合式顺序”和“改变枚举变量”来化简。

数论分块

具体就不写了,简要记录证明思路,算了,也不写了。

题目及更多技巧等会写

原文地址:https://www.cnblogs.com/Hikigaya/p/11616447.html

时间: 2024-10-03 23:06:47

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

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

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

狄利克雷卷积 && 莫比乌斯反演 狄利克雷卷积 数论函数及其运算 数论函数是指定义域是正整数,值域是一个数集的函数. 加法,逐项相加,即\((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