进阶版找素数：埃氏筛法、区间筛法——《挑战程序设计竞赛》

进阶版找素数：埃氏筛法、区间筛法

接着上一篇初级版找素数，书上讲了一种进阶版找素数，我给这两段代码加了一些注释，方便初学者理解。
上一篇：初级版找素数 : 素数判定、约数枚举、整数分解

埃氏筛法：

给定整数n，请问n以内有多少个素数？

限制条件：
n <= 10⁶

思路：
将2~n放入bool数组中，初始化假设所有数都是素数（值置为true），从 i = 2 开始遍历，如果 i 是素数，将数组中 i 的所有倍数划去（值置为false），然后继续遍历，对于当前遍历的数，如果是素数，就把它的倍数都划去，通过这种方法就能一次枚举n以内的素数。

区间筛法

给定整数a和b，请问区间 [ a , b ) 内有多少个素数？

限制条件：
a < b <= 10¹²
b - a <= 10⁶

思路：
基于上述埃氏筛法，建立两个数组，一个数组在 [2,√b) 中筛素数，另一个数组在 [a,b) 中筛素数，从数组**[2,√b)** 中筛得素数的同时，把它的倍数也从数组 [a,b) 中划去。最后，[a,b) 数组中为true的就是素数。

埃氏筛法代码实现：

#include<iostream>

using namespace std;

const int MAX_N = 1000001;

int n;

int prime[MAX_N];//第i个素数
bool is_prime[MAX_N];//is_prime[i]为true表示i为素数

int sieve(int n) {
	int p = 0;//记录素数的个数
	for (int i = 0; i <= n; i++) {
		is_prime[i] = true;//将每个数都初始化当做素数
	}
	is_prime[0] = is_prime[1] = false;//已知0和1不是素数
	for (int i = 2; i <= n; i++) {//开始遍历
		if (is_prime[i]) {//如果遍历到i的时候is_prime[i]还没有被变成false,就说明i并不是前面任何一个数（除去1以外）的倍数
			prime[p++] = i;//表明第p个素数是i，记录在prime数组中
			for (int j = 2 * i; j <= n; j += i) {//数组后面的数中是i倍数的数都不是素数，直接去掉
				is_prime[j] = false;
			}
		}
	}
	return p;//返回素数的个数p
}

int main() {
	cin >> n;
	int p = sieve(n);
	cout << p << endl;
	/*for (int i = 0; i < p; i++) {//列举出n以内所有素数
		cout << prime[i] << endl;
	}*/
	return 0;
}

/*
1000000
*/

区间筛法代码实现：

#include<iostream>

using namespace std;

const int MAX_L = 1000001;//区间[a,b)的最大长度
const int MAX_SQRT_B = 1000001;//根号b的最大值
typedef long long ll;

ll a, b;
bool is_prime[MAX_L];
bool is_prime_small[MAX_SQRT_B];

void segment_sieve(ll a, ll b) {
	//先初始化两个数组，假设都是素数
	for (int i = 0; (ll)i * i < b; i++)
		is_prime_small[i] = true;
	for (int i = 0; i < b - a; i++)
		is_prime[i] = true;

	//开始筛
	for (int i = 2; (ll)i * i < b; i++) {//循环遍历[2,√b)
		if (is_prime_small[i]) {//如果i是素数
			//先筛[2,根号b)
			for (int j = 2 * i; (ll)j * j < b; j += i)
				is_prime_small[j] = false;//把is_prime_small数组中所有【i的倍数】 都划去

			//再筛[a,b)
			for (ll j = max(2LL, (a + i - 1) / i) * i; j < b; j += i)//这里2LL就是把2转换成long long型
				is_prime[j - a] = false;//把is_prime数组中所有 【i的倍数】 也划去
		}
	}

	int p = 0;//记录素数个数
	for (ll i = 0; i < b-a; i++) {
		if (is_prime[i])//is_prime[i-a] = true表示i是素数
		{
			p++;
			cout << i + a << endl;//遍历输出素数
		}
	}
	cout << p << endl;
}

int main() {
	cin >> a >> b;
	segment_sieve(a, b);
	return 0;
}

/*
22 37
*/