有关素数的基础算法

所谓素数，是指恰好有2个约数的整数。因为n的约束都不超过n，所以只要检查2~n-1的所有整数是否整出n就能判定n是不是素数。在此，如果d是n的约数，那么n/d也是n的约数。由n=d×n/d可知min(d,n/d)<=，所以只要检查2~的所有整数就够了。同理可知，整数分解和约数枚举都可以在时间完成，虽然还有更高效的算法（费马测试，算法，数域筛法等），不过多数情况下这已经足够了。——挑战程序设计竞赛p117

素数测试

给定一个整数n，请判断n是不是素数。（1<=n<=10的9次方）

思路：根据上述推导，只需依次检查从1到根号n之间的数i如果n%i==0说明n能被i整除，则n不是素数，返回false。如果所有的i都不整除n则返回true;

代码：

#include<iostream>
#include<vector>
#include<map>
#include<algorithm>
using namespace std;
//假设输入都是正数
//素性测试
bool is_prime(int n) {
	for (int i = 2; i*i <= n; i++) {
		if (n%i == 0)
			return false;
	}
	return n != 1;//1是例外
}
//约数枚举
vector<int> divisor(int n) {
	vector<int> res;
	for (int i = 1; i*i <= n; i++) {
		if (n%i == 0) {
			res.push_back(i);
			if (i != n / i)
				res.push_back(n / i);
		}
	}
	return res;
}
//整数分解，res中值为1的数下标相乘可得n
map<int, int> prime_factor(int n) {
	map<int, int> res;
	for (int i = 2; i*i <= n; i++) {
		while (n%i == 0) {
			++res[i];
			n /= i;
		}
	}
	if (n != 1)
		res[n] = 1;
	return res;
}

埃氏筛法

如果要对许多整数进行素性测试，有更高效的算法，例如当我们要求n以内有多少个素数。

当i是素数时，2*i，3*i，，，k*i 就不是素数，所以一旦检测到i是素数就提前把i的倍数是否是素数设为假

int prime[MAX_N];        //用来保存1到n之间哪些数是素数
bool is_prime[MAX_N + 1];//is_prime[i]为true表示i为素数
//返回n以内素数的个数
int sieve(int n) {
	int p = 0;  //p用以记录素数的个数
	for (int i = 0; i <= n; i++)
		is_prime[i] = true;  //先开始所有的数都为真
	is_prime[0] = is_prime[1] = false; 
	for (int i = 2; i <= n; i++) {
		if (is_prime[i]) {   //最开始2是素数
			prime[p++] = i;//将素数存入prime
			for (int j = 2 * i; j <= n; j += i)//i是素数，则将i的倍数都剔除掉
				is_prime[j] = false;
		}
	}
	return p;//返回素数个数
}

区间筛法

b以内的合数的最小质因子一定不超过根号b，如果有b以内的质数表的话，就可以把埃氏筛法运用在[a,b)上了。也就是说，先分别做好[a,)的表和[a,b)的表，然后从[2,)的表中筛得素数的同时，也将其倍数从[a,b)的表中划去，最后剩下的就是区间[a,b)的素数了。

typedef long long ll;
bool is_prime[MAX_L];
bool is_prime_small[MAX_SQRT_B];
//对区间[a,b)内的整数执行筛法，is_prime[i-a]=true等价于i是素数
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++) {
		if (is_prime_small[i]) {
			for (int j = 2 * i; (ll)j*j < b; j += i)
				is_prime_small[i] = false;  //筛[2,根号b]
			for (ll j = max(2LL, (a + i - 1) / i)*i; j < b; j += i)
				//2LL指将2转换为long long类型，(a+i-1)/i*i可以得到离a最近的i的整数倍
				is_prime[j - a] = false;  //筛[a,b)
		}
	}
}