普通法筛选素数
--o(nlognlogn)
1.存在一个数X>0,则X的倍数不为素数;
例如X=3,则6,9,12不为素数,因为存在因子X=3;
大体思路为: 1.默认所有数为素数
2.从2(因为2为最小的素数)开始循环,如果X之前没有被筛选过,则X为素数,反之反之;
eg:X=4(判断(0,4]的素数)
2.从2开始判断是否为素数--因为2为最小的素数
3.从2开始循环,2的倍数被筛选掉;
4.判断3是否为素数,这里因为3不是2的倍数,所以之前没有被筛选掉,进而得3为素数;
5.从3开始循环,3的倍数被筛选掉;
6.判断4是否为素数,这里4因为是2的倍数,所以被步骤2筛选掉,进而判断出4不是素数并记录;
测试:输出(0,X]的素数
源码:
#include <cstdio> #include <iostream> using namespace std; const int MAX=10000100; //之所以定一百万的长度,是因为一百万以内的数据,普通筛选法比线性筛选法(近似=0(n))快 int prime_table[MAX]={0}; //建立素数表 0代表素数,1代表非素数 void prime(int x) { int i,j; for(i=2;i*i<=x;i++) //从最小的素数2开始判断 { if(!prime_table[i]) //判断是否在之前被筛选过 { for(j=2*i;j<=x;j+=i) //j+=i为j的倍数,因为j=2i; prime_table[j]=1; //将j的整数倍标记为非素数 } } } /* int prime(int x) //原来常用判断是否为素数,效率偏低,但一般够用 { int i; if(x==1) return 0; for(i=2;i*i<=x;i++) if(x%i==0) return 0; return 1; } */ int main() { int x; int i; cin>>x; prime(x); for(i=2;i<=x;i++) { if(!prime_table[i]) //看对应素数表数组元素是否为0,是0则为素数进而则输出 cout<<i<<' '; } return 0; }
之后写线性筛选法。
//之所以定一百万的长度,是因为一百万以内的数据,普通筛选法比线性筛选法(近似=0(n))快
每天比昨天更好一些