这里最近看到了一个方法Sieve of Eratosthenes的方法,突然觉得以前的方法太坑了,我以前都是直接使用穷举法。
先说一下我理解的原理:
首先任何一个大于2的合数都等于一个素数乘以一个正整数,根据这个就可以开始推下面的方法。
问题:求1-n的所有素数
步骤:
1、创建一个n大小的数组A,将1标记为非素数,将2标记为素数,然后将3-n中所有的2的倍数标记为合数,同时将2赋值给K;
2、判断A中大于K的数是否已经全部被标记或者K等于n,如果是空跳到第4步;
3、求得A中最小的未标记的数m,将m标记为素数,将B中所有m的倍数标记为合数,同时将m赋值给K,返回第2步;
4、输出数组A中所有素数。
下面贴出代码:
#include
#include
#define UNLABELED 0 //未标记
#define PRINE 1 //素数
#define COMPOSITE 2 //合数
#define NOPRINE 3 //非素数
#define SIZE 100
int main()
{
int i = 0; //i表示整数和对应的下标
int j = 0; //j表示正要处理的素数j之前的已处理j之后的未处理
int k = 0; //k表示正在处理的j的倍数从2开始到 j * k < SIZE
int *p = NULL; //控制循环
int a[SIZE + 1] = { 0 }; //下标表示整数内容判断是否为素数
for (p = a; p < a + SIZE; ++p) { //初始化数组全是UNLABELED
*p = UNLABELED;
}
p = NULL;
a[0] = a[1] = NOPRINE; //设置前面两个不是素数的数的状态为NOPRINE
i = 2;
while (i < SIZE + 1) {
//判断是否未标记,找到下一个素数
if (a[i] == UNLABELED) {
j = i;
a[i++] = PRINE;
}
else {
//已经标记了继续i++向下走
i++;
continue;
}
for (k = 2; j * k < SIZE; ++k) { //处理素数的倍数
a[j * k] = COMPOSITE;
}
}
for (p = a; p < a + SIZE; ++p) { //打印出素数
if (*p == PRINE) {
printf("%8d", p - a);
}
}
p = NULL;
printf("\n");
system("pause");
return 0;
}
我把运行结果复制一遍吧,之前没仔细看了,计算到一千,但是999没有访问到,所以1000的时候999打出来了,100的时候99打出来了,多谢指点了。
一个是1-100
这个是1-1000
-_-.