数学知识：质数

试除法判定质数

题解
Y总

给定 n 个正整数 ai，判定每个数是否是质数。

输入格式
第一行包含整数 n。

接下来 n 行，每行包含一个正整数 ai。

输出格式
共 n 行，其中第 i 行输出第 i 个正整数 ai 是否为质数，是则输出 Yes，否则输出 No。

数据范围
1 ≤ n ≤ 100,
1 ≤ ai ≤ 231−1
输入样例：
2
2
6

输出样例：
Yes
No

#include <iostream>
#include <algorithm>
using namespace std;
bool is_prime(int a)
{
	if (a < 2) return false;
	for (int i = 2; i <= a / i; i ++ ) 
		if (a % i == 0)
			return false;
	return true;
}
int main()
{
	int n;
	cin >> n;
	while (n -- ) {
		int a;
		cin >> a;
		if(is_prime(a)) puts("Yes");
		else puts("No");
	}
	return 0;
} 

分解质因数

题解
Y总

给定 n 个正整数 ai，将每个数分解质因数，
并按照质因数从小到大的顺序输出每个质因数的底数和指数。

输入格式
第一行包含整数 n。

接下来 n 行，每行包含一个正整数 ai。

输出格式
对于每个正整数 ai，按照从小到大的顺序输出其分解质因数后，
每个质因数的底数和指数，每个底数和指数占一行。

每个正整数的质因数全部输出完毕后，输出一个空行。

数据范围
1 ≤ n ≤ 100,
1 ≤ ai ≤ 2×109
输入样例：
2
6
8

输出样例：
2 1
3 1

2 3

#include <iostream>
#include <algorithm>
using namespace std;
void divide(int a) 
{
	for (int i = 2; i <= a / i; i ++ ) 
		if (a % i == 0) {
			int s = 0;
			while (a % i == 0) a /= i, s ++ ;
			cout << i << " " << s << endl;
		}
	if (a > 1) cout << a << " " << 1 << endl;
	cout << endl;
}
int main()
{
	int n;
	cin >> n;
	while (n -- ) {
		int a;
		cin >> a;
		divide(a);
	}
	return 0;
}

筛质数

题解
Y总

给定一个正整数 n，请你求出 1∼n 中质数的个数。

输入格式
共一行，包含整数 n。

输出格式
共一行，包含一个整数，表示 1∼n 中质数的个数。

数据范围
1 ≤ n ≤ 106

输入样例：
8

输出样例：
4

#include <iostream>
#include <algorithm>
using namespace std;
const int N = 1000010;
int primes[N], cnt;
bool st[N];

void get_primes(int n)
{
	// 普通 
	for (int i = 2; i <= n; i ++ ) {
		if (!st[i]) primes[cnt ++ ] = i;
		for (int j = i; j <= n; j += i ) 
			st[j] = true;
	}
	
	// 诶氏筛法 
	for (int i = 2; i <= n; i ++ ) {
		if (st[i]) continue;
		primes[cnt ++ ] = i;
		for (int j = i; j <= n; j += i) 
			st[j] = true;
	} 
	
	// 线性筛法
	for (int i = 2; i <= n; i ++ ) {
		if (!st[i]) primes[cnt ++ ] = i;
		// primes[j]<=n/i：变形一下得到 primes[j]*i<=n,把大于n的合数都筛了就没啥意义了
		for (int j = 0; primes[j] <= n / i; j ++ ) {
			// 用最小质因子去筛合数
			st[primes[j] * i] = true;
			// 1)当 i%primes[j]!=0 时，说明此时遍历到的 primes[j] 不是 i 的质因子，
			// 那么只可能是此时的 primes[j]<i 的最小质因子，
			// 所以 primes[j]*i 的最小质因子就是 primes[j]； 
            // 2)当有 i%primes[j]==0 时，说明 i 的最小质因子是 primes[j]，
			// 因此 primes[j]*i 的最小质因子也就应该是
            // prime[j]，之后接着用 st[primes[j+1]*i]=true 去筛合数时，就不是用最小质因子去更新了， 
	     	// 	因为 i 有最小质因子 primes[j]<primes[j+1]，
			 // 此时的 primes[j+1]不是 primes[j+1]*i 的最小质因子， 
			// 此时就应该退出循环，避免之后重复进行筛选。
			if (i % primes[j] == 0) break;
		}
	} 
}
int main()
{
	int n;
	cin >> n;
	get_primes(n);
	cout << cnt;
	return 0;
}