前述:有关于数论这一部分也看了不少了,正好今天去看病,路上总结一下有关于素数的一些相关知识,汇总一下我所见过的素数相关判别法,当然有可能不怎么官方的说法,按照自己的理解写的。
素数的定义:对于一个数n,若n除了1和它本身没有其他因子,那么则称n是一个素数。
素数的判别法:
1.定义法:这也就是为什么要再提一遍这个简单定义的原因,定义在数学中是极为重要的,所有的定理都是由公理和定义引出来的。这个方法也就是司空见惯的,从2到sqrt(n)找有没有被n整除的数。那么为什么不搜索超过sqrt(n)的数,若n的其他因子有小于sqrt(n),且大于1的,那么就已经可以判断这个数不是一个素数了,现在假设取x>sqrt(n)能被n整除,由于x<n,且有其其他因子也大于sqrt(n),且大于等于x,假设为y,那么x*y>sqrt(n)*sqrt(n)=n,与二者的乘积也是n的因子矛盾,故搜索的时候只用搜索到sqrt(n)。这里和欧拉函数的求法区别开,欧拉函数的筛法需要判断最后剩下的那个数是否为大于1的素数,是因为欧拉函数的求法是在找他的所有素因子,而不是简单判断它是不是一个素数。当然,最后一个数一定是素数的理由和这里为什么不搜到超过sqrt(n)道理相同。
伪代码:
bool isprime(long long n){
for (int i=2;i<=sqrt(n);i++){
if (n%i==0)return 0;
}
return 1;
}
这里,算法的复杂度为x*sqrt(n),x为要判断的数的个数,显然,不是一个高效率的算法,当需要判断的数个数少而且范围小于10^12可以使用。
2.线性筛法:这个算法可以在较短时间内判断大量的数是不是素数,也是一个被广泛应用的算法。算法想法就是用一个bool数组存这个数是否是素数,从2开始判断,若为素数,那么他的倍数都不是素数,bool取和是素数的值的相反值,想法很简单,应用的是定义的逆否命题,若不是素数,则有大于等于2且不为本身的因子。
伪代码:
void getprime(){
int i,j;
memset (visit,0,sizeof(visit));
for (i=2;i<maxn;i++){
if (!visit[i]){
for (j=2*i;j<maxn;j+=i){
visit[i]=1;
}
}
}
}
3.miller-rabin算法:
很明显,上面的两个算法都是比较稳定的算法,分别用时间或者空间换取了相应的效率,但是,如何判断一个long long范围内的数是否是一个素数,这两个算法都望尘莫及,那么怎么办呢,当然就是这个相对“不靠谱”的算法。这个算法的基础是费马小定理,根据素数的某些特性,随机验证多个数是否符合条件,最后得到这个数是否是素数的结论。