来自:【数据结构与算法分析——C语言描述】练习2.14
问题描述:Eratosthenes筛是一种用于计算小于N的所有素数的方法。我们从制作整数2到N的表开始。我们找出最小的未被删除的整数i,打印i,然后删除i, 2i, 3i, …, 当i > √N时,算法终止。
首先,没必要做2到N的表,在一个循环内遍历2到N即可。
其次,所谓最小也没必要判断,依次遍历时整数i自然是它到最后一个数之间的最小值。
最后,整数i是否被删除等价于整数i是否素数flag[i]==1或0表示,1表示素数,0表示非素数。
思路:如果整数i是素数,打印它,然后删除它的倍数。i==2时删除2的倍数,等于3时删除3的倍数,等于5时删除5的倍数,… ,直到N为止。
void Eratosthenes(int N) { char * flag; int i, j; flag = (char *)malloc(sizeof(char)*N); memset(flag, '1', sizeof(char)*N); for (i = 2; i < N; i++) { if (flag[i] == '1') { printf("%d ", i); for (j = 2; j*i < N; j++) //删除素数i的倍数 flag[j*i] = '0'; } } free(flag); }
测试结果:
其中数据IO花费了不少的时间。
这是筛选法最简单的思路之一,还可以继续优化,比如:因为除了2所有的偶数都不是素数,所以可以排除一半的数据量。也可以对外层循环只遍历到根号N,可以减小循环的规模。
最新评论