厄拉多塞筛(C语言)
若i为素数,则设置a[i]为1;反之则设置为0.
首先,将所有数组的元素设置为1,表示没有已知的非素数。然后将已知为非素数(即为已知素数的倍数)的索引对应的数组元素设置为0。如果将所有较小的素数的倍数都设置为0之后,a[i]仍然保持为1,则可判断它是所找的素数。
#include<stdio.h>
#define N 10000
int main()
{ int i, j, a[N];
for (i = 2; i < N; i++) a[i] = 1;
for (i = 2; i < N; i++)
if (a[i])
for (j = i; j < N/i; j++) a[i*j] = 0;
for (i = 2; i < N; i++)
if (a[i]) printf("%4d ", i);
printf("\n");
return 0;
}
int main()
{ int i, j, a[N];
for (i = 2; i < N; i++) a[i] = 1;
for (i = 2; i < N; i++)
if (a[i])
for (j = i; j < N/i; j++) a[i*j] = 0;
for (i = 2; i < N; i++)
if (a[i]) printf("%4d ", i);
printf("\n");
return 0;
}
因为程序使用一个数组来包含最简单元素类型,0和1两个值,如果我们使用位的数组,则可以获得更高的空间有效性。
而且,如果N庞大,一些编程环境可能要求数组为全局,或是可以动态分配它。
#include <stdlib.h>
main(int argc, char *argv[])
{ long int i, j, N = atoi(argv[1]);
int *a = malloc(N*sizeof(int));
if (a == NULL)
{ printf("Insufficient memory.\n"); return; }
...
{ long int i, j, N = atoi(argv[1]);
int *a = malloc(N*sizeof(int));
if (a == NULL)
{ printf("Insufficient memory.\n"); return; }
...
作者:cpoint
本文版权归作者和博客园共有,欢迎转载,但未经作者同意必须保留此段声明,且在文章页面明显位置给出原文连接,否则保留追究法律责任的权利.