【学习笔记】Min25筛

Min25筛是一种求解积性函数前缀和的高效算法,适用于求解∑i=1N[i是质数]∗f(i)。算法时间复杂度为O(LogNN43),空间复杂度为O(N)。它还能辅助计算每个x=⌊iN⌋处的特定值。通过线性筛得到质数及其前缀k次方和,定义并计算g(N,i)和s(N,i),分别用于处理合数和质数的贡献。最后通过模板题【LOJ6053】展示实际应用。" 122968931,9387022,Python二级考试复习:创建10字符字符串并获取指定位置字符,"['Python', '开发语言', '后端']
摘要由CSDN通过智能技术生成

【定理简介】

  • 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 以内的质因数,因此

  • 11
    点赞
  • 30
    收藏
    觉得还不错? 一键收藏
  • 5
    评论
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值