引入:
这道题目大家都不陌生吧,可以先动手试试然后在这个提交链接Hcode OnlineJudge提交自己的代码试试,再来往后看。
。
。
。
。
。
如何?如何AC了,恭喜你你掌握了后面的方法可以尝试更高阶的东西,如果失败了请往下看:
你的代码应该与这大同小异,为什么结果会超时呢?看看题目给出的n的数据,1e7,如果使用上面的代码我们循环的次数将近是n*n次,相当大,这时就要用到我们的方法--埃氏筛了。
什么是埃氏筛呢?
筛法,简称埃氏筛或爱氏筛,是一种由希腊数学家埃拉托斯特尼所提出的一种简单检定素数的算法。要得到自然数n以内的全部素数,必须把不大于根号n的所有素数的倍数剔除,剩下的就是素数。
代码如下:从2-sqrt(n)范围内筛选,1、如果是素数则进入第二层循环,将此数的倍数全部置为false即非素数;2、如果是非素数,则代表自己是前面的数的倍数,代表此数的倍数也是前面的的某个数的倍数,所以不用进入第二层循环。IsPrime函数结束后,0-n之间的素数已经筛选完毕,isprime[i]=true代表是素数,false代表不是素数。
#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
bool isprime[1000001];
void IsPrime();
int main(void)
{
memset(isprime, true, sizeof(isprime));
IsPrime();
int n;
cin >> n;
int res = 0;
for (int i = 2; i <= n; i++)
if (isprime[i])
res++;
cout << res;
return 0;
}
void IsPrime()//预处理
{
isprime[0] = isprime[1] = false;
for (int i = 2; i <= sqrt(1e6); i++)
{
if(isprime[i])
for (int j = 2; j * i <= 1e6; j++)
isprime[i * j] = false;
}
}