质数筛
- 1不是质数
- 一个数只能被1和自己整除的数叫做质数
- 数N,如果不是质数,那么N一定有一个$<=sqrt(N)$的约数
试除法,一个一个试
E氏筛
任意一个数N的倍数一定不是质数,2N,3N……
求1~N
之间的质数1
2
3
4
5
6
7
8
9
10
11const int maxn =1e6;
int prime[maxn];
for (int i=2;i<=sqrt(n);i++) {
//无论数i是不是质数,但是他的倍数一定是质数
for (int j=i*2;j<=n;j+=i)
prime[j]=1;
}
for (int i=2;i<=n;i++) {
if ( !prime[i] ) cout << "YES" << endl;
else cout << "NO" << endl;
}线性筛
e氏筛 会重复的标记一些合数,每个合数都会被他的最小质因数筛去,所以不会重复筛,复杂度降低
通过从小到大累计质数因子来标记合数
个人理解:通过 当前
数*已推出的质数
来标记合数的,索引为合数,存的是合数的最小质因子
- 扫描不大于v[i]的每一个质数p,在数i的基础上累计一个质因子,应为
p<v[i]
,所以p是合数p*i的最小质因子
v[i]已经是求出的最小质因子的数,那么我们累加此数上累加一个质因子,只要是小于v[i]的质数,都能维持了合数所对应的最小质因子 - 当
i % prime[j] == 0
那么i*prime[j+1]
这个合数肯定会被prime[j]
筛掉,i = k*prime[j]
->i * prime[j+1] = k* prime[j]* prime[j+1]
-> i * prime[j+1] = k1 * prime[j]
1 |
|
本博客所有文章除特别声明外,均采用 CC BY-SA 4.0 协议 ,转载请注明出处!