一般的求解n以内素数的算法
#include <iostream>
#include <cmath>
#define N 10000
int prime[N+1];
int main()
{
int i, j, num = 0;
for(i=2; i<N; i++)
{
for(j=2; j<=sqrt(i); j++)
if( j%i==0 ) break;
if( j>sqrt(i) ) prime[num++] = i;
}
for(i=2; i<100; i++) //由于输出将占用太多io时间,所以只输出2-100内的素数。可以把100改为N
if( prime[i] )
cout << prime[i];
return 0;
}
复杂度是o(n*sqrt(n)),如果n很小的话,这种算法不会耗时很多。
当n很大时,会花费很多时间。所以在程序设计竞赛中就必须要设计出一种更好的算法要求能在几秒钟甚至一秒钟之内找出n以内的所有素数。于是就有了素数筛法。
素数筛选法
1.开一个大的bool型数组prime[],大小就是n+1就可以了.先把所有的下标为奇数的标为true,下标为偶数的标为false.
2.然后用奇数的下标i把是它的倍数的下标对应值设为false
for( i=3; i<=sqrt(n); i+=2 )
{
if(prime[i])
for( j=i+i; j<=n; j+=i )prime[j]=false;
}
3.最后输出bool数组中的值为true的单元的下标,就是所求的n以内的素数了。