试除法
试除法是我们平时用的最多的方法,也是最简单的方法。顾名思义,试除法就是通过穷举每所有小于欲判断的数n的数字的方法来判断n是否为素数。其代码如下
//判断n是否为素数
#include<iostream>
#include<cmath>
using namespace std;
bool prime(int n){
for (int i=2;i<=sqrt(n);i++)
if ((n % i) == 0) return false;
return true;
}
int main(){
int n;
cin>>n;
for(int i=2;i<=n;i++) if( prime(i))cout<<i<<' ';
return 0;
}
埃氏筛法
原理:如果一个数是素数,则删去它的倍数。
如素数2,删去4,6,8,10…
其代码如下
//求0-n之间的所有质数
#include<iostream>
#include<cmath>
using namespace std;
bool E_sieve(int n){
bool num[n+1];
for(int i=2;i<=n;i++)
num[i]=1;
for (int i=2;i<=sqrt(n);i++)
if (num[i]){//如果num[i]是素数
for(int j=2*i;j<=n;j+=i)
num[j]=false;//num[i]的倍数一定不是素数
}
for(int i=2;i<=n;i++)
if(num[i])
cout<<i<<' ';//打印0-n之间的素数
}
int main(){
int n;
cin>>n;
E_sieve(n); ;
return 0;
}
埃氏筛虽好,但它也不是最完美的,也有瑕疵。对于一个如12,18这样的数而言,它会被多次删除。什么意思呢,如12,可以分解为2 * 6也可以分解为3*4那么在i=2,3的时候各会删除一次,这样就会循环删除。那么有没有什么办法让一个数只被筛选一次呢?有,那就是接下来介绍的方法——欧拉筛
欧拉筛法
原理:保证在 2~n 范围中的每一个合数都能被唯一分解成它的最小质因数与除自己外最大的因数相乘的形式。通俗易懂的来说就是把一个非质数尽可能分解成最小的质因子与最大的因子的乘积,然后再删去它。如12=2*6,删去12。
#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;
const int INF=1e6;
bool num[INF];
int prime[INF],cnt=1;
void Oula(int n){//欧拉筛
for(int i=2;i<=n;i++){//i是用来表示倍数的
if(num[i]) prime[cnt++]=i;
for(int j=1;j<cnt;j++)
{
if(i*prime[j]>n)break;
num[i*prime[j]]=0;
if(i%prime[j]==0)break;//此处难懂请跳转底部
} }
for(int i=1;i<cnt;i++)
cout<<prime[i]<<' '; }
int main(){
int n;
cin>>n;//输出0-n之间的所有素数
memset(num,1,sizeof(num));
Oula(n);
return 0;
}
对于任何i%prime[j]==0需要break可以这么理解,i%prime[j]==0时,说明i能分解成prime[j]*k,那么后面的循环中会有一个数能被分解为i * prime[j+1]或prime[j]*k *prime[j+1],这样一来,最小的质因数显然是prime[j]而不是prime[j+1],就让后面的循环中的某一个倍数值(最大因子)i=prime[j+1]*k
相信这样一来就很清晰了吧。