质数-判断、分解质因子、筛质数

1、质数

（1）>=2的正整数；
（2）只被1和本身整除。

2、判断是否为素数

bool is_prime(int x){
	if(x <= 2) return false;
	for(int i = 2; i * i <= x; ++i){
		if(x % i == 0) return false;
	}
	return true;
}

3、分解质因子

void divide(int x){
	for(int i = 2; i * i <= x; ++i){
		if(x % i == 0){
			int cnt = 0;
			while(x % i == 0){
				cnt++;
				x /= i;
			}
			// 输出质因子个数和相应的质因子
			cout << cnt << ' ' << i << endl;
		}
	}
	// 如果还有大于根号x的质因子，那么输出
	if(x > 1) cout << 1 << ' ' << x << endl;
}

4、筛素数

4.1 埃式筛法

思想：例如对于2，将4， 6， 8. . .进行筛选出来；对于3，将6， 9， 12. . .进行筛选。
时间复杂度O(NloglogN)。
可以看出来6被筛选过两次了，所以就有了线性筛。

const int N = 100005;
int primes[N];
bool st[N]; // st[i]==false表示i是素数
int cnt = 0;

void get_primes(int x){
	for(int i = 2; i <= x; ++i){
		// 用素数来进行筛选
		if(!st[i]){
			primes[cnt++] = i;
			for(int j = 2 * i; j <= x; j += i){
				st[j] = true;
		}
	}
}

4.2 线性筛

线性筛考虑的是用最小质因子来进行筛选，所以不会重复进行筛选，时间复杂度O(N)。
代码中st[i * primes[j]] = true，筛选时就是用最小质因子。
当i % primes[j] == 0，因为primes存储的素数是从小到大进行排序的，所以primes[j]是i的最小质因子，primes[j] 也是 i * primes[j] 的最小质因子，满足了筛选时用最小质因子。
当i % primes[j] != 0，primes[j]比i的所有质因子都要小，原因也是上面所述，那么primes[j]也是i * primes[j] 的最小质因子，也满足筛选时用最小质因子。

代码中if(i % primes[j] == 0) break，当primes[j]是i的最小质因子，如果我们不进行break，继续遍历primes，当遍历到将 i * primes[k] 进行筛选，因为此刻i == primes[j] * p，
元素primes[j] * p * primes[k]被筛选，那么当i变化时，
当i == primes[k] * p时，primes[k] * p * primes[j]再次被筛选，不满足每次用最小质因子进行筛选。

const int N = 100005;
int primes[N];
bool st[N];
int cnt = 0; // cnt记录primes的素数个数

void get_primes(int x){
	for(int i = 2; i <= x; ++i){
		if(!st[i]) primes[cnt++] = i;
		for(int j = 0; j < cnt && i * primes[j] <= x; ++j){
			st[i * primes[j]] = true;
			if(i % primes[j] == 0) break;
		}
	}
}