筛选法,是求不超越自然数N(N>1)的全部质数的一种方法。
具体做法是:
- 把N个自然数按顺序排列起来。
- 1不是质数,要划去。
- 第二个数2是质数留下来,而把2后面2的倍数都划去。
- 2后面第一个没划去的数是3,把3留下,再把3后面3的倍数都划去。
- 3后面第一个没划去的数是5,把5留下,再把5后面5的倍数都划去。
- 这样一直做下去,留下的就是不超越N的全部质数。
第一种代码:
#include <iostream>
using namespace std;
const int m = 100;
int main()
{
int a[m], b[m], n = 0;
//初始化
for (int i = 0; i < m; i++)
{
b[i] = 1;
}
//判断
for (int i = 2; i < m; i++)
{
for (int j = 2; i * j <= m; j++)
{
b[i * j] = 0;
}
if (b[i] != 0)
{
a[n] = i;
n++;
}
}
//输出
cout << n << endl;
for (int i = 0; i < n; i++)
{
cout << a[i] << " ";
}
return 0;
}
代码的解释:我用到两个数组来实现,首先我将b数组全部初始化为0,最关键的一步:
for (int i = 2; i < m; i++)
{
for (int j = 2; i * j <= m; j++)
{
b[i * j] = 0;
}
if (b[i] != 0)
{
a[n] = i;
n++;
}
}
i从2开始,直到m结束,j从2开始,2的倍数分别是: 2 * 3, 2 * 4,2 * 5…
相当于: i * j
所以里面的第一个 for 是用来找出2的倍数,并将下标为2的倍数的b数组对应的元素赋为 0 ;
if语句:此时 i=2 ;将 2 赋给 a[0] ;
第一次只将一个数赋给 a数组 ;
一共执行 m-2 次
第二种代码:
#include <iostream>
using namespace std;
const int m = 100;
int main()
{
int a[m], n = 0;
bool b[m];
//初始化
for (int i = 0; i < m; i++)
{
b[i] = 1;
}
//判断
for (int i = 2; i < m; i++)
{
for (int j = 2; i * j <= m; j++)
{
b[i * j] = 0;
}
if (b[i])
{
a[n] = i;
n++;
}
}
//输出
cout << n << endl;
for (int i = 0; i < n; i++)
{
cout << a[i] << " ";
}
return 0;
}
第二种代码也没什么特别,只是将 int b数组 改为 bool b数组,即布尔数组, 0 是 true ,其他是 false ,所以将 b数组 初始化为 1 时,实际上是将 b数组 初始化为 true
结果:
其他的博文用到 memset() 函数,其实就是初始化的意思,结束。