素数判断/以及素数相关编程题

• 素数判断

bool is_prime(int n) {
	if (n < 2)return false;
	if (n == 2 || n == 3)return true;
	if ((n + 1) % 6 != 0 && (n - 1) % 6 != 0) return false;//5以后的素数肯定在6的倍数旁边，
                                                           // 先排除掉不在6倍数旁边的数字
	else{
		for (int i = 5; i <= sqrt(n); i += 2) {    //从5开始，偶数肯定不是素数，排除。另外一个 
                                                   //数的最大质因子<=n开方
			if (n%i == 0)return false;
		}
	}
	return true;
}

• 找出1-n内的所有素数

void find_all_primes(int n) {
	if (n < 2)return;
	cout << 2 << " ";
	for (int i = 3; i <= n; i += 2) {
		if (is_prime(i)) {
			cout << i << " ";
		}
	}
}

• 给定数字n，找出它的质因子

list<int> get_all_prime_divisors(int n)
{
	list<int> divisor;
	if (n < 2)return divisor;
	if (n == 2 || n == 3) {
		divisor.push_back(n);
		return divisor;
	}
	for (int i = n/2; i > 1; i--) {
		while (1) {
			if (n%i == 0 && is_prime(i)) {
                //这一步把重复的质因子剔除，如果想要完整的分解质因子可以不要这一步
				if (divisor.size() == 0 || i != divisor.back()) {
					divisor.push_back(i);
				}
				n = n / i;
				if (is_prime(n)) {
					divisor.push_back(n);
					return divisor;
				}
			}
			else {
				break;
			}
		}
	}
	return divisor;
}