初试线性筛

今天在做题时遇到一题求最小质因子的,作为菜鸟的我,一上来就将下面的代码敲了上去,

#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
long long ans=0;
int f[10000100];
int main()
{
 int n,j,i,flag;
 cin>>n;
 for(i=2; i<=n; i++)
 {
  flag=1;
  for(j=2; j<=sqrt(i); j++)
  {
   if(i%j==0)//j是i的最小质因子 
   {
    flag=0;
    break;
   }
  }
  if(flag==1)
  ans+=i;//素数的最小质因子就是它本身 
  else
  ans+=j; 
 }
 cout<<ans;
 return 0;
}

果不其然,超了,,正当我苦苦思索时,一位大佬过来瞟了一眼说,这用个线性筛就行了,再不行用欧拉筛。刚开始还不知道,然后在网上找了找。了解了一下。果然大佬就是大佬。以下是针对那道题目而言的欧拉筛的解法:题目连接:https://ac.nowcoder.com/acm/contest/3620/A

#include<iostream>
#include<algorithm>
using namespace std;

int n, i, j,cnt,tmp;
int a[30000001]//初始为0,代表是素数
int r[30000000];
long long ans;//保存答案
int main()
{
 cin >> n;
 for(i=2; i<=n; i++)
 {
  if (!a[i]) //是质数,因为它里面还没有保存过数
        { 
            r[++cnt]=i;  //r[]数组里存的是质数
            a[i] = i;//i就是它本身的最小质因子
        }
        for (j = 1; (tmp=i * r[j]) <= n; j++)//对于i的质数倍
        {
            if (!a[tmp])//是质数
                a[tmp] = r[j];//tmp为质数倍数,不是质数,赋值为tmp的最小质因子,
                              //tmp的最小质因子为min(i的最小质因子,r[j]的最小质因子)
            if (i % r[j] == 0) //r[j]是i的最小因子
                break;
        }
 }
 for(i=3;i<=n;i+=2) ans+=a[i];//将结果加起来,因为a[i]中存的都是最小质因子
 cout << ans+n/2*2;//稍微进行了一下优化,因为上一行中i是从3开始的
 return 0;
}

其实线性筛就是避免了普通筛选中重复的情况,关键代码为(tmp=i * r[j]) 。比如算1~6的最小质因子之和,应该为2+3+2(4的最小质因子为2)+5(5是质数,最小质因子就是它本身)+2(6的最小质因子为2)=14.在普通筛选中,到4和6的时候还要重新算,而线性筛则在算2和3时就已经将这些数的最小质因子算过了。
//--------------------------------------------------------------------------------------------------------------
如果有什么问题请各位大佬及时指正

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值