筛选法求素数

      如果整数 n ≠ 0，±1。如果除了显然因数 ±1 和 ±n 以外，n没有其他因数，那么 n 叫做素数。如：2，3，5，7是素数，4，6，10，15是合数。

　　根据素数的性质，有定理：1、设 n 是一个正合数，p 是 n 的一个大于1的最小正因数，则 p 一定是素数。2、设 n 是一个正整数，如果对所有的素数 p ≤ n^1/2，都有 p 不整除 n，则 n 一定是素数。

　　根据两条定理，可以得到寻找素数的群定性方法，通常叫做“埃拉托斯散筛法”。
　　筛法具体描述如下：

      对于任意给定的整数 N ，要求出所有不超过 N 的素数。列出 N 个整数，从中删除小于等于N^1/2的所有素数p1, …, pk的倍数。然后依次删除余下的整数（不包括1）就是所要求的不超过 N 的素数。
　　如果让你用筛法求出2 到某个数之间的所有素数，如何用代码实现呢？

　　思路：用一个数组x[0], x[1], …来存储3,5,7,9,11…這些奇数，因此，x[i]中所有存储的数就是2i+3。将所有元素初始状态设置为未筛除，然后依次去掉其中的合数。

　　现在考虑如何删除合数。从x[0],x[1]…如果x[i]未被筛出，则其为质数，删除其所有倍数，依次往下走。但x[i]中储存的值是2i+3,而2i+3的倍数在数组中的什么地方呢？因为数组中全部是奇数，所以2i+3的倍数可以表示为(2n+1)(2i+3),根据(2n+1)(2i+3) = 2n(2i+3)+2i+3 = 2[n(2i+3)+i]+3，知道x[i]对应的数(即2i+3)的倍数的位置为n(2i+3)+i。当n=1是，值为(2i+3)+i,当n=2时，值为(2i+3)+i+(2i+3), 依次类推，用(2i+3)+i作为初值，删掉对应的x[]，再加上2i+3，删掉对应的数，直到无数可删。因为最后的一个素数可以表示为2N+3,所以筛法求出的为2到N之间的素数。

      代码如下：

 #include <stdio .h> #define MAXSIZE 200 #define DELETED 1 #define KEPT 0 void main(void) { int sieve[MAXSIZE+1]; /* 待筛数组 */ int count = 1; /* 素数总数 */ int prime; int i, k; printf("/n埃拉托斯散筛法"); printf("/n========================="); printf("/n/n2 到 %d 之间的素数为：/n", MAXSIZE*2+3); for (i = 0; i < = MAXSIZE; i++) sieve[i] = KEPT; /* 初始状态未筛出 */ for (i = 0; i <= MAXSIZE; i++) if (sieve[i] == KEPT) { prime = i + i + 3; count++; for (k = prime + i; k <= MAXSIZE; k += prime) sieve[k] = DELETED; //筛出 } printf("/n%6d", 2); for (i = 0, k = 2; i <= MAXSIZE; i++) if (sieve[i] == KEPT) { if (k > 10) { printf("/n"); k = 1; } printf("%6d", 2*i+3); k++; } printf("/n/n总共有 %d 个素数。", count); }