筛选法判断素数

                                      普通法筛选素数 

                                                            --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))快

每天比昨天更好一些微笑

猜你喜欢

转载自blog.csdn.net/wen2471145970/article/details/80480545