积性函数的整理

定义

如果一个数论函数\(f(n)\)满足

\[f(pq)=f(p)f(q),p\perp q\]

则称\(f(n)\)是一个积性函数。

特别的,如果不要求\(p\perp q\)且依然满足上述式子的话,则称\(f(n)\)是一个完全积性函数。

常见积性函数

\[e(n)=[n=1]\]

\[1(n)=1\]

\[\mu(n)=\begin{cases}(-1)^k&n=p_1p_2p_3\dots p_k\\0&n=p^2q\end{cases}\]

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

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

\[id(n)=n\]

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

扫描二维码关注公众号,回复: 6663510 查看本文章

至于它们为什么是积性函数,,,,

我也管不了那么多了。

狄利克雷卷积

定义

\[f*g(n)=\sum_{d|n}f(d)g(\frac nd)\]

要计算的话可以把枚举约数换成枚举倍数(下面会讲到),以调和级数\(O(nlogn)\)的复杂度求出\(f*g\)的前\(n\)项。

交换律

\[f*g(n)=\sum_{d|n}f(d)g(\frac nd)=\sum_{d|n}g(d)f(\frac nd)=g*f(n)\]

结合律

\[f*g*h(n)=\sum_{d|n}f(d)\sum_{t|\frac nd}g(t)h(\frac n{dt})=\sum_{d_1d_2d_3=n}f(d_1)g(d_2)h(d_3)=f*(g*h)(n)\]

常见积性函数的卷积

\[\forall f(n),e*f(n)=f(n)\]

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

\[id*1(n)=\sum_{d|n}d=\sigma(n)\]

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

这个需要特别说明一下。

假设\(n=p_1^{k_1}p_2^{k_2}\dots p_m^{k_m}\),那么上式可以改写成:

\[\mu*1(n)=\sum_{c_1=0}^{k_1}\sum_{c_2=0}^{k_2}\dots\sum_{c_m=0}^{k_m}\mu(p_1^{c_1}p_2^{c_2}\dots p_m^{c_m})\]

观察\(\mu(n)\)的定义,可以发现它与\(n\)的质因子个数有关,而且当某个质因子出现不止一次时,\(\mu(n)=0\)。这样的话,只要\(c_1\)\(c_n\)中有任意一个大于\(1\),后面那个\(\mu\)值就为\(0\)

这样的话,我们就可以大大降低枚举范围:

\[\mu*1(n)=\sum_{c_1=0}^1\sum_{c_2=0}^1\dots\sum_{c_m=0}^1\mu(p_1^{c_1}p_2^{c_2}\dots p_m^{c_m})\]

\[=\sum_{c_1=0}^1\sum_{c_2=0}^1\dots\sum_{c_m=0}^1(-1)^{\sum_{i=1}^mc_i}\]

\[=\sum_{i=0}^m(-1)^i\dbinom mi\]

\[=[m=0]=[n=1]=e(n)\]

至于最后那个式子为什么是\([m=0]\),证法多种多样,这里不再赘述。

那么继续:

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

这个又是为什么?

直接证比较麻烦,我们利用\(\mu*1=e\)证一个反过来的式子:

\[\varphi(n)=id*\mu(n)\]

直接暴力拆式子即可。要用到下面讲到的一些技巧,可以回头再来看。

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

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

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

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

\[=id*\mu(n)\]

这样,我们就完成了上述积性函数的大一统:

\[\mu\xrightarrow{*1}e\xrightarrow{*1}1\xrightarrow{*1}d\]

\[\varphi\xrightarrow{*1}id\xrightarrow{*1}\sigma\]

反过来也可以:

\[\mu\xleftarrow{*\mu}e\xleftarrow{*\mu}1\xleftarrow{*\mu}d\]

\[\varphi\xleftarrow{*\mu}id\xleftarrow{*\mu}\sigma\]

莫比乌斯反演

类似于二项式反演和斯特林反演,有一个这样的反演式子:

\[f(n)=\sum_{d|n}g(d)\iff g(n)=\sum_{d|n}\mu(\frac nd)f(d)\]

其实没毛用,因为它的本质就是这样的:

\[f=g*1\iff g=\mu*f\]

但这样看的话,似乎就是废话了。

常用技巧

无关项提前

本质是分配律。

\[\sum_{i=1}^n\sum_{j=1}^ma_ib_j=\sum_{i=1}^na_i\sum_{j=1}^mb_j\]

交换枚举顺序

\[\sum_{i=1}^na_i\sum_{j=1}^mb_j=\sum_{j=1}^mb_j\sum_{i=1}^na_i\]

这个看上去还是很\(naive\)

比较重要的是枚举约数变成枚举倍数:

\[\sum_{i=1}^na_i\sum_{d|i}b_d=\sum_{d=1}^nb_d\sum_{i=1}^{\lfloor\frac nd\rfloor}a_{id}\]

反演

就是活用上面那几个常用积性函数的狄利克雷卷积关系式,尤其是\(\mu*1=e\)

数论分块

这个居然在我很\(naive\)的时候自己\(yy\)出来了。

举个例子,求\(\sum_{i=0}^n\lfloor\frac ni\rfloor,n\leq10^9\)

不会的话打个表,然后就会了

\(n=100\)时,可以发现\(\lfloor\frac n{100}\rfloor,\lfloor\frac n{99}\rfloor,\dots,\lfloor\frac n{51}\rfloor\)都是\(1\)\(\lfloor\frac n{50}\rfloor,\lfloor\frac n{49}\rfloor,\dots,\lfloor\frac n{34}\rfloor\)都是\(2\),这样的话,我们得到了一个可靠的结论:

\(\lfloor\frac n{\lfloor\frac ni\rfloor+1}\rfloor,\lfloor\frac n{\lfloor\frac ni\rfloor+2}\rfloor,\dots,\lfloor\frac n{\lfloor\frac n{i+1}\rfloor}\rfloor\)的结果都是\(i+1\)

这样我们就得到了一个快速的根号算法:

对于大于\(\sqrt n\)的数,不同的\(\lfloor\frac ni\rfloor\)只有\(\sqrt n\)个;

小于\(\sqrt n\)的数只有\(\sqrt n\)个,直接暴力算即可。

……

猜你喜欢

转载自www.cnblogs.com/star-city/p/11101991.html
今日推荐