【1】直观判断
根据质数的定义,直接遍历2~n-1,看其中是否有n的因子。
bool is_prim1(int x){
for(int i=2;i<x;i++){
if(x%i==0)return 0;
}
return 1;
}
【2】改进
对于一个数,假设有n=x*y;不难发现,x,y中必有一个数>=sqrt(n),另一个数<=sqrt(n);所以,我们其实只需要遍历2~sqrt(n),当其中不含n的因子时,那么n就是质数。
bool is_prim2(int x){
int tmp=sqrt(x);
for(int i=2;i<=tmp;i++){
if(x%i==0)return 0;
}
return 1;
}
【3】找规律
观察下面一组质数:
可以发现,这些数都是6x+1
或6x+5
形式;因此可以得出一个结论
质数除了2,3之外,其余都是6x+1
或6x+5
形式。
证明
首先,可以确定,任意数都可以表示成6x,6x+1,6x+2,6x+3,6x+4,6x+5的形式。而6x+2必有因子2,6x+3必有因子3,6x+4必有因子2;
所以,质数必定在6x+1
或6x+5
中。
于是可以得出下面的代码:
bool is prim3(int x){
if(x==2||x==3)return 1;
//1
if(x%6!=1||x%6!=5)return 0;
//2.非6x+1、6x+5形式的都是非质数
int tmp=sqrt(x);
for(int i=2;i<=tmp;i++){
if(x%i==0)return 0;
}
return 1;
}
再进一步优化,6x+2,6x+3,6x+4形式都不用判断,必定是非质数,所以循环的步长可以改成6.
bool is prim3(int x){
if(x==2||x==3)return 1;
//1
if(x%6!=1||x%6!=5)return 0;
//2.非6x+1、6x+5形式的都是非质数
int tmp=sqrt(x);
for(int i=5;i<=tmp;i+=6){
if(x%i==0||x%(i+2)==0)return 0;
}
//3.又在6x+1、6x+5形式中,6x+2,+3,+4都必定是偶数,所以不用判断
return 1;
}
over over