求解一个算法,我们首先要知道它的数学含义。依据这个原则,首先我们要知道什么是素数。; 素数是这样的整数,它除了表示为它自己和1的乘积以外,无论他表示为任何两个整数的乘积。
找素数的方法多种多样。
1:是从2开始用"是则留下,不是则去掉"的方法把所有的数列出来(一直列到你不想再往下列为止,比方说,一直列到10,000)。个数是2,它是一个素数,所以应当把它留下来,然后继续往下数,每隔一个数删去一个数,这样就能把所有能被2整除、因而不是素数的数都去掉。在留下的小的数当中,排在2后面的是3,这是第二个素数,因此应该把它留下,然后从它开始往后数,每隔两个数删去一个,这样就能把所有能被3整除的数全都去掉。下一个未去掉的数是5,然后往后每隔4个数删去一个,以除去所有能被5整除的数。再下一个数是7,往后每隔6个数删去一个;再下一个数是11,往后每隔10个数删一个;再下一个是13,往后每隔12个数删一个。就这样依法做下去。
但是编程我们一般不采用上面的方法,并不说这中方法计算机实现不了,或者说实现算法比较复杂。因为它更像一个数学推理。我们也给一个算法。
下面我们介绍几种长用的编程方法。
2:遍历2以上N的平方根以下的每一个整数,是不是能整除N;(这是基本的方法)
3:遍历2以上N的平方根以下的每一个素数,是不是能整除N;(这个方法是上面方法的改进,但要求N平方根以下的素数已全部知道)
4:采用Rabin-Miller算法进行验算;
例如:N=2^127-1是一个38位数,要验证它是否为素数,用上面几个不同的方法:
验算结果,假设计算机能每秒钟计算1亿次除法,那么
算法2要用4136年,算法3要用93年,算法4只要不到1秒钟!(这些数据是通过计算得到)
另外印度有人宣称素数测试是P问题,我一直没有找到那篇论文,听说里面有很多数学理论。如果那位大人有这篇论文,麻烦转发一份。
下面我们分别实现上面的三种算法:
以下算法我们不涉及内存溢出,以及大数字的问题。如果测试数字超过2^32,发生内存溢出,你需要自己使用策略解决这个问题,在这里只讨论32位机有效数字算法。
1: 算法0:是从2开始用"是则留下,不是则去掉"的方法把所有的数列出来
数组中不为0的数字就是要查找的素数。
void PrimeNumber0()
{
int time GetTickCount();
cout start time time endl;
int Max[MAX_NUMBER]; 在栈上分配,栈上空间要求一般都在2M之间,
如果你需要更大空间,请在堆上申请空间(就是通过malloc,new来申请)。
memset(Max,0,MAX_NUMBER);
for(int i = 0 ; i MAX_NUMBER; ++i)