前言:一般情况下,判断一个数是否为素数,直接对这个数进行简单素性测试,如下
bool is_prime(int n){
for(int i=2;i*i<=n;i++){
if(n%i==0) return false;
}
return n!=1; //1是例外
若要对许多整数进行素性测试,可以用埃氏筛法
- 1既不是素数也不是合数
- 最小的素数为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