求解素数几种方法
2006-02-16 15:28
387 查看
[align=center]转贴文章请注明：逸学堂[/align]
求解一个算法,我们首先要知道它的数学含义.依据这个原则,首先我们要知道什么是素数.; 素数是这样的整数,它除了表示为它自己和1的乘积以外,无论他表示为任何两个整数的乘积。
找素数的方法多种多样。
1：是从2开始用“是则留下，不是则去掉”的方法把所有的数列出来（一直列到你不想再往下列为止，比方说，一直列到10，000）。第一个数是2，它是一个素数，所以应当把它留下来，然后继续往下数，每隔一个数删去一个数，这样就能把所有能被2整除、因而不是素数的数都去掉。在留下的最小的数当中，排在2后面的是3，这是第二个素数，因此应该把它留下，然后从它开始往后数，每隔两个数删去一个，这样就能把所有能被3整除的数全都去掉。下一个未去掉的数是5，然后往后每隔4个数删去一个，以除去所有能被5整除的数。再下一个数是7，往后每隔6个数删去一个；再下一个数是11，往后每隔10个数删一个；再下一个是13，往后每隔12个数删一个。就这样依法做下去。
但是编程我们一般不采用上面的方法，并不说这中方法计算机实现不了，或者说实现算法比较复杂。因为它更像一个数学推理。最后我们也给一个算法。
下面我们介绍几种长用的编程方法。
2：遍历２以上Ｎ的平方根以下的每一个整数，是不是能整除Ｎ;(这是最基本的方法）
3：遍历２以上Ｎ的平方根以下的每一个素数，是不是能整除Ｎ;(这个方法是上面方法的改进，但要求N平方根以下的素数已全部知道)
4：采用Rabin-Miller算法进行验算；
例如：Ｎ=２^１２７-１是一个３８位数，要验证它是否为素数，用上面几个不同的方法：
验算结果，假设计算机能每秒钟计算１亿次除法，那么
算法2要用４１３６年，算法3要用９３年，算法4只要不到１秒钟！(这些数据是通过计算得到)
另外印度有人宣称素数测试是P问题，我一直没有找到那篇论文，听说里面有很多数学理论。如果那位大人有这篇论文，麻烦转发一份（E-Mail:exuetang@126.com）有任何问题也可以写信与我讨论。
[align=left] [/align]
[align=left]下面我们分别实现上面的三种算法：[/align]
[align=left]以下算法我们不涉及内存溢出，以及大数字的问题。如果测试数字超过2^32，发生内存溢出，你需要自己使用策略解决这个问题，在这里只讨论32位机有效数字算法。[/align]
[align=left]1：// 算法0：是从2开始用“是则留下，不是则去掉”的方法把所有的数列出来[/align]
[align=left]// 最后数组中不为0的数字就是要查找的素数。[/align]
[align=left]void PrimeNumber0()[/align]
[align=left]{[/align]
[align=left] // int time ::GetTickCount();[/align]
[align=left] // cout << "start time:" << time << endl;[/align]
[align=left] [/align]
[align=left] int Max[MAX_NUMBER]; // 在栈上分配，栈上空间要求一般都在2M之间，[/align]
[align=left] // 如果你需要更大空间，请在堆上申请空间(就是通过malloc，new来申请).[/align]
[align=left] memset(Max,0,MAX_NUMBER);[/align]
[align=left] for(int i = 0 ; i < MAX_NUMBER; ++i)[/align]
[align=left] {[/align]
[align=left] Max[i] = i;[/align]
[align=left] }[/align]
[align=left] int cout = 0;// 记录当前i的位置[/align]
[align=left] [/align]
[align=left] // 遍历整个数组[/align]
[align=left] for(i = 1; i < MAX_NUMBER; ++i)[/align]
[align=left] {[/align]
[align=left] if(Max[i] != 0 )// 如果数据不为0，说明是一个素数[/align]
[align=left] {[/align]
[align=left] int iCout = i;[/align]
[align=left] int j = Max[i];// 记录数组中数组位的数字，以便设置[/align]
[align=left] while((iCout+=j) < MAX_NUMBER)[/align]
[align=left] {[/align]
[align=left] // 把不是素数的数位在数组中置为0[/align]
[align=left] Max[iCout] = 0;[/align]
[align=left] }[/align]
[align=left] ++cout;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] [/align]
[align=left] // int time ::GetTickCount();[/align]
[align=left] // cout << "end time:" << time << endl;[/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]2：这个算法可以修改成为，验证一个给定数字是否是一个素数。[/align]
[align=left]// 因为我们讨论多个算法，所以我们把每个算法都单独[/align]
[align=left]// 写在一个或多个函数内。这些函数并不要求输入值和返回值[/align]
[align=left]// 如果你需要这些结果，可以自己修改。[/align]
[align=left] [/align]
[align=left]// 算法1：遍历２以上Ｎ的平方根以下的每一个整数，是不是能整除Ｎ;[/align]
[align=left]void PrimeNumber1()[/align]
[align=left]{[/align]
[align=left]// int time ::GetTickCount();[/align]
[align=left]// cout << "start time:" << time << endl;[/align]
[align=left] [/align]
[align=left] int Max[MAX_NUMBER/2]; // 在栈上分配，栈上空间要求一般都在2M之间，[/align]
[align=left]// 如果你需要更大空间，请在堆上申请空间(就是通过malloc，new来申请).素数的个数很少[/align]
[align=left]// 所以没有必要申请和所求数字同样大小的空间。[/align]
[align=left] memset(Max,0,MAX_NUMBER);[/align]
[align=left] Max[0] = 2;// 放入第一个素数，有人说2不是素数，如果你是其中一员，就改成3吧[/align]
[align=left] int cout = 1;// 记录素数个数[/align]
[align=left] [/align]
[align=left] // 挨个数进行验证[/align]
[align=left] bool bflag = true;[/align]
[align=left] for(int i = 3; i < MAX_NUMBER; ++i)[/align]
[align=left] {[/align]
[align=left] bflag = true;[/align]
[align=left] // 需要是使用数学库（math.h）中sqrt[/align]
[align=left] int iTemp = (int)sqrt((float)i);// 强制转换成int类型,有的人在这里使用i＋1就是为了增加sqrt的精度[/align]
[align=left] // 没有特殊函数，你也可以使用int iTemp = (int)sqrt(i)＋1;来提高进度[/align]
[align=left] for (int j = 2; j < iTemp; ++j)[/align]
[align=left] {[/align]
[align=left] if(i%j == 0)// 求余，如果为0说明，可以整除，不是素数。[/align]
[align=left] {[/align]
[align=left] bflag = false;[/align]
[align=left] break;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] // 经过验证是素数，放入数组。[/align]
[align=left] if(bflag)[/align]
[align=left] {[/align]
[align=left] Max[cout++] = i;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] [/align]
[align=left]// int time ::GetTickCount();[/align]
[align=left]// cout << "end time:" << time << endl;[/align]
[align=left] [/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]3：这个方法是上面方法的改进，但要求N平方根以下的素数已全部知道[/align]
[align=left]// 算法2：遍历２以上Ｎ的平方根以下的每一个素数，是不是能整除Ｎ;[/align]
[align=left]// (这个方法是上面方法的改进，但要求N平方根以下的素数已全部知道)[/align]
[align=left]void PrimeNumber2()[/align]
[align=left]{[/align]
[align=left] // int time ::GetTickCount();[/align]
[align=left] // cout << "start time:" << time << endl;[/align]
[align=left] [/align]
[align=left] int Max[MAX_NUMBER/2]; // 在栈上分配，栈上空间要求一般都在2M之间，[/align]
[align=left] // 如果你需要更大空间，请在堆上申请空间(就是通过malloc，new来申请).素数的个数很少[/align]
[align=left] // 所以没有必要申请和所求数字同样大小的空间。[/align]
[align=left] memset(Max,0,MAX_NUMBER);[/align]
[align=left] Max[0] = 2;// 放入第一个素数，有人说2不是素数，如果你是其中一员，就改成3吧[/align]
[align=left] int cout = 1;// 记录素数个数[/align]
[align=left] [/align]
[align=left] // 挨个数进行验证[/align]
[align=left] bool bflag = true;[/align]
[align=left] for(int i = 3; i < MAX_NUMBER; ++i)[/align]
[align=left] {[/align]
[align=left] bflag = true;[/align]
[align=left] // 需要是使用数学库（math.h）中sqrt[/align]
[align=left] int iTemp = (int)sqrt((float)i);// 强制转换成int类型,有的人在这里使用i＋1就是为了增加sqrt的精度[/align]
[align=left] // 没有特殊函数，你也可以使用int iTemp = (int)sqrt(i)＋1;来提高进度[/align]
[align=left]/////////////////////////////////////////////////////////////////////////////////[/align]
[align=left] // 修改的是这里以下的部分[/align]
[align=left] for (int j = 0; j < cout; ++j)[/align]
[align=left] {[/align]
[align=left] if(i%Max[j] == 0)// 求余，如果为0说明，可以整除，不是素数。[/align]
[align=left] {[/align]
[align=left] bflag = false;[/align]
[align=left] break;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] // 修改的是这里以上的部分[/align]
[align=left]//////////////////////////////////////////////////////////////////////////////////[/align]
[align=left] // 经过验证是素数，放入数组。[/align]
[align=left] if(bflag)[/align]
[align=left] {[/align]
[align=left] Max[cout++] = i;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] [/align]
[align=left] // int time ::GetTickCount();[/align]
[align=left] // cout << "end time:" << time << endl;[/align]
[align=left] [/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]4：采用Rabin-Miller算法进行验算，Rabin-Miller算法是典型的验证一个数字是否为素数的方法。判断素数的方法是Rabin-Miller概率测试，那么他具体的流程是什么呢。假设我们要判断n是不是素数，首先我们必须保证n 是个奇数，那么我们就可以把n 表示为 n = (2^r)*s+1,注意s 也必须是一个奇数。然后我们就要选择一个随机的整数a (1<=a<=n-1)，接下来我们就是要判断 a^s=1 (mod n) 或a^((2^j)*s)= -1（mod n）(0<=j<R),如果任意一式成立，我们就说n通过了测试，但是有可能不是素数也能通过测试。所以我们通常要做多次这样的测试，以确保我们得到的是一个素数。（DDS的标准是要经过50次测试） [/align]
[align=left] [/align]
[align=left]// 算法3：采用Rabin-Miller算法进行验算[/align]
[align=left]//首先选择一个代测的随机数p，计算b，b是2整除p-1的次数。然后计算m，使得n=1+(2^b)m。[/align]
[align=left] [/align]
[align=left]//（1） 选择一个小于p的随机数a。[/align]
[align=left]//（2） 设j=0且z=a^m mod p[/align]
[align=left]//（3） 如果z=1或z=p-1，那麽p通过测试，可能使素数[/align]
[align=left]//（4） 如果j>0且z=1, 那麽p不是素数[/align]
[align=left]//（5） 设j=j+1。如果j<b且z<>p-1,设z=z^2 mod p，然后回到(4)。如果z=p-1，那麽p通过测试，可能为素数。[/align]
[align=left]//（6） 如果j=b 且z<>p-1，不是素数[/align]
[align=left] [/align]
[align=left]// 判定是否存在 a^s=1 (mod n) 或a^((2^j)*s)= -1（mod n）(0<=j<R),[/align]
[align=left] [/align]
[align=left]bool Witness(int a,int n)[/align]
[align=left]{[/align]
[align=left] // 解释一下数学词汇：[/align]
[align=left] // ceil求不小于x的最小整数，函数原型extern float ceil(float x);求得i的最大值[/align]
[align=left] // log计算x的自然对数,函数原型extern float log(float x);[/align]
[align=left] long i,d=1,x;[/align]
[align=left] for (i=(int)ceil(log((double)n-1)/log(2.0))-1;i>=0;--i) [/align]
[align=left] {[/align]
[align=left] x=d;[/align]
[align=left] d=(d*d)%n;[/align]
[align=left] if ((1==d) && (x!=1) && (x!=n-1))[/align]
[align=left] {[/align]
[align=left] return 1;[/align]
[align=left] }[/align]
[align=left] if ((n-1)&(1<0)) [/align]
[align=left] {[/align]
[align=left] d=(d*a)%n;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] return (d!=1);[/align]
[align=left] [/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]// 参数n，是要测定的数字，s是要内部测试的次数。[/align]
[align=left]bool Rabin_Miller(int n,int s)[/align]
[align=left]{[/align]
[align=left] for (int j = 0;j < s; ++j)[/align]
[align=left] {[/align]
[align=left] int a = rand()*(n-2)/RAND_MAX + 1;// 获得一个随机数1<=a<=n-1[/align]
[align=left] if (Witness(a,n))// 利用这个随即数和n进行判断对比，只要有一次返回true，就说明n不是一个素数[/align]
[align=left] {[/align]
[align=left] return false;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] return true;// 通过验证是一个素数[/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]// 算法3：采用Rabin-Miller算法进行验算[/align]
[align=left]// 这个算法是求大素数使用的。所以你的必须想办法支持大数字运算，[/align]
[align=left]// 不然极易造成内存访问失效,我在我的机子上，MAX_NUMBER＝10000时就会出现问题，1000就没有问题[/align]
[align=left]void PrimeNumber3()[/align]
[align=left]{[/align]
[align=left] [/align]
[align=left] int Max[MAX_NUMBER/2];// 在栈上分配，栈上空间要求一般都在2M之间，[/align]
[align=left] // 如果你需要更大空间，请在堆上申请空间(就是通过malloc，new来申请).素数的个数很少[/align]
[align=left] // 所以没有必要申请和所求数字同样大小的空间。[/align]
[align=left] int cout = 0;// 记录素数个数[/align]
[align=left] memset(Max,0,MAX_NUMBER/2);[/align]
[align=left] [/align]
[align=left] for(int i = 2; i < 1000; ++i)[/align]
[align=left] {[/align]
[align=left] if(Rabin_Miller(i,20))[/align]
[align=left] {[/align]
[align=left] Max[cout++] = i;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]以上程序都经过测试，测试环境Window 2003＋VC7.1[/align]
[align=left]小结：以上只是简单介绍一下求素数的几种常用方法，实际方法远远不知这些。对于普通用户来说算法2就可以方便应用了。Rabin－miller算法，确实很高效，但是如果你对数学，编程知识不是很熟悉，建议不要使用。[/align]
[align=left]如果使用过程中有任何疑问请与我联系。[/align]
[align=left]QQ：35091551[/align]
[align=left]MSN：ugg_xchj@163.com[/align]
[align=left]下载源文件http://www.exuetang.net/article/NewsViewAdmin.aspx?NewsId=191[/align]
求解一个算法,我们首先要知道它的数学含义.依据这个原则,首先我们要知道什么是素数.; 素数是这样的整数,它除了表示为它自己和1的乘积以外,无论他表示为任何两个整数的乘积。
找素数的方法多种多样。
1：是从2开始用“是则留下，不是则去掉”的方法把所有的数列出来（一直列到你不想再往下列为止，比方说，一直列到10，000）。第一个数是2，它是一个素数，所以应当把它留下来，然后继续往下数，每隔一个数删去一个数，这样就能把所有能被2整除、因而不是素数的数都去掉。在留下的最小的数当中，排在2后面的是3，这是第二个素数，因此应该把它留下，然后从它开始往后数，每隔两个数删去一个，这样就能把所有能被3整除的数全都去掉。下一个未去掉的数是5，然后往后每隔4个数删去一个，以除去所有能被5整除的数。再下一个数是7，往后每隔6个数删去一个；再下一个数是11，往后每隔10个数删一个；再下一个是13，往后每隔12个数删一个。就这样依法做下去。
但是编程我们一般不采用上面的方法，并不说这中方法计算机实现不了，或者说实现算法比较复杂。因为它更像一个数学推理。最后我们也给一个算法。
下面我们介绍几种长用的编程方法。
2：遍历２以上Ｎ的平方根以下的每一个整数，是不是能整除Ｎ;(这是最基本的方法）
3：遍历２以上Ｎ的平方根以下的每一个素数，是不是能整除Ｎ;(这个方法是上面方法的改进，但要求N平方根以下的素数已全部知道)
4：采用Rabin-Miller算法进行验算；
例如：Ｎ=２^１２７-１是一个３８位数，要验证它是否为素数，用上面几个不同的方法：
验算结果，假设计算机能每秒钟计算１亿次除法，那么
算法2要用４１３６年，算法3要用９３年，算法4只要不到１秒钟！(这些数据是通过计算得到)
另外印度有人宣称素数测试是P问题，我一直没有找到那篇论文，听说里面有很多数学理论。如果那位大人有这篇论文，麻烦转发一份（E-Mail:exuetang@126.com）有任何问题也可以写信与我讨论。
[align=left] [/align]
[align=left]下面我们分别实现上面的三种算法：[/align]
[align=left]以下算法我们不涉及内存溢出，以及大数字的问题。如果测试数字超过2^32，发生内存溢出，你需要自己使用策略解决这个问题，在这里只讨论32位机有效数字算法。[/align]
[align=left]1：// 算法0：是从2开始用“是则留下，不是则去掉”的方法把所有的数列出来[/align]
[align=left]// 最后数组中不为0的数字就是要查找的素数。[/align]
[align=left]void PrimeNumber0()[/align]
[align=left]{[/align]
[align=left] // int time ::GetTickCount();[/align]
[align=left] // cout << "start time:" << time << endl;[/align]
[align=left] [/align]
[align=left] int Max[MAX_NUMBER]; // 在栈上分配，栈上空间要求一般都在2M之间，[/align]
[align=left] // 如果你需要更大空间，请在堆上申请空间(就是通过malloc，new来申请).[/align]
[align=left] memset(Max,0,MAX_NUMBER);[/align]
[align=left] for(int i = 0 ; i < MAX_NUMBER; ++i)[/align]
[align=left] {[/align]
[align=left] Max[i] = i;[/align]
[align=left] }[/align]
[align=left] int cout = 0;// 记录当前i的位置[/align]
[align=left] [/align]
[align=left] // 遍历整个数组[/align]
[align=left] for(i = 1; i < MAX_NUMBER; ++i)[/align]
[align=left] {[/align]
[align=left] if(Max[i] != 0 )// 如果数据不为0，说明是一个素数[/align]
[align=left] {[/align]
[align=left] int iCout = i;[/align]
[align=left] int j = Max[i];// 记录数组中数组位的数字，以便设置[/align]
[align=left] while((iCout+=j) < MAX_NUMBER)[/align]
[align=left] {[/align]
[align=left] // 把不是素数的数位在数组中置为0[/align]
[align=left] Max[iCout] = 0;[/align]
[align=left] }[/align]
[align=left] ++cout;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] [/align]
[align=left] // int time ::GetTickCount();[/align]
[align=left] // cout << "end time:" << time << endl;[/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]2：这个算法可以修改成为，验证一个给定数字是否是一个素数。[/align]
[align=left]// 因为我们讨论多个算法，所以我们把每个算法都单独[/align]
[align=left]// 写在一个或多个函数内。这些函数并不要求输入值和返回值[/align]
[align=left]// 如果你需要这些结果，可以自己修改。[/align]
[align=left] [/align]
[align=left]// 算法1：遍历２以上Ｎ的平方根以下的每一个整数，是不是能整除Ｎ;[/align]
[align=left]void PrimeNumber1()[/align]
[align=left]{[/align]
[align=left]// int time ::GetTickCount();[/align]
[align=left]// cout << "start time:" << time << endl;[/align]
[align=left] [/align]
[align=left] int Max[MAX_NUMBER/2]; // 在栈上分配，栈上空间要求一般都在2M之间，[/align]
[align=left]// 如果你需要更大空间，请在堆上申请空间(就是通过malloc，new来申请).素数的个数很少[/align]
[align=left]// 所以没有必要申请和所求数字同样大小的空间。[/align]
[align=left] memset(Max,0,MAX_NUMBER);[/align]
[align=left] Max[0] = 2;// 放入第一个素数，有人说2不是素数，如果你是其中一员，就改成3吧[/align]
[align=left] int cout = 1;// 记录素数个数[/align]
[align=left] [/align]
[align=left] // 挨个数进行验证[/align]
[align=left] bool bflag = true;[/align]
[align=left] for(int i = 3; i < MAX_NUMBER; ++i)[/align]
[align=left] {[/align]
[align=left] bflag = true;[/align]
[align=left] // 需要是使用数学库（math.h）中sqrt[/align]
[align=left] int iTemp = (int)sqrt((float)i);// 强制转换成int类型,有的人在这里使用i＋1就是为了增加sqrt的精度[/align]
[align=left] // 没有特殊函数，你也可以使用int iTemp = (int)sqrt(i)＋1;来提高进度[/align]
[align=left] for (int j = 2; j < iTemp; ++j)[/align]
[align=left] {[/align]
[align=left] if(i%j == 0)// 求余，如果为0说明，可以整除，不是素数。[/align]
[align=left] {[/align]
[align=left] bflag = false;[/align]
[align=left] break;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] // 经过验证是素数，放入数组。[/align]
[align=left] if(bflag)[/align]
[align=left] {[/align]
[align=left] Max[cout++] = i;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] [/align]
[align=left]// int time ::GetTickCount();[/align]
[align=left]// cout << "end time:" << time << endl;[/align]
[align=left] [/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]3：这个方法是上面方法的改进，但要求N平方根以下的素数已全部知道[/align]
[align=left]// 算法2：遍历２以上Ｎ的平方根以下的每一个素数，是不是能整除Ｎ;[/align]
[align=left]// (这个方法是上面方法的改进，但要求N平方根以下的素数已全部知道)[/align]
[align=left]void PrimeNumber2()[/align]
[align=left]{[/align]
[align=left] // int time ::GetTickCount();[/align]
[align=left] // cout << "start time:" << time << endl;[/align]
[align=left] [/align]
[align=left] int Max[MAX_NUMBER/2]; // 在栈上分配，栈上空间要求一般都在2M之间，[/align]
[align=left] // 如果你需要更大空间，请在堆上申请空间(就是通过malloc，new来申请).素数的个数很少[/align]
[align=left] // 所以没有必要申请和所求数字同样大小的空间。[/align]
[align=left] memset(Max,0,MAX_NUMBER);[/align]
[align=left] Max[0] = 2;// 放入第一个素数，有人说2不是素数，如果你是其中一员，就改成3吧[/align]
[align=left] int cout = 1;// 记录素数个数[/align]
[align=left] [/align]
[align=left] // 挨个数进行验证[/align]
[align=left] bool bflag = true;[/align]
[align=left] for(int i = 3; i < MAX_NUMBER; ++i)[/align]
[align=left] {[/align]
[align=left] bflag = true;[/align]
[align=left] // 需要是使用数学库（math.h）中sqrt[/align]
[align=left] int iTemp = (int)sqrt((float)i);// 强制转换成int类型,有的人在这里使用i＋1就是为了增加sqrt的精度[/align]
[align=left] // 没有特殊函数，你也可以使用int iTemp = (int)sqrt(i)＋1;来提高进度[/align]
[align=left]/////////////////////////////////////////////////////////////////////////////////[/align]
[align=left] // 修改的是这里以下的部分[/align]
[align=left] for (int j = 0; j < cout; ++j)[/align]
[align=left] {[/align]
[align=left] if(i%Max[j] == 0)// 求余，如果为0说明，可以整除，不是素数。[/align]
[align=left] {[/align]
[align=left] bflag = false;[/align]
[align=left] break;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] // 修改的是这里以上的部分[/align]
[align=left]//////////////////////////////////////////////////////////////////////////////////[/align]
[align=left] // 经过验证是素数，放入数组。[/align]
[align=left] if(bflag)[/align]
[align=left] {[/align]
[align=left] Max[cout++] = i;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] [/align]
[align=left] // int time ::GetTickCount();[/align]
[align=left] // cout << "end time:" << time << endl;[/align]
[align=left] [/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]4：采用Rabin-Miller算法进行验算，Rabin-Miller算法是典型的验证一个数字是否为素数的方法。判断素数的方法是Rabin-Miller概率测试，那么他具体的流程是什么呢。假设我们要判断n是不是素数，首先我们必须保证n 是个奇数，那么我们就可以把n 表示为 n = (2^r)*s+1,注意s 也必须是一个奇数。然后我们就要选择一个随机的整数a (1<=a<=n-1)，接下来我们就是要判断 a^s=1 (mod n) 或a^((2^j)*s)= -1（mod n）(0<=j<R),如果任意一式成立，我们就说n通过了测试，但是有可能不是素数也能通过测试。所以我们通常要做多次这样的测试，以确保我们得到的是一个素数。（DDS的标准是要经过50次测试） [/align]
[align=left] [/align]
[align=left]// 算法3：采用Rabin-Miller算法进行验算[/align]
[align=left]//首先选择一个代测的随机数p，计算b，b是2整除p-1的次数。然后计算m，使得n=1+(2^b)m。[/align]
[align=left] [/align]
[align=left]//（1） 选择一个小于p的随机数a。[/align]
[align=left]//（2） 设j=0且z=a^m mod p[/align]
[align=left]//（3） 如果z=1或z=p-1，那麽p通过测试，可能使素数[/align]
[align=left]//（4） 如果j>0且z=1, 那麽p不是素数[/align]
[align=left]//（5） 设j=j+1。如果j<b且z<>p-1,设z=z^2 mod p，然后回到(4)。如果z=p-1，那麽p通过测试，可能为素数。[/align]
[align=left]//（6） 如果j=b 且z<>p-1，不是素数[/align]
[align=left] [/align]
[align=left]// 判定是否存在 a^s=1 (mod n) 或a^((2^j)*s)= -1（mod n）(0<=j<R),[/align]
[align=left] [/align]
[align=left]bool Witness(int a,int n)[/align]
[align=left]{[/align]
[align=left] // 解释一下数学词汇：[/align]
[align=left] // ceil求不小于x的最小整数，函数原型extern float ceil(float x);求得i的最大值[/align]
[align=left] // log计算x的自然对数,函数原型extern float log(float x);[/align]
[align=left] long i,d=1,x;[/align]
[align=left] for (i=(int)ceil(log((double)n-1)/log(2.0))-1;i>=0;--i) [/align]
[align=left] {[/align]
[align=left] x=d;[/align]
[align=left] d=(d*d)%n;[/align]
[align=left] if ((1==d) && (x!=1) && (x!=n-1))[/align]
[align=left] {[/align]
[align=left] return 1;[/align]
[align=left] }[/align]
[align=left] if ((n-1)&(1<0)) [/align]
[align=left] {[/align]
[align=left] d=(d*a)%n;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] return (d!=1);[/align]
[align=left] [/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]// 参数n，是要测定的数字，s是要内部测试的次数。[/align]
[align=left]bool Rabin_Miller(int n,int s)[/align]
[align=left]{[/align]
[align=left] for (int j = 0;j < s; ++j)[/align]
[align=left] {[/align]
[align=left] int a = rand()*(n-2)/RAND_MAX + 1;// 获得一个随机数1<=a<=n-1[/align]
[align=left] if (Witness(a,n))// 利用这个随即数和n进行判断对比，只要有一次返回true，就说明n不是一个素数[/align]
[align=left] {[/align]
[align=left] return false;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left] return true;// 通过验证是一个素数[/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]// 算法3：采用Rabin-Miller算法进行验算[/align]
[align=left]// 这个算法是求大素数使用的。所以你的必须想办法支持大数字运算，[/align]
[align=left]// 不然极易造成内存访问失效,我在我的机子上，MAX_NUMBER＝10000时就会出现问题，1000就没有问题[/align]
[align=left]void PrimeNumber3()[/align]
[align=left]{[/align]
[align=left] [/align]
[align=left] int Max[MAX_NUMBER/2];// 在栈上分配，栈上空间要求一般都在2M之间，[/align]
[align=left] // 如果你需要更大空间，请在堆上申请空间(就是通过malloc，new来申请).素数的个数很少[/align]
[align=left] // 所以没有必要申请和所求数字同样大小的空间。[/align]
[align=left] int cout = 0;// 记录素数个数[/align]
[align=left] memset(Max,0,MAX_NUMBER/2);[/align]
[align=left] [/align]
[align=left] for(int i = 2; i < 1000; ++i)[/align]
[align=left] {[/align]
[align=left] if(Rabin_Miller(i,20))[/align]
[align=left] {[/align]
[align=left] Max[cout++] = i;[/align]
[align=left] }[/align]
[align=left] }[/align]
[align=left]}[/align]
[align=left] [/align]
[align=left]以上程序都经过测试，测试环境Window 2003＋VC7.1[/align]
[align=left]小结：以上只是简单介绍一下求素数的几种常用方法，实际方法远远不知这些。对于普通用户来说算法2就可以方便应用了。Rabin－miller算法，确实很高效，但是如果你对数学，编程知识不是很熟悉，建议不要使用。[/align]
[align=left]如果使用过程中有任何疑问请与我联系。[/align]
[align=left]QQ：35091551[/align]
[align=left]MSN：ugg_xchj@163.com[/align]
[align=left]下载源文件http://www.exuetang.net/article/NewsViewAdmin.aspx?NewsId=191[/align]
内容来自用户分享和网络整理，不保证内容的准确性，如有侵权内容，可联系管理员处理 