ACM_数论a.pptVIP

下面是中学生的做法： bool IsPrime(unsigned n) { for(unsigned i=2;in/2+1;++i) { // 和比它的一半小数相除，如果都除不尽，证明素数 if ( 0 == n % i ) { // 除尽了，合数 return false; } } return true; // 都没除尽，素数 }//复杂度是o(n*sqrt(n))， 中学生的做法 但是当n很大的时候，比如n，n*sqrt(n)30000000000,数量级相当大。在一般的机子它不是一秒钟跑不出结果，它是好几分钟都跑不出结果 -_-！，想锻炼耐心的同学不妨试一试 在程序设计竞赛中就必须要设计出一种更好的算法要求能在几秒钟甚至一秒钟之内找出n以内的所有素数。于是就有了素数筛法。 筛法求素数 void Prime(int max) { memset(flag, true, sizeof(flag)); for(int i = 2; i = max / 2; i++) { if(flag[i]) for(int j = i 1 ; j = max; j += i) flag[j] = false; } for(int i = 2 ; i = max; i++) if(flag[i]) prime[size++] = i; } 筛法求素数 可以再进一步化为二次筛选法，就是欲求n以内的素数，就先把sqrt(n)内的素数求出来，用已经求得的素数来筛出后面的合数 一个合数必然可以由两个或多个质数相乘而得到。那么如果一个数不能被比它的一半小的所有的质数整除，那么比它一半小的所有的合数也一样不可能整除它。 高斯猜测，n以内的素数个数大约与n/ln(n)相当，或者说，当n很大时，两者数量级相同。这就是著名的素数定理。 二分法in乘方 为了讲清这个算法，举一个例子2^7：2*2*2*2*2*2*2 两两分开：(2*2)*(2*2)*(2*2)*2 如果用2*2来计算，那么指数就可以除以2了，不过剩了一个，稍后再单独乘上它。 再次两两分开，指数除以2： ((2*2)*(2*2))*(2*2)*2 实际上最后一个括号里的2 * 2是这回又剩下的，那么，稍后再单独乘上它 现在指数已经为1了，可以计算最终结果了：16*4*2=128 unsigned Power(unsigned n, unsigned p) { // 计算n的p次方 unsigned odd = 1; //oddk用来计算“剩下的”乘积 while (p 1) { // 一直计算到指数小于或等于1 if (( p 1 )!=0) { // 判断p是否奇数，偶数的最低位必为0 odd *= n; // 若odd为奇数，则把“剩下的”乘起来 } n *= n; // 主体乘方 p /= 2; // 指数除以2 } return n * odd; // 最后把主体和“剩下的”乘起来作为结果 } 蒙格马利快速幂模算法 后面我们会用到这样一种运算：(X^Y)%Z 问题是当X和Y很大时，只有32位的整型变量如何能够有效的计算出结果？ 考虑上面那份最终的优化代码和再上面提到过的积模分解公式，我想你也许会猛拍一下脑门，吸口气说：“哦，我懂了！”。 蒙格马利快速幂模算法 下面的讲解是给尚没有做出这样动作的同学们准备的。X^Y可以看作Y个X相乘，即然有积模分解公式，那么我们就可以把Y个X相乘再取模的过程分解开来， 比如：(17^25)%29则可分解为： ( ( 17 * 17 ) % 29 * ( 17 * 17 ) % 29 * …… 如果用上面的代码将这个过程优化，那么我们就得到了著名的“蒙格马利”快速幂模算法： unsigned Montgomery(unsigned n, unsigned p, unsigned m) { // 快速计算 (n ^ e) % m 的值，与power算法极类似 unsigned r = n % m; // 这里的r可不能省 unsigned k = 1; while (p 1) { if ((p 1)!=0) { k = (k * r)