先上代码
- 埃氏筛法
#include<bits/stdc++.h>
using namespace std;
//时间复杂度(nloglogn)
//适用10^6下
const int MAX = 1000000;
void FindPrime(vector<int>& prime);
int main() {
vector<int>prime;
FindPrime(prime);
for (auto i : prime)
printf("%d ", i);
return 0;
}
void FindPrime(vector<int>& prime) {
array<int, MAX>flag;
int j = 0;
fill(flag.begin(), flag.end(), 0);
for (int i = 2; i < MAX; i++) {
if (flag[i] == 0) {
prime.push_back(i);
}
for (int k = i + i; k < MAX; k += i) {
flag[k] = 1;
}
}
}
- 欧式筛法
#include<bits/stdc++.h>
using namespace std;
//时间复杂度(n)
//适用低于10^7
const int MAX = 10000000;
void FindPrime(vector<int>& prime);
int main() {
vector<int>prime;
FindPrime(prime);
for (auto i : prime)
printf("%d ", i);
return 0;
}
void FindPrime(vector<int>& prime) {
array<int, MAX>flag;
fill(flag.begin(), flag.end(), 0);
for (int i = 2; i < MAX; i++) {
if (flag[i]==0) {
prime.push_back(i);
}
for (int j = 0; j < prime.size() && i * prime[j] < MAX; j++) {
flag[i * prime[j]] = 1;
if (i % prime[j] == 0) {
break;
}
}
}
}
一些细节: 埃氏和欧式的共同点都是从已知的最小素数 2起手,并将2放入素数数组,然后开始针对2的倍数进行筛选。
其他: 假如是vs编译的话,这么大的数据量需要修改堆栈区的大小,否则一定会抛出异常。
我的修改堆栈区大小为:99999999——》100m
分析
在开始分析之前,请忘了数论的那些术语,质因数等等
知道什么是素数即可
进行手动计算过程
原因是:2分解成乘法,可有12=2×6;12=3×4
而其中最小的质数组合为2×6
也就是说欧式筛法,每次筛的是最小素数的组合,这样保证了一定不会多筛重复的数
从代码来分析
欧式筛法的脱出条件共有3条
1. j < prime.size()
内部脱出:遍历整个素数数组,筛掉i与素数的乘法组合
2. i * prime[j] < MAX
外部进入判断:如果第一次筛选已经超过上界,则脱出
3. if (i % prime[j] == 0)
内部脱出:保证了筛掉的为最小素数的乘法组合(专业术语叫质因数)
其中1和3的组合构成素数筛选的核心判断
2其实可有可无,只为了控制边界而已
后记
如果只是想学习算法,不想做科研,大可不必深入去研究数论的那些术语,大部分算法都可通过手工模拟去理解核心原理,再去理解算法。
但是 假如有数论的基础,相对看起代码也会更快,不至于都去手工模拟。
所以如何学习,根据自己的目标量力而行即可。