挑战程序设计竞赛-求素数算法.pptVIP

挑战程序设计竞赛-求素数算法 素数（prime number） ，又称质数，指在大于1的自然数中，除了1和此整数自身外，无法被其他自然数整除的数（也可定义为只有1和本身两个因数的数）。 比1大但不是素数的数称为合数。1和0既非素数也非合数。素数在数论中有着非常重要的地位。 最小的素数是2，也是素数中唯一的偶数；其他素数都是奇数。素数有无限多个，所以不存在最大的素数。 围绕著素数存在很多问题、猜想和定理。著名的有“孪生素数猜想”和“哥德巴赫猜想”。 关于素数的算法是信息学竞赛和程序设计竞赛中经常出现的基础算法，虽然题目各不相同，但都要涉及验证一个自然数是否为素数的问题。下面探讨寻找一定范围内素数的几个算法。 根据以上思路，可编写以下的试除法函豹： 其实，可以对素数的定义进行进一步的分析。要判断数n是否为素数，不需要用n一直除到n-1才能确认，而只需要除到√n 而就可以了。例如，n=15，则能被15整除的除数有1、3、5，对于除数5就不用判断，因为n被3整除时其商就是5，也就表示n能被5整除了。 验证一个自然数是否为素数，这个问题早在中世纪就引起人们的注意，当时人们还在寻找一个公式可以一劳永逸地解决求素数的问题，直到高斯提出了素数产生的不确定性后，人们才认识到没有一个公式可以简单确认素数，从而人们才转去寻找各种验证素数的方法。 一、试除法 试除法是根据素数的定义得出的一种方法，用需要验证的数n逐个除以从2开始至n-1中的所有数，若能被某一个数整除，表示有一个因数，说明数n不是素数：若直到n-1都不能被整除，则说明该数是素数。 根据以上思路，可编写以下的试除法函数： 试除法判断是否为素数的函数: int is_prime(int n) { int i; if (n=1) return 0; //n不是素数，返回零 for(i=2; i=n; i++) { if( n%i==0 ) return 0; //判断n能否被i整除 } return 1; } 其实，可以对素数的定义进行进一步的分析。要判断数n是否为素数，不需要用n一直除到n-1才能确认，而只需要除到√n 而就可以了。例如，n=15，则能被15整除的除数有1、3、5，对于除数5就不用判断，因为n被3整除时其商就是5，也就表示n能被5整除了。 改进后的是否为素数的函数: int is_prime(int n) { int i; if (n=1) return 0; //n不是素数，返回零 for(i=2; i*i=n; i++) //for(int i = 2; i = sqrt(n); ++i) { if( n%i==0 ) return 0; //判断n能否被i整除 } return 1; } 改进后的程序中，在for循环中以i*i既i的平方与n值进行比较，就可以显著地减少循环的次数，提高验证的效率。 这里用了 i*i = n 来代替 sqrt(n)，可以避免调用函数sqrt()，其消耗时间较多，特别是在大量数据测试的时候消耗很明显。 求出1000以内的所有素数: #include cstdio int is_prime(int n) { for(int i=2;i*i=n;i++) if(n%i==0) return 0; return 1; } int main() { int i, n=1000, num = 0; for(i=2; i=n; i++) { if(is_prime(i)) printf(%3d ,i); } return 0; } 在上面的代码中，通过is_prime()函数来验证指定区间(2～1000)中的每一个数是否为素数，而is_prime()函数中又通过循环进行验证。这种双循环会导致程序执行效率不高。 试除法求解n以内素数的算法。