【定理简介】
- M i n 25 Min25 Min25筛是一种能够求解积性函数 f ( x ) f(x) f(x)的前缀和 ∑ i = 1 N f ( i ) \sum_{i=1}^{N}f(i) ∑i=1Nf(i)的筛法,其前提条件为 ∑ i = 1 N [ i i s a p r i m e ] ∗ f ( i ) \sum_{i=1}^{N}[i\ is\ a\ prime]*f(i) ∑i=1N[i is a prime]∗f(i)可以由 ∑ i = 1 x [ i i s a p r i m e ] ∗ i k \sum_{i=1}^{x}[i\ is\ a\ prime]*i^k ∑i=1x[i is a prime]∗ik简单表示,或能够直接快速计算,其中 k k k为一个较小的常数。
- 其时间复杂度为 O ( N 3 4 L o g N ) O(\frac{N^{\frac{3}{4}}}{LogN}) O(LogNN43),空间复杂度为 O ( N ) O(\sqrt{N}) O(N)。
- 在求解过程中,我们还可以顺带求出对于每一个 x = ⌊ N i ⌋ x=\lfloor\frac{N}{i}\rfloor x=⌊iN⌋, ∑ i = 1 x [ i i s a p r i m e ] ∗ i k \sum_{i=1}^{x}[i\ is\ a\ prime]*i^k ∑i=1x[i is a prime]∗ik的值。在一些题目中,上述功能可以直接帮助我们解决问题。
【算法流程】
首先,我们先来解决对于每一个 x = ⌊ N i ⌋ x=\lfloor\frac{N}{i}\rfloor x=⌊iN⌋,求解 ∑ i = 1 x [ i i s a p r i m e ] ∗ i k \sum_{i=1}^{x}[i\ is\ a\ prime]*i^k ∑i=1x[i is a prime]∗ik的值。
我们先通过线性筛得到 N \sqrt{N} N以内所有的质数 p r i m e i prime_i primei以及它们的前缀 k k k次方和 s u m i = ∑ j = 1 i p r i m e j k sum_{i}=\sum_{j=1}^{i}prime_j^k sumi=∑j=1iprimejk。
定义 g ( N , i ) = ∑ j = 1 N [ j i s a p r i m e o r M i n j > p r i m e i ] ∗ j k g(N,i)=\sum_{j=1}^{N}[j\ is\ a\ prime\ or\ Min_j>prime_i]*j^k g(N,i)=∑j=1N[j is a prime or Minj>primei]∗jk,其中 M i n i Min_i Mini表示 i i i最小的质因数。
直观地来说, g ( N , i ) g(N,i) g(N,i)表示的就是 N N N以内的在埃拉特斯特尼筛算法进行第 i i i轮后尚未被筛去的数的 k k k次方和。
一个合数 X X X一定存在一个 X \sqrt{X} X以内的质因数,因此