埃氏筛法(对许多整数进行筛出)

前言:一般情况下,判断一个数是否为素数,直接对这个数进行简单素性测试,如下

bool is_prime(int n){
	for(int i=2;i*i<=n;i++){
		if(n%i==0)	return false;
		}
		return n!=1;			//1是例外

若要对许多整数进行素性测试,可以用埃氏筛法

  1. 1既不是素数也不是合数
  2. 最小的素数为2,进行如下步骤:从2到给定的n的范围内,将2的倍数的数标记(即不是素数),再从下一个未被标记的数开始(即3),重复进行该标记过程。最终未被标记的数都未素数。

代码如下:

#include<iostream>
using namespace std;
int prime[10000];
bool is_prime[10001];
int sieve(int n){
		int p=0;
		for(int i=0;i<=n;i++)	is_prime[i]=true;
		is_prime[0]=is_prime[1]=false;
		for(int i=2;i<=n;i++){
			if(is_prime[i]){
				prime[p++]=i;			//存取素数
				for(int j=2*i;j<=n;j+=i)		//将素数的倍数标记
				is_prime[j]=false;
			}
		}
		return p;
}
int main(){
	int a;
	cin>>a;
	cout<<sieve(a)<<endl;
	for(int i=0;i<sieve(a);i++)
	cout<<prime[i]<<" ";
}

输入11
输出结果为
5
2 3 5 7 11

猜你喜欢

转载自blog.csdn.net/weixin_43556295/article/details/86635630
今日推荐