今天在做题时遇到一题求最小质因子的,作为菜鸟的我,一上来就将下面的代码敲了上去,
#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时就已经将这些数的最小质因子算过了。
//--------------------------------------------------------------------------------------------------------------
如果有什么问题请各位大佬及时指正