一般来说超大素数生成都是基于概率的素性检验算法实现的。
这里的基于概率是指这个数可能极大概率是素数,也可能极小概率是合数;但是只要将后者的概率控制在一个非常小的范围内,就能将这个数当作大素数来使用。
一般的大素数生成流程简单的超乎想象,说白了就是挨个数试。试某个数是否是合数,如果是的话就试下一个数;如果不是的话就认定这个数是素数。值得一提的是,这种素数一般被称为伪素数;因为这些被算法判定为非合数的数不一定是真的素数。
说的更清楚一点儿就是:这些算法能够快速的确定大部分数是否是合数,但不能确定所有的数是否是合数;也就是不能确定某个数是否是真的素数。
有很多算法可以实现快速的合数确认;有趣的是,这些算法并不能确定合数的分解因子,但是就是能确认这些数是不是合数。
简单说一下比较常用的几个算法。
第一个算法是基于费马小定理的费马算法。费马小定理是一个很有意思的定理,大意就是随便选一个素数p,再选一个和p不成倍数关系的整数β,必然满足β的p次幂和β对p同余。
那么这个命题的逆否命题就是,不满足β的p次幂和β对p同余的,p必然不是素数。
当然这个算法局限性还是很大的,因为你只能证明某个数不是素数,而无法证明某个数一定是素数。一个经典的伪素数就是561,还有1105和1729,这些“蒙混过关”通过费马算法的伪素数被称为Carmichael整数。
很不幸的是,这一类可以通过费马算法的伪素数有无限个。所以费马算法只适用于初步的简单筛选。
第二个算法是基于欧拉准则和勒让德符号的Solovay-Strassen算法。欧拉准则涉及一个平方剩余的概念,不是很好描述。馒馒来。
首先先说平