数学问题04 素数

一般判断方法

• 证明一个数是否为素数，可用这种方法
• 如果一个数n不能被2~sqrt(n)的任意一个数整除，那这个数就一定是质数

int isPrime(int n) {
	for (int i = 2; i <= (int)sqrt(1.0 * n); ++i) {
		//(int)1.9 = 1;(int)(-1.9) = -1;
		if (n % i == 0) {
			//能被整除
			return 0;
		}
	}
	return 1;
}

埃拉托斯特尼筛法（埃氏筛法）

• https://baike.baidu.com/item/%E5%9F%83%E6%8B%89%E6%89%98%E6%96%AF%E7%89%B9%E5%B0%BC%E7%AD%9B%E6%B3%95/374984?fr=aladdin
• 适合求得自然数n以内的所有素数

example1:

• https://pintia.cn/problem-sets/994805260223102976/problems/994805309963354112
  

#include <iostream>
#include <algorithm>
#define maxSize 1000000
using namespace std;

void getAllPrime(int n, int allPrime[]) {
	//获取前n个素数，保存在allPrime[]中
	int count = 0;
	int* temp = new int[maxSize];
	for (int i = 0; i < maxSize; ++i) {
		temp[i] = 0;
	}
//	int temp[maxSize] = {};//init to zero
	for (int i = 2; i < maxSize; ++i) {
		if (temp[i] == 0) {//0代表是素数
			allPrime[count++] = i;
		}
		if (count == n) {
			return;//已经够了
		}
		for (int j = i + i; j < maxSize; j += i) {
			//j = 2*i,j = 3*i....
			temp[j] = 1;//代表不是素数
		}
	}
}
int main() {
	int *allPrime = new int[10000];//前10000个素数
	int m, n;
	cin >> m >> n;
	getAllPrime(n, allPrime);
	int flag = 0;
	for (int i = m - 1; i < n; ++i) {
		cout << allPrime[i];
		if (i == n - 1) {
			break;
		}
		++flag;
		if (flag == 10) {
			flag = 0;
			cout << endl;
		}
		else {
			cout << " ";
		}
	}
	return 0;

}