一、代码
#include<iostream>
#include<vector>
#include<cmath>
using namespace std;
int main() {
int n; //求n以内的质数
int count=0; //统计质数的个数
cin >> n;
vector<int> prime;
for (int i = 0; i <n + 1; i++) {// 假设给定范围内的所有数都是潜在的质数,标记为1
prime.push_back(1);
}
for (int i = 2; i <= sqrt(n); i++) {
if (prime[i]) {
for (int j = i*i; j <= n; j += i) {//从 i*i 开始剔除倍数是因为所有小于 i*i 的 i 的倍数必定已经被更小的质数筛选过了
prime[j] = 0;
}
}
}
for (int i = 2; i <= n; i++) {
if (prime[i]) {
count++;
cout << i << " ";
}
}
cout << endl;
cout << "质数的个数 = " <<count<< endl;
return 0;
}
二、什么是埃氏筛法?
-
目的:埃氏筛法的主要目的是找出一个给定范围内的所有质数。
-
开始:从最小的质数2开始,假设2是质数。
-
筛选过程:将2的所有倍数(4, 6, 8, ...)标记为非质数,因为它们至少有一个因数2。
-
寻找下一个质数:在完成2的倍数的标记后,寻找列表中下一个未被标记的数。这个数是3,因为除了2之外的所有偶数已经被标记了。
-
重复筛选:将3标记为质数,然后继续筛选,剔除3的所有倍数(9, 15, 21, ...)。
-
继续寻找下一个未标记的数:完成3的筛选后,再次寻找下一个未被标记的数,这次是5,并将5标记为质数。
-
剔除5的倍数:剔除5的所有倍数(25, 35, ...),这些数不是质数。
-
过程重复:重复上述过程,每次筛选后,都寻找下一个未被标记的数,将其标记为质数,并剔除它的所有倍数。
-
终止条件:当当前筛选的质数的下一个数的平方大于你考虑的范围时,停止筛选。例如,如果你正在筛选100以内的质数,当你筛选到13时,13的下一个数是14,而13的平方是169,大于100,所以筛选结束。
-
结果:筛选结束后,所有未被标记的数都是质数。这些数就是给定范围内的所有质数。