为了表述方便,下用 P \mathbb P P 表示素数集,并用 p 1 < p 2 < ⋯ < p c p_1<p_2<\cdots<p_c p1<p2<⋯<pc 表示 1 1 1 到 n n n 内的所有质数, t t t 表示第一个大于 n \sqrt n n 的质数下标.
特别的,我们令 p 0 = 1 p_0=1 p0=1.
算法介绍
一、简介
Min_25
筛是一种能快速求出某些满足特殊条件的积性函数前缀和的亚线性筛法,在通常的递归实现中复杂度为 Θ ( n 1 − ϵ ) \Theta(n^{1-\epsilon}) Θ(n1−ϵ),而经过略复杂的后缀和优化后复杂度甚至能达到 O ( n 3 4 log 2 n ) \Omicron(\frac{n^{\frac34}}{\log_2 n}) O(log2nn43) 级别.
由于后者实现较困难且前者在数据范围较小时效率更高,这里只介绍前者.
二、适用范围
要求所求函数 f ( x ) f(x) f(x) 为积性函数,同时对任意 p ∈ P , e ∈ N + p\in \mathbb{P},e\in\mathbb{N}^+ p∈P,e∈N+ 要能够快速计算 f ( p e ) f(p^e) f(pe) 的值且 f ( p e ) f(p^e) f(pe) 最好能表示成关于 p p p 的低阶多项式.
三、算法流程
为了能快速转移出所有数的和,我们设
S ( m , k ) = ∑ i = 2 m f ( i ) [ min p ∈ P , p ∣ i { p } > p k ] S(m,k)=\sum_{i = 2}^mf(i)\Big[\min_{p\in \mathbb{P},p|i}\{p\} > p_k\Big] S(m,k)=i=2∑mf(i)[p∈P,p∣imin{
p}>pk]
则所求答案为 S ( n , 0 ) S(n,0) S(n,0). 注意到 i i i 是从 2 2 2 开始循环的,这是因为 1 1 1 没有最小质因子,计算时容易引起麻烦,所以我们最后再加上 f ( 1 ) f(1) f(1) 即可.
我们考虑对其递归计算,将其分为质数与合数分别计算贡献。对于合数,我们枚举剩下的数的最小质因子与其次数进行递归转移,且因为我们只算合数,只枚举枚举小于等于 n \sqrt n n 的质数即可;而质数部分似乎无法直接转移,故我们考虑再设出一个类似的函数 h ( m ) h(m) h(m) 表示 f f f 在小于等于 m m m 的所有质数处的和,那么 S S S 的转移式就非常显然了:
S ( m , k ) = h ( m ) − h ( p k ) + ∑ k < j < t ∑ p j e ⩽ m